00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023 #ifndef O2SCL_GSL_SERIES_H
00024 #define O2SCL_GSL_SERIES_H
00025
00026 #include <iostream>
00027 #include <o2scl/collection.h>
00028 #include <gsl/gsl_sum.h>
00029
00030 #ifndef DOXYGENP
00031 namespace o2scl {
00032 #endif
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043 class gsl_series {
00044 public:
00045
00046
00047
00048 gsl_series(int size=1);
00049
00050 virtual ~gsl_series();
00051
00052
00053
00054
00055 double series_accel(double *x, double &err);
00056
00057
00058
00059
00060 double series_accel_err(double *x, double &err);
00061
00062
00063 int set_size(int new_size);
00064
00065 #ifndef DOXYGENP
00066
00067 friend class io_tlate<gsl_series>;
00068
00069 protected:
00070
00071
00072 gsl_sum_levin_u_workspace *w;
00073
00074
00075 gsl_sum_levin_utrunc_workspace *wt;
00076
00077
00078 int size;
00079
00080 #endif
00081
00082 };
00083
00084 template<> int io_tlate<gsl_series>::input
00085 (cinput *co, in_file_format *ins, gsl_series *gs);
00086 template<> int io_tlate<gsl_series>::output
00087 (coutput *co, out_file_format *outs, gsl_series *gs);
00088 template<> const char *io_tlate<gsl_series>::type();
00089
00090 typedef io_tlate<gsl_series> gsl_series_io_type;
00091
00092 #ifndef DOXYGENP
00093 }
00094 #endif
00095
00096 #endif