solvetteq.cpp

Go to the documentation of this file.

#include <iostream>
#include <cmath>
using namespace std;
int main()
{
//      const double tps = 4.5;
//      const double tppss = 20.9;
//      const double tppps = 15.8;
//      const double p = 0.064;
        const double tps = 4.5;
        const double tppss = 19.6;
        const double tppps = 15.1;
        const double p = 0.07;
        const double psq = p * p;
        double h,vs,vp;
        const double starth = 20;
        const double endh = 45;
        const double deltah = 0.1;
        const double startvs = 3;
        const double endvs = 4;
        const double deltavs = 0.01;
        const double startvp = startvs * sqrt(3);
        const double endvp = endvs * sqrt(3);
        const double deltavp = 0.01;
        const double pickerror = 0.1;
        double deltatps, deltatppss, deltatppps;
        double besth,bestvp, bestvs;
        double bestdelta = 1e16;
        double currdelta;
        for (double currvs = startvs; currvs < endvs; currvs+=deltavs)
        {
                const double vsterm = sqrt(1./(currvs*currvs) - psq);
                //const double currvp = sqrt(3) * currvs;
                for (double currvp = startvp; currvp < endvp; currvp+=deltavp)
                {
                        const double vpterm = sqrt(1./(currvp*currvp) -psq);
                        for (double currh = starth; currh < endh; currh+=deltah)
                        {
                                deltatps = tps - currh *(vsterm - vpterm);
                                deltatppps = tppps - currh *(vsterm + vpterm);
                                deltatppss = tppss - 2*currh * vsterm;
                                currdelta = deltatps * deltatps + deltatppps* deltatppps + deltatppss*deltatppss;
                                if (currdelta < bestdelta)
                                {
                                        besth = currh;
                                        bestvp = currvp;
                                        bestvs = currvs;
                                        bestdelta = currdelta;
                                }
                        }
                }
        }
        cout << "H: " << besth << " Vp: " << bestvp << " Vs: " <<  bestvs << " Delta: " << sqrt(bestdelta/3.) << " Vp/Vs: " << bestvp/bestvs << endl; 
}

---