diff options
Diffstat (limited to 'externals/grill/vasp/source/ops_cplx.cpp')
| -rw-r--r-- | externals/grill/vasp/source/ops_cplx.cpp | 30 |
1 files changed, 0 insertions, 30 deletions
diff --git a/externals/grill/vasp/source/ops_cplx.cpp b/externals/grill/vasp/source/ops_cplx.cpp index b2f01788..013571ef 100644 --- a/externals/grill/vasp/source/ops_cplx.cpp +++ b/externals/grill/vasp/source/ops_cplx.cpp @@ -16,13 +16,6 @@ WARRANTIES, see the file, "license.txt," in this distribution. // ----------------------------------------------------- -template<class T> V f_polar(T &rv,T &iv,T ra,T ia) { rv = sqrt(sqabs(ra,ia)),iv = arg(ra,ia); } -template<class T> V f_rect(T &rv,T &iv,T ra,T ia) { rv = ra*cos(ia),iv = ra*sin(ia); } - -BL VecOp::d_polar(OpParam &p) { d__cun(f_polar<S>,p); } -BL VecOp::d_rect(OpParam &p) { d__cun(f_rect<S>,p); } - - VASP_UNARY("vasp.polar",polar,true,"convert complex vector pair from rectangular to polar coordinates") VASP_UNARY("vasp.rect",rect,true,"convert complex vector pair from polar to rectangular coordinates") @@ -30,16 +23,6 @@ VASP_UNARY("vasp.rect",rect,true,"convert complex vector pair from polar to rect // ----------------------------------------------------- -template<class T> V f_radd(T &rv,T &iv,T ra,T ia,T rb,T) -{ - register const R _abs = sqrt(sqabs(ra,ia))+rb; - register const R _phi = arg(ra,ia); - - rv = _abs*cos(_phi),iv = _abs*sin(_phi); -} - -BL VecOp::d_radd(OpParam &p) { d__cbin(f_radd<S>,p); } - Vasp *VaspOp::m_radd(OpParam &p,CVasp &src,const Argument &arg,CVasp *dst) { Vasp *ret = NULL; @@ -65,22 +48,9 @@ VASP_ANYOP("vasp.r+",radd,0,true,VASP_ARG_R(0),"add offset to complex radius (of // ----------------------------------------------------- -template<class T> V f_cnorm(T &rv,T &iv,T ra,T ia) -{ - register T f = sqabs(ra,ia); - if(f) { f = 1./sqrt(f); rv = ra*f,iv = ia*f; } - else rv = iv = 0; -} - -BL VecOp::d_cnorm(OpParam &p) { d__cun(f_cnorm<S>,p); } - VASP_UNARY("vasp.cnorm",cnorm,true,"normalize complex radius to 1 (but preserve angle)") // ----------------------------------------------------- -template<class T> inline V f_cconj(T &,T &iv,T,T ia) { iv = -ia; } - -BL VecOp::d_cconj(OpParam &p) { D__cun(f_cconj<S>,p); } - VASP_UNARY("vasp.cconj",cconj,true,"complex conjugate: multiply imaginary part with -1") // should be replaced by an abstraction |