by Jessica Ramella-Roman

In this page I will post any upgrade to the programs

  • 6-23-2005, "nan" error

    There was an error for the remote case of theta=nangles, at initializatio the terms s11, s12, s33, and s34 are initialized so that

    for (i=0; i < nangles; ++i){
        s11[i]=....
    

    this was replaced with

    for(i=0;  i <= nangles; ++i) {
        s11[i]=....
    

    This solves the problem of nan results. All the programs were updated to solve this issue

  • 7-14-2005, roundoff error

    Round off error in mcMeridian program, in the previous version I had taken out two lines

    if (cosi > 1) cosi=1;
    if (cosi < -1) cosi=-1; 
    

    immediately after rejection method. Taking those line out gave wrong nan results, lines were re-added to programs