1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
|
/*
VASP modular - vector assembling signal processor / objects for Max/MSP and PD
Copyright (c) 2002 Thomas Grill (xovo@gmx.net)
For information on usage and redistribution, and for a DISCLAIMER OF ALL
WARRANTIES, see the file, "license.txt," in this distribution.
*/
#include "main.h"
#include "ops_cmp.h"
#include "opdefs.h"
#include "util.h"
#include <math.h>
// --------------------------------------------------------------
template<class T> V f_gate(T &rv,T ra,T rb) { rv = fabs(ra) >= rb?ra:0; }
template<class T> V f_igate(T &rv,T ra,T rb) { rv = fabs(ra) <= rb?ra:0; }
template<class T> V f_rgate(T &rv,T &iv,T ra,T ia,T rb,T)
{
register const T _abs = sqabs(ra,ia);
if(_abs >= rb*rb) rv = ra,iv = ia;
else rv = iv = 0;
}
template<class T> V f_rigate(T &rv,T &iv,T ra,T ia,T rb,T)
{
register const T _abs = sqabs(ra,ia);
if(_abs <= rb*rb) rv = ra,iv = ia;
else rv = iv = 0;
}
BL VecOp::d_gate(OpParam &p) { D__rbin(f_gate<S>,p); }
BL VecOp::d_igate(OpParam &p) { d__rbin(f_igate<S>,p); }
BL VecOp::d_rgate(OpParam &p) { d__cbin(f_rgate<S>,p); }
BL VecOp::d_rigate(OpParam &p) { d__cbin(f_rigate<S>,p); }
Vasp *VaspOp::m_gate(OpParam &p,Vasp &src,const Argument &arg,Vasp *dst)
{
Vasp *ret = NULL;
RVecBlock *vecs = GetRVecs(p.opname,src,dst);
if(vecs) {
if(arg.IsList() && arg.GetList().Count() >= 1 && flext::CanbeFloat(arg.GetList()[0]))
p.rbin.arg = flext::GetAFloat(arg.GetList()[0]);
else {
post("%s - argument is invalid -> set to 1",p.opname);
p.rbin.arg = 1;
}
ret = DoOp(vecs,VecOp::d_gate,p);
delete vecs;
}
return ret;
}
Vasp *VaspOp::m_rgate(OpParam &p,Vasp &src,const Argument &arg,Vasp *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 1",p.opname);
p.cbin.rarg = 1;
}
p.cbin.iarg = 0; // not used anyway
ret = DoOp(vecs,VecOp::d_rgate,p);
delete vecs;
}
return ret;
}
VASP_ANYOP("vasp.gate",gate,1,true,VASP_ARG_R(1),"set destination to 0 if source < argument")
VASP_ANYOP("vasp.rgate",rgate,1,true,VASP_ARG_R(1),"complex radius gate: set destination to 0 if rad(complex source) < rad(complex argument)")
// --------------------------------------------------------------
|
