Revision: | 3640 |
Committed: | Thu Aug 12 20:25:00 2010 UTC (14 years ago) by gezelter |
File size: | 389 byte(s) |
Log Message: | new import of Langevin Hull paper |
# | User | Rev | Content |
---|---|---|---|
1 | gezelter | 3640 | vert0 = {x0, y0, z0} |
2 | vert1 = {x1, y1, z1} | ||
3 | vert2 = {x2, y2, z2} | ||
4 | |||
5 | centroid = (vert0 + vert1 + vert2) / 3 | ||
6 | |||
7 | a = vert0 - vert1 | ||
8 | b = vert0 - vert2 | ||
9 | c = vert1 - vert2 | ||
10 | |||
11 | s3 = Sqrt[Cross[a,b] . Cross[a,b]] / 2 | ||
12 | dr3 = centroid - vert0 | ||
13 | l3 = 1 / (dr3.dr3) | ||
14 | |||
15 | G3 = (IdentityMatrix[3] + Outer[Times, dr3, dr3]*l3)*Sqrt[l3]*s3*8/(Pi*eta) | ||
16 | |||
17 | Hmat = G3 | ||
18 | |||
19 | |||
20 | |||
21 | |||
22 | normal = Cross[a,b] | ||
23 | |||
24 | nhat = normal / Sqrt[normal.normal] |