home *** CD-ROM | disk | FTP | other *** search
/ InfoMagic Source Code 1993 July / THE_SOURCE_CODE_CD_ROM.iso / gnu / nihcl-3.0 / ex / AutoDeriv.h < prev    next >
Encoding:
C/C++ Source or Header  |  1990-05-15  |  2.2 KB  |  78 lines
  1. #ifndef AutoDerivH
  2. #define AutoDerivH
  3.  
  4. // AutoDeriv.h -- one-dimensional automatic derivative
  5.  
  6. // $Header: /afs/alw.nih.gov/unix/sun4_40c/usr/local/src/nihcl-3.0/share/ex/RCS/AutoDeriv.h,v 3.0 90/05/15 22:43:25 kgorlen Rel $
  7.  
  8. #include <math.h>
  9.  
  10. class ostream;
  11.  
  12. class AutoDeriv {
  13. private:
  14.     double u,du;
  15. public:
  16.     AutoDeriv(double v = 0.0,double dv = 0.0) { u=v; du=dv; }
  17.  
  18.     friend double F(const AutoDeriv& x) { return x.u; }
  19.     friend double dF(const AutoDeriv& x) { return x.du; }
  20.  
  21.     void operator=(const AutoDeriv& x);
  22.     int operator==(const AutoDeriv& x) const;
  23.  
  24.     AutoDeriv operator-() const;
  25.     AutoDeriv pow(int e) const;
  26.  
  27.     friend AutoDeriv operator+(const AutoDeriv&,
  28.                                const AutoDeriv&);
  29.     friend AutoDeriv operator-(const AutoDeriv&,
  30.                                const AutoDeriv&);
  31.     friend AutoDeriv operator*(const AutoDeriv&,
  32.                                const AutoDeriv&);
  33.     friend AutoDeriv operator/(const AutoDeriv&,
  34.                                const AutoDeriv&);
  35.     void printOn(ostream& strm) const;
  36. };
  37. inline AutoDeriv operator+(const AutoDeriv& x,
  38.                            const AutoDeriv& y)
  39. { 
  40.     return AutoDeriv(F(x)+F(y),dF(x)+dF(y)); 
  41. }
  42. inline AutoDeriv operator-(const AutoDeriv& x,
  43.                            const AutoDeriv& y)
  44. { 
  45.     return AutoDeriv(F(x)-F(y),dF(x)-dF(y)); 
  46. }
  47. inline AutoDeriv operator*(const AutoDeriv& x,
  48.                            const AutoDeriv& y)
  49. {
  50.     return AutoDeriv(F(x)*F(y),
  51.                      F(y)*dF(x)+F(x)*dF(y));
  52. }
  53. inline AutoDeriv operator/(const AutoDeriv& x,
  54.                            const AutoDeriv& y)
  55. { 
  56.     return AutoDeriv(F(x)/F(y),
  57.                      (F(y)*dF(x)-F(x)*dF(y))/F(y)*F(y)); 
  58. }
  59. inline AutoDeriv Exp(const AutoDeriv& x) 
  60. { 
  61.     return AutoDeriv(exp(F(x)),dF(x)*exp(F(x))); 
  62. }
  63. inline AutoDeriv Log(const AutoDeriv& x) 
  64. { 
  65.     return AutoDeriv(log(F(x)),dF(x)/F(x)); 
  66. }
  67. inline AutoDeriv Sin(const AutoDeriv& x) 
  68. { 
  69.     return AutoDeriv(sin(F(x)),dF(x)*cos(F(x))); 
  70. }
  71. inline AutoDeriv Cos(const AutoDeriv& x) 
  72. { 
  73.     return AutoDeriv(cos(F(x)),-dF(x)*sin(F(x))); 
  74. }
  75. extern ostream& operator<<(ostream& strm, 
  76.                            const AutoDeriv& x);
  77. #endif/* AutoDerivH */
  78.