ParametricGeoMagChange.cc

Go to the documentation of this file.

#include "ParametricGeoMagChange.h"

namespace fwk {

  ParametricGeoMagChange::ParametricGeoMagChange(size_t order, double g[], double h[]) :
    fOrder(order), fG(std::valarray<double> (g,order*(order+1)/2)), fH(std::valarray<double> (h, order*(order+1) / 2)) 
  {}

  
  ParametricGeoMagneticField 
  ParametricGeoMagChange::Extrapolate (double years, 
                                       const ParametricGeoMagneticField &baseModel)
  {
    // new_coefficients = old_coefficients + dt * change_coefficients
    //
    // If change and base model differ in size, resize accordingly.
    // That should be a corner case, but its implementation is trivial.

    ParametricGeoMagneticField  model = baseModel; //copy constructor
    
    // we are friends ...
    if (fOrder > model.fOrder)
      {
        model.Grow(fOrder);
      }

    if (model.fOrder == fOrder)
      {
        model.fG += years * fG;
        model.fH += years * fH;
      }
    else  // unlikely and maybe senseless case
      {
        size_t order = model.fOrder;
        std::valarray<double> g = fG;
        std::valarray<double> h = fH;
        g.resize(order*(order+1) / 2, 0.0);
        h.resize(order*(order+1) / 2, 0.0);

        model.fG += years * g;
        model.fH += years * h;
      }

    return model;    
  }


  /* Meant for testing 
     Parameters should rather be read from file/db.
  */
  ParametricGeoMagChange
  ParametricGeoMagChange::SV2010_2015 ()
  {
    double g[] ={ 
      11.40,
      16.70,
      -11.30,
      -3.90,
      2.70,
      1.30,
      -3.90,
      -2.90,
      -8.10,
      -1.40,
      2.00,
      -8.90,
      4.40,
      -2.30,
      -0.50,
      0.50,
      -1.50,
      -0.70,
      1.30,
      1.40,
      -0.30,
      -0.30,
      -0.30,
      1.90,
      -1.60,
      -0.20,
      1.80,
      0.20,
      -0.10,
      -0.60,
      1.40,
      0.30,
      0.10,
      -0.80,
      0.40,
      -0.10,
      0.10,
      -0.50,
      0.30,
      -0.30,
      0.30,
      0.20,
      -0.50,
      0.20,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00
    };

    double h[] = {
      0.00,
      -28.80,
      0.00,
      -23.00,
      -12.90,
      0.00,
      8.60,
      -2.90,
      -2.10,
      0.00,
      0.40,
      3.20,
      3.60,
      -0.80,
      0.00,
      0.50,
      1.50,
      0.90,
      3.70,
      -0.60,
      0.00,
      -0.10,
      -2.10,
      -0.40,
      -0.50,
      0.80,
      0.50,
      0.00,
      0.60,
      0.30,
      -0.20,
      -0.10,
      -0.80,
      -0.30,
      0.20,
      0.00,
      0.00,
      0.20,
      0.50,
      0.40,
      0.10,
      -0.10,
      0.40,
      0.40,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00,
      0.00
    };

    return ParametricGeoMagChange (14, g, h);
  }

 
}