| Revision: | 3640 |
| Committed: | Thu Aug 12 20:25:00 2010 UTC (13 years, 11 months 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] |