interpolate_CF2()

template<typename T >

void cntr::interpolate_CF2	(	int	tstp,
		cntr::function< T > &	H,
		cdmatrix &	Hinter,
		int	kt,
		bool	fixHam = false
	)

Interpolation for the second order propagator

Purpose

Interpolates H(t-dt/2) from H(t+dt), H(t), H(t-dt)

where H(t+dt) is obtained from the extrapolation if not yet known H is original function If fixham=true we assume that the hamiltonian is known for all times and there is no extrapolation

Parameters

tstp	timestep at which propagator is calculated
H	Time dependent Hamiltonian
Hinter	Interpolated Hamiltonian
kt	Order of integrator used for extrapolation and interpolation
fixHam	Hamiltonian is known for all times and no extrapolation is needed for the predictor/corrector

Definition at line 648 of file cntr_equilibrium_impl.hpp.

References cntr::function< T >::get_value(), integration::I< double >(), and cntr::function< T >::size1_.

Referenced by propagator_exp().

                                                                                           {
    int size=H.size1_;
    cdmatrix tmp1(size,size),tmp2(size,size);
    int ktextrap=(tstp<=kt ? tstp : kt);
    int n1=(tstp<=kt ? kt : tstp); //Once you are at kt the next point needs to be included
    // Extrapolate to get H(t + dt)
    if(!fixHam && tstp>kt){
        extrapolate_timestep(tstp-1,H,integration::I<double>(ktextrap));
    }
    // Interpolate to get H(t + dt/2)
    H.get_value(tstp-1,tmp1);
    H.get_value(tstp,tmp2);

    Hinter=interpolation(n1,(tstp-0.5),H,integration::I<double>(kt)); //ktextrap+1 since we added one more point
}

Here is the call graph for this function:

Here is the caller graph for this function: