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 > > > >
