Commit 5b93b2ef authored by Doak, Peter W.'s avatar Doak, Peter W.
Browse files

more work on rashba test and inconsistent cplx division

parent 8a0306e2

include/dca/phys/dca_step/cluster_solver/ctaux/ctaux_cluster_solver.hpp

+2 −0
Original line number Diff line number Diff line
@@ -112,6 +112,8 @@ public: 


  void setSampleConfiguration(const io::Buffer&) {}



  /** used for testing */

  auto& getG0() { return g0_; };

protected:

  void warmUp(Walker& walker);

include/dca/phys/dca_step/cluster_solver/ctaux/ctaux_walker.hpp

+0 −1
Original line number Diff line number Diff line
@@ -1658,7 +1658,6 @@ void CtauxWalker<device_t, Parameters, Data>::recomputeMCWeight() { 
      return;



    const auto [log_det, phase] = linalg::matrixop::logDeterminant(m);

    std::cout << "phase: " << phase.getSign() << '\n';

    mc_log_weight_ -= log_det;  // MC weight is proportional to det(N^-1)

    phase_.divide(phase);

  };

include/dca/phys/dca_step/cluster_solver/ctaux/walker/ct_aux_walker_tools.hpp

+2 −0
Original line number Diff line number Diff line
@@ -100,6 +100,8 @@ private: 
  bool test_max_min(int n, dca::linalg::Matrix<Scalar, dca::linalg::CPU>& Gamma_LU, Real max,

                    Real min);



  static Scalar consistentScalarInv(Scalar gamma_k);



private:

  dca::linalg::Vector<Scalar, dca::linalg::CPU> r;

  dca::linalg::Vector<Scalar, dca::linalg::CPU> c;

include/dca/phys/dca_step/cluster_solver/ctaux/walker/tools/g0_interpolation/g0_interpolation_cpu.hpp

+2 −0
Original line number Diff line number Diff line
@@ -63,6 +63,8 @@ public: 
    return 0;

  }



  /// For testing

  auto& getAkimaCoefficients() { return Base::akima_coefficients; }

private:

  auto interpolate(int nu_0, int nu_1, int delta_r, Real delta_time) -> Scalar const;

include/dca/phys/dca_step/cluster_solver/ctaux/walker/tools/g0_interpolation/g0_interpolation_gpu.hpp

+4 −1
Original line number Diff line number Diff line
@@ -64,6 +64,8 @@ public: 
           r1_minus_r0.deviceFingerprint();

  }



  /// for testing

  auto& getAkimaCoefficients() { return akima_coefficients_CPU; }

 private:

  template <class Configuration>

  void uploadConfiguration(const Configuration& configuration);
@@ -156,6 +158,7 @@ void G0Interpolation<dca::linalg::GPU, Parameters>::initialize(MOMS_type& MOMS) 
  G0_r_t_GPU = G0_r_t_CPU;

  grad_G0_r_t_GPU = grad_G0_r_t_CPU;



  // this is the rearrangement for the GPU

  for (int t_ind = 0; t_ind < shifted_t::dmn_size(); t_ind++)

    for (int r_ind = 0; r_ind < r_dmn_t::dmn_size(); r_ind++)

      for (int nu1_ind = 0; nu1_ind < b::dmn_size(); nu1_ind++)