On Mon, Apr 4, 2016 at 1:27 PM, Jose E. Roman wrote:
>
> > El 4 abr 2016, a las 19:13, David Knezevic
> escribió:
> >
> > OK, thanks, I'll have a look at that paper.
> >
> > And just to confirm that I've understood properly: You're saying that
> the GNHEP solver should converge more robustly tha
> El 4 abr 2016, a las 19:13, David Knezevic
> escribió:
>
> OK, thanks, I'll have a look at that paper.
>
> And just to confirm that I've understood properly: You're saying that the
> GNHEP solver should converge more robustly than the GHIEP solver for
> symmetric indefinite problems? Are t
On Mon, Apr 4, 2016 at 12:47 PM, Jose E. Roman wrote:
>
> > El 4 abr 2016, a las 18:32, David Knezevic
> escribió:
> >
> > Hi Jose,
> >
> > I'm interested to know more about this comment from your email below:
> > "You can also try the symmetric-indefinite solver, but this is not
> recommended s
> El 4 abr 2016, a las 18:43, Manav Bhatia escribió:
>
> Thanks!
> I upgraded to 3.6.3 and ran the code, and it went through all the load steps
> without problem.
>
> A followup question:
>
> It seems like the eigenvalues are sorted based on abs(lambda). If I am
> interested in only the l
> El 4 abr 2016, a las 18:32, David Knezevic
> escribió:
>
> Hi Jose,
>
> I'm interested to know more about this comment from your email below:
> "You can also try the symmetric-indefinite solver, but this is not
> recommended since it is not numerically stable."
>
> I frequently use SLEPc's
Thanks!
I upgraded to 3.6.3 and ran the code, and it went through all the load steps
without problem.
A followup question:
It seems like the eigenvalues are sorted based on abs(lambda). If I am
interested in only the lowest N eigenvalues greater than 0, is there a way to
tell slepc to ignor
Hi Jose,
I'm interested to know more about this comment from your email below:
"You can also try the symmetric-indefinite solver, but this is not
recommended since it is not numerically stable."
I frequently use SLEPc's symmetric-indefinite solver for buckling
eigenvalue problems (which are natur
> El 4 abr 2016, a las 18:23, Manav Bhatia escribió:
>
> Thanks, Jose!
>
> I am currently running 3.6.2, and will update to 3.6.3.
>
> Is there a recommended strategy to automatically switch from GHEP to GNHEP
> for some subset of problems? Or should I choose to run all my eigenprobelms
> w
Thanks, Jose!
I am currently running 3.6.2, and will update to 3.6.3.
Is there a recommended strategy to automatically switch from GHEP to GNHEP for
some subset of problems? Or should I choose to run all my eigenprobelms with
GNHEP?
Regards,
Manav
> On Apr 4, 2016, at 11:18 AM, Jose E. Rom
> El 4 abr 2016, a las 18:06, Manav Bhatia escribió:
>
> Hi Jose,
>
>I also read these matrices into matlab and found the eigenvalues as
>
> >>A = PetscBinaryRead('A.petsc’);
> >>B = PetscBinaryRead(‘B.petsc’);
> >> [v,d] = eigs(A,B)
> (*** got a lot of output about poor-conditioning ***)
Hi Jose,
I also read these matrices into matlab and found the eigenvalues as
>>A = PetscBinaryRead('A.petsc’);
>>B = PetscBinaryRead(‘B.petsc’);
>> [v,d] = eigs(A,B)
(*** got a lot of output about poor-conditioning ***)
>> diag(d)
ans =
1.0e-05 *
-0.2219
0.0229
0.0229
0.0
I just sent you the updated matrices on slepc-main.
Sorry about sending the wrong matrices earlier.
Regards,
Manav
> On Apr 4, 2016, at 9:42 AM, Manav Bhatia wrote:
>
> Ok. So, I ran ex7 with the same command-line options in your email, and got a
> result. This is on my Mac OS X, without a
It turns out that the matrices that I sent you was from a load step before the
one that failed.
I am calling eps_solve once for each load step, and I get the error at the 14th
load step. So, the matrices that I sent were from step 13.
The eigenvalues you got match the 1/lambda that I get from
On Mon, Apr 4, 2016 at 9:42 AM, Manav Bhatia wrote:
> Ok. So, I ran ex7 with the same command-line options in your email, and
> got a result. This is on my Mac OS X, without any changes to the
> lapack/blas/slepc library.
>
> I also ran my code on a linux machine with lapack/blas build from sourc
Following is the output of -eps_view from ex7 and from my code. I have included
only what is different.
There are differences in nev, ncv, mpd, etc.
-Manav
output from ex7:
EPS Object: 1 MPI processes
type: arnoldi
problem type: generalized symmetric eigenvalue problem
selected po
> El 4 abr 2016, a las 16:42, Manav Bhatia escribió:
>
> Ok. So, I ran ex7 with the same command-line options in your email, and got a
> result. This is on my Mac OS X, without any changes to the lapack/blas/slepc
> library.
>
> I also ran my code on a linux machine with lapack/blas build fr
Ok. So, I ran ex7 with the same command-line options in your email, and got a
result. This is on my Mac OS X, without any changes to the lapack/blas/slepc
library.
I also ran my code on a linux machine with lapack/blas build from source
(obtained from netlib), and got the same error as on my m
Thanks for looking into this, Jose.
I will try with the source LAPACK. I am currently running on MAC OS X with the
system provided lapack and bias. I have heard of some issues with these
libraries, but have not had any issues until now.
Hopefully that will sort the issue. I will let you know w
> El 3 abr 2016, a las 22:17, Manav Bhatia escribió:
>
> I just send you the matrices.
>
> Thanks,
> Manav
I cannot reproduce the problem. I was able to compute the eigenvalues without
problems with Krylov-Schur (note that I scaled matrix A by 1e7 because the
eigenvalues are tiny):
$ ./ex7
I just send you the matrices.
Thanks,
Manav
> On Apr 3, 2016, at 3:03 PM, Jose E. Roman wrote:
>
>
>> El 3 abr 2016, a las 21:45, Manav Bhatia escribió:
>>
>> Hi Jose,
>>
>> I did not specify Arnoldi. I am running on a single CPU, so maybe Arnoldi
>> is the default for 1 cpu runs?
>
> El 3 abr 2016, a las 21:45, Manav Bhatia escribió:
>
> Hi Jose,
>
>I did not specify Arnoldi. I am running on a single CPU, so maybe Arnoldi
> is the default for 1 cpu runs?
No. Maybe you set it in code or maybe the option is being picked from
PETSC_OPTIONS or petscrc file.
>
>
Hi Jose,
I did not specify Arnoldi. I am running on a single CPU, so maybe Arnoldi
is the default for 1 cpu runs?
How would you like me to send you the matrices? Is there a specific format
you would like?
Thanks,
Manav
> On Apr 3, 2016, at 2:39 PM, Jose E. Roman wrote:
>
>
>> E
> El 3 abr 2016, a las 20:47, Manav Bhatia escribió:
>
> Hi,
>
>I am using slepc to solve for the natural frequencies of a
> small-disturbance modal eigensolution of a plate with compressive stresses.
>
> M x = (1/omega^2) K(lambda) x
>
>lambda is the loading parameter.
Hi,
I am using slepc to solve for the natural frequencies of a small-disturbance
modal eigensolution of a plate with compressive stresses.
M x = (1/omega^2) K(lambda) x
lambda is the loading parameter. I first solve for the nonlinear static
solution for a given load, which
24 matches
Mail list logo