Thanks, Barry! Yes, I do have my own command line arguments -1 and 1.
Good to know that it's safe to ignore the warning.
Xiangdong
On Wed, Nov 23, 2011 at 11:30 PM, Barry Smith wrote:
>
> On Nov 23, 2011, at 10:26 PM, Xiangdong Liang wrote:
>
>> Hello everyone,
>>
>> When I compile my program wi
Hello everyone,
When I compile my program with debug mode and ran it, I got
WARNING! There are options you set that were not used!
WARNING! could be spelling mistake, etc!
Option left: name:-1 value: -1
I even got the warning if I do not provide any options. What's the option left?
However, su
On Nov 23, 2011, at 10:26 PM, Xiangdong Liang wrote:
> Hello everyone,
>
> When I compile my program with debug mode and ran it, I got
>
> WARNING! There are options you set that were not used!
> WARNING! could be spelling mistake, etc!
> Option left: name:-1 value: -1
You must have some of
t part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/1462159e/attachment-0001.htm>
they used to do. Make sure you
>> are
>> > > assembling correctly. If it's still confusing, run the old and new
>> code with
>> > >
>> > > -ksp_monitor_true_residual -ksp_converged_reason -ksp_view -pc_type
>> > > bjacobi -sub_pc_type sor
>> > >
>> > > and send the output of both for us to look at.
>> > >
>> > > Also note that you can use -pc_type sor even in parallel. There are
>> > > options for local iterations and full iterations.
>> > >
>> >
>>
>
>
--
Eng. Manuel Ricardo Perez Cerquera. MSc. Ph.D student
Antenna and EMC Lab (LACE)
Istituto Superiore Mario Boella (ISMB)
Politecnico di Torino
Via Pier Carlo Boggio 61, Torino 10138, Italy
Email: manuel.perezcerquera at polito.it
Phone: +39 0112276704
Fax: +39 011 2276 299
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/73dce149/attachment.htm>
to
> > > > > MatSetOption().
> > > > >
> > > > > And also, with the previous version of PETSc I've used, the were
> no
> > > > >> problem using -pc_type bjacobi and -sub_pc_type sor, juste to
> solve
> > > linear
> > > > >> system doing parallel computations and because SOR is not
> > > parallelized. But
> > > > >> now, when I use -pc_type bjacobi and -sub_pc_type sor, with 3
> rank, I
> > > > >> experiment some convergence problem during my simulation.
> > > > >>
> > > > >
> > > > > These options should do the same thing they used to do. Make sure
> you
> > > are
> > > > > assembling correctly. If it's still confusing, run the old and new
> > > code with
> > > > >
> > > > > -ksp_monitor_true_residual -ksp_converged_reason -ksp_view -pc_type
> > > > > bjacobi -sub_pc_type sor
> > > > >
> > > > > and send the output of both for us to look at.
> > > > >
> > > > > Also note that you can use -pc_type sor even in parallel. There are
> > > > > options for local iterations and full iterations.
> > > > >
> > > >
> > >
> >
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/391adccc/attachment.htm>
El 23/11/2011, a las 17:12, Satish Balay escribi?:
> On Wed, 23 Nov 2011, Jose E. Roman wrote:
>
>> I have a question related to this. We tried today (with petsc-dev)
>> -mat_mumps_icntl_28 2 -mat_mumps_icntl_29 2
>> and it works, but
>> -mat_mumps_icntl_28 2 -mat_mumps_icntl_29 1
>> gives an e
But
> > >> now, when I use -pc_type bjacobi and -sub_pc_type sor, with 3 rank, I
> > >> experiment some convergence problem during my simulation.
> > >>
> > >
> > > These options should do the same thing they used to do. Make sure you
> are
&
On Nov 23, 2011, at 1:58 PM, Abraham Taicher wrote:
> Hi,
>
> I'm trying to add a Dirichlet boundary condition to my finite element
> code. I want the vector with the boundary data to split across
> processors. I have a rectangular grid so I could just make 4 different
> 1D DA's and split them
On Nov 23, 2011, at 1:31 PM, Shiyuan wrote:
> I've check it in Matlab. the matrix A is symmetric and has no negative
> eigenvalues( only a zero eigen value). I've also checked that the nullspace
> is correct( norm(A*phi)<1e-11);
Hmm. You can run with -ksp_view_binary and email to petsc-ma
_converged_reason -ksp_view -pc_type
> bjacobi -sub_pc_type sor
>
> and send the output of both for us to look at.
>
> Also note that you can use -pc_type sor even in parallel. There are
> options for local iterations and full iterations.
>
-- next part --
An
El 23/11/2011, a las 16:25, Barry Smith escribi?:
>
> On Nov 23, 2011, at 7:18 AM, Jed Brown wrote:
>
>> On Wed, Nov 23, 2011 at 05:43, Dominik Szczerba
>> wrote:
>> In my procedure considerable time is spent to partition the domain.
>> When using MUMPS as a solver for my matrix I see the mes
On Wed, Nov 23, 2011 at 3:09 PM, Jed Brown
wrote:> On Wed, Nov 23, 2011 at 08:06, Dominik Szczerba
> wrote: I looked through the help index
page for AIJ, matlab and similar>> keywords, but can not find a
suitable function.>> As I said in the other mail, just view to a PETSc
binary file and rea
You are right, thanks!
On Wed, Nov 23, 2011 at 2:25 PM, Jed Brown wrote:
> On Wed, Nov 23, 2011 at 07:22, Dominik Szczerba
> wrote:
>>
>> Would not MatConvert do the same, but easier?
>
> No, how would MatConvert() know on which communicator to put the result?
> What would all the other processe
, each time I do the simulation with
one more processor, I lost time). Any idea of what's wrong ?
I would appreciate some help in that purpose
Best Regards.
John
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/pet
I have converted my MATMPIAIJ to MATSEQAIJ using MatGetSubMatrices.
I need to dump it to a hdf5 file in form of three vectors: si, sj
(int) and sa (double), with their sizes equal to the number of
nonzeros.
I looked through the help index page for AIJ, matlab and similar
keywords, but can not find
I am a bit stuck trying to inspect some elements in a big matrix.
I need to dump a matrix and a vector to a file, that is given, but I
need to make sure some of their values are exactly as I expect them.
ASCII is reportedly not an option, the matrix be too big.
HDF5 Viewer complains it does not sup
Would not MatConvert do the same, but easier?
Thanks
Dominik
On Mon, Nov 21, 2011 at 4:41 PM, Jed Brown wrote:
> On Mon, Nov 21, 2011 at 09:36, Dominik Szczerba
> wrote:
>>
>> I am thinking of the simplest but not necessarily efficient way to
>> pull a distributed matrix onto one CPU... I would
Hi,
I'm trying to add a Dirichlet boundary condition to my finite element
code. I want the vector with the boundary data to split across
processors. I have a rectangular grid so I could just make 4 different
1D DA's and split them up according to the way the 2D DA is split. Is
that what DMDA_BO
;CHKERRV(ierr);
> > ierr=PCSetType(prec,PCNONE);CHKERRV(ierr);
> > ierr=KSPSetTolerances(ksp,1e-5,1e-20,1e5,1);CHKERRV(ierr);
> > ierr=KSPSetFromOptions(ksp);CHKERRV(ierr);
> > ierr=KSPSetUp(ksp);CHKERRV(ierr);
> >
> > Shiyuan
>
>
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/940092a4/attachment.htm>
Print out the matrix for a small grid then check in Matlab that it is
symmetric and has no negative eigenvalues.
Barry
On Nov 23, 2011, at 9:50 AM, Shiyuan wrote:
> Hi,
> I want to solve a singular system with a known nullspace. However, I the
> KSP solve diverges with KSP_INDEFIN
Hong,
Thank you for your help!
I tried the example code of petsc-3.2/src/mat/examples/tests/ex1.c.
But the result is the same. I attach the modified code and the
results here. I don't know what is the problem exactly. By the way I am
using PETSc on a 64bit ubuntu 11.04. I complied MPICH2 myself,
Hong,
Thank you for your help!
I tried the example code of petsc-3.2/src/mat/examples/tests/ex1.c.
But the result is the same. I attach the modified code and the
results here. I don't know what is the problem exactly. By the way I am
using PETSc on a 64bit ubuntu 11.04. I complied MPICH2 myself,
Hong,
Thank you for your help!
I tried the example code of petsc-3.2/src/mat/examples/tests/ex1.c.
But the result is the same. I attach the modified code and the
results here. I don't know what is the problem exactly. By the way I am
using PETSc on a 64bit ubuntu 11.04. I complied MPICH2 myself,
In my procedure considerable time is spent to partition the domain.
When using MUMPS as a solver for my matrix I see the message:
"Ordering based on METIS"
and this seems to take a lot of time. Is it not possible to take over
the existing ordering, or is this ordering something else than matrix
p
> [1]PETSC ERROR: d_nnz cannot be less than 0: local row 78 value -750763693!
There is something wrong in your code. You'll have to verify the
prototypes of all functions - with the petsc-3.2 documentation.
And valgrind is an easy way to pinpoint to the problem source. [also a
debugger should sho
In addition I set FieldCtx* FC = (FieldCtx*) ctx;
in the _petsc_resiadualVector before any calculations
Thanks a lot
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111123/7b600894/attachment.htm>
ian
6- set opitions
7- solve with SNESSolve()
Is it look right with SNES ???
Thanks again for your attention.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/5d65352e/attachment-0001.htm>
?
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111123/2340cc81/attachment.htm>
.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111123/b312f980/attachment.htm>
;
FC.f = f;
FC.flw = flw;
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111123/98f6c28e/attachment.htm>
Since you get a SEGV - I would suggest running the code in the
debugger - to check where its crashing.
Also run with valgrind to see where problems start.. Mostlikely the
issues would be change in prototypes for PETSc functions - between
releases.
Satish
On Wed, 23 Nov 2011, jean-frederic thebau
On Wed, 23 Nov 2011, Jose E. Roman wrote:
> I have a question related to this. We tried today (with petsc-dev)
> -mat_mumps_icntl_28 2 -mat_mumps_icntl_29 2
> and it works, but
> -mat_mumps_icntl_28 2 -mat_mumps_icntl_29 1
> gives an error of MUMPS complaining that PTScotch was not enabled.
>
>
d...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/e2045a09/attachment.htm>
On Nov 23, 2011, at 7:18 AM, Jed Brown wrote:
> On Wed, Nov 23, 2011 at 05:43, Dominik Szczerba
> wrote:
> In my procedure considerable time is spent to partition the domain.
> When using MUMPS as a solver for my matrix I see the message:
>
> "Ordering based on METIS"
>
> This is an ordering
can pull from
http://petsc.cs.iit.edu/petsc/releases/petsc-3.2or apply the attached
patch.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/d75e35c5/attachment-0001.htm>
-- next part
send the output of both for us to look at.
Also note that you can use -pc_type sor even in parallel. There are options
for local iterations and full iterations.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/3d7ef281/attachment.htm>
aryRead.m
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/ff2ea452/attachment-0001.htm>
Values()
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/0a75f816/attachment.htm>
ngs.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/f852cdd6/attachment.htm>
MUMPS developers
or use a different package if you don't like these properties.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/40e31ad9/attachment.htm>
HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/e21877e8/attachment.htm>
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/b6737352/attachment.htm>
the state passed into your SNES residual function, not based on whatever is
in this other struct.
-- next part --
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/2023/001a3e3c/attachment.htm>
44 matches
Mail list logo