00001
00010
00011
00012
00013
00014 #ifndef NUM
00015 # define NUM double
00016 #endif
00017
00018 #ifdef PRAGMA_I
00019
00020 #endif
00021 #include "my_nr.h"
00022
00023
00024
00025 #define T NUM
00026 #include "my_nr_inst.h"
00027
00028 #if 0
00029 #define TEMPLATE template
00030
00031
00032
00033 TEMPLATE int basis_trafo <NUM> (long int value, long int basis, Vector<NUM>& index);
00034
00035 TEMPLATE NUM partial_del <NUM> (NUM (*func)(const Vector<NUM>&, const Vector<NUM>&),
00036 const Vector<NUM>& x, const Vector<NUM>& p,
00037 int mu, NUM h, NUM& err);
00038
00039 TEMPLATE NUM chisq <NUM> (Vector<NUM>& data, Vector<NUM>& func, Vector<NUM>& sig);
00040
00041 TEMPLATE NUM chisquare <NUM> (NUM (*func)(const Vector<NUM>&, const Vector<NUM>&),
00042 Vector<NUM>& param, Vector<NUM>& x,
00043 Vector<NUM>& y, Vector<NUM>& sigma);
00044
00045
00046
00047
00048
00049
00050 TEMPLATE TVector<NUM> grid_min <NUM> (NUM (*func)(const Vector<NUM>&, const Vector<NUM>&),
00051 Vector<NUM>& x, Vector<NUM>& y, Vector<NUM>& sig,
00052 Vector<NUM>& pmin, Vector<NUM>& pmax, Vector<NUM>& lambda,
00053 long int res, int max_it, char v);
00054
00055 TEMPLATE NUM lev_mar <NUM> (Vector<NUM>& x, Vector<NUM>& y, Vector<NUM>& sig,
00056 Vector<NUM>& param, Vector<NUM>& plow, Vector<NUM>& phigh,
00057 NUM (*func)(const Vector<NUM>&, const Vector<NUM>&),
00058 Vector<NUM> lambda, Vector<NUM>& lscale, NUM tol,
00059 int iter, char ver);
00060
00061 TEMPLATE void hpsort <NUM> (Matrix<NUM>& x, int pos);
00062
00063
00064
00065
00066
00067
00068 #endif
00069
00070