GreatCircle_Xn.h

Go to the documentation of this file.

//- ****************************************************************************
//-
//- Copyright 2009 National Technology & Engineering Solutions of Sandia, LLC
//- (NTESS). Under the terms of Contract DE-NA0003525 with NTESS, the U.S.
//- Government retains certain rights in this software.
//-
//- BSD Open Source License
//- All rights reserved.
//-
//- Redistribution and use in source and binary forms, with or without
//- modification, are permitted provided that the following conditions are met:
//-
//-   1. Redistributions of source code must retain the above copyright notice,
//-      this list of conditions and the following disclaimer.
//-
//-   2. Redistributions in binary form must reproduce the above copyright
//-      notice, this list of conditions and the following disclaimer in the
//-      documentation and/or other materials provided with the distribution.
//-
//-   3. Neither the name of the copyright holder nor the names of its
//-      contributors may be used to endorse or promote products derived from
//-      this software without specific prior written permission.
//-
//- THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
//- AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
//- IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
//- ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
//- LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
//- CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
//- SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
//- INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
//- CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
//- ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
//- POSSIBILITY OF SUCH DAMAGE.
//-
//- ****************************************************************************
 
#ifndef GreatCircle_Xn_H
#define GreatCircle_Xn_H
 
// **** _SYSTEM INCLUDES_ ******************************************************
 
#include <vector>
#include <string>
#include <iostream>
 
using namespace std;
 
// **** _LOCAL INCLUDES_ *******************************************************
 
#include "SLBMGlobals.h"
#include "GreatCircle.h"
#include "Grid.h"
#include "Location.h"
#include "CrustalProfile.h"
#include "LayerProfile.h"
 
// **** _BEGIN SLBM NAMESPACE_ **************************************************
 
namespace slbm {
 
// **** _CLASS DEFINITION_ *****************************************************
//
class SLBM_EXP_IMP GreatCircle_Xn : public GreatCircle
{
 
public:
 
    GreatCircle_Xn(
        const int& _phase,
        Grid& _grid,
        const double& latSource,
        const double& lonSource,
        const double& depthSource,
        const double& latReceiver,
        const double& lonReceiver,
        const double& depthReceiver,
        const double& chMax);
 
    ~GreatCircle_Xn();
 
    GreatCircle_Xn(const GreatCircle_Xn &other);
 
    GreatCircle_Xn& operator=(const GreatCircle_Xn& other);
 
    void getZhaoParameters(double& Vm, double& Gm, double& H, double& C, double& Cm, int& udSign);
 
    string toString(const int& verbosity);
 
private:
 
    double ch_max;
 
    double Vm;
 
    double Gm;
 
    double H;
 
    double rMoho;
    double rZm;
 
    int udSign;
 
    double V0;
 
    void computeTravelTime();
 
    bool inCrust;
 
    void computeTravelTimeCrust();
 
    void computeTravelTimeMantle();
 
    double c, cm, zm, zhao_r, cz, cmz;
 
    const double cmin;
 
    double func(const double& h);
 
    void mnbrak(double &ax, double &bx, double &cx,
                double &fa, double &fb, double &fc);
 
    double brent(const double ax, const double bx, const double cx,
        const double tol, double &xmin);
 
    void SWAP(double &a, double &b) {double dum=a; a=b; b=dum;}
 
    double SIGN(const double &a, const double &b)
        {return b >= 0 ? (a >= 0 ? a : -a) : (a >= 0 ? -a : a);};
 
    void shft3(double &a, double &b, double &cc, const double d)
        { a=b; b=cc; cc=d;};
 
};
 
inline void GreatCircle_Xn::getZhaoParameters(double& Vm_, double& Gm_, double& H_, double& C_, double& Cm_, int& udSign_)
{
    Vm_ = Vm;
    Gm_ = Gm;
    H_ = H;
    C_ = c;
    Cm_ = cm;
    udSign_ = udSign;
}
 
} // end slbm namespace
 
#endif // GreatCircle_Xn.h