reconstruct_ADMM.cpp File Reference

#include "reconstruct_ADMM.h"
#include "symmetrize.h"
#include <core/metadata_extension.h>

Include dependency graph for reconstruct_ADMM.cpp:

Go to the source code of this file.

Macros
#define	SYMMETRIZE_PROJECTIONS

Functions
void	addGradientTerm (double mu, AdmmKernel &kernel, MultidimArray< double > &L, FourierTransformer &transformer, MultidimArray< std::complex< double > > &fourierKernelV, char direction)

Macro Definition Documentation

SYMMETRIZE_PROJECTIONS

#define SYMMETRIZE_PROJECTIONS

Definition at line 37 of file reconstruct_ADMM.cpp.

Function Documentation

addGradientTerm()

void addGradientTerm	(	double	mu,
		AdmmKernel &	kernel,
		MultidimArray< double > &	L,
		FourierTransformer &	transformer,
		MultidimArray< std::complex< double > > &	fourierKernelV,
		char	direction
	)

Definition at line 444 of file reconstruct_ADMM.cpp.

{
    kernel.computeGradient(L,direction);

    transformer.FourierTransform();
    double K=mu*MULTIDIM_SIZE(L);
    auto xdim_2=(double)(XSIZE(L)/2);
    double Kargument=2*PI*xdim_2/XSIZE(L);
    FOR_ALL_ELEMENTS_IN_ARRAY3D(fourierKernelV)
    {
        // Calculation of L^t*L
        // -L^2(r)*exp(-i*2*pi*(N/2)/N*(k+i+j)): the last term is a phase shift due to FFT
        double argument=Kargument*(k+i+j);
        double s, c;
        //sincos(argument,&s,&c);
        s = sin(argument);
        c = cos(argument);
        A3D_ELEM(fourierKernelV,k,i,j)+=-A3D_ELEM(transformer.fFourier,k,i,j)*
                A3D_ELEM(transformer.fFourier,k,i,j)*std::complex<double>(K*c,K*s);
    }
}