[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: On entry to DGEMV parameter number 6 had an illegal value




KSPSolve_BCGSL() in src/ksp/ksp/impls/bcgsl/bcgsl.c uses calls to gemv() directly.


Could you please try the following. Edit the file and for alll calls to BLAS... and LAPACK... replace the character string by its first character.
For example replace "noTr" with "N". Next before each BLASgemv() put in
a printf("ldMZ %d\n",ldMZ); run make lib shared in that directory
then relink you code and run it. Let us know what happens.


   Barry




On Wed, 14 Jun 2006, Jordi Marcé Nogué wrote:

En/na owner-petsc-users@xxxxxxxxxxx ha escrit:

This an error message from BLAS.

  BLAS has a terrible model for error handling; it just prints
a message to screen that you cannot control and does not
return an error code :-(

   Are you
1) using a BAIJ matrix with a block size bigger than 8?

No. I'm using a SeqAIJ matrix


2) the bcgsl KSP solver or
Yes. And the error happens when I'm calling the "KSPsolve" function

3) dense matrices (how do you create the dense matrix?)

No. I'm using a sparse matrix

4) are you using 64bit integers (very unlikely)

No. I'm using 32bit integers


The best thing to do is run in the debugger to see when/where/why this is happening.

  Barry