Re: [petsc-users] unsorted local columns in 3.8?

2017-11-02 Thread Randy Michael Churchill
rror Message >> ------ >> [0]PETSC ERROR: No support for this operation for this object type >> [0]PETSC ERROR: unsorted iscol_local is not implemented yet >> [1]PETSC ERROR: - Error Mess

Re: [petsc-users] unsorted local columns in 3.8?

2017-11-02 Thread Randy Michael Churchill
Thanks Mark, both your suggestions of using pc hypre or turning off repartitioning does indeed make the error go away. On Thu, Nov 2, 2017 at 8:51 AM, Mark Adams <mfad...@lbl.gov> wrote: > > > On Wed, Nov 1, 2017 at 9:36 PM, Randy Michael Churchill <rchur...@pppl.gov > >

Re: [petsc-users] unsorted local columns in 3.8?

2017-11-01 Thread Randy Michael Churchill
Doing some additional testing, the issue goes away when removing the gamg preconditioner line from the petsc.rc: -pc_type gamg On Wed, Nov 1, 2017 at 8:23 PM, Hong wrote: > Randy: > Thanks, I'll check it tomorrow. > Hong > > OK, this might not be completely satisfactory,

[petsc-users] unsorted local columns in 3.8?

2017-10-30 Thread Randy Michael Churchill
I'm running a Fortran code that was just changed over to using petsc 3.8 (previously petsc 3.7.6). An error was thrown during a KSPSetUp() call. The error is "unsorted iscol_local is not implemented yet" (see full error below). I tried to trace down the difference in the source files, but where

Re: [petsc-users] Updating Fortran code to petsc 3.8

2017-10-30 Thread Randy Michael Churchill
e you verify it behaves as the previous > version. Note that on any computer including those at NERSc if 3.8 is not > installed by default you may install it yourself. > >Barry > > > On Oct 30, 2017, at 12:18 PM, Randy Michael Churchill <rchur...@pppl.gov> > wro

[petsc-users] Updating Fortran code to petsc 3.8

2017-10-30 Thread Randy Michael Churchill
I'm updating my Fortran code to petsc 3.8. I have several modules with types, and often when a type is initialized in the code, there are certain Petsc variables that aren't ready to be initialized, so they are set to 0. Later, in several parts of the code, these are checked with and if statments,

Re: [petsc-users] TAO setup with modules in Fortran 90

2017-10-10 Thread Randy Michael Churchill
rds that. On Fri, Oct 6, 2017 at 5:59 PM, Barry Smith <bsm...@mcs.anl.gov> wrote: > > > On Oct 6, 2017, at 11:18 PM, Randy Michael Churchill <rchur...@pppl.gov> > wrote: > > > > So if I'm limited to petsc 3.7.6 for reasons of eventually using within > an exist

Re: [petsc-users] TAO setup with modules in Fortran 90

2017-10-06 Thread Randy Michael Churchill
pages/ > Sys/UsingFortran.html > >Matt > > >> >>Barry >> >> >> > On Oct 5, 2017, at 11:48 PM, Randy Michael Churchill <rchur...@pppl.gov> >> wrote: >> > >> > A simple setup question with TAO: if I were to convert

[petsc-users] TAO setup with modules in Fortran 90

2017-10-05 Thread Randy Michael Churchill
A simple setup question with TAO: if I were to convert the rosenbrock1f.F90 example to use a module instead of common structures, how would I setup the include statements? I've tried various combinations (using petscXXXdef.h, petscXXX.h, petscXXX.h90, along with use petscXXX), but seem to get