Changeset 19504
- Timestamp:
- Sep 11, 2008, 2:59:00 PM (18 years ago)
- File:
-
- 1 edited
-
trunk/psLib/src/astro/psCoord.c (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/psLib/src/astro/psCoord.c
r15254 r19504 10 10 * @author GLG, MHPCC 11 11 * 12 * @version $Revision: 1.14 1$ $Name: not supported by cvs2svn $13 * @date $Date: 200 7-10-09 19:25:44$12 * @version $Revision: 1.142 $ $Name: not supported by cvs2svn $ 13 * @date $Date: 2008-09-12 00:59:00 $ 14 14 * 15 15 * Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii … … 87 87 88 88 psPlaneTransform *out = psPlaneTransformAlloc(1, 1); 89 90 if (0) {91 // XXX: Get rid of this code.92 psF64 A = transform->x->coeff[1][0];93 psF64 B = transform->x->coeff[0][1];94 psF64 C = transform->y->coeff[1][0];95 psF64 D = transform->y->coeff[0][1];96 psF64 E = transform->x->coeff[0][0];97 psF64 F = transform->y->coeff[0][0];98 psF64 invDet = 1.0 / (A * D - B * C);99 100 out->x->coeff[1][0] = invDet * D;101 out->x->coeff[0][1] = invDet * -B;102 out->y->coeff[1][0] = invDet * -C;103 out->y->coeff[0][1] = invDet * A;104 105 //106 // This comes from substituting the (x_2 = x_1 - e) and (y_2 = y_1 - f) equations:107 //108 out->x->coeff[0][0] = invDet * (B * F - D * E);109 out->y->coeff[0][0] = invDet * (E * C - A * F);110 111 }112 89 113 90 psF64 r12 = 0.0;
Note:
See TracChangeset
for help on using the changeset viewer.
