public static double RadiusVector(double JD, bool bHighPrecision) { if (bHighPrecision) { return(AASVSOP87D_Mars.R(JD)); } double rho = (JD - 2451545) / 365250; double rhosquared = rho * rho; double rhocubed = rhosquared * rho; double rho4 = rhocubed * rho; //Calculate R0 int nR0Coefficients = g_R0MarsCoefficients.Length; double R0 = 0; int i; for (i = 0; i < nR0Coefficients; i++) { R0 += g_R0MarsCoefficients[i].A * Math.Cos(g_R0MarsCoefficients[i].B + g_R0MarsCoefficients[i].C * rho); } //Calculate R1 int nR1Coefficients = g_R1MarsCoefficients.Length; double R1 = 0; for (i = 0; i < nR1Coefficients; i++) { R1 += g_R1MarsCoefficients[i].A * Math.Cos(g_R1MarsCoefficients[i].B + g_R1MarsCoefficients[i].C * rho); } //Calculate R2 int nR2Coefficients = g_R2MarsCoefficients.Length; double R2 = 0; for (i = 0; i < nR2Coefficients; i++) { R2 += g_R2MarsCoefficients[i].A * Math.Cos(g_R2MarsCoefficients[i].B + g_R2MarsCoefficients[i].C * rho); } //Calculate R3 int nR3Coefficients = g_R3MarsCoefficients.Length; double R3 = 0; for (i = 0; i < nR3Coefficients; i++) { R3 += g_R3MarsCoefficients[i].A * Math.Cos(g_R3MarsCoefficients[i].B + g_R3MarsCoefficients[i].C * rho); } //Calculate R4 int nR4Coefficients = g_R4MarsCoefficients.Length; double R4 = 0; for (i = 0; i < nR4Coefficients; i++) { R4 += g_R4MarsCoefficients[i].A * Math.Cos(g_R4MarsCoefficients[i].B + g_R4MarsCoefficients[i].C * rho); } return((R0 + R1 * rho + R2 * rhosquared + R3 * rhocubed + R4 * rho4) / 100000000); }
public void RTest(double jd, double expectedResult) { double result = AASVSOP87D_Mars.R(jd); Assert.Equal(expectedResult, result); }