[petsc-users] How to create a quadrature object?

[Matthew Knepley]

On Thu, Jul 28, 2016 at 7:08 PM, Andrew Ho <andrewh0 at uw.edu> wrote:

>
> On Thu, Jul 28, 2016 at 4:35 PM, Matthew Knepley <knepley at gmail.com>
> wrote:
>
>>
>> This is strange. First, you should only have uninitialized classids if
>> you build without dynamics libraries. Did you? If you
>> include the entire error output, it would show us.
>>
>>
> This is pretty much the entire error output; I trimmed off the ending "MPI
> Abort called" message (it's the standard message printed by OpenMPI when
> MPIAbort is called). The code snippet I had is a complete working example
> which is able to give me this error.
>
>
>> If so, then PetscInitialize should set the value of this classid, unless
>> you built this with PETSC_HAVE_DYNAMIC_LIBRARIES,
>> but linked against static libraries, as the error message says. Do you
>> have multiple versions of PETSc on your machine?
>>
>
> I have a single install of PETSc at any one time. I've tried builds with
> --with-shared-libraries on and off with the same results (I made sure to
> purge any previous install first). Here's a diff if you want to test the
> same configs I tried running.
>

Crap. We changed the way initialization works, but this was left behind.
You can make this change

diff --git a/src/dm/dt/interface/dt.c b/src/dm/dt/interface/dt.c
index 5d12959..d6f3454 100644
--- a/src/dm/dt/interface/dt.c
+++ b/src/dm/dt/interface/dt.c
@@ -50,7 +50,7 @@ PetscErrorCode PetscQuadratureCreate(MPI_Comm comm,
PetscQuadrature *q)

   PetscFunctionBegin;
   PetscValidPointer(q, 2);
-  ierr = DMInitializePackage();CHKERRQ(ierr);
+  ierr = PetscSysInitializePackage();CHKERRQ(ierr);
   ierr =
PetscHeaderCreate(*q,PETSC_OBJECT_CLASSID,"PetscQuadrature","Quadrature","DT",comm,PetscQuadratureDestroy,PetscQuadratureView);CHKERRQ(ierr);
   (*q)->dim       = -1;
   (*q)->order     = -1;

and then rebuild

  cd $PETSC_DIR
  make -f ./gmakefile

and try running your example again.

I will get this change in.

  Thanks,

     Matt


> Just run *make ex4* in the src/dm/dt/examples/tests folder
>
> Configure command:
>
> ./configure --with-debugging=0 COPTFLAGS="-O3 -march=native"
> CXXOPTFLAGS="-O3 -march=native" FOPTFLAGS="-O3 -march=native"
>
> --
> Andrew Ho
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160728/24da1d50/attachment-0001.html>