Package sfc :: Package geometry :: Module geometric_quantities
[hide private]
[frames] | no frames]

Source Code for Module sfc.geometry.geometric_quantities

 1   
2 -def gen_vx_names(d):
3 return ["vx[%d]" % i for i in range(d)]
4 -def gen_vx_code(d): # FIXME
5 num_vertices = d+1 # FIXME 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
15 -def gen_x0_names(d):
16 return ["x0[%d]" % i for i in range(d)]
17 -def gen_x0_code(d): # FIXME
18 expr = ", ".join("1.0" 19 for i in range(d)) 20 code = "const double x0[%d] = { %s };" % (d, expr) 21 return [code] 22
23 -def gen_G_names(d):
24 return ["G[%d]" % i for i in range(d*d)]
25 -def gen_G_code(d): # FIXME
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
32 -def gen_detG_names(d):
33 return ["detG"]
34 -def gen_detG_code(d):
35 expr = "1.0" # FIXME 36 code = "const double detG = %s;" % expr 37 return [code]
38
39 -def gen_Ginv_names(d):
40 return ["Ginv[%d]" % i for i in range(d*d)]
41 -def gen_Ginv_code(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
48 -def gen_detGsign_names(d):
49 return ["detGsign"]
50 -def gen_detGsign_code(d):
51 expr = "1.0" # FIXME 52 code = "const double detGsign = %s;" % expr 53 return [code]
54
55 -def gen_D_names(d):
56 return ["D"]
57 -def gen_D_code(d):
58 expr = "1.0" # FIXME 59 code = "const double D = %s;" % expr 60 return [code]
61
62 -def gen_n_names(d):
63 return ["n[%d]" % i for i in range(d)]
64 -def gen_n_code(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
70 -def gen_K_vol_names(d):
71 return ["K_vol"]
72 -def gen_K_vol_code(d):
73 expr = "1.0" # FIXME 74 code = "const double K_vol = %s;" % expr 75 return [code]
76
77 -def gen_K_rad_names(d):
78 return ["K_rad"]
79 -def gen_K_rad_code(d):
80 expr = "1.0" # FIXME 81 code = "const double K_rad = %s;" % expr 82 return [code]
83
84 -def gen_K_surf_names(d):
85 return ["K_surf"]
86 -def gen_K_surf_code(d):
87 expr = "1.0" # FIXME 88 code = "const double K_surf = %s;" % expr 89 return [code]
90
91 -def gen_F_area_names(d):
92 return ["F_area"]
93 -def gen_F_area_code(d):
94 expr = "1.0" # FIXME 95 code = "const double F_area = %s;" % expr 96 return [code]
97