← Back to team overview

cbc.block team mailing list archive

[Bug 1063737] [NEW] Crash with explicit matrix multiplies in trilinos

 

Public bug reported:

Code such as

M = collapse(C*B)
prec = ML(M)

crashes with a memory access error.

It seems like EpetraExt.Multiply returns a malformed Epetra matrix when
B and C are rectangular matrices, which then ML chokes on. This used to
work with older Trilinos releases.

I will investigate further as time allows, but for now let's call it a
known problem.

Example code: demo/navierstokes.py

** Affects: cbc.block
     Importance: Undecided
     Assignee: Joachim Haga (jobh)
         Status: New

-- 
You received this bug notification because you are a member of cbc.block
maintainers, which is subscribed to CBC.Block.
Matching subscriptions: cbc.block bugs
https://bugs.launchpad.net/bugs/1063737

Title:
  Crash with explicit matrix multiplies in trilinos

Status in Block operators and algebra for DOLFIN:
  New

Bug description:
  Code such as

  M = collapse(C*B)
  prec = ML(M)

  crashes with a memory access error.

  It seems like EpetraExt.Multiply returns a malformed Epetra matrix
  when B and C are rectangular matrices, which then ML chokes on. This
  used to work with older Trilinos releases.

  I will investigate further as time allows, but for now let's call it a
  known problem.

  Example code: demo/navierstokes.py

To manage notifications about this bug go to:
https://bugs.launchpad.net/cbc.block/+bug/1063737/+subscriptions


Follow ups

References