# Diff of /trunk/base/generic/solver/calc.h

revision 33 by ben.allan, Sun Dec 26 20:06:01 2004 UTC revision 34 by ben.allan, Sun Dec 26 20:35:09 2004 UTC
# Line 69  extern boolean calc_print_errors; Line 69  extern boolean calc_print_errors;
69
70  #define calc_ERF_COEF    1.1283791670955130   /**<  = 2 / sqrt(PI) */  #define calc_ERF_COEF    1.1283791670955130   /**<  = 2 / sqrt(PI) */
71  #define calc_LOG10_COEF  0.4342944819032518   /**< = log10(e) = 1/ln(10) */  #define calc_LOG10_COEF  0.4342944819032518   /**< = log10(e) = 1/ln(10) */
72  #define calc_PI          3.1415926535897930 /* really should be using M_PI */  #define calc_PI          3.1415926535897930 /**< really should be using M_PI */
73
74  extern real64 calc_upower(real64, unsigned);  extern real64 calc_upower(real64, unsigned);
75  extern real64 calc_factorial(unsigned);  extern real64 calc_factorial(unsigned);
# Line 81  extern real64 calc_erf_inv(); Line 81  extern real64 calc_erf_inv();
81  #endif /* HAVE_ERF */  #endif /* HAVE_ERF */
82  extern real64 calc_lnm_inv();  extern real64 calc_lnm_inv();
83  extern int calc_is_int();  extern int calc_is_int();
84  /**  /**<
85   ***  value = calc_upower(x,n)   Compute x^n.   ***  value = calc_upower(x,n)   Compute x^n.
86   ***  value = calc_factorial(n)  Compute n!.   ***  value = calc_factorial(n)  Compute n!.
87   ***  value = calc_rec(x)        Compute 1/x.   ***  value = calc_rec(x)        Compute 1/x.
# Line 151  extern real64 calc_log_D1(double); Line 151  extern real64 calc_log_D1(double);
151  extern real64 calc_sqrt_D1(double);  extern real64 calc_sqrt_D1(double);
152  extern real64 calc_cbrt_D1(double);  extern real64 calc_cbrt_D1(double);
153  extern real64 calc_fabs_D1(double);  extern real64 calc_fabs_D1(double);
154  /**  /**<
155   ***  value = calc_<uop>_D1(x)   ***  value = calc_<uop>_D1(x)
156   ***  real64 value,x;   ***  real64 value,x;
157   ***   ***
# Line 182  extern real64 calc_log_D2(); Line 182  extern real64 calc_log_D2();
182  extern real64 calc_sqrt_D2();  extern real64 calc_sqrt_D2();
183  extern real64 calc_cbrt_D2();  extern real64 calc_cbrt_D2();
184  extern real64 calc_fabs_D2();  extern real64 calc_fabs_D2();
185  /**  /**<
186   ***  value = calc_<uop>_D2(x)   ***  value = calc_<uop>_D2(x)
187   ***  real64 value,x;   ***  real64 value,x;
188   ***   ***
# Line 204  extern real64 calc_sqr_Dn(); Line 204  extern real64 calc_sqr_Dn();
204  extern real64 calc_sqrt_Dn();  extern real64 calc_sqrt_Dn();
205  extern real64 calc_tan_Dn();  extern real64 calc_tan_Dn();
206  extern real64 calc_fabs_Dn();  extern real64 calc_fabs_Dn();
207  /**  /**<
208   ***  value = calc_<uop>_Dn(x,n)   ***  value = calc_<uop>_Dn(x,n)
209   ***  real64 value,x;   ***  real64 value,x;
210   ***  int n;   n >= 0   ***  int n;   n >= 0
# Line 218  extern real64 calc_fabs_Dn(); Line 218  extern real64 calc_fabs_Dn();
218  #define calc_div_D0(x,y) calc_mul_D0(x,calc_rec(y))  #define calc_div_D0(x,y) calc_mul_D0(x,calc_rec(y))
219  #define calc_ipow_D0(x,y) asc_ipow(x,(int)y)  #define calc_ipow_D0(x,y) asc_ipow(x,(int)y)
220  extern real64 calc_pow_D0();  extern real64 calc_pow_D0();
221  /**  /**<
222   ***  value = calc_<binop>_D0(x,y)   ***  value = calc_<binop>_D0(x,y)
223   ***  real64 value,x,y;   ***  real64 value,x,y;
224   ***   ***
# Line 231  extern real64 calc_pow_D0(); Line 231  extern real64 calc_pow_D0();
231  extern real64 calc_div_D1();  extern real64 calc_div_D1();
232  #define calc_ipow_D1(x,y) asc_d1ipow(x,(int)y)      #define calc_ipow_D1(x,y) asc_d1ipow(x,(int)y)
233  extern real64 calc_pow_D1();  extern real64 calc_pow_D1();
234  /**  /**<
235   ***  diff = calc_<binop>_D1(x,y,wrt)   ***  diff = calc_<binop>_D1(x,y,wrt)
236   ***  real64 diff,x,y;   ***  real64 diff,x,y;
237   ***  int wrt;   0 ==> w.r.t 1st arg, 1 ==> w.r.t. 2nd arg   ***  int wrt;   0 ==> w.r.t 1st arg, 1 ==> w.r.t. 2nd arg
# Line 259  extern real64 calc_sub_Dn(); Line 259  extern real64 calc_sub_Dn();
259  extern real64 calc_mul_Dn();  extern real64 calc_mul_Dn();
260  extern real64 calc_div_Dn();  extern real64 calc_div_Dn();
261  extern real64 calc_pow_Dn();  extern real64 calc_pow_Dn();
262  /**  /**<
263   ***  diffn = calc_<binop>_Dn(x,y,nwrt0,nwrt1)   ***  diffn = calc_<binop>_Dn(x,y,nwrt0,nwrt1)
264   ***  real64 diffn,x,y;   ***  real64 diffn,x,y;
265   ***  int nwrt0,nwrt1;   ***  int nwrt0,nwrt1;

Legend:
 Removed from v.33 changed lines Added in v.34