On Tue, Feb 16, 2016 at 10:19 PM, Harshad Sahasrabudhe
wrote:
> Hi John,
>
> I was running an older version of LibMesh when I got the previous segfault
> in libMesh::DofMap::clear_sparsity(), which happened after changing the
> assemble for CondensedEigenSystem to only set diagonal dofs in B. The
Hi John,
I was running an older version of LibMesh when I got the previous segfault
in libMesh::DofMap::clear_sparsity(), which happened after changing the
assemble for CondensedEigenSystem to only set diagonal dofs in B. Then I
upgraded to libMesh 0.9.4 and that segfault went away. Now the code r
Hi John and David,
Thanks for your comments. I'm working on rebuilding our environment to get
the full profiling. I have done some preliminary profiling, the scaling
that I posted above was only for the call "_system->solve()" on the
CondensedEigenSystem.
I suggest you consult the SLEPc manual, a
On Tue, Feb 16, 2016 at 10:12 AM, Harshad Sahasrabudhe
wrote:
> Hi John,
>
> Thanks for the response.
>
>
>> I'm using KrylovSchur eigensolver
>>> in SLEPc and I find that the eigenvalue computation anti-scales when I
>>> use
>>>
>>
>> anti-scales?
>>
>
> Yes, I get the following scaling
>
> Proc
On Tue, Feb 16, 2016 at 12:18 PM, Harshad Sahasrabudhe
wrote:
> Hi David,
>
>
>> One interesting thing that I didn't realize for a long time is that PETSc
>> condenses the sparsity pattern when you call matrix->close(), i.e.:
>>
>> http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/M
Hi David,
> One interesting thing that I didn't realize for a long time is that PETSc
> condenses the sparsity pattern when you call matrix->close(), i.e.:
>
> http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatAssemblyBegin.html
> "Space for preallocated nonzeros that is not fill
Hi John,
Thanks for the response.
> I'm using KrylovSchur eigensolver
>> in SLEPc and I find that the eigenvalue computation anti-scales when I use
>>
>
> anti-scales?
>
Yes, I get the following scaling
Processes Eigensolver time(s)
1 68.3037
2 51.960
On Tue, Feb 16, 2016 at 12:00 PM, John Peterson
wrote:
> On Tue, Feb 16, 2016 at 9:51 AM, Harshad Sahasrabudhe >
> wrote:
>
> > Hi,
> >
> > I'm using QTrap to integrate a differential equation for a generalized
> > eigenvalue problem. I'm using CondensedEigenSystem to solve the equation.
> > Whe
On Tue, Feb 16, 2016 at 9:51 AM, Harshad Sahasrabudhe
wrote:
> Hi,
>
> I'm using QTrap to integrate a differential equation for a generalized
> eigenvalue problem. I'm using CondensedEigenSystem to solve the equation.
> When I use QTrap, I get a diagonal B matrix in
>
> A \psi = \lambda B \ps
Hi,
I'm using QTrap to integrate a differential equation for a generalized
eigenvalue problem. I'm using CondensedEigenSystem to solve the equation.
When I use QTrap, I get a diagonal B matrix in
A \psi = \lambda B \psi
So my question is: Is the sparsity pattern of B set according to the
qua
10 matches
Mail list logo
