On entry to DGEMV parameter number 6 had an illegal value
Barry Smith
bsmith at mcs.anl.gov
Wed Jun 14 13:14:32 CDT 2006
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 at mcs.anl.gov 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
>>
>>
>
>
More information about the petsc-users
mailing list