/*! \brief Finds peaks or valleys by radius in a complex vasp. \param arg argument list \param arg.rep repetition count \param dst destination vasp (NULL for in-place operation) \param inv true for valley operation \return normalized destination vasp */ Vasp *VaspOp::m_rpeaks(OpParam &p,CVasp &src,CVasp *dst,BL inv) { Vasp *ret = NULL; CVecBlock *vecs = GetCVecs(p.opname,src,dst); if(vecs) { p.peaks.cx = true; ret = DoOp(vecs,inv?d_rvalleys:d_rpeaks,p); delete vecs; } return ret; }
Vasp *VaspOp::m_ccopy(OpParam &p,CVasp &src,CVasp &arg) { Vasp *s = NULL,*d = NULL; CVecBlock *vecs = GetCVecs(p.opname,src,&arg); if(vecs) { d = DoOp(vecs,VecOp::d_ccopy,p); s = vecs->SrcVasp(); if(d) arg = *d; else { arg.Clear(); delete d; } delete vecs; } return s; }
Vasp *VaspOp::m_radd(OpParam &p,CVasp &src,const Argument &arg,CVasp *dst) { Vasp *ret = NULL; CVecBlock *vecs = GetCVecs(p.opname,src,dst); if(vecs) { if(arg.IsList() && arg.GetList().Count() >= 1 && flext::CanbeFloat(arg.GetList()[0])) p.cbin.rarg = flext::GetAFloat(arg.GetList()[0]); else { post("%s - argument is invalid -> set to 0",p.opname); p.cbin.rarg = 0; } p.cbin.iarg = 0; // not used anyway ret = DoOp(vecs,VecOp::d_radd,p); delete vecs; } return ret; }
Vasp *VaspOp::m_cpowi(OpParam &p,CVasp &src,const Argument &arg,CVasp *dst) { Vasp *ret = NULL; CVecBlock *vecs = GetCVecs(p.opname,src,dst); if(vecs) { I powi = 1; if(arg.IsList() && arg.GetList().Count() >= 1 && flext::CanbeInt(arg.GetList()[0])) powi = flext::GetAInt(arg.GetList()[0]); else post("%s - power arg is invalid -> set to 1",p.opname); if(powi < 0) { post("%s - negative integer power is not allowed",p.opname); } else { switch(powi) { case 0: { p.cbin.rarg = 1,p.cbin.iarg = 0; ret = DoOp(vecs,VecOp::d_cset,p); break; } case 1: { // set arg to src ret = DoOp(vecs,VecOp::d_ccopy,p); break; } case 2: { ret = DoOp(vecs,VecOp::d_csqr,p); break; } default: { p.ibin.arg = powi; ret = DoOp(vecs,VecOp::d_cpowi,p); break; } } } delete vecs; } return ret; }