コード例 #1
0
ファイル: almabench.c プロジェクト: ultimate-pa/ultimate
//---------------------------------------------------------------------------
// Entry point
// Calculate RA and Dec for noon on every day in 1900-2100
int main(int argc, char ** argv)
{
    int i, n, p;
    double jd[2];
    double pv[2][3];
    double position[8][3];
    bool   ga_testing = false;

    // do we have verbose output?
    if (argc > 1)
    {
        for (i = 1; i < argc; ++i)
        {
            if (!strcmp(argv[1],"-ga"))
            {
                ga_testing = true;
                break;
            }
        }
    }

    // get starting time

    // main loop
    for (i = 0; i < TEST_LOOPS; ++i)
    {
        jd[0] = J2000;
        jd[1] = 0.0;

        for (n = 0; n < TEST_LENGTH; ++n)
        {
            jd[0] += 1.0;

            for (p = 0; p < 8; ++p)
            {
                planetpv(jd,p,pv);
                radecdist(pv,position[p]);
            }
        }
    }

    for (p = 0; p < 8; ++p)
      printf("%f %f %f\n", position[p][0], position[p][1], position[p][2]);


    // get final time

    // report runtime

    fflush(stdout);

    return 0;
}
コード例 #2
0
ファイル: almabench.c プロジェクト: uwplse/peek
static void test(void)
{
    int p;
    double jd[2];
    double pv[2][3];
    double position[3];
    
    jd[0] = J2000;
    jd[1] = 0.0;
    for (p = 0; p < 8; ++p)
      {
        planetpv(jd,p,pv);
        radecdist(pv,position);
        printf("p = %d  position[0] = %g  position[1] = %g\n",
               p, position[0], position[1]);
      }
}
コード例 #3
0
ファイル: almabench.c プロジェクト: uwplse/peek
static void bench(int nloops)
{
    int i, n, p;
    double jd[2];
    double pv[2][3];
    double position[3];
    
    for (i = 0; i < nloops; ++i)
    {
        jd[0] = J2000;
        jd[1] = 0.0;

        for (n = 0; n < TEST_LENGTH; ++n)
        {
            jd[0] += 1.0;
            
            for (p = 0; p < 8; ++p)
            {
                planetpv(jd,p,pv);
                radecdist(pv,position);
            }
        }
    }
}