cctbx/eltbx/xray_scattering/form_factor.h

#ifndef CCTBX_ELTBX_XRAY_SCATTERING_FORM_FACTOR_H
#define CCTBX_ELTBX_XRAY_SCATTERING_FORM_FACTOR_H

#include <scitbx/array_family/shared.h>
#include <cctbx/import_scitbx_af.h>

namespace cctbx { namespace eltbx {

namespace xray_scattering {

  template<class Derived>
  class isotropic_form_factor_mixin
  {
    Derived const & heir() const {
      return static_cast<Derived const &>(*this);
    }

    public:
      double
      at_stol_sq(double stol_sq) const
      {
        return heir().at_x_sq(stol_sq);
      }

      double
      at_stol(double stol) const
      {
        return heir().at_x_sq(stol * stol);
      }

      double
      at_d_star_sq(double d_star_sq) const
      {
        return heir().at_x_sq(d_star_sq / 4);
      }

      af::shared<double>
      at_d_star_sq(af::const_ref<double> const& d_star_sq) const
      {
        af::shared<double>
          result(d_star_sq.size(), af::init_functor_null<double>());
        for(std::size_t i=0;i<d_star_sq.size();i++) {
          result[i] = at_d_star_sq(d_star_sq[i]);
        }
        return result;
      }

      double
      at_d_star(double d_star) const
      {
        return heir().at_x_sq(d_star * d_star / 4);
      }
  };

}}} // cctbx::eltbx::xray_scattering

#endif // GUARD

---