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