MatMatMult_MPIDense_MPIDense() works currently?

[Barry Smith]

    This mail belongs in petsc-maint at mcs.anl.gov You will need to run  
in the debugger to see what it is crashing or even better run in the  
debugger with valgrind (see valgrind.org).

    Barry

On Sep 18, 2009, at 6:17 PM, Fábio Leite wrote:

> Thanks Hong, I tested this examples (ex103.c and ex107.c) and worked  
> fine !
> Now I want to multiply two dense matrices using BLAS3 routines. I  
> try to use the MatMatMultNumeric_MPIDense_MPIDense() function  
> implemented in src/mat/impls/dense/mpi/mpidense.c but I see this  
> error in console:
>
>
> [...]
> [0]PETSC ERROR:  
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation  
> Violation, probably memory access out of range
> [0]PETSC ERROR: Try option -start_in_debugger or - 
> on_error_attach_debugger
> [0]PETSC ERROR: or see http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#Signal 
> [0]PETSC ERROR: or try http://valgrind.org on linux or man  
> libgmalloc on Apple to find memory corruption errors
> [0]PETSC ERROR: likely location of problem given in stack below
> [0]PETSC ERROR: ---------------------  Stack Frames  
> ------------------------------------
> [0]PETSC ERROR: Note: The EXACT line numbers in the stack are not  
> available,
> [0]PETSC ERROR:       INSTEAD the line number of the start of the  
> function
> [0]PETSC ERROR:       is given.
> [0]PETSC ERROR: [0] MatMPIDenseCopyToPlapack line 1028 /home/hpcin/ 
> soft/petsc-3.0.0-p8/src/mat/impls/dense/mpi/mpidense.c
> [0]PETSC ERROR: [0] MatMatMultNumeric_MPIDense_MPIDense line 1078 / 
> home/hpcin/soft/petsc-3.0.0-p8/src/mat/impls/dense/mpi/mpidense.c
> [0]PETSC ERROR: --------------------- Error Message  
> ------------------------------------
> [0]PETSC ERROR: Signal received!
> [0]PETSC ERROR:  
> ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Release Version 3.0.0, Patch 8, Fri Aug 21  
> 14:02:12 CDT 2009
> [0]PETSC ERROR: See docs/changes/index.html for recent updates.
> [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> [0]PETSC ERROR: See docs/index.html for manual pages.
> [0]PETSC ERROR:  
> ------------------------------------------------------------------------
> [0]PETSC ERROR: ./mult on a linux-gnu named hpcin08 by hpcin Fri Sep  
> 18 19:14:55 2009
> [0]PETSC ERROR: Libraries linked from /home/hpcin/soft/petsc-3.0.0- 
> p8/linux-gnu-c-debug/lib
> [0]PETSC ERROR: Configure run at Fri Sep 18 16:03:03 2009
> [0]PETSC ERROR: Configure options --download-f-blas-lapack=1 -- 
> download-plapack --with-mpi-dir=/usr/local/bin/mpich2-1.1.1p1 --with- 
> scalar-type=real --with-precision=d[1]
> [...]
>
>
>
> Can I really use this function to multiply two dense matrices ? If  
> not, how can I multiply these matrices using BLAS3 routines ?
> I am attaching my code if you want to see how I am doing.
>
> Regards.
>
>
>
> --
> Fábio Leite Soares
> Undergraduate Student of Computing Engineering
> Centro de Informática - UFPE - BRAZIL
> <matrix.txt><mult.c>