Hello,
I'trying to solve an ODE system using PETSc TS through its wrapper
petsc4py, and I'm having problems with the setting. The problem is of the
kind: u_t = G(t, u) with u being a vector with DoF from a FEM problem and G
is computed by a series of matrix operations and manipulations, which are
i
Hello, I am trying to save some field stored in a vector which has values
associated with vertexes, edges and cells in a DMPlex.
This vector was created (using petsc4py) from a petcs section, setting this
as default section and then creating the vector from the DMPlex.
The thing is that when I save
;>
>
> Justin is correct. However, I will note that this strategy is not
> scalable. Usually the data is much larger than the mesh,
> so this runs out of memory much sooner. Are you sure that this is what you
> want?
>
> Thanks,
>
> Matt
>
>
>> On Thu,
2
> (and maybe others?) the first step is handled through the mesh generation
> functions and step 3 is handled through step 4
>
> Thanks,
> Justin
>
> On Wednesday, December 16, 2015, Alejandro D Otero
> wrote:
>
>> Hi, I need some help understanding how to
Hi, I need some help understanding how to distribute data together with a
dmplex representing a FE mesh.
At the beginning I define the structure of the dmplex assigning certain
number of DoF to cells, edges and vertexes, in one process (the dmplex in
the rest is empty)
I create a petscsecton and I
Ok, I'll think about it.
Thank you!
On Wed, Nov 18, 2015 at 4:29 PM, Jed Brown wrote:
> Alejandro D Otero writes:
>
> > Hi, I am trying to define a spectral element mesh in the frame of DMPlex
> > structure.
> > To begin with I am thinking of quadrilateral 2d
Hi, I am trying to define a spectral element mesh in the frame of DMPlex
structure.
To begin with I am thinking of quadrilateral 2d elements. Lets say I have a
16 node element with 2 nodes on each edge besides the 2 corners and 4 nodes
inside the element. I'd like to treat all nodes as vertex in th
Hi all,
I am trying to nest to square matrices to form a rectangular one with dims
2n x n (with n the size of the square matrices)
I want the two matrices to have interlaced rows so I passed the following 2
list of indices to MatCreateNest:
[[ 0 2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 3
Hi all, is it posible to build a rectangular matrix that relates two fields
with (posible) different number of DoF per point defining the nonzero
structure from 2 petscsections representing one field each?
As an example, I want to build an operator that given a 2D velocity
computes the correspondi
Hi, I seems to be having problems with PETSc vector VTK viewer. The output
file I get has some values which are different from what I actually have.
I have a scalar field represented by a vector. If I print it on screen I
got the following output.
Vec Object: 2 MPI processes
type: mpi
Process [
Thanks Matt.
It is much clearer for me now.
Regards,
Alejandro
On Fri, Jul 10, 2015 at 6:14 PM, Matthew Knepley wrote:
> On Fri, Jul 10, 2015 at 6:48 AM, Alejandro D Otero
> wrote:
>
>> Hi Matt, thanks for your answer.
>>
>> I got a vector from getCoordinates(
ks for your help,
A.
On Thu, Jul 9, 2015 at 5:18 PM, Matthew Knepley wrote:
> On Thu, Jul 9, 2015 at 7:42 AM, Alejandro D Otero
> wrote:
>
>> Hi, sorry if this is an obvious question, but I cannot figure out how to
>> recover finite element nodes coordinates once I have distr
Hi, sorry if this is an obvious question, but I cannot figure out how to
recover finite element nodes coordinates once I have distributed a mesh
stored as a dmplex. I am using petsc4py as interface to petsc rutines.
I first created a dmplex using:
dm.createFromCellList()
In a sequential run I got
13 matches
Mail list logo