[petsc-users] Fwd: Re: Question on zero sized DAs

[Shri]

Barry,
      Sorry i forgot to mention that i'm using petsc-3.1 and not petsc-dev. I'll add your change to petsc-3.1.

Thanks,
Shri

----- Barry Smith <bsmith at mcs.anl.gov> wrote:
> 
>   Pull and recompile in src/dm/da/src/ then try running again and report any problems to petsc-maint
   
>    Do you have a stencil width greater than 0?

>  If you are using a stencil width of zero we can try to bypass this error check and see if things go through.

  BTW: this is a petsc-maint question :-) or petsc-users


   Barry
 
> 
> On Sep 8, 2010, at 8:56 PM, Shri wrote:
> 
> > I have three DAs da1,da2, and da3 (1-dimensional DAs) for my power system project where da1 is for the generator subsystem, da2 is for the network subsystem, and da3 is for the load subsystem. Here,the load subsystem size is zero and hence i need to create da3 with zero nodes. However,DACreate1D does not allow zero nodes or zero degrees of freedom and gives an error
> > 
> > [0]PETSC ERROR: --------------------- Error Message ------------------------------------
> > [0]PETSC ERROR: Argument out of range!
> > [0]PETSC ERROR: More processors than data points! 1 0!
> > [0]PETSC ERROR: ------------------------------------------------------------------------
> > [0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 4, unknown
> > [0]PETSC ERROR: See docs/changes/index.html for recent updates.
> > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
> > [0]PETSC ERROR: See docs/index.html for manual pages.
> > [0]PETSC ERROR: ------------------------------------------------------------------------
> > [0]PETSC ERROR: /Users/Shri/TS-dev/TS1 on a osx-debug named dhcp55.swwn2.iit.edu by Shri Wed Sep  8 19:56:14 2010
> > [0]PETSC ERROR: Libraries linked from /Users/Shri/petsc-3.1/osx-debug/lib
> > [0]PETSC ERROR: Configure run at Fri Jul 23 16:06:56 2010
> > [0]PETSC ERROR: Configure options --with-mpi-dir=/Users/Shri/packages/install/mpich2 --with-clanguage=cxx
> > [0]PETSC ERROR: ------------------------------------------------------------------------
> > [0]PETSC ERROR: DACreate_1D() line 152 in src/dm/da/src/da1.c
> > [0]PETSC ERROR: DASetType() line 48 in src/dm/da/src/dareg.c
> > [0]PETSC ERROR: DASetTypeFromOptions_Private() line 65 in src/dm/da/src/dacreate.c
> > [0]PETSC ERROR: DASetFromOptions() line 131 in src/dm/da/src/dacreate.c
> > [0]PETSC ERROR: DACreate1d() line 393 in src/dm/da/src/da1.c
> > [0]PETSC ERROR: CreateLoadSysDA() line 80 in createsubsysda.c
> > [0]PETSC ERROR: main() line 64 in TS1.c
> > application called MPI_Abort(MPI_COMM_WORLD, 63) - process 0
> > 
> > 1) Is it possible to create 1d DAs with zero nodes?? 
> > 2) I can try using vectors instead of DAs for the three subsystems since zero sized vectors can be created. However, i want to pack these three subsystem vectors later on similar to what the DMComposite object does. Did PETSc have a vector packer object similar to DMComposite in any of the previous releases?
> > 
> > Thanks,
> > Shri
>