dolfin team mailing list archive

[Johan Hake <johan.hake@xxxxxxxxx>]

This error should have been picked up by a buildbot as:

  demo/undomcumented/trilinos/python

also fails. I think we need to sort out the PyTrilinos linking bug first as 
Garth suggested and then look for other sources.

Johan

On Friday November 18 2011 07:20:23 Martin Sandve Alnæs wrote:
> The problem is that dolfin.cpp.EpetraMatrix.mat() returns a
> SwigPyObject instead of a PyTrilinos.Epetra.FECrsMatrix.
> 
> With your cbc.block patch, I'm running the following code:
> 
>         print type(Puu)
> 
>         A = dolfin.down_cast(Puu)
>         print type(A)
> 
>         B = A.mat()
>         print type(B)
> 
>         C = PyTrilinos.ML.MultiLevelPreconditioner(B, 0)
>         print type(C)
> 
> With dolfin/1.0.x -r6439 I get:
> 
> <class 'dolfin.cpp.Matrix'>
> <class 'dolfin.cpp.EpetraMatrix'>
> <class 'PyTrilinos.Epetra.FECrsMatrix'>
> <class 'PyTrilinos.ML.MultiLevelPreconditioner'>
> 
> With dolfin/1.0.x -r6447 I get:
> 
> <class 'dolfin.cpp.Matrix'>
> <class 'dolfin.cpp.EpetraMatrix'>
> <type 'SwigPyObject'>
> Traceback (most recent call last):
>   File "/home/martinal/dev/cbc/invbench/scripts/invbench", line 5, in
> <module> simfra.introrun.run(locals())
>   File "/home/martinal/dev/cbc/invbench/src/simfra/introrun.py", line 79,
> in run return func(**kwargs)
>   File "/home/martinal/dev/cbc/invbench/src/invbench/oneshot.py", line
> 160, in case_fruitfly_oneshotblock
>     C = PyTrilinos.ML.MultiLevelPreconditioner(B, 0)
>   File
> "/home/martinal/opt/fenics/dorsal-unstable/lib/python2.7/site-packages/PyT
> rilinos/ML.py", line 1442, in __init__
>     this = _ML.new_MultiLevelPreconditioner(*args)
> NotImplementedError: Wrong number or type of arguments for overloaded
> function 'new_MultiLevelPreconditioner'.
>   Possible C/C++ prototypes are:
>    
> ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_RowMa
> trix const &,bool const)
>    
> ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_RowMa
> trix const &,Teuchos::ParameterList const &,bool const)
>    
> ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(ML_Operator
> *,Teuchos::ParameterList const &,bool const)
>    
> ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(ML_Operator
> *,Teuchos::ParameterList const &,Epetra_RowMatrix
> **,Teuchos::ParameterList *,int,bool const)
>    
> ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_RowMa
> trix const &,Epetra_RowMatrix const &,Epetra_RowMatrix const
> &,Teuchos::ParameterList const &,bool const,bool const)
>    
> ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_RowMa
> trix const &,Epetra_RowMatrix const &,Epetra_RowMatrix const
> &,Epetra_RowMatrix const &,Teuchos::ParameterList const &,bool const)
>    
> ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_MsrMa
> trix const &,ML_Operator *,AZ_MATRIX *,int *,Teuchos::ParameterList const
> &,bool const)
> 
> 
> Martin
> 
> On 18 November 2011 16:03, Joachim Berdal Haga <jobh@xxxxxxxxx> wrote:
> >>   self.ml_prec = MultiLevelPreconditioner(A.down_cast().mat(), 0)
> > 
> > I assume Puu is a dolfin.GenericMatrix, so this should be equivalent to
> >   PyTrilinos.ML.MultiLevelPreconditioner(dolfin.down_cast(Puu).mat(), 0)
> > and the only reason I can see for that to fail is if the returned matrix
> > from dolfin isn't (recognised as) an Epetra matrix.
> > Does the explicit call above give the same error? If so, what is the type
> > of the down_cast(Puu).mat()?
> > -j.
> > 
> > On 18 November 2011 15:47, Martin Sandve Alnæs <martinal@xxxxxxxxx> wrote:
> >> On 18 November 2011 13:52, Joachim Berdal Haga <jobh@xxxxxxxxx> wrote:
> >> > Yep, sorry, forgot the list. Here's a patch for cbc.block which
> >> > removes the
> >> > workaround, can you see if it helps?
> >> > -j.
> >> 
> >> Now I get:
> >> 
> >> Traceback (most recent call last):
> >>  File "/home/martinal/dev/cbc/invbench/scripts/invbench", line 5, in
> >> <module>
> >>    simfra.introrun.run(locals())
> >>  File "/home/martinal/dev/cbc/invbench/src/simfra/introrun.py", line 79,
> >> in run
> >>    return func(**kwargs)
> >>  File "/home/martinal/dev/cbc/invbench/src/invbench/oneshot.py", line
> >> 168, in case_fruitfly_oneshotblock
> >>    P = block_mat([[ML(Puu), 0, 0],
> >>  File
> >> "/home/martinal/opt/fenics/dorsal-unstable/lib/python2.7/site-packages/b
> >> lock/algebraic/trilinos/MLPrec.py", line 25, in __init__
> >>    self.ml_prec = MultiLevelPreconditioner(A.down_cast().mat(), 0)
> >>  File
> >> "/home/martinal/opt/fenics/dorsal-unstable/lib/python2.7/site-packages/P
> >> yTrilinos/ML.py", line 1442, in __init__
> >>    this = _ML.new_MultiLevelPreconditioner(*args)
> >> NotImplementedError: Wrong number or type of arguments for overloaded
> >> function 'new_MultiLevelPreconditioner'.
> >>  Possible C/C++ prototypes are:
> >> 
> >>  ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_Ro
> >> wMatrix const &,bool const)
> >> 
> >>  ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_Ro
> >> wMatrix const &,Teuchos::ParameterList const &,bool const)
> >> 
> >>  ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(ML_Operat
> >> or *,Teuchos::ParameterList const &,bool const)
> >> 
> >>  ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(ML_Operat
> >> or *,Teuchos::ParameterList const &,Epetra_RowMatrix
> >> **,Teuchos::ParameterList *,int,bool const)
> >> 
> >>  ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_Ro
> >> wMatrix const &,Epetra_RowMatrix const &,Epetra_RowMatrix const
> >> &,Teuchos::ParameterList const &,bool const,bool const)
> >> 
> >>  ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_Ro
> >> wMatrix const &,Epetra_RowMatrix const &,Epetra_RowMatrix const
> >> &,Epetra_RowMatrix const &,Teuchos::ParameterList const &,bool const)
> >> 
> >>  ML_Epetra::MultiLevelPreconditioner::MultiLevelPreconditioner(Epetra_Ms
> >> rMatrix const &,ML_Operator *,AZ_MATRIX *,int *,Teuchos::ParameterList
> >> const &,bool const)
> >> 
> >> 
> >> --
> >> Martin
> 
> _______________________________________________
> Mailing list: https://launchpad.net/~dolfin
> Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin
> More help   : https://help.launchpad.net/ListHelp