Ejemplos de Epetra_MultiVector::Dot en C++ (Cpp)

Ejemplo n.º 1

Archivo: MxGeoMultigridPrec.cpp Proyecto: achellies/maxwell

void MxGeoMultigridPrec::removeConstField(Epetra_MultiVector & x) const {
  //x.Comm().Barrier();
  Epetra_MultiVector ones(x);
  double dotProd, norm;
  ones.PutScalar(1);
  ones.Norm2(&norm);
  x.Dot(ones, &dotProd);
  //std::cout << "norm = " << norm << ", dotProd = " << dotProd << "\n";
  x.Update(-dotProd / norm / norm, ones, 1.);
  x.Dot(ones, &dotProd);
  //std::cout << "const vec part: " << dotProd << "\n";
}

Ejemplo n.º 2

Archivo: MxGeoMultigridPrec.cpp Proyecto: achellies/maxwell

double MxGeoMultigridPrec::constPart(const Epetra_MultiVector & x) const {
  //x.Comm().Barrier();
  Epetra_MultiVector ones(x);
  double dotProd;
  ones.PutScalar(1);
  x.Dot(ones, &dotProd);
  return dotProd;
}