diff options
Diffstat (limited to 'mass_profile/dump_fit_qdp.cpp')
| -rw-r--r-- | mass_profile/dump_fit_qdp.cpp | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/mass_profile/dump_fit_qdp.cpp b/mass_profile/dump_fit_qdp.cpp index 556edea..85307d1 100644 --- a/mass_profile/dump_fit_qdp.cpp +++ b/mass_profile/dump_fit_qdp.cpp @@ -10,13 +10,13 @@ namespace opt_utilities os<<"la x \"radius (kpc)\""<<std::endl; os<<"la y \"surface brightness (cts s\\u-1\\d pixel\\u-2\\d)\""<<std::endl; os<<"li on 2"<<std::endl; - for(int i=1;i<r.size();++i) + for(size_t i=1;i<r.size();++i) { os<<(r[i]+r[i-1])/2*cm_per_pixel/kpc<<"\t"<<(r[i]-r[i-1])/2*cm_per_pixel/kpc<<"\t"<<y[i-1]<<"\t"<<ye[i-1]<<std::endl; } os<<"no no no"<<std::endl; std::vector<double> p=f.get_all_params(); - for(int i=1;i<r.size();++i) + for(size_t i=1;i<r.size();++i) { double x=(r[i]+r[i-1])/2; os<<x*cm_per_pixel/kpc<<"\t"<<0<<"\t"<<f.eval_model_raw(x,p)<<"\t"<<0<<std::endl; @@ -30,26 +30,26 @@ namespace opt_utilities os<<"la x \"radius (kpc)\""<<std::endl; os<<"la y \"density (cm\\u-3\\d)\""<<std::endl; os<<"li on 2"<<std::endl; - - for(int i=1;i<r.size();++i) + + for(size_t i=1;i<r.size();++i) { double x=(r[i]+r[i-1])/2; double y=sbps[i-1]; double ye=sbpe[i-1]; os<<x*cm_per_pixel/kpc<<"\t0\t"<<y<<"\t"<<ye<<std::endl; } - + os<<"no no no"<<std::endl; std::vector<double> p=f.get_all_params(); std::vector<double> mv=f.eval_model_raw(r,p); - for(int i=1;i<r.size();++i) + for(size_t i=1;i<r.size();++i) { double x=(r[i]+r[i-1])/2; double y=mv[i-1]; double ye=0; os<<x*cm_per_pixel/kpc<<"\t0\t"<<y<<"\t"<<ye<<std::endl; } - + } }; |