Hi,
I have been PETSc together with my CFD code. There seems to be a bug
with the Intel compiler such that when I call some DM routines such as
DMLocalToLocalBegin, a segmentation violation will occur if full
optimization is used. I had posted this question a while back. So the
current soluti
The resulting matrix has something like
>>> 11808*11808*1.e-6
14,399,953,920.036863
nonzero entries. It is possible that some integer operations are overflowing
since C int can only go up to about 4 billion before overflowing.
You can building with a different PETSC_ARCH value us
> On Jun 5, 2017, at 11:20 AM, Francesco Migliorini
> wrote:
>
> Dear Stefano,
> Thank you for your answer. I tried to use VecScatterCreateToAll as you
> suggested but it does not work since the first processor can only view its
> part of the vector. Here's how I managed the code:
>
> Vec
On 6 June 2017 at 17:45, Franck Houssen wrote:
> How to VecScatter from global to local vector, and then, VecGather back ?
>
> This is a very simple use case: I need to split a global vector in local
> (possibly overlapping) pieces, then I need to modify each local piece (x2),
> and finally I nee
How to VecScatter from global to local vector, and then, VecGather back ?
This is a very simple use case: I need to split a global vector in local
(possibly overlapping) pieces, then I need to modify each local piece (x2), and
finally I need to assemble (+=) back local parts into a global vecto
> On 5 Jun 2017, at 23:01, Matthew Knepley wrote:
>
> To do a FEM integral on a facet f I need:
>
> i) to evaluate coefficients at quadrature points (on the facet)
> ii) to evaluate basis functions at quadrature points (on the facet)
>
> for (i), I need all the dofs in closure(support(f)).
>