[petsc-users] How to debug problem with MatSetValues

Jed Brown wrote: On Mon, 05 Jul 2010 15:22:50 +0200, Thomas Witkowski thomas.witkowski at tu-dresden.de wrote: Hi, I've some trouble with matrix values that are set by MatSetValues, but are missing in the final matrix. I can reduce the problem to a 75x75 matrix that is created on

[petsc-users] Slepc - SVD routines

Hi, First off I apologize for the slepc question, they don't have a users lists so I'm hoping someone on here might be able to offer some guidance. Problem: I'm working on a graph partitioning problem. Basically I have the laplacian of a large graph and am interested in extracting its

[petsc-users] Slepc - SVD routines

On 07/07/2010, Luke Bloy wrote: Hi, First off I apologize for the slepc question, they don't have a users lists so I'm hoping someone on here might be able to offer some guidance. Problem: I'm working on a graph partitioning problem. Basically I have the laplacian of a large graph

[petsc-users] my code runs too slow

total: nonzeros=1830 mallocs used during MatSetValues calls =1830 Looks like you are zero-ing out the non-zero structure - before assembling the matrix. Are you calling MatZeroRows() or MatZeroEntries() or something else - before assembling the matrix? Satish On Wed, 7 Jul 2010, Xuan YU

[petsc-users] my code runs too slow

at psu.edu Xuan YU (??) xxy113 at psu.edu -- next part -- An HTML attachment was scrubbed... URL: http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100707/7cbd2b82/attachment-0001.htm

[petsc-users] my code runs too slow

On Wed, 7 Jul 2010, Xuan YU wrote: ierr = MatCreateSeqAIJ(PETSC_COMM_SELF,N,N,10,PETSC_NULL,J);CHKERRQ(ierr); ierr = MatAssemblyBegin(J,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr); ierr = MatAssemblyEnd(J,MAT_FINAL_ASSEMBLY);CHKERRQ(ierr); This assembly removes the unused space here. Since no values

[petsc-users] my code runs too slow

/petsc-users/attachments/20100707/b68eb0b5/attachment-0001.htm

[petsc-users] my code runs too slow

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/20100707/ac7b65c1/attachment.htm

[petsc-users] my code runs too slow

ierr = SNESComputeJacobian(ts_snes,CV_Y,J,J,flag);CHKERRQ(ierr); Perhaps there is a function thats set to compute jacobinan thats not assembling the matrix properly? Also - please dont' truncate error messages - when you send them. Incomplete info is not useful. Satish On Wed, 7 Jul 2010,

[petsc-users] my code runs too slow

(??) xxy113 at psu.edu Xuan YU (??) xxy113 at psu.edu Xuan YU (??) xxy113 at psu.edu Xuan YU (??) xxy113 at psu.edu -- next part -- An HTML attachment was scrubbed... URL: http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100707/aec3bc29

[petsc-users] my code runs too slow

What function are you using to compute jacobian? Perhaps you are missing MatAssembly() calls in that function? Or reordering the code as follows should work - but the correct location for MatAssembly() calls are in the routine you are using for assembling the jacobian. Satish ierr =

[petsc-users] my code runs too slow

/pipermail/petsc-users/attachments/20100707/e44e44a4/attachment-0001.htm

[petsc-users] my code runs too slow

Xuan YU (??) xxy113 at psu.edu Xuan YU ( ?? ) xxy113 at psu.edu -- next part -- An HTML attachment was scrubbed... URL: http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100707

[petsc-users] my code runs too slow

Xuan YU (??) xxy113 at psu.edu Xuan YU ( ?? ) xxy113 at psu.edu -- next part -- An HTML attachment was scrubbed... URL: http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100707/19a1f9cd/attachment