[petsc-users] ex52_integrateElement.cu

Matthew Knepley knepley at gmail.com
Tue Mar 27 20:37:18 CDT 2012


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.
>

David and Blaise,  I have converted this function, now
DMComplexCreateExodus(). Its not tested, but I think
Blaise has some stuff we can use to test it.

  Thanks,

     Matt


> 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/0d75e6a6/attachment.htm>


More information about the petsc-users mailing list