MatView_SeqAIJ_Draw_Zoom() in src/mat/impls/aij/seq/aij.c
To find any function, I would use the debugger.
Matt
On Jan 23, 2008 8:38 PM, Yujie wrote:
> To generate the sparsity picture, could you tell me how to find the code
> used by "-mat_view_draw"? I have searched the whole PETSc package.
HYPRE_PARCSR,IFAIL)
> ! write(6,*)'HYPRE_IJSetObjectType has returned ifail = ',
> !&ifail
> CALL HYPRE_IJVectorInitialize(ij,IFAIL)
> ! write(6,*)'HYPRE_IJVectorInitialize has returned ifail = ',
> !&ifail
> CALL VecGetLocalSize(x,n,IFAIL)
> if(n.GT.maxcols)STOP 'Must increase maxcols'
> CALL VecGetArray(x,x_array,i_x,IFAIL)
> k = 0
> do i = rstart,rend-1
> k = k + 1
> cols(k) = i
> enddo
> CALL HYPRE_IJVectorSetValues(ij,n,cols,x_array(i_x + 1),IFAIL)
> ! write(6,*)'HYPRE_IJVectorSetValues has returned ifail = ',
> !&ifail
> CALL VecRestoreArray(x,x_array,i_x,IFAIL)
> CALL HYPRE_IJVectorAssemble(ij,IFAIL)
> CALL HYPRE_IJVectorPrint(ij,solfile,IFAIL)
> ! write(6,*)'HYPRE_IJVectorPrint has returned ifail = ',
> !&ifail
> CALL HYPRE_IJVectorDestroy(ij,IFAIL)
> RETURN
> END
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080123/e0e0d7c3/attachment.htm>
iles do not make sense for large matrices. You should use
> > binary
> > > > files.
> > > > >
> > > > > Matt
> > > > >
> > > > > > Regards,
> > > > > > Yujie
> > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > >
> > > > > --
> > > > > What most experimenters take for granted before they begin their
> > > > > experiments is infinitely more interesting than any results to
> which
> > > > > their experiments lead.
> > > > > -- Norbert Wiener
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > >
> > > --
> > > What most experimenters take for granted before they begin their
> > > experiments is infinitely more interesting than any results to which
> > > their experiments lead.
> > > -- Norbert Wiener
> > >
> > >
> >
> >
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which
> their experiments lead.
> -- Norbert Wiener
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080123/6b0b41df/attachment.htm>
gt; >Matt
> > > > >
> > > > > > Regards,
> > > > > > Yujie
> > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > >
> > > > > --
> > > > > What most experimenters take for granted before they begin their
> > > > > experiments is infinitely more interesting than any results to
> which
> > > > > their experiments lead.
> > > > > -- Norbert Wiener
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> > >
> > > --
> > > What most experimenters take for granted before they begin their
> > > experiments is infinitely more interesting than any results to which
> > > their experiments lead.
> > > -- Norbert Wiener
> > >
> > >
> >
> >
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which
> their experiments lead.
> -- Norbert Wiener
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080123/5c7a840b/attachment.htm>
On Jan 23, 2008 2:18 PM, Yujie wrote:
> Thank you for your further explanation. I just want to use this data in
> other packages. I think that ASCII file is likely better. Because I don't
> know the format of the binary file? how to find it?
Look at MatView_SeqAIJ_Binary() in src/mat/impls/aij/se
st experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which
> their experiments lead.
> -- Norbert Wiener
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080123/3a6051ae/attachment.htm>
> > > Yujie
> > > >
> > > >
> > > >
> > >
> > >
> > >
> > > --
> > > What most experimenters take for granted before they begin their
> > > experiments is infinitely more interesting than any results to which
> > > their experiments lead.
> > > -- Norbert Wiener
> > >
> > >
> >
> >
>
>
>
> --
> What most experimenters take for granted before they begin their
> experiments is infinitely more interesting than any results to which
> their experiments lead.
> -- Norbert Wiener
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080123/f68a27c7/attachment.htm>
alidate() line 139 in src/sys/memory/mtr.c
> > > [0]PETSC ERROR: SNESComputeFunction() line 950 in
> src/snes/interface/snes.c
> > > [0]PETSC ERROR: SNESSolve_LS() line 158 in src/snes/impls/ls/ls.c
> > > [0]PETSC ERROR: SNESSolve() line 1871 in src/snes/interface/snes.c
> > > [0]PETSC ERROR: main() line 395 in unknowndirectory/myProg.C
> > >
> > > If you need some more information, feel free to let me know. Thanks in
> advance !
> > >
> > > Vijay
> > >
> > >
> >
> >
> >
> > --
> > What most experimenters take for granted before they begin their
> > experiments is infinitely more interesting than any results to which
> > their experiments lead.
> > -- Norbert Wiener
> >
> >
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080123/558c8b00/attachment.htm>
-unspecified text was scrubbed...
Name: printmat.F
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20080123/8fe41d16/attachment.diff>
On Jan 22, 2008 11:01 PM, Yujie wrote:
> Dear Matt:
>
> thank you for your reply. Do you have any method to generate an ascii file
> of the huge sparse matrix? thanks
I think you miss my point. The PETSc function is not a bad way to generate
ASCII matrices. ASCII matrices make "no sense" for larg
10 matches
Mail list logo
