Skip to content
Snippets Groups Projects
Commit e9e4e7cc authored by Wuttke, Joachim's avatar Wuttke, Joachim
Browse files

same name for member and argument variables

parent 81ed04f8
No related branches found
No related tags found
No related merge requests found
...@@ -35,7 +35,7 @@ DecouplingApproximationStrategy::scalarCalculation(const SimulationElement& sim_ ...@@ -35,7 +35,7 @@ DecouplingApproximationStrategy::scalarCalculation(const SimulationElement& sim_
{ {
double intensity = 0.0; double intensity = 0.0;
complex_t amplitude = complex_t(0.0, 0.0); complex_t amplitude = complex_t(0.0, 0.0);
for (auto& ffw : m_formfactor_wrappers) { for (auto& ffw : m_weighted_formfactors) {
complex_t ff = ffw.evaluate(sim_element); complex_t ff = ffw.evaluate(sim_element);
if (std::isnan(ff.real())) if (std::isnan(ff.real()))
throw Exceptions::RuntimeErrorException( throw Exceptions::RuntimeErrorException(
...@@ -57,7 +57,7 @@ DecouplingApproximationStrategy::polarizedCalculation(const SimulationElement& s ...@@ -57,7 +57,7 @@ DecouplingApproximationStrategy::polarizedCalculation(const SimulationElement& s
Eigen::Matrix2cd mean_amplitude = Eigen::Matrix2cd::Zero(); Eigen::Matrix2cd mean_amplitude = Eigen::Matrix2cd::Zero();
const auto& polarization_handler = sim_element.polarizationHandler(); const auto& polarization_handler = sim_element.polarizationHandler();
for (auto& ffw : m_formfactor_wrappers) { for (auto& ffw : m_weighted_formfactors) {
Eigen::Matrix2cd ff = ffw.evaluatePol(sim_element); Eigen::Matrix2cd ff = ffw.evaluatePol(sim_element);
if (!ff.allFinite()) if (!ff.allFinite())
throw Exceptions::RuntimeErrorException( throw Exceptions::RuntimeErrorException(
......
...@@ -23,14 +23,14 @@ ...@@ -23,14 +23,14 @@
IInterferenceFunctionStrategy::IInterferenceFunctionStrategy( IInterferenceFunctionStrategy::IInterferenceFunctionStrategy(
const std::vector<FormFactorCoherentSum>& weighted_formfactors, const std::vector<FormFactorCoherentSum>& weighted_formfactors,
const IInterferenceFunction* p_iff, const SimulationOptions& sim_params, bool polarized) const IInterferenceFunction* p_iff, const SimulationOptions& sim_params, bool polarized)
: m_formfactor_wrappers(weighted_formfactors) : m_weighted_formfactors(weighted_formfactors)
, m_iff(p_iff ? p_iff->clone() : new InterferenceFunctionNone()) , m_iff(p_iff ? p_iff->clone() : new InterferenceFunctionNone())
, m_options(sim_params) , m_options(sim_params)
, m_polarized(polarized) , m_polarized(polarized)
, m_integrator( , m_integrator(
make_integrator_miser(this, &IInterferenceFunctionStrategy::evaluate_for_fixed_angles, 2)) make_integrator_miser(this, &IInterferenceFunctionStrategy::evaluate_for_fixed_angles, 2))
{ {
ASSERT(!m_formfactor_wrappers.empty()); ASSERT(!m_weighted_formfactors.empty());
} }
IInterferenceFunctionStrategy::~IInterferenceFunctionStrategy() = default; IInterferenceFunctionStrategy::~IInterferenceFunctionStrategy() = default;
......
...@@ -50,7 +50,7 @@ public: ...@@ -50,7 +50,7 @@ public:
double evaluate(const SimulationElement& sim_element) const; double evaluate(const SimulationElement& sim_element) const;
protected: protected:
std::vector<FormFactorCoherentSum> m_formfactor_wrappers; std::vector<FormFactorCoherentSum> m_weighted_formfactors;
std::unique_ptr<IInterferenceFunction> m_iff; std::unique_ptr<IInterferenceFunction> m_iff;
const SimulationOptions m_options; const SimulationOptions m_options;
......
...@@ -55,7 +55,7 @@ SSCApproximationStrategy::SSCApproximationStrategy( ...@@ -55,7 +55,7 @@ SSCApproximationStrategy::SSCApproximationStrategy(
, m_kappa(kappa) , m_kappa(kappa)
{ {
m_mean_radius = 0.0; m_mean_radius = 0.0;
for (auto& ffw : m_formfactor_wrappers) for (auto& ffw : m_weighted_formfactors)
m_mean_radius += ffw.relativeAbundance() * ffw.radialExtension(); m_mean_radius += ffw.relativeAbundance() * ffw.radialExtension();
} }
...@@ -67,7 +67,7 @@ double SSCApproximationStrategy::scalarCalculation(const SimulationElement& sim_ ...@@ -67,7 +67,7 @@ double SSCApproximationStrategy::scalarCalculation(const SimulationElement& sim_
double qp = sim_element.getMeanQ().magxy(); double qp = sim_element.getMeanQ().magxy();
double diffuse_intensity = 0.0; double diffuse_intensity = 0.0;
complex_t ff_orig = 0., ff_conj = 0.; // original and conjugated mean formfactor complex_t ff_orig = 0., ff_conj = 0.; // original and conjugated mean formfactor
for (auto& ffw : m_formfactor_wrappers) { for (auto& ffw : m_weighted_formfactors) {
complex_t ff = ffw.evaluate(sim_element); complex_t ff = ffw.evaluate(sim_element);
double fraction = ffw.relativeAbundance(); double fraction = ffw.relativeAbundance();
diffuse_intensity += fraction * std::norm(ff); diffuse_intensity += fraction * std::norm(ff);
...@@ -78,7 +78,7 @@ double SSCApproximationStrategy::scalarCalculation(const SimulationElement& sim_ ...@@ -78,7 +78,7 @@ double SSCApproximationStrategy::scalarCalculation(const SimulationElement& sim_
ff_conj += prefac * std::conj(ff); ff_conj += prefac * std::conj(ff);
} }
complex_t mean_ff_norm = ff_orig * ff_conj; complex_t mean_ff_norm = ff_orig * ff_conj;
complex_t p2kappa = getCharacteristicSizeCoupling(qp, m_formfactor_wrappers); complex_t p2kappa = getCharacteristicSizeCoupling(qp, m_weighted_formfactors);
complex_t omega = getCharacteristicDistribution(qp, m_iff.get()); complex_t omega = getCharacteristicDistribution(qp, m_iff.get());
double iff = 2.0 * (mean_ff_norm * omega / (1.0 - p2kappa * omega)).real(); double iff = 2.0 * (mean_ff_norm * omega / (1.0 - p2kappa * omega)).real();
double dw_factor = m_iff->DWfactor(sim_element.getMeanQ()); double dw_factor = m_iff->DWfactor(sim_element.getMeanQ());
...@@ -93,7 +93,7 @@ double SSCApproximationStrategy::polarizedCalculation(const SimulationElement& s ...@@ -93,7 +93,7 @@ double SSCApproximationStrategy::polarizedCalculation(const SimulationElement& s
const auto& polarization_handler = sim_element.polarizationHandler(); const auto& polarization_handler = sim_element.polarizationHandler();
Eigen::Matrix2cd ff_orig = Eigen::Matrix2cd::Zero(); Eigen::Matrix2cd ff_orig = Eigen::Matrix2cd::Zero();
Eigen::Matrix2cd ff_conj = Eigen::Matrix2cd::Zero(); Eigen::Matrix2cd ff_conj = Eigen::Matrix2cd::Zero();
for (auto& ffw : m_formfactor_wrappers) { for (auto& ffw : m_weighted_formfactors) {
Eigen::Matrix2cd ff = ffw.evaluatePol(sim_element); Eigen::Matrix2cd ff = ffw.evaluatePol(sim_element);
double fraction = ffw.relativeAbundance(); double fraction = ffw.relativeAbundance();
diffuse_matrix += fraction * (ff * polarization_handler.getPolarization() * ff.adjoint()); diffuse_matrix += fraction * (ff * polarization_handler.getPolarization() * ff.adjoint());
...@@ -103,7 +103,7 @@ double SSCApproximationStrategy::polarizedCalculation(const SimulationElement& s ...@@ -103,7 +103,7 @@ double SSCApproximationStrategy::polarizedCalculation(const SimulationElement& s
ff_orig += prefac * ff; ff_orig += prefac * ff;
ff_conj += prefac * ff.adjoint(); ff_conj += prefac * ff.adjoint();
} }
complex_t p2kappa = getCharacteristicSizeCoupling(qp, m_formfactor_wrappers); complex_t p2kappa = getCharacteristicSizeCoupling(qp, m_weighted_formfactors);
complex_t omega = getCharacteristicDistribution(qp, m_iff.get()); complex_t omega = getCharacteristicDistribution(qp, m_iff.get());
Eigen::Matrix2cd interference_matrix = (2.0 * omega / (1.0 - p2kappa * omega)) Eigen::Matrix2cd interference_matrix = (2.0 * omega / (1.0 - p2kappa * omega))
* polarization_handler.getAnalyzerOperator() * ff_orig * polarization_handler.getAnalyzerOperator() * ff_orig
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment