dolfin team mailing list archive

["Johannes Ring" <johannr@xxxxxxxxx>]

On Tue, April 29, 2008 DOLFIN wrote:
> One or more new changesets pushed to the primary dolfin repository.
> A short summary of the last three changesets is included below.
>
> changeset:   4117:6e2f799b6fe1a300e804f6253021f1b4c578009d
> tag:         tip
> user:        Kent-Andre Mardal <kent-and@xxxxxxxxx>
> date:        Tue Apr 29 09:36:46 2008 +0200
> files:       demo/la/trilinos/python/demo.py
> description:
> clean up by Ilmar

This update broke the demo on my computer:

johannr@simula-x61:python$ pwd
/home/johannr/src/hg/dolfin/demo/la/trilinos/python
johannr@simula-x61:python$ python demo.py
[simula-x61:13319] mca: base: component_find: unable to open osc pt2pt:
file not found (ignored)
libibverbs: Fatal: couldn't read uverbs ABI version.
--------------------------------------------------------------------------
[0,0,0]: OpenIB on host simula-x61 was unable to find any HCAs.
Another transport will be used instead, although this may result in
lower performance.
--------------------------------------------------------------------------
Assembling matrix over cells (finished).
Assembling vector over cells (finished).
Assembling vector over exterior facets (finished).
Creating sub domain markers for boundary condition.
Computing Dirichlet boundary values (finished).
Applying boundary conditions to linear system.

*********************************************************
ML failed to compute the multigrid preconditioner. The
most common problem is an incorrect  data type in ML's
parameter list (e.g. 'int' instead of 'bool').

Note: List.set("ML print parameter list",true) might help
figure out the bad one.
*********************************************************

ML::ERROR:: -1, ./Utils/ml_MultiLevelPreconditioner.cpp, line 1797

		*******************************************************
		***** Problem: Epetra::CrsMatrix
		***** Preconditioned CG solution
		***** User-defined preconditioner
		***** No scaling
		*******************************************************

ML::ERROR:: -10, ./Utils/ml_MultiLevelPreconditioner.cpp, line 2148

Error in Epetra Object with label:  Epetra_Aztec_precond X Vector
Epetra Error:  Error in call to Epetra_Operator for preconditioner  Error
Code:  -10
terminate called after throwing an instance of 'int'
[simula-x61:13319] *** Process received signal ***
[simula-x61:13319] Signal: Aborted (6)
[simula-x61:13319] Signal code:  (-6)
[simula-x61:13319] [ 0] [0xb7fb7440]
[simula-x61:13319] [ 1] /lib/tls/i686/cmov/libc.so.6(abort+0x101)
[0xb7e32a01]
[simula-x61:13319] [ 2]
/usr/lib/libstdc++.so.6(_ZN9__gnu_cxx27__verbose_terminate_handlerEv+0x150)
[0xb5872480]
[simula-x61:13319] [ 3] /usr/lib/libstdc++.so.6 [0xb586fd05]
[simula-x61:13319] [ 4] /usr/lib/libstdc++.so.6 [0xb586fd42]
[simula-x61:13319] [ 5] /usr/lib/libstdc++.so.6 [0xb586fe6a]
[simula-x61:13319] [ 6]
/home/johannr/src/trilinos-8.0.3/local/lib/libaztecoo.so(Epetra_Aztec_precond+0x1fe)
[0xb5f18e4e]
[simula-x61:13319] [ 7]
/home/johannr/src/trilinos-8.0.3/local/lib/libaztecoo.so(AZ_pcg_f+0xd8e)
[0xb5f2d36e]
[simula-x61:13319] [ 8]
/home/johannr/src/trilinos-8.0.3/local/lib/libaztecoo.so(AZ_oldsolve+0x500)
[0xb5f53cc0]
[simula-x61:13319] [ 9]
/home/johannr/src/trilinos-8.0.3/local/lib/libaztecoo.so(AZ_iterate+0x10e)
[0xb5f548ce]
[simula-x61:13319] [10]
/home/johannr/src/trilinos-8.0.3/local/lib/libaztecoo.so(_ZN7AztecOO7IterateEid+0xcb)
[0xb5f16ccb]
[simula-x61:13319] [11]
/home/johannr/src/trilinos-8.0.3/local/lib/python2.5/site-packages/PyTrilinos/_AztecOO.so
[0xb2978d4d]
[simula-x61:13319] [12] python(PyObject_Call+0x27) [0x805cb37]
[simula-x61:13319] [13] python(PyEval_EvalFrameEx+0x4064) [0x80c7ce4]
[simula-x61:13319] [14] python(PyEval_EvalCodeEx+0x6e7) [0x80cb0d7]
[simula-x61:13319] [15] python(PyEval_EvalFrameEx+0x565e) [0x80c92de]
[simula-x61:13319] [16] python(PyEval_EvalCodeEx+0x6e7) [0x80cb0d7]
[simula-x61:13319] [17] python(PyEval_EvalCode+0x57) [0x80cb227]
[simula-x61:13319] [18] python(PyRun_FileExFlags+0xf8) [0x80ea6d8]
[simula-x61:13319] [19] python(PyRun_SimpleFileExFlags+0x199) [0x80ea979]
[simula-x61:13319] [20] python(Py_Main+0xa35) [0x8059335]
[simula-x61:13319] [21] python(main+0x22) [0x80587f2]
[simula-x61:13319] [22]
/lib/tls/i686/cmov/libc.so.6(__libc_start_main+0xe0) [0xb7e1c450]
[simula-x61:13319] [23] python [0x8058761]
[simula-x61:13319] *** End of error message ***
Aborted
johannr@simula-x61:python$


The fix was as simple as importing the DOLFIN stuff after the Trilinos
stuff, that is, calling

  from dolfin import *

after

  from PyTrilinos import Epetra, AztecOO, TriUtils, ML

(patch attached)

Johannes


>
> changeset:   4116:e860012bf5050538f42ed023597710c28bc1856e
> user:        Kent-Andre Mardal <kent-and@xxxxxxxxx>
> date:        Mon Apr 28 22:19:04 2008 +0200
> files:       dolfin/la/EpetraSparsityPattern.cpp
> dolfin/la/EpetraVector.cpp sandbox/la/poisson/BlockLinearAlgebra.py
> description:
> merge + updated EpetraSparsityPattern for rectangular matrices
>
>
> changeset:   4115:8b8c80ebc18e678814cd16f44d8e451af3a1ba77
> user:        "Martin Sandve Alnæs <martinal@xxxxxxxxx>"
> date:        Mon Apr 28 19:56:39 2008 +0200
> files:       dolfin/la/uBlasVector.cpp
> description:
> Removed printouts in uBlasVector.
>
> ----------------------------------------------------------------------
> For more details, visit http://www.fenics.org/hg/dolfin
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev
>

diff -r 6e2f799b6fe1 demo/la/trilinos/python/demo.py
--- a/demo/la/trilinos/python/demo.py	Tue Apr 29 09:36:46 2008 +0200
+++ b/demo/la/trilinos/python/demo.py	Tue Apr 29 10:11:13 2008 +0200
@@ -8,6 +8,13 @@ __date__ = "2008-04-24"
 __date__ = "2008-04-24"
 __copyright__ = "Copyright (C) 2008 Kent-Andre Mardal"
 
+try:
+    from PyTrilinos import Epetra, AztecOO, TriUtils, ML 
+except:
+    print "You Need to have PyTrilinos with Epetra, AztecOO, TriUtils and ML installed for this demo to run",
+    print "Exiting."
+    exit()
+
 from dolfin import *
 
 # Test for Trilinos:
@@ -15,13 +22,6 @@ try:
     dolfin.EpetraMatrix
 except:
     print "PyDOLFIN has not been configured with Trilinos. Exiting."
-    exit()
-
-try:
-    from PyTrilinos import Epetra, AztecOO, TriUtils, ML 
-except:
-    print "You Need to have PyTrilinos with Epetra, AztecOO, TriUtils and ML installed for this demo to run",
-    print "Exiting."
     exit()
 
 # Create mesh and finite element