dolfin team mailing list archive

[Johan Hake <706406@xxxxxxxxxxxxxxxxxx>]

My comment relates to bug: #70540, not this one...

-- 
You received this bug notification because you are a member of DOLFIN
Team, which is subscribed to DOLFIN.
https://bugs.launchpad.net/bugs/706406

Title:
  Failure creating PyTrilinos object when imported after dolfin

Status in DOLFIN:
  New

Bug description:
  The following test program fails unless the import order is switched.
  The raw_input at the end is to confirm that the crash happens before
  exit.

  -------
  from dolfin import *
  from PyTrilinos import ML
  dolfin.parameters["linear_algebra_backend"] = "Epetra"
  mesh = UnitSquare(8,8)
  V = FunctionSpace(mesh, "CG", 1)
  u,v = TrialFunction(V), TestFunction(V)
  a = u*v*dx
  A = assemble(a)
  MLList = { "aggregation: type"         : "Uncoupled"  }
  ml_prec = ML.MultiLevelPreconditioner(down_cast(A).mat(), 0)
  ml_prec.SetParameterList(MLList)
  ml_agg = ml_prec.GetML_Aggregate()
  ml_prec.ComputePreconditioner()
  raw_input('Finished, press enter to continue')
  ----------

  See bug 705401 (https://bugs.launchpad.net/dolfin/+bug/705401) which
  is about another failure when PyTrilinos is imported after DOLFIN.

  It works fine when PyTrilinos is imported before DOLFIN, so the
  workaround is simple. The problem is that the reported error messages
  do not point obviously to the solution. Some suggestions for dealing
  with the problem are given in the original bug report.