[petsc-users] DMPlex - Section.

[Afanasiev Michael]

Hi,

I’d like to define vectors on a subset of points defined by a DM (created by DMPlex). For example, if I’m modelling wave propagation in a model with both solid and fluid regions, I’d like to define ux, uy, uz vectors in the solid part, and p in the acoustic part. Basically, just a varying number of dofs per integration point. Looking through the documentation I’ve found a pretty comprehensive discussion of setting up the PetscSection object. Is an inroad here to define a custom PetscSection on a subset of elements, faces, edges, vertices, etc….? Or is there a better way to approach this problem? Or is this functionality supported at all? We’re pretty familiar with PLEX by now, just looking for ways to define vectors on subsets of the whole domain.

Thanks,
Mike.
--
Michael Afanasiev
Ph.D. Candidate
Computational Seismology
Institut für Geophysik
ETH Zürich

Sonneggstrasse 5, NO H 39.2
CH 8092 Zürich
michael.afanasiev at erdw.ethz.ch