- Timestamp:
- Oct 12, 2015, 4:59:00 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/ephemeris.cpp
r7487 r7501 1481 1481 double doty = dotr*sin(u) + r*cos(u)*dotu; 1482 1482 1483 1484 1483 const double iMaxGEO = 10.0 / 180.0 * M_PI; 1485 1484 … … 1550 1549 - yp*sini*cosom*doti; 1551 1550 1552 //double vz = sini *doty + yp*cosi *doti; 1553 1554 ColumnVector V(3); 1555 V[0] = cosom *vx - cosi*sinom *vy // dX / dr 1556 - xp*sinom*dotom - yp*cosi*cosom*dotom // dX / dOMEGA 1557 + yp*sini*sinom*doti; // dX / di 1558 1559 V[1] = sinom *vx + cosi*cosom *vy 1560 + xp*cosom*dotom - yp*cosi*sinom*dotom 1561 - yp*sini*cosom*doti; 1562 1563 V[2] = sini *vy + yp*cosi *doti; 1551 double vz = sini *doty + yp*cosi *doti; 1552 1553 ColumnVector V(3); V << vx << vy << vz; 1564 1554 1565 1555 Matrix RdotZ(3,3);
Note:
See TracChangeset
for help on using the changeset viewer.