VecGetValue?

2006-02-22 Thread Matthew Knepley
billy at dem.uminho.pt writes: > Hi, > > Why does PETSc have a VecSetValue and not a VecGetValue? We could introduce one, but its not usually what you want since it would require a collective Assemble() call before you could use the value. You can accomplish the same thing with a scatter.

Matrix operation

2006-02-22 Thread Matthew Knepley
billy at dem.uminho.pt writes: > Thank you very much for your answers. > > I was thinking why doesn't PETSc allow functions inside functions (like > LASPack)? > > h += (Diag(A) - A) * u <=> AddAsgn_VV (h, Mul_QV (Sub_QQ (Diag_Q (A), A), u)); Its an error reporting discipline since C does not a

New changes for CG set type in 2.3.1-p7 ?

2006-02-22 Thread Satish Balay
On Wed, 22 Feb 2006, abdul-rahman at tu-harburg.de wrote: > Satish, > > On Tue, 21 Feb 2006, Satish Balay wrote: > > > petsc-2.3.0 changelog has the following entry: > > > > >>> > > -ksp_cg_Hermitian and -ksp_cg_symmetric have been changed > > to -ksp_cg_type Hermitian or symmetric > > oops, th

New changes for CG set type in 2.3.1-p7 ?

2006-02-22 Thread abdul-rah...@tu-harburg.de
Satish, On Tue, 21 Feb 2006, Satish Balay wrote: > petsc-2.3.0 changelog has the following entry: > > >>> > -ksp_cg_Hermitian and -ksp_cg_symmetric have been changed > to -ksp_cg_type Hermitian or symmetric oops, thanks for pointing it out. It works as above. However, I still have problem with t

Matrix operation

2006-02-22 Thread bi...@dem.uminho.pt
Thank you very much for your answers. I was thinking why doesn't PETSc allow functions inside functions (like LASPack)? h += (Diag(A) - A) * u <=> AddAsgn_VV (h, Mul_QV (Sub_QQ (Diag_Q (A), A), u)); Billy Quoting Matthew Knepley : > Hong Zhang writes: > > I would recommend the later