[petsc-users] ex52_integrateElement.cu

[Matthew Knepley]

On Tue, Mar 27, 2012 at 7:18 PM, David Fuentes <fuentesdt at gmail.com> wrote:

> thanks! when generating the header file ex52.h
>
> what should the input arguments be to the PetscGenerateFEMQuadrature.py
> file ?
>

Andy Terrel already busted me about this today. I pushed a new tarball to
the ftp site, so this will work

  rm externalpackages/Generator
  ./$PETSC_ARCH/conf/reconfigure-${PETSC_ARCH}.py

or if you are more sophisticated with Mercurial

  cd externalpackages/Generator
  hg revert --all
  hg pull -u https://bitbucket.org/knepley/code-generator

Then it should run.

  Thanks,

     Matt


> SCRGP2$ $PETSC_DIR/bin/pythonscripts/PetscGenerateFEMQuadrature.py 3 1 8 1
> laplacian ex52.h
> [{(-1.0, -1.0, -1.0): [(1.0, ())]}, {(1.0, -1.0, -1.0): [(1.0, ())]},
> {(-1.0, 1.0, -1.0): [(1.0, ())]}, {(-1.0, -1.0, 1.0): [(1.0, ())]}]
> {0: {0: [0], 1: [1], 2: [2], 3: [3]}, 1: {0: [], 1: [], 2: [], 3: [], 4:
> [], 5: []}, 2: {0: [], 1: [], 2: [], 3: []}, 3: {0: []}}
> Perm: [0, 1, 2, 3]
> Creating /home/fuentes/snestutorials/ex52.h
> Traceback (most recent call last):
>   File
> "/opt/apps/PETSC/petsc-dev/bin/pythonscripts/PetscGenerateFEMQuadrature.py",
> line 33, in <module>
>     generator.run(elements, numBlocks, operator, filename)
>   File "/opt/apps/PETSC/petsc-dev/config/PETSc/FEM.py", line 967, in run
>     self.outputElementSource(self.getElementSource(elements, numBlocks,
> operator, sourceType = 'GPU'),
> os.path.splitext(filename)[0]+'_gpu'+os.path.splitext(filename)[1])
>   File "/opt/apps/PETSC/petsc-dev/config/PETSc/FEM.py", line 926, in
> getElementSource
>     defns.extend(self.getComputationTypes(element, n))
>   File "/opt/apps/PETSC/petsc-dev/config/PETSc/FEM.py", line 440, in
> getComputationTypes
>     vecType.type       = self.Cxx.typeMap['float'+str(dim)]
> KeyError: 'float3'
>
>
>
>
> On Tue, Mar 27, 2012 at 2:10 PM, Blaise Bourdin <bourdin at lsu.edu> wrote:
>
>>
>> On Mar 27, 2012, at 1:23 PM, Matthew Knepley wrote:
>>
>> On Tue, Mar 27, 2012 at 12:58 PM, David Fuentes <fuentesdt at gmail.com>wrote:
>>
>>> Hi,
>>>
>>> I had a question about the status of example 52.
>>>
>>>
>>> http://petsc.cs.iit.edu/petsc/petsc-dev/file/a8e2f2c19319/src/snes/examples/tutorials/ex52.c
>>>
>>> http://petsc.cs.iit.edu/petsc/petsc-dev/file/a8e2f2c19319/src/snes/examples/tutorials/ex52_integrateElement.cu
>>>
>>>
>>> Can this example be used with a DM object created from an unstructured
>>> exodusII mesh, DMMeshCreateExodus, And the FEM assembly done on GPU ?
>>>
>>
>> 1) I have pushed many more tests for it now. They can be run using the
>> Python build system
>>
>>   ./config/builder2.py check src/snes/examples/tutorials/ex52.c
>>
>>   in fact, you can build any set of files this way.
>>
>> 2) The Exodus creation has to be converted to DMComplex from DMMesh. That
>> should not take me very long. Blaise maintains that
>>      so maybe there will be help :) You will just replace
>> DMComplexCreateBoxMesh() with DMComplexCreateExodus(). If you request
>>      it, I will bump it up the list.
>>
>>
>> DMMeshCreateExodusNG is much more flexible than DMMeshCreateExodus in
>> that it can read meshes with multiple element types and should have a much
>> lower memory footprint. The code should be fairly easy to read. you can
>> email me directly if you have specific questions. I had looked at creating
>> a DMComplex and it did not look too difficult, as long as interpolation is
>> not needed. I have plans to write DMComplexCreateExodus, but haven't had
>> time too so far. Updating the Vec viewers and readers may be a bit more
>> involved. In perfect world, one would write an EXODUS viewer following the
>> lines of the VTK and HDF5 ones.
>>
>> Blaise
>>
>>
>>
>> Let me know if you can run the tests.
>>
>>   Thanks
>>
>>      Matt
>>
>>
>>> Thanks,
>>> David
>>>
>> --
>> 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
>>
>>
>>  --
>> Department of Mathematics and Center for Computation & Technology
>> Louisiana State University, Baton Rouge, LA 70803, USA
>> Tel. +1 (225) 578 1612, Fax  +1 (225) 578 4276
>> http://www.math.lsu.edu/~bourdin
>>
>>
>>
>>
>>
>>
>>
>>
>


-- 
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/20120327/00f7749f/attachment-0001.htm>