[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Some basic questions about PETSC
- To: petsc-users@xxxxxxxxxxx
- Subject: Some basic questions about PETSC
- From: "Ben Tay" <zonexo@xxxxxxxxx>
- Date: Mon, 4 Dec 2006 11:41:29 +0800
- Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=beta; d=gmail.com; h=received:message-id:date:from:to:subject:mime-version:content-type; b=VMAeAKmN0IzHN1ak3RXuEBUjXMSjioPVxbAwMPIcUPWEbtv7Lwf1rND25IKI0YB5beIyGorfJPPVI6Ia2hhgAo+fHvs3Xy9qbas+02wkVUGu/2IxAt/JabCNyhjHs1IlIdVVijPXYjOC+PlraMK0MP/TJtMUvWhL378zAMmRlNc=
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
Hi,
I'm still relatively new to PETSC but I've finally managed to solve my own poisson eqn using it in about a week's time.Thank you for writing such a good solver.
However, I have some questions still about PETSC:
1. Are petscscaler the same as real nos.? I am writing in double precision or real(8). So are these 2 interchangeable? (same as petscint and integer).
2. while linking in compaq visual fortran, I get the warning message about conflict library with LIBCMT. Is it important or can I just ignore it?
3. I am writing in f90 free form format. When I modify from fixed f77 format to f90, it gives me a lot of error, which seems to be the include file error. So, is it possible to write in f90 format?
4. I am writing a Navier stokes solver and I'm using PETSC to solve the poisson eqn. If only the RHS changes with each time step, do I need to call
VecCreateSeq
<insert values in RHS b>
VecAssemblyBegin
VecAssemblyEnd
VecDestroy
at each time step? I don't think I can (or I should) change values of b after VecAssemblyBegin,VecAssemblyEnd, or can I?
Thank you very much and have a nice day
Regards,
Ben