CMTDataComp.cpp

Go to the documentation of this file.

#include "CMTDataComp.h"

#ifndef _STDLIB_H_INCLUDED_ 
#define _STDLIB_H_INCLUDED_
#include "stdlib.h"
#endif

#include <cmath>
#include <algorithm>

CMTDataComp::CMTDataComp(const int nfreq):
name("unknown"),
rhoa(nfreq),
phi(nfreq),
drhoa(nfreq),
dphi(nfreq),
Z(nfreq),
dZ(nfreq) 
{ 
        
} 

CMTDataComp::~CMTDataComp()
{
}

void CMTDataComp::Assign(const int nfreq)
{
        rhoa.assign(nfreq,0);
        phi.assign(nfreq,0);
        drhoa.assign(nfreq,0);
        dphi.assign(nfreq,0);
        Z.assign(nfreq,0);
        dZ.assign(nfreq,0);
}

CMTDataComp& CMTDataComp::operator= (const CMTDataComp& source)
{
        this->rhoa.assign(source.Z.size(),0);
        this->phi.assign(source.Z.size(),0);
        this->dphi.assign(source.Z.size(),0);
        this->drhoa.assign(source.Z.size(),0);
        this->Z.assign(source.Z.size(),0);
        this->dZ.assign(source.Z.size(),0);
        this->name = source.name;
        copy(source.rhoa.begin(),source.rhoa.end(),this->rhoa.begin());
        copy(source.phi.begin(),source.phi.end(),this->phi.begin());
        copy(source.drhoa.begin(),source.drhoa.end(),this->drhoa.begin());
        copy(source.dphi.begin(),source.dphi.end(),this->dphi.begin());
        copy(source.Z.begin(),source.Z.end(),this->Z.begin());
        copy(source.dZ.begin(),source.dZ.end(),this->dZ.begin());
        return *this;
}

/*! Given a vector with frequencies
 * that correspond to the impedance values, Update calculates apparent
 * resistivity and phase
 */
void CMTDataComp::Update(trealdata frequency)
{
   //! If we calculate for the first time, the derived vectors have size 0
   if (rhoa.empty())
   {
                //! So we have to allocate memory and set them to 0
                rhoa.assign(frequency.size(),0);
                phi.assign(frequency.size(),0);
                drhoa.assign(frequency.size(),0);
                dphi.assign(frequency.size(),0);
   }
   //! Then we can calculate the values for each frequency
   for (unsigned int k = 0; k < frequency.size(); ++k)
   {
                rhoa.at(k) = 1./(2 * PI * frequency.at(k)) * mu * (pow(Z.at(k).real(),2)+pow(Z.at(k).imag(),2)) * pow(1000.0,2);
                phi.at(k)  = atan(Z.at(k).imag()/Z.at(k).real())/ PI *180;
                drhoa.at(k) = 1./(PI * frequency.at(k)) * mu * abs(Z.at(k)) * dZ.at(k)* pow(1000.0,2);
                dphi.at(k) = dZ.at(k) * abs(Z.at(k)) * 1./(1.+ pow(Z.at(k).imag()/Z.at(k).real(),2)) * 1./pow(Z.at(k).real(),2)/ PI *180;  
   }
}

/*! If we want to erase one entry, we have to erase all the corresponding derived quantities as well*/
void CMTDataComp::erase(const int index)
{
        rhoa.erase(rhoa.begin()+index);
        phi.erase(phi.begin()+index);
        drhoa.erase(drhoa.begin()+index);
        dphi.erase(dphi.begin()+index);
        Z.erase(Z.begin()+index);
        dZ.erase(dZ.begin()+index);
}

---