Dear all,
I am in the process of upgrading from 3.2 to 3.3.
I am a little bit puzzled by the following change:
> VecSetBlockSize() cannot be called after VecCreateSeq() or
> VecCreateMPI() and must be called before VecSetUp() or
> VecSetFromOptions() or before either VecSetType() or VecSetSizes()
With the earlier release I used to do the following:
CALL VecCreateSeq(PETSC_COMM_SELF,NPOIN*NOFVAR,DT,IFAIL)
C
C
IF(NOFVAR.GT.1) CALL VecSetBlockSize(DT,NOFVAR,IFAIL)
with 3.3 it looks like the following is required :
CALL VecCreate(PETSC_COMM_SELF,DT,IFAIL)
CALL VecSetType(DT,VECSEQ,IFAIL)
CALL VecSetBlockSize(DT,NOFVAR,IFAIL)
CALL VecSetSizes(DT,NPOIN*NOFVAR,PETSC_DECIDE,IFAIL)
Is there a simpler (i.e. less library calls) way to achieve the same result?
Regards,
Aldo
--
Dr. Aldo Bonfiglioli
Associate professor of Fluid Flow Machinery
Scuola di Ingegneria
Universita' della Basilicata
V.le dell'Ateneo lucano, 10 85100 Potenza ITALY
tel:+39.0971.205203 fax:+39.0971.205215
Publications list <http://publicationslist.org/aldo.bonfiglioli>
