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 }
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 }