You are right, of course, Wolfgang.
I should offer to contribute in the development of this method. I will.
Thank you
Alberto
*Alberto Salvadori* Dipartimento di Ingegneria Meccanica e Industriale
(DIMI)
Universita` di Brescia, via Branze 43, 25123 Brescia
Italy
tel 030 3715426
e-mail:
alber
On 12/20/20 2:58 AM, Alberto Salvadori wrote:
Point centroid = cell->face(face_number)->barycenter() ;
provides always centroid = {0,0,0}, whereas the following
Point centroid = cell->face(face_number)->center() ;
returns the correct amount.
Am I misunderstanding something or may it be a bug
Zachary,
I haven't tried to debug your code snippets, but would suggest to take a look
how step-40 sets up and builds the sparsity pattern. There, we call
make_sparsity_pattern(), which really just adds certain entries to the
sparsity pattern like you do, so if you understand how the set up w
I should mention how I am creating my dynamic sparsity pattern. The
structure is like this:
dsp.reinit(local_number_of_rows, number_of_columns); // the number of rows
is divided among processors
for (int i =0; i < local_number_of_rows; ++i )
for (int j = 0; j < number_of_columns; ++j )
Hi there,
I am trying to build a parallel PETSc sparse matrix by building a local
sparsity pattern and subsequently distributing the pattern to all processes so
I then have the full sparsity pattern for the PETSc sparse matrix (I read this
needs to be done but please correct me if I am wrong).
Dear community
I observed that this instruction:
Point centroid = cell->face(face_number)->barycenter() ;
provides always centroid = {0,0,0}, whereas the following
Point centroid = cell->face(face_number)->center() ;
returns the correct amount.
Am I misunderstanding something or may it be a b