EMM_OperatorsTest.hpp

Go to the documentation of this file.

#ifndef EMM_OperatorsTest_HPP
#define EMM_OperatorsTest_HPP

#include "EMM_Array.h"

template <class T>
tReal EMM_OperatorsTest::compareDiscretizedData(const tUSInt& d,const tString& fn,const EMM_MagneticExcitationOperator& op,tUIndex& indexMax) {
    EMM_Array<T> array;
    EMM_Array<T> arrayRef;
    tString dataName;
    tUIndex i,n;
    tDimension dim;
    const T * values;
    tBoolean ok=true;
    tReal errMax=0,err;
    if (!op.getDataField(d,dataName,n,dim,values)) return CORE_Object::getRealInfinity();
    array.setValues(n*dim,values);
    array.saveToFile(fn+".txt");
    indexMax=0;
    if (arrayRef.loadFromFile(fn+"_ref.txt")) {
        ok = ok && (array.getSize()==arrayRef.getSize());
        if (!ok) cout << "reference array & array "<<fn<<" have not same size.\n";
        if (ok) {
            errMax=0;
            for (i=0;i<array.getSize();i++) {
                err=(fabs(arrayRef[i])>1)?fabs(1.-array[i]/arrayRef[i]):fabs(array[i]-arrayRef[i]);
                if (errMax<err) {
                    errMax=err;
                    indexMax=i;
                    if (errMax>1.e-4) cout <<"Warning "<<array[indexMax]<<" != "<<arrayRef[indexMax]<<" err="<<errMax<<"\n";
                }
            }
            
        }
    }else {
        cout << "impossible to read "<<fn<<"_ref.txt"<<"\n";
        ok=false;   
    }
    
    if (!ok) errMax=CORE_Object::getRealInfinity();
    return errMax;
}

#endif