newsolver.cc File Reference

#include <madness/mra/mra.h>
#include <madness/mra/operator.h>
#include "nonlinsol.h"

Include dependency graph for newsolver.cc:

Functions
double	rifunction (const coord_3d &r)
double	iterate_ground (World &world, NonlinearSolver &solver, real_function_3d &V, real_function_3d &psi, double &eps)
double	iterate_excite (World &world, NonlinearSolver &solver, real_function_3d &V, real_function_3d &psi, real_function_3d &dpsi, double &eps, real_function_3d &ri)
double	inner (const F &a, const F &b)
template<class solverT >
double	iterate_xy (World &world, solverT &solver, const real_function_3d &V, const real_function_3d &psi, double &eps, const real_function_3d &ri, real_function_3d &x, real_function_3d &y, const double omega)
	solve the dynamic response equations More...
int	main (int argc, char **argv)

Function Documentation

double inner	(	const F &	a,
		const F &	b
	)

References madness::inner().

double iterate_excite	(	World &	world,
		NonlinearSolver &	solver,
		real_function_3d &	V,
		real_function_3d &	psi,
		real_function_3d &	dpsi,
		double &	eps,
		real_function_3d &	ri
	)

References madness::Function< T, NDIM >::norm2(), op(), madness::print(), madness::World::rank(), madness::Function< T, NDIM >::scale(), sqrt(), madness::Function< T, NDIM >::truncate(), and madness::NonlinearSolverND< NDIM >::update().

Referenced by main().

double iterate_ground	(	World &	world,
		NonlinearSolver &	solver,
		real_function_3d &	V,
		real_function_3d &	psi,
		double &	eps
	)

References madness::inner(), madness::norm(), madness::Function< T, NDIM >::norm2(), op(), madness::print(), psi(), madness::World::rank(), madness::Function< T, NDIM >::scale(), sqrt(), and madness::NonlinearSolverND< NDIM >::update().

Referenced by main().

template<class solverT >

double iterate_xy	(	World &	world,
		solverT &	solver,
		const real_function_3d &	V,
		const real_function_3d &	psi,
		double &	eps,
		const real_function_3d &	ri,
		real_function_3d &	x,
		real_function_3d &	y,
		const double	omega
	)

solve the dynamic response equations

Parameters

[in]	world	the world
[in]	solver	the KAIN solver
[in]	V	the local potential (here: V_nuc + 1/2 J )
[in]	psi	the orbital
[in]	eps	orbital energy for psi
[in]	ri	the external perturbation (here: z)
[in,out]	x	the x part of the response vector
[in,out]	y	the y part of the response vector
[in]	omega	the frequency of the external perturbation

Returns

the current error in the residual of the response equations

References madness::inner(), madness::norm2(), madness::print(), psi(), madness::Function< T, NDIM >::scale(), sqrt(), and madness::XNonlinearSolver< T, C, Alloc >::update().

Referenced by main().

int main	(	int	argc,
		char **	argv
	)

References SafeMPI::COMM_WORLD, e1(), madness::WorldGopInterface::fence(), madness::finalize(), madness::World::gop, guess(), madness::initialize(), madness::inner(), iterate_excite(), iterate_ground(), iterate_xy(), k, L, madness::Function< T, NDIM >::norm2(), op(), madness::plot_line(), potential(), madness::print(), psi(), R, madness::World::rank(), rifunction(), madness::Function< T, NDIM >::scale(), madness::FunctionDefaults< NDIM >::set_cubic_cell(), madness::FunctionDefaults< NDIM >::set_initial_level(), madness::FunctionDefaults< NDIM >::set_k(), madness::FunctionDefaults< NDIM >::set_refine(), madness::FunctionDefaults< NDIM >::set_thresh(), madness::FunctionDefaults< NDIM >::set_truncate_mode(), madness::square(), madness::startup(), thresh, madness::Function< T, NDIM >::truncate(), V(), and madness::vec().

double rifunction

(

const coord_3d &

r

)

Referenced by main().