← Back to branch summary

~chaffra/+junk/trilinos

« back to all changes in this revision

Viewing changes to packages/Sundance/tests-std-framework/Problem/HighOrderPoisson1D.cpp

  • Committer: Bazaar Package Importer
  • Author(s): Christophe Prud'homme, Christophe Prud'homme, Johannes Ring
  • Date: 2009-12-13 12:53:22 UTC
  • mfrom: (5.1.2 sid)
  • Revision ID: james.westby@ubuntu.com-20091213125322-in0nrdjc55deqsw9
Tags: 10.0.3.dfsg-1
[Christophe Prud'homme]
* New upstream release

[Johannes Ring]
* debian/patches/libname.patch: Add prefix 'libtrilinos_' to all
  libraries. 
* debian/patches/soname.patch: Add soversion to libraries.
* debian/watch: Update download URL.
* debian/control:
  - Remove python-numeric from Build-Depends (virtual package).
  - Remove automake and autotools from Build-Depends and add cmake to
    reflect switch to CMake.
  - Add python-support to Build-Depends.
* debian/rules: 
  - Cleanup and updates for switch to CMake.

Show diffs side-by-side

added added

removed removed

Lines of Context:
packages/Sundance/tests-std-framework/Problem/HighOrderPoisson1D.cpp
44
44
  try
45
45
                {
46
46
      Sundance::init(&argc, &argv);
47
 
      int np = MPIComm::world().getNProc();
48
 
 
49
 
      //   cerr << "pid = " << getpid() << endl;
50
 
      //int zippy = 1;
51
 
      //while (zippy) ;
52
47
 
53
48
      /* We will do our linear algebra using Epetra */
54
49
      VectorType<double> vecType = new EpetraVectorType();
55
50
 
56
51
      /* Create a mesh. It will be of type BasisSimplicialMesh, and will
57
52
       * be built using a PartitionedLineMesher. */
58
 
      int nx = 10;
 
53
      int nx = 32;
59
54
      MeshType meshType = new BasicSimplicialMeshType();
60
55
      MeshSource mesher = new PartitionedLineMesher(0.0, 1.0, nx, meshType);
61
56
      Mesh mesh = mesher.getMesh();
97
92
      cerr << "matrix = " << endl << prob.getOperator() << endl;
98
93
      cerr << "rhs = " << endl << prob.getRHS() << endl;
99
94
 
 
95
 
 
96
#ifdef HAVE_CONFIG_H
100
97
      ParameterXMLFileReader reader(searchForFile("SolverParameters/bicgstab.xml"));
 
98
#else
 
99
      ParameterXMLFileReader reader("bicgstab.xml");
 
100
#endif
101
101
      ParameterList solverParams = reader.getParameters();
102
102
      cerr << "params = " << solverParams << endl;
103
103
 
109
109
      Expr soln = prob.solve(solver);
110
110
 
111
111
 
112
 
      TEST_FOR_EXCEPTION(!(prob.solveStatus().finalState() == SolveConverged),
113
 
                         RuntimeError,
114
 
                         "solve failed");
115
 
 
116
112
      Expr errExpr = Integral(interior, 
117
113
                              pow(soln-exactSoln, 2),
118
114
                              quad);