I was looking into the sparsity pattern code a bit, and I was wondering
if the !need_full_sparsity_pattern branch in
SparsityPattern::Build::parallel_sync() over-estimates n_nz and n_oz? It
doesn't do the equivalent of the my_row.erase operation that we have in
the need_full_sparsity_pattern br
That is indeed the case. Petsc allows you to overestimate with no major issue,
so it's faster to build the graph that way. Others require the full graph,
which is more expensive to build.
> On Oct 13, 2014, at 5:17 PM, David Knezevic
> wrote:
>
> I was looking into the sparsity pattern cod
Makes sense, thanks.
On 10/13/2014 09:33 PM, Kirk, Benjamin (JSC-EG311) wrote:
> That is indeed the case. Petsc allows you to overestimate with no major
> issue, so it's faster to build the graph that way. Others require the full
> graph, which is more expensive to build.
>
>
>
>> On Oct 13, 20
