43 Teuchos::RCP<Epetra_CrsMatrix> A = Teuchos::rcp(
new Epetra_CrsMatrix(Copy,map,3));
46 double values[] = {-1,2,-1};
52 for(
int i=0;i<4;i++) {
54 ind[0] = map.GID(i-1);
56 ind[2] = map.GID(i+1);
69 A->InsertGlobalValues(gid,entries,valuesPtr,indPtr);
76int main(
int argc,
char *argv[])
79 MPI_Init(&argc, &argv);
80 Epetra_MpiComm Comm(MPI_COMM_WORLD);
82 Epetra_SerialComm Comm;
85 Teuchos::ParameterList GaleriList;
87 Teuchos::RCP<Epetra_Map> Map =
buildMap(Comm);
88 Teuchos::RCP<Epetra_CrsMatrix> Matrix =
buildMatrix(*Map);
91 Epetra_MultiVector x(*Map,NumVectors);
92 Epetra_MultiVector x_exact(*Map,NumVectors);
93 Epetra_MultiVector b(*Map,NumVectors);
95 Matrix->Apply(x_exact,b);
101 Epetra_LinearProblem Problem(&*Matrix, &x, &b);
104 Teuchos::ParameterList List;
105 List.set(
"Reindex",
true);
107 Solver.SetParameters(List);
112 double norm = ComputeNorm(&*Matrix, &x_exact, &b);
113 if (Comm.MyPID() == 0)
114 std::cout <<
"norm = " << norm << std::endl;
123 return(EXIT_SUCCESS);
int main(int argc, char *argv[])
Teuchos::RCP< Epetra_Map > buildMap(Epetra_Comm &Comm)
Teuchos::RCP< Epetra_CrsMatrix > buildMatrix(Epetra_Map &map)