diff options
Diffstat (limited to 'externals/grill/vasp/source/ops_arith.cpp')
-rw-r--r-- | externals/grill/vasp/source/ops_arith.cpp | 30 |
1 files changed, 0 insertions, 30 deletions
diff --git a/externals/grill/vasp/source/ops_arith.cpp b/externals/grill/vasp/source/ops_arith.cpp index 6faa1aad..5eb11e7d 100644 --- a/externals/grill/vasp/source/ops_arith.cpp +++ b/externals/grill/vasp/source/ops_arith.cpp @@ -14,22 +14,6 @@ WARRANTIES, see the file, "license.txt," in this distribution. #include "util.h" #include <math.h> -template<class T> inline V f_radd(T &v,T a,T b) { v = a+b; } -template<class T> inline V f_rsub(T &v,T a,T b) { v = a-b; } -template<class T> inline V f_rsubr(T &v,T a,T b) { v = b-a; } -template<class T> inline V f_rmul(T &v,T a,T b) { v = a*b; } -template<class T> inline V f_rdiv(T &v,T a,T b) { v = a/b; } -template<class T> inline V f_rdivr(T &v,T a,T b) { v = b/a; } -template<class T> inline V f_rmod(T &v,T a,T b) { v = fmod(a,b); } - -BL VecOp::d_add(OpParam &p) { D__rbin(f_radd<S>,p); } -BL VecOp::d_sub(OpParam &p) { D__rbin(f_rsub<S>,p); } -BL VecOp::d_subr(OpParam &p) { D__rbin(f_rsubr<S>,p); } -BL VecOp::d_mul(OpParam &p) { D__rbin(f_rmul<S>,p); } -BL VecOp::d_div(OpParam &p) { D__rbin(f_rdiv<S>,p); } -BL VecOp::d_divr(OpParam &p) { D__rbin(f_rdivr<S>,p); } -BL VecOp::d_mod(OpParam &p) { D__rbin(f_rmod<S>,p); } - VASP_BINARY("vasp.+",add,true,VASP_ARG_R(0),"Adds a value, envelope or vasp") VASP_BINARY("vasp.-",sub,true,VASP_ARG_R(0),"Subtracts a value, envelope or vasp") @@ -41,25 +25,11 @@ VASP_BINARY("vasp.%",mod,true,VASP_ARG_R(0),"Calculates the remainder of the div // ----------------------------------------------------- -template<class T> inline V f_rsqr(T &v,T a) { v = a*a; } -template<class T> inline V f_rssqr(T &v,T a) { v = a*fabs(a); } - -BL VecOp::d_sqr(OpParam &p) { D__run(f_rsqr<S>,p); } -BL VecOp::d_ssqr(OpParam &p) { d__run(f_rssqr<S>,p); } - VASP_UNARY("vasp.sqr",sqr,true,"Calculates the square") VASP_UNARY("vasp.ssqr",ssqr,true,"Calculates the square with preservation of the sign") - // ----------------------------------------------------- -template<class T> inline V f_rsign(T &v,T a) { v = (a == 0?0:(a < 0?-1.:1.)); } -template<class T> inline V f_rabs(T &v,T a) { v = fabs(a); } - -BL VecOp::d_sign(OpParam &p) { D__run(f_rsign<S>,p); } -BL VecOp::d_abs(OpParam &p) { D__run(f_rabs<S>,p); } - - VASP_UNARY("vasp.sign",sign,true,"Calculates the sign (signum function)") VASP_UNARY("vasp.abs",abs,true,"Calulates the absolute value") |