| Revision: | 3640 |
| Committed: | Thu Aug 12 20:25:00 2010 UTC (15 years, 1 month ago) by gezelter |
| File size: | 389 byte(s) |
| Log Message: | new import of Langevin Hull paper |
| # | Content |
|---|---|
| 1 | 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] |