Ejemplo n.º 1
0
double SpinAdapted::dotproduct(const ColumnVector& a, const ColumnVector& b)
{
  assert(a.Nrows() == b.Nrows());
#ifdef BLAS
  return DDOT(a.Storage(), a.Store(), 1, b.Store(), 1);
#else
  return a.t() * b;
#endif
}
Ejemplo n.º 2
0
static void print_vector2(const ColumnVector &xNames) {
#ifdef DEBUG_VC
    int count = xNames.Storage();
    String s = "->";
    for (int i=1; i<=count; i++ ) {
	s += ", ";
	s += m.format(xNames(i),0);
    }
    DEBUG_LOG(s);
#endif
}