[petsc-users] Changing the diagonal of a matrix via a vector

2011-04-27 Thread Jed Brown
/lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110427/742187be/attachment.htm>

[petsc-users] Changing the diagonal of a matrix via a vector

2011-04-27 Thread Jed Brown
x27;t propagate up correctly). http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/Sys/CHKERRQ.html -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110427/1a638b5f/attachment.htm>

[petsc-users] create a new large vector via combining existing small vectors

2011-04-27 Thread Xiangdong Liang
with Petsc's built in functions? Thanks. Best, Xiangdong -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110427/bf0be569/attachment.htm>

[petsc-users] create a new large vector via combining existing small vectors

2011-04-27 Thread Barry Smith
On Apr 27, 2011, at 3:35 PM, Xiangdong Liang wrote: > Hello everyone, > > I have a problem with creating a new large vector via combining existing > small vectors. Suppose I have two vectors v1 and v2 (size n-by-1) already. I > want to have a new vector vout (size 2n-by-1) with vout(1:n) = v1

[petsc-users] Changing the diagonal of a matrix via a vector

2011-04-27 Thread Xiangdong Liang
cCreate and MatCreate. they are supposed to return PetscErrorCode. However, I still get "warning: return makes pointer from integer without a cast" for these built-in functions. > > http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/Sys/CHKERRQ.html > -- next part -- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110427/595b387c/attachment.htm>

[petsc-users] Changing the diagonal of a matrix via a vector

2011-04-27 Thread Xiangdong Liang
permail/petsc-users/attachments/20110427/b1c63a14/attachment.htm>

[petsc-users] KSP solver increases the solution time

2011-04-27 Thread Jose E. Roman
El 27/04/2011, a las 09:17, Alejandro Marcos Arag?n escribi?: > I understand that, but I'm trying to provide default behavior for the solver > because the default one (no parameters) works very bad in my case. > However, I'm stuck because I can't set the same parameters that I obtain with > com

[petsc-users] KSP solver increases the solution time

2011-04-27 Thread Alejandro Marcos Aragón
odes, limit used is 5 >>> linear system matrix = precond matrix: >>> Matrix Object: >>> type=seqaij, rows=2020, cols=2020 >>> total: nonzeros=119396, allocated nonzeros=163620 >>> using I-node routines: found 676 nodes, limit used is 5 >>> linear system matrix = precond matrix: >>> Matrix Object: >>> type=mpiaij, rows=6058, cols=6058 >>> total: nonzeros=365026, allocated nonzeros=509941 >>> using I-node (on process 0) routines: found 676 nodes, limit used is 5 >>> >>> [1] System solved in 3664 iterations... 42.683 s >>> >>> >>> >>> As you can see, the second iteration takes more than 40 seconds to solve. >>> Could some explain why this is happening and why he number of iterations is >>> increasing dramatically between solves? Thank you all, >>> >>> Alejandro M. Arag?n, Ph.D. >> > > > > > -- > 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/20110427/a143dc5f/attachment-0001.htm>