CDPP / AMDA_Kernel

Commit e959cf9ce45a2077351b1a94ef2e51d9580be667

Authored by Hacene SI HADJ MOHAND
1 parent dbf061ca
Exists in master and in 80 other branches 10038, 10098, 5390, 5770, 6728, 7616, 7811, 7875, 8041, 8045_8044, 8295, 8459, 8528, 8928, 9016, 9016Bis, 9017, 9017-new, 9049, 9084, 9259, 9326, 9347, 9383, 9719, 9719.new, AllLegends, CDF_ISTP, CDF_ISTP_HOPEFULLY_ALL_WORKING, FER_10400, FER_10412, FER_10663, FER_10765, FER_11000, FER_11100, FER_11242, FER_11277, FER_11500, FER_11537, FER_11588, FER_11651, FER_19836, FER_5344, FER_6031, FER_6598, FER_6966, FER_7318, FER_7618, FER_7857, FER_8455, FER_9629, FER_CR41, FER_CR42, FER_CR43, FER_TFCat, Histo1D2D, JUNO_JEDI, MAZ_10936, MAZ_6018, MAZ_9016_New, MAZ_FIX_9016, SpeasyGet, amdadev, amdadev_test, cu2022, das2, density_maps, develop, evol_9042, feature/9259, frames-transformation, pickle, rm_6912, rm_7926, rm_8050, rm_8932, shue, shue97, test-merge, ubuntu

extention libraire 1945-2000

Showing 1 changed file with 18 additions and 0 deletions   Show diff stats
Inline Side-by-side
src/ExternLib/Maglib/Posmag.hh
  Diff comments   View file @e959cf9
... ... @@ -75,6 +75,8 @@ namespace AMDA {
75 75 double rrmag = std::sqrt(in[0]*in[0]+in[1]*in[1]+in[2]*in[2]);
76 76 double thetr = std::acos(in[3]/rrmag);
77 77 double phir = std::atan2(in[1] , in[0]);
  78 +
  79 +
78 80 /*
79 81 PAR magout (I) : type de champ magnetique externe
80 82  
... ... @@ -139,11 +141,27 @@ namespace AMDA {
139 141 double tgl = 1.0;
140 142 int ifail=10;
141 143  
  144 +
142 145 maglib::maglibWarpper::getPosmag(1900 + tmp->tm_year, 1 + tmp->tm_mon, tmp->tm_mday, tmp->tm_hour, tmp->tm_min, tmp->tm_sec,
143 146 rrmag, thetr, phir,
144 147 _isatex, _magout,
145 148 tgl, flg, xlamb,ifail);
146 149  
  150 +
  151 + // test results shud be as
  152 + /**
  153 + const double PI =3.141592653589793238463;
  154 + rrmag = 7.0; thetr=30*PI/180.0; phir=25*PI/180;
  155 + maglib::maglibWarpper::getPosmag(2017,02,15,11,30,0,
  156 + rrmag, thetr, phir,
  157 + _isatex, _magout,
  158 + tgl, flg, xlamb,ifail);
  159 + * expected results
  160 + * tgl 14.08
  161 + flg 18.76
  162 + xlamb 76.65
  163 + tglc 15.24
  164 + */
147 165 DataType tgl_res = tgl;
148 166 DataType xlamb_res = (DataType) xlamb*180/ (2*std::asin(1.0));
149 167 DataType flg_res = (DataType) flg;
... ...