Ejemplo n.º 1
0
void construct_blas_matrix(MAT & M, const operator_container_t &creation_operators,
  const operator_container_t &annihilation_operators, double BETA, const HYB & F)
{

  int N = creation_operators.size();
  M.resize(N, N);
  int row = -1;
  int column = -1;
  for (operator_container_t::iterator ita = annihilation_operators.begin(); ita != annihilation_operators.end(); ita++) {
    row++;
    for (operator_container_t::iterator itc = creation_operators.begin(); itc != creation_operators.end(); itc++) {
      column++;

      double argument = ita->time() - itc->time();
      double sign = 1;
      if (argument < 0) {
        argument += BETA;
        sign = -1;
      }
      M(row, column) = interpolate_F(argument, BETA, F[ita->flavor()][itc->flavor()]) * sign;
    }
    column = -1;
  }
}
Ejemplo n.º 2
0
Archivo: st2010.hpp Proyecto: cmsi/bcl
 static void generate_transition_matrix_resize(VEC const& weights, MAT& tm) {
   std::size_t n = weights.size();
   tm.resize(n);
   for (std::size_t i = 0; i < n; ++i) tm[i].resize(n);
   generate_transition_matrix(weights, tm);
 }