LMU/musrfit
LMU/musrfit
5
0
Fork 0
Code Issues 6 Pull Requests Actions Packages Projects Releases 3 Activity

fixed bug with the call to averaged function

Browse Source
This commit is contained in:
suter_a 2014-01-27 09:55:11 +00:00
parent 7665696e38
commit 135ca8f99b
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/external/libPhotoMeissner/test/photo_meissner.cpp vendored
View File

@ -409,7 +409,7 @@ void calcFields(const double zStart, const double zEnd, const PhotoParam param,
bb.push_back(b); bb.push_back(b);
bbL.push_back(bL); bbL.push_back(bL);
z += dz; z += dz;
} while (z <= zEnd); } while (z <= zEnd+dz);
} else { // film } else { // film
double beta = exp(-(param.filmThickness/2.0-param.deadLayer)/param.z0); double beta = exp(-(param.filmThickness/2.0-param.deadLayer)/param.z0);
double NN = InuMinus(nuPlus, nuPhoto*beta)*gsl_sf_bessel_Inu(nuPlus, nuPhoto) - double NN = InuMinus(nuPlus, nuPhoto*beta)*gsl_sf_bessel_Inu(nuPlus, nuPhoto) -
@ -428,7 +428,7 @@ void calcFields(const double zStart, const double zEnd, const PhotoParam param,
bb.push_back(b); bb.push_back(b);
bbL.push_back(bL); bbL.push_back(bL);
z += dz; z += dz;
} while (z <= zEnd); } while (z <= zEnd+dz);
} }
} }
@ -842,8 +842,8 @@ int main(int argc, char *argv[])
nn.clear(); nn.clear();
matched_nz(rgeZ[0], rgeN[0], zz, nn); matched_nz(rgeZ[0], rgeN[0], zz, nn);
zMean.push_back(averaged(zz, nn)); zMean.push_back(averaged(zz, nn));
BMean.push_back(averaged(zz, bb)); BMean.push_back(averaged(bb, nn));
BLMean.push_back(averaged(zz, bbL)); BLMean.push_back(averaged(bbL, nn));
rgeFlnList.push_back(rgeFln); rgeFlnList.push_back(rgeFln);
if (!strcmp(fln, "")) if (!strcmp(fln, ""))
cout << "<z>=" << zMean[0] << " (nm), <B>/Bext=" << BMean[0] << ", <BL>/Bext=" << BLMean[0] << " for rge-file: '" << rgeFln << "'." << endl; cout << "<z>=" << zMean[0] << " (nm), <B>/Bext=" << BMean[0] << ", <BL>/Bext=" << BLMean[0] << " for rge-file: '" << rgeFln << "'." << endl;