[petsc-users] PETSc DAE solver

[Barry Smith]

On Feb 20, 2013, at 2:08 PM, "Jin, Shuangshuang" <Shuangshuang.Jin at pnnl.gov> wrote:

> Hi, Shri, I was able to run the ex9bus.c code successfully with ‘mpirun –n 1 ex9bus’, but failed with ‘mpirun –n 2 ex9bus’ with the following message:

   This particular example is not set up to run in parallel. The PETSc solvers are all parallel but to use the ODE integrators in parallel one must provide the function evaluation (and Jacobian optionally) also in parallel. For this example we did not provide the parallel function evaluation.

   Many of the examples in the ts/examples/tutorials are parallel as are in the subdirectory advection-diffusion-reaction

   Barry

>  
> [0]PETSC ERROR: --------------------- Error Message ------------------------------------
> [0]PETSC ERROR: No support for this operation for this object type!
> [0]PETSC ERROR: Only for sequential runs!
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Development HG revision: fab2442df69c95ef24bf6631a38c7dc028d7ee34  HG Date: Tue Feb 19 16:16:15 2013 -0600
> [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: ex9bus on a arch-linux2-c-debug named olympus.local by d3m956 Wed Feb 20 11:36:02 2013
> [0]PETSC ERROR: Libraries linked from /pic/projects/ds/petsc-dev/arch-linux2-c-debug/lib
> [0]PETSC ERROR: Configure run at Wed Feb 20 11:01:15 2013
> [0]PETSC ERROR: Configure options
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: main() line 807 in src/ts/examples/tutorials/power_grid/stability_9busex9bus.c
>  
> Does that mean  the parallelized version of DAE solver in PETSc is not available yet? Or I missed some PETSc options to enable the function?
>  
> Another question is I also found two other examples in previous installed non-dev PETSc version:
> petsc-3.3-p6/src/ts/examples/tutorials/ex8.c and ex19.c. I felt that they are somewhat similar to my problem. Do you think those two are useful examples to me as well? However, it seems like they are uniprocessor examples only too. Please correct me if I understand it wrong.
>  
> Thanks,
> Shuangshuang
>  
>  
> From: petsc-users-bounces at mcs.anl.gov [mailto:petsc-users-bounces at mcs.anl.gov] On Behalf Of Shri
> Sent: Wednesday, February 20, 2013 10:49 AM
> To: PETSc users list
> Subject: Re: [petsc-users] PETSc DAE solver
>  
> http://www.mcsanl.gov/petsc/developers/index.html
>  
> On Feb 20, 2013, at 12:39 PM, Jin, Shuangshuang wrote:
> 
> 
> Shri, thank you! That’s very helpful information. I’m going to obtain the development version of PETSc now.
>  
> Thanks,
> Shuangshuang
>  
> From: petsc-users-bounces at mcs.anl.gov [mailto:petsc-users-bounces at mcs.anl.gov] On Behalf Of Shri
> Sent: Wednesday, February 20, 2013 10:34 AM
> To: PETSc users list
> Subject: Re: [petsc-users] PETSc DAE solver
>  
> Jin,
>    If you are working on integrating power grid DAE equations  then you should use petsc-dev instead as we've added an example for it. petsc-dev/src/ts/examples/tutorials/power_grid/stabiltiy_9bus/ex9bus.c. There are other examples in the power_grid directory.
>    
> Shri
>  
> On Feb 20, 2013, at 12:19 PM, Jin, Shuangshuang wrote:
> 
> 
> 
> Hi, I’m trying to use PETSc to solve a set of first-order differential and algebraic equations such as:
> x' =  f(x,y)      
> 0 = g(x,y)
>  
> The resources I can find online so far are
> 	• PETSc Manual: Chapter 6 TS: Scalable ODE and DAE Solvers
> 	• petsc-3.3-p6/src/ts/examples/tutorials
> which still seems  vague to me at this moment.
>  
> Can anyone tell me if there’s any more specific references I can refer to as a starting point to tackle this problem?
>  
> Thanks,
> Shuangshuang
>  
>