volume_halves_restoration.cpp File Reference

#include "volume_halves_restoration.h"
#include <data/numerical_tools.h>

Include dependency graph for volume_halves_restoration.cpp:

Go to the source code of this file.

Functions
double	restorationSigmaCost (double *x, void *_prm)

Function Documentation

restorationSigmaCost()

double restorationSigmaCost	(	double *	x,
		void *	_prm
	)

Definition at line 292 of file volume_halves_restoration.cpp.

{
    auto *prm=(ProgVolumeHalvesRestoration *) _prm;
    double sigma1=x[1];
    double sigma2=x[2];
    if (sigma1<0 || sigma2<0 || sigma1>2 || sigma2>2)
        return 1e38;
    double K1=-0.5/(sigma1*sigma1);
    double K2=-0.5/(sigma2*sigma2);
    double error=0;
    double N=0;
    const MultidimArray< std::complex<double> > &fV=prm->fVol;
    const MultidimArray< std::complex<double> > &fV1r=prm->fV1r;
    const MultidimArray< std::complex<double> > &fV2r=prm->fV2r;
    const MultidimArray<double> &R2=prm->R2;
    FOR_ALL_DIRECT_ELEMENTS_IN_MULTIDIMARRAY(fV)
    {
        double R2n=DIRECT_MULTIDIM_ELEM(R2,n);
        if (R2n<=0.25)
        {
            double H1=exp(K1*R2n);
            double H2=exp(K2*R2n);
            error+=abs(DIRECT_MULTIDIM_ELEM(fV,n)*H1-DIRECT_MULTIDIM_ELEM(fV1r,n))+
                   abs(DIRECT_MULTIDIM_ELEM(fV,n)*H2-DIRECT_MULTIDIM_ELEM(fV2r,n));
            N++;
        }
    }
    return error;
}