testcosine.cc File Reference

#include <madness/mra/mra.h>
#include <madness/tensor/solvers.h>

Include dependency graph for testcosine.cc:

Classes
class	CosPotentialFunctor
class	ExpFunctor< Q >
class	ExpFunctor3d< Q >
class	PWFunctor

Typedefs
typedef SeparatedConvolution< double_complex, 3 >	coperatorT
typedef SeparatedConvolution< double, 3 >	operatorT
typedef std::shared_ptr< operatorT>	poperatorT

Functions
complex_function_3d	apply_periodic_bsh (World &world, const complex_function_3d &f, const double &kx, const double &ky, const double &kz, const double &energy, const double &L)
vector_complex_function_3d	apply_potential (const real_function_3d &potential, const vector_complex_function_3d &psi)
std::vector< coord_3d >	get_coeffs_pw (double maxKlen)
static const double_complex	I (0, 1)
int	main (int argc, char **argv)
vector_complex_function_3d	make_basis (World &world, double maxKlen)
tensor_complex	make_kinetic_matrix (World &world, const vector_complex_function_3d &v)
real_tensor	make_pw_matrix (const std::vector< coord_3d > &gvecs)
void	orthogonalize (World &world, vector_complex_function_3d &psi)
vector_complex_function_3d	update (World &world, const vector_complex_function_3d &psi, vector_complex_function_3d &vpsi, const tensor_real &e, int iter)

Variables
static const double	alpha = 2.5
static const double	kwavelet = 4
static const double	kx =0.5*twopi/L
static const double	ky =0.5*twopi/L
static const double	kz =0.5*twopi/L
static const double	L = 5.0
static const double	pi = constants::pi
static const double	thresh = 1e-2
static const int	truncate_mode = 0
static const double	twopi = 2.0*constants::pi

Typedef Documentation

coperatorT

typedef SeparatedConvolution<double_complex,3> coperatorT

operatorT

typedef SeparatedConvolution<double,3> operatorT

poperatorT

typedef std::shared_ptr<operatorT> poperatorT

Function Documentation

apply_periodic_bsh()

complex_function_3d apply_periodic_bsh	(	World &	world,
		const complex_function_3d &	f,
		const double &	kx,
		const double &	ky,
		const double &	kz,
		const double &	energy,
		const double &	L )

References madness::apply(), e(), energy, f, madness::FunctionFactory< T, NDIM >::functor(), madness::g, madness::FunctionDefaults< NDIM >::get_thresh(), I, kx, ky, kz, L, op(), madness::PeriodicBSHOperator3D(), madness::FunctionFactory< T, NDIM >::truncate_mode(), and madness::FunctionFactory< T, NDIM >::truncate_on_project().

Referenced by update().

apply_potential()

vector_complex_function_3d apply_potential	(	const real_function_3d &	potential,
		const vector_complex_function_3d &	psi )

References potential(), and psi().

Referenced by main().

get_coeffs_pw()

std::vector< coord_3d > get_coeffs_pw

(

double

maxKlen

)

References a, b, L, and pi.

I()

const double_complex I ( 0 , 1  )

static

main()

int main	(	int	argc,
		char **	argv )

References madness::_(), std::abs(), apply_potential(), madness::BC_PERIODIC, c, SafeMPI::COMM_WORLD, madness::copy(), e(), madness::FunctionFactory< T, NDIM >::functor(), madness::initialize(), kwavelet, L, make_basis(), make_kinetic_matrix(), madness::matrix_inner(), madness::print(), madness::project(), psi(), madness::World::rank(), madness::FunctionDefaults< NDIM >::set_bc(), madness::FunctionDefaults< NDIM >::set_cubic_cell(), madness::FunctionDefaults< NDIM >::set_k(), madness::FunctionDefaults< NDIM >::set_thresh(), madness::FunctionDefaults< NDIM >::set_truncate_mode(), madness::startup(), madness::sygv(), thresh, madness::transform(), madness::FunctionFactory< T, NDIM >::truncate_mode(), truncate_mode, madness::FunctionFactory< T, NDIM >::truncate_on_project(), update(), and v.

make_basis()

vector_complex_function_3d make_basis	(	World &	world,
		double	maxKlen )

References madness::FunctionFactory< T, NDIM >::functor(), I, L, madness::normalize(), pi, madness::FunctionFactory< T, NDIM >::truncate_on_project(), and twopi.

Referenced by main().

make_kinetic_matrix()

tensor_complex make_kinetic_matrix	(	World &	world,
		const vector_complex_function_3d &	v )

References madness::apply(), f1, f2, f3, I, kx, ky, kz, madness::matrix_inner(), and v.

Referenced by main().

make_pw_matrix()

real_tensor make_pw_matrix

(

const std::vector< coord_3d > &

gvecs

)

References madness::abs(), alpha, e(), L, madness::Vector< T, N >::normf(), pi, and madness::BaseTensor::size().

orthogonalize()

void orthogonalize	(	World &	world,
		vector_complex_function_3d &	psi )

References madness::compress(), madness::Function< T, NDIM >::gaxpy(), madness::inner(), madness::Function< T, NDIM >::norm2(), psi(), and madness::Function< T, NDIM >::scale().

Referenced by update().

update()

vector_complex_function_3d update	(	World &	world,
		const vector_complex_function_3d &	psi,
		vector_complex_function_3d &	vpsi,
		const tensor_real &	e,
		int	iter )

References apply_periodic_bsh(), e(), madness::gaxpy(), k, kx, ky, kz, L, madness::norm2(), madness::normalize(), orthogonalize(), madness::print(), psi(), madness::World::rank(), madness::scale(), shift, and madness::truncate().

Referenced by madness::Exchange< T, NDIM >::ExchangeImpl< T, NDIM >::MacroTaskExchangeSimple::compute_diagonal_batch_in_symmetric_matrix(), main(), and madness::Nemo::solve().

Variable Documentation

alpha

const double alpha = 2.5

static

Referenced by GaussianGuess< T, NDIM >::GaussianGuess(), GaussianGuess< T, NDIM >::GaussianGuess(), refpotfunctor::refpotfunctor(), madness::FunctionImpl< T, LDIM >::accumulate_trees(), madness::CCPotentials::apply_KffK(), madness::population::atomic_populations(), madness::cblas::axpy(), madness::cblas::axpy(), madness::cblas::axpy(), madness::cblas::axpy(), bind_function_type(), madness::DiagnosticMatrix< T, NDIM >::build_Gab_bra(), cggev_(), compute_madelung_energy(), compute_madelung_energy_PWSCF(), ctrsm_(), madness::Solver< T, NDIM >::do_rhs(), madness::Solver< T, NDIM >::do_rhs_simple(), dtrsm_(), exchange_anchor_test(), madness::exchange_anchor_test(), function(), madness::GFit< T, NDIM >::gaussian_spherical_moments(), madness::Function< double, 3 >::gaxpy(), madness::FunctionImpl< T, LDIM >::gaxpy(), madness::gaxpy(), madness::gaxpy(), madness::gaxpy(), madness::GenTensor< T >::gaxpy(), madness::GenTensor< T >::gaxpy(), madness::SVDTensor< T >::gaxpy(), madness::Tensor< double >::gaxpy(), madness::TensorTrain< T >::gaxpy(), gaxpy1(), madness::Function< double, 3 >::gaxpy_ext(), madness::FunctionImpl< T, LDIM >::gaxpy_ext(), madness::FunctionImpl< T, LDIM >::gaxpy_ext_node(), madness::FunctionImpl< T, LDIM >::gaxpy_ext_recursive(), madness::FunctionImpl< T, LDIM >::gaxpy_inplace(), madness::FunctionNode< T, NDIM >::gaxpy_inplace(), madness::FunctionImpl< T, LDIM >::gaxpy_inplace_reconstructed(), madness::Function< double, 3 >::gaxpy_oop(), madness::gaxpy_oop(), madness::gaxpy_oop(), madness::gaxpy_oop(), madness::FunctionImpl< T, NDIM >::gaxpy_oop_reconstructed(), madness::FunctionImpl< T, LDIM >::gaxpy_oop_reconstructed(), madness::gaxpy_oop_reconstructed(), madness::cblas::gemm(), madness::cblas::gemm(), madness::cblas::gemm(), madness::cblas::gemm(), madness::cblas::gemv(), madness::cblas::gemv(), madness::cblas::gemv(), madness::cblas::gemv(), generalized_laguerre(), madness::cblas::ger(), madness::cblas::ger(), madness::cblas::ger(), madness::cblas::ger(), gt(), gt1(), guess(), madness::SRConf< T >::inplace_add(), madness::GFit< T, NDIM >::InverseRSqFit(), main(), main(), Ansatz0::make_guess(), make_operands(), make_pw_matrix(), make_sp(), madness::mul(), madness::mul(), CosPotentialFunctor::operator()(), madness::operator*(), madness::operator*(), madness::cblas::scal(), madness::cblas::scal(), madness::cblas::scal(), madness::cblas::scal(), madness::cblas::scal(), madness::cblas::scal(), slater(), madness::Solver< T, NDIM >::solve(), strsm_(), test_coulomb(), madness::test_exchange(), test_exchange(), test_inv_rsq_gfit(), test_XCOperator(), uinitial(), uinitial(), V_cosine1D(), V_cosine3D(), V_cosine3D_sep(), zggev_(), and ztrsm_().

kwavelet

const double kwavelet = 4

static

Referenced by MiniDFT::MiniDFT(), MiniDFT::doit(), and main().

kx

const double kx =0.5*twopi/L

static

Referenced by madness::CCPotentials::apply_Kfxy(), apply_periodic_bsh(), main(), madness::BasisFunctions::make_cartesian_guess(), madness::Solver< T, NDIM >::make_kinetic_matrix(), make_kinetic_matrix(), make_kinetic_matrix(), makeao_slow(), update(), and update().

ky

const double ky =0.5*twopi/L

static

Referenced by apply_periodic_bsh(), main(), madness::BasisFunctions::make_cartesian_guess(), madness::Solver< T, NDIM >::make_kinetic_matrix(), make_kinetic_matrix(), make_kinetic_matrix(), makeao_slow(), update(), and update().

kz

const double kz =0.5*twopi/L

static

Referenced by apply_periodic_bsh(), main(), madness::BasisFunctions::make_cartesian_guess(), madness::Solver< T, NDIM >::make_kinetic_matrix(), make_kinetic_matrix(), make_kinetic_matrix(), makeao_slow(), update(), and update().

L

const double L = 5.0

static

Referenced by PWFunctor::PWFunctor(), and CosPotentialFunctor::operator()().

pi

const double pi = constants::pi

static

Referenced by DSphere::DSphere(), Guess::Guess(), MolecularNuclearPotentialFunctor2::MolecularNuclearPotentialFunctor2(), PWFunctor::PWFunctor(), Sphere::Sphere(), madness::HartreeFock< T, NDIM >::calculate_tot_coulomb_energy(), madness::HartreeFock< T, NDIM >::calculate_tot_exchange_energy(), madness::DFT< T, NDIM >::calculate_tot_ke_sp(), madness::HartreeFock< T, NDIM >::calculate_tot_ke_sp(), madness::HartreeFock< T, NDIM >::calculate_tot_pe_sp(), madness::Nemo::compute_IR_intensities(), erfaroverr(), madness::detail::adqtest::exact(), exact_cosine(), exact_spheropole(), f_cosine(), f_dipole(), f_quadrupole(), f_spheropole(), madness::gauss_legendre_numeric(), get_coeffs_pw(), madness::DistributedMatrixDistribution::get_range(), madness::GradBSHOperator(), GradBSHOperator_Joel(), madness::GradCoulombOperator(), gradPbarA(), madness::GradSlaterOperator(), madness::SeparatedConvolution< double, 3 >::initialize_images_only(), main(), make_basis(), make_fit(), make_grad_operator(), make_pw_matrix(), madness::EigSolverOp< T, NDIM >::multi_op_o(), madness::EigSolverOp< T, NDIM >::multi_op_r(), madness::EigSolver< T, NDIM >::multi_solve(), CosPotentialFunctor::operator()(), opposing_field_potential(), CFFT::Perform(), slymer::PrimitiveGaussian(), madness::EigSolver< T, NDIM >::solve(), test(), test_coulomb_gfit_discard_large_alpha(), test_in(), test_out(), madness::test_rnlij_rangelimited(), madness::test_rnlp(), madness::test_rnlp_rangelimited(), madness::test_rnlp_rangelimited_erf(), and testNavierStokes().

thresh

const double thresh = 1e-2

static

truncate_mode

const int truncate_mode = 0

static

Referenced by main().

twopi

const double twopi = 2.0*constants::pi

static