Yup, this comes from trying to patch more and more flexibility on
a system
that wasn't originally designed for that particular flexibility. The
matrix and vector
are in a strange state at this point that is not defined by our
original design.
One way we can mangle our minds around
I've just noticed that now SNES does not create the inner KSP at the
time SNESCreate() is called. Then I could easily manage to break
things (basically, if SNESSetFromOptions() is never called)
Does it make sense to add:
if (!snes-ksp) { SNESGetKSP(snes, snes-ksp);}
near the begining of
before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20081202/b284bade
I still fu**ing hate hg. There is no trivial way I know to
determine
when this happened. i.e. when I was so f**ing stupid!
This is seriously, seriously, seriously bad! We do not/cannot
require
that SNESSetFromOptions() be called.
I will fix this inconceivably stupid design change
Matt,
You are wrong!
I will fix this today.
Barry
On Dec 2, 2008, at 8:29 AM, Matthew Knepley wrote:
I am not against requiring SetFromOptions() to create default
objects. We
already require it to set a default type.
Matt
On Tue, Dec 2, 2008 at 8:08 AM, Lisandro
their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20081202/3ba0456a/attachment.html
On Dec 2, 2008, at 8:43 AM, Matthew Knepley wrote:
On Tue, Dec 2, 2008 at 8:41 AM, Barry Smith bsmith at mcs.anl.gov
wrote:
Matt,
You are wrong!
I will fix this today.
Okay, thats fine. For hg, you can do annotate on the file to see
what lines
come from what version,
Other (really nice!) change I've noticed is that now we can call
{Vec|Mat}Create() and next {Vec|Mat}SetType() and it works!
However, see this.
In [1]: from petsc4py import PETSc
In [2]: x = PETSc.Vec().create()
In [3]: x.setType('seq')
In [4]: print x.getType()
None
In [5]: A =
their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20081202/7173c780/attachment.html
On Tue, Dec 2, 2008 at 2:06 PM, Matthew Knepley knepley at gmail.com wrote:
This is strange. I just tried KSP ex2 and after MatSetType() the type_name
is
set. What does *.setType() call?
It's a plain call to {Vec|Mat}SetType() . Matt, look at the actual
functions in core PETSc. If the sizes
On Mon, 1 Dec 2008, Lisandro Dalcin wrote:
I'm testing configure passing '--prefix=', and after building, I do
$ su -c 'make install'
But then, some files ('configure.log', 'RDict.db') are left with
'root:root' owner and group. If you want to build/rebuild again, this
fails bad.
And
they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-- next part --
An HTML attachment was scrubbed...
URL:
http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20081202/d0e1aeec/attachment.html
On Tue, 2 Dec 2008, Matthew Knepley wrote:
As a hack - I can add os.remove('RDict.db') in config/installer.py
I am for this.
ok - I pushed the hack [deletes more than RDict.db - as I see
different files created in different situations..]
Satish
On Tue, Dec 2, 2008 at 4:46 PM, Barry Smith bsmith at mcs.anl.gov wrote:
Yup, this comes from trying to patch more and more flexibility on a system
that wasn't originally designed for that particular flexibility. The matrix
and vector
are in a strange state at this point that is not defined
Ugg, I have added this check.
Barry
This bug came from another design extension that was slightly
half-assed.
Originally SNESCreate() always created a KSP and KSPCreate() always
created
a PC and one could not change them. But I realized that PCCreate()
did not
automatically
However, the argument (at source code) that PCGetOperators should
return new Mat's if they were never set never, never convinced me. For
me, Mat's are to KSP/PC like Function/Jacobian to SNES. Stuff that the
user has to provide. Supose you call PCGetOperators(), and then that
call returns a just
On Dec 2, 2008, at 4:02 PM, Lisandro Dalcin wrote:
However, the argument (at source code) that PCGetOperators should
return new Mat's if they were never set never, never convinced me. For
me, Mat's are to KSP/PC like Function/Jacobian to SNES. Stuff that the
user has to provide. Supose you
dded tests_DATAFILESPATH and the variable testexamples_DATAFILESPATH
to handle all the examples that use the matrix files
HG: Enter commit message. Lines beginning with 'HG:' are removed.
HG: --
HG: user: bsmith at barry-smiths-macbook-pro-17.local
HG: branch 'default'
HG: changed
18 matches
Mail list logo