Matrix input file format

Barry Smith bsmith at mcs.anl.gov
Wed Jun 7 11:35:36 CDT 2006


http://www-unix.mcs.anl.gov/web-mail-archive/lists/petsc-users/2006/06/msg00012.html

which can be reached from
http://www-unix.mcs.anl.gov/web-mail-archive/lists/petsc-users/2006/06/threads.html

   Barry


On Wed, 7 Jun 2006, Evrim Dizemen wrote:

> Hi all,
>
> First thanks for your comments. I still have the same problem but i agree 
> with you all that reading a matrix from a file is not parallel programming. 
> So i wonder if you can recommend an efficient and faster way to read a 
> 30000x30000 sparse matrix where i got those values from my fortran code on 
> aeroacoustics similation. I tried calling MatSetValues as i attached below. I 
> found that way from the ex11f.F in the matrix examples and not sure if it is 
> appropriate for my situation.
>
> !      do, i=1,n
> !       do, j=1,n
> !        value=Aval(i,j)
> !       if(value .ne. 0)then
> !       l=i-1
> !       m=j-1
> !        call MatSetValues(A,1,l,1,m,value,INSERT_VALUES,ierr)
> !       endif
> !       enddo
> !       enddo
>
> Thanks soooooo much
>
> EVRIM
>
>
>
>




More information about the petsc-users mailing list