1
3 return ["vx[%d]" % i for i in range(d)]
5 num_vertices = d+1
6 pre = "const double vx[%d] = { " % (num_vertices*d,)
7 lines = [", ".join("c.coordinates[%d][%d]" % (i,j) for j in range(d))
8 for i in range(num_vertices)]
9 sep = "\n" + (" "*len(pre))
10 expr = sep.join(lines)
11 post = "\n" + (" "*len(pre)) + "};"
12 code = "%s%s%s" % (pre, expr, post)
13 return [code]
14
16 return ["x0[%d]" % i for i in range(d)]
18 expr = ", ".join("1.0"
19 for i in range(d))
20 code = "const double x0[%d] = { %s };" % (d, expr)
21 return [code]
22
24 return ["G[%d]" % i for i in range(d*d)]
26 expr = ", ".join("1.0" if i == j else "0.0"
27 for i in range(d)
28 for j in range(d))
29 code = "const double G[%d] = { %s };" % (d*d, expr)
30 return [code]
31
35 expr = "1.0"
36 code = "const double detG = %s;" % expr
37 return [code]
38
40 return ["Ginv[%d]" % i for i in range(d*d)]
42 expr = ", ".join("1.0" if i == j else "0.0"
43 for i in range(d)
44 for j in range(d))
45 code = "const double Ginv[%d] = { %s };" % (d*d, expr)
46 return [code]
47
51 expr = "1.0"
52 code = "const double detGsign = %s;" % expr
53 return [code]
54
58 expr = "1.0"
59 code = "const double D = %s;" % expr
60 return [code]
61
63 return ["n[%d]" % i for i in range(d)]
65 expr = ", ".join("1.0" if i == 0 else "0.0"
66 for i in range(d))
67 code = "const double n[%d] = { %s };" % (d, expr)
68 return [code]
69
73 expr = "1.0"
74 code = "const double K_vol = %s;" % expr
75 return [code]
76
80 expr = "1.0"
81 code = "const double K_rad = %s;" % expr
82 return [code]
83
87 expr = "1.0"
88 code = "const double K_surf = %s;" % expr
89 return [code]
90
94 expr = "1.0"
95 code = "const double F_area = %s;" % expr
96 return [code]
97