luke
/
ngspice
mirror of https://git.code.sf.net/p/ngspice/ngspice
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

add the atan2 derivatives

This commit is contained in:
dwarning 2017-09-21 15:28:18 +02:00 committed by Holger Vogt
parent a4f696720c
commit 67d52e8cba
2 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/spicelib/devices/adms/admst/analogfunction.xml
View File

@ -76,7 +76,8 @@ inline double _d0_hypot(double x,double y) { return (x)/sqrt((x)*(x)+(y)*(y)); }
inline double _d1_hypot(double x,double y) { return (y)/sqrt((x)*(x)+(y)*(y)); }
inline double _atan2(double x,double y) { return atan2(x,y); }
// TODO atan2 derivatives?
inline double _d0_atan2(double x,double y) { return (-y)/((x)*(x)+(y)*(y)); }
inline double _d1_atan2(double x,double y) { return (x)/((x)*(x)+(y)*(y)); }
inline double _max(double x,double y) { return ((x)>(y))?(x):(y); }
inline double _d0_max(double x,double y) { return ((x)>(y))?1.0:0.0; }

3
src/spicelib/devices/adms/admst/ngspiceMODULE.hxx.xml
View File

@ -65,7 +65,8 @@
#define m00_limexp(v00,x) v00 = ((x)<90.0?exp(x):EXP90*(x-89.0));
#define m10_limexp(v10,v00,x) v10 = ((x)<90.0?(v00):EXP90);
#define m00_atan2(v00,x,y) v00 = atan2(x,y);
// TODO atan2 derivatives ?
#define m10_atan2(v10,x,y) v10 = (-y/(x*x+y*y));
#define m11_atan2(v11,x,y) v11 = (x/(x*x+y*y));
#define m00_acosh(v00,x) v00 = acosh(x);
#define m10_acosh(v10,v00,x) v10 = (1.0/(sqrt(x-1)*sqrt(x+1)));
#define m00_asinh(v00,x) v00 = asinh(x);