trilinos / Trilinos

Primary repository for the Trilinos Project
https://trilinos.org/
Other
1.21k stars 564 forks source link

Amesos-Pardiso: input inconsistent error #12965

Open ewu9 opened 6 months ago

ewu9 commented 6 months ago

Bug Report

@srajama1

I've compiled Trilinos with Intel MKL support and trying to use the Pardiso solver in Xyce. I am trying to run the Pardiso solver with mpi and threads. Both OpenMPI and OpenMP are installed. I configured trilinos and xyce using the link line advisor.

When I try to solve my circuit with the pardiso solver, an error occurs during the symbolic factorization step in Pardiso Error message: Input inconsistent AMESOS ERROR -1

I've tried several input matrices for the solver, all which are simple circuits and are solved with other Trilinos solvers like KLU, Aztecoo. None are solvable with Pardiso.

Steps to Reproduce

  1. Configure script: attached both Trilinos and Xyce autoconf scripts (attachment 1 and 2)
  2. Intel specific environment variables set
  3. OMP_NUM_THREADS=2
  4. call inside the spice file. .options linsol tr_singleton_filter=1 tr_amd=0 type=pardiso tr_partition=1 (attachment 3)
  5. function call: mpirun -np 4 Xyce -hspice-ext all (output in attachment 4)

Any help would be greatly appreciated.

files changed to .txt extension so they can be uploaded reconfigure-trilinos.txt reconfigure-xyce.txt or2_3v3_tb.spice.txt xyce_terminal_output.txt

ccober6 commented 6 months ago

@hkthorn

hkthorn commented 6 months ago

@ewu9 That code often indicates that Pardiso thinks the matrix structure is singular. In practice, neither distributed-memory nor shared-memory parallelism will provide any speedup for a linear system that has been reduced to a 3x3 through singleton filtering. The circuit and resulting linear systems need to be substantially larger before parallelism provides any benefit.