Re: [petsc-users] DM question.

2015-10-23 Thread Matthew Knepley
On Fri, Oct 23, 2015 at 3:01 AM, Afanasiev Michael < michael.afanas...@erdw.ethz.ch> wrote: > Hi Matthew, > > So I’m discretizing via a tensor product of Lagrange polynomials > co-located at the Gauss-Lobatto-Legendre (GLL) points. The polynomial order > might range from 4 to 10 or something like

[petsc-users] DMPlex natural to global mappings

2015-10-23 Thread Lawrence Mitchell
Dear all, I'm trying to understand how to use the recent natural-global mappings that have made it into dmplex. The goal is to do checkpoint restart (onto, perhaps, differing numbers of processes) with DMs with a non-zero overlap (I can see that this is currently unsupported, but let's try

Re: [petsc-users] DMCreateMatrix with some dense row

2015-10-23 Thread Sang pham van
Hi Barry, The routine DMCreateMatrix_DA_3d_MPIAIJ has 2 input arguments (DM da, Mat J), the function pointer in DMDASetGetMatrix() only accept function with that two arguments. As you suggested, I am writing a routine (based on DMCreateMatrix_DA_3d_MPIAIJ()) to preallocate the matrix in the way I

[petsc-users] Question - read text file

2015-10-23 Thread Marcelo Xavier Guterres
Dear friends, I need to read a matrix from a text file (txt). This matrix has a high dimension (2000 X 2000) . Example: date.txt - provision in the text file: 1 2 3 5 5 5 7 8 5 6 8 9 What's the proper way to perform the action? thanks for help -- -- -- -- -- -- -- -- -- -- -- -- -- -- --

Re: [petsc-users] petsc-users你好:(AD)

2015-10-23 Thread シ?
petsc-users 淘宝2015双十一5000万红包24号0点首发,不管大小,人人有份,晚了无。 点�艚�入红包专区

Re: [petsc-users] DM question.

2015-10-23 Thread Afanasiev Michael
Hi Matthew, So I’m discretizing via a tensor product of Lagrange polynomials co-located at the Gauss-Lobatto-Legendre (GLL) points. The polynomial order might range from 4 to 10 or something like that. The current problem is solved on 2D hexes. I had found DMPlexCreateSection, and followed