ErrMin.cpp

Go to the documentation of this file.

 
#include "SD.h"
 
 
namespace HepLib::SD {
 
IntegratorBase *ErrMin::Integrator = NULL;
MinimizeBase *ErrMin::miner = NULL;
qREAL *ErrMin::paras = NULL;
dREAL ErrMin::err_max;
dREAL ErrMin::err_min = -0.001;
size_t ErrMin::MaxRND = 50;
size_t ErrMin::RunRND = 0;
dREAL *ErrMin::lambda = NULL;
dREAL ErrMin::hjRHO = 0.75;
ex ErrMin::lastResErr;
 
dREAL ErrMin::IntError(int nvars, dREAL *las, dREAL *n1, dREAL *n2) {
    RunRND++;
    dREAL dlas[nvars];
    qREAL qlas[nvars];
    mpREAL mplas[nvars];
    for(int i=0; i<nvars; i++) {
        dlas[i] = las[i];
        qlas[i] = las[i];
        mplas[i] = las[i];
    }
    Integrator->qLambda = qlas;
    Integrator->dLambda = dlas;
    Integrator->mpLambda = mplas;
    auto res = Integrator->Integrate();
    if(res.has(NaN)) return 1.E100L;
 
    auto err = res.subs(VE(w0, w1)==w1);
    numeric nerr = numeric(1.E100);
    try {
        nerr = ex_to<numeric>(NN(abs(err)));
        if(nerr > numeric(1.E100)) nerr = numeric(1.E100);
        if(err_max > ex2q(nerr)) {
            auto diff = VESimplify(lastResErr - res);
            diff = diff.subs(VE(0,0)==0);
            exset ves;
            diff.find(VE(w0, w1), ves);
            for(auto ve : ves) {
                if(abs(ve.op(0)) > ve.op(1)) return 1.E100L;
            }
            if(Verbose > 3) {
                cout << "\r                             \r";
                cout << Color_HighLight << "     " << RunRND << ": " << RESET;
                for(int i=0; i<nvars; i++) cout << las[i] << " ";
                cout << endl << "     " << res.subs(VE(0,0)==0).subs(VE(w1,w2)==VEO(w1,w2)) << endl;
            }
            err_max = (dREAL)ex2q(nerr);
            for(int i=0; i<nvars; i++) lambda[i] = las[i];
            if(err_max<=err_min) {
                cout << "\r                             \r";
                cout << "     ------------------------------" << endl;
                miner->ForceStop();
                return 0.;
            }
        } else {
            if(Verbose > 3) {
                cout << "\r                             \r";
                cout << Color_HighLight << "     [ " << RunRND << " / " << MaxRND << " ] ..." << RESET << flush;
            }
        }
    } catch(domain_error &ex) {
        throw ex;
    } catch(...) { }
    if(RunRND>=MaxRND) {
        cout << "\r                             \r";
        cout << "     ------------------------------" << endl;
        miner->ForceStop();
        return 0.;
    }
    
    auto pid = getpid();
    ostringstream fn;
    fn << pid << ".las.done";
    if(file_exists(fn.str().c_str())) {
        ostringstream cmd;
        cmd << "rm " << fn.str();
        auto rc = system(cmd.str().c_str());
        cout << "\r                             \r";
        if(Verbose>3) cout << "     Exit: " << fn.str() << endl;
        cout << "     ------------------------------" << endl;
        miner->ForceStop();
        return 0.;
    }
    dREAL ret = (dREAL)ex2q(nerr);
    return ret;
}
 
 
}