Re: strict static pivoting vs pliable static pivoting

Florin Dobrian <[email protected]> Wed, 13 Sep 2006 11:37:48 -0400 (EDT)
Newsgroups gmane.comp.mathematics.csc
Message-ID <[email protected]>
i like "restricted  pivoting" in order to describe procedures that do not
change supernodes.

one problem with "static pivoting" is the granularity of the term. does
one imply a particular pivot or the whole factorization? if just one pivot
then "static pivoting" is really obvious. but for the whole factorization
it could be ambiguous. if someone tells me that the whole factorization
uses static pivoting then i would immediately assume that every single
pivot was handled statically. one probably has to be very specific in such
a case. otherwise, "restricted pivoting" does the job.

as for the word dynamic, i know it is redundant but i always prefer
redundancy over ambiguity. no harm in saying dynamic pivoting when there
is some static pivoting around as well. i could say that during a
particular factorization i used static pivoting in column x and dynamic
pivoting in column y. instead of saying that i used static pivoting in
column x and pivoting in column y. the former sounds much better, in my
opinion. of course, if there is no static pivoting involved then i would
never use the term "dynamic pivoting".

the term "delayed pivots" seems to be the right one for describing columns
that are moved from on supernode to its parent, changing thus supernodal
data structures (something not allowed in restricted pivoting). would
"delayed pivoting" be the opposite of "restricted pivoting"? or shall
we use "unrestricted pivoting"?

as for preprocessing, as i said, i would never consider that part of
static pivoting since it could be used without any static pivoting at
all, just to reduce the number of swaps and delayed pivots. i like to
say that there is preprocessing for structural reasons (like fill reducing
orderings) and preprocessing for numerical reasons (like the matching
based codes such as mc64).

the word preassigned does not seem to be a good choice. it is ok if every
single pivot is handled statically. but one can have a scenario where
preassigned does not really describe what i going on. for example column j
might be in position x, within a particular supernode, preassigned by
the preprocessing procedure. let's say that the factorization uses
restricted pivoting. now, when we get to position x we may actually decide
to use dynamic pivoting and bring some other column k, from the same
supernode, into position x. then column j would end up in some other
position y, still within the same supernode, but no longer preassigned.
later, when we get to position y, we may no longer have a choice and
use static pivoting. i don't know if any code does that but it is a
possibility.

florin

On Wed, 13 Sep 2006, Tim Davis wrote:

> I'm pliable enough to accept the rejection of "pliable static pivoting"
> in favor of "restricted pivoting"  :-)
>
> So we could have "static pivoting" to apply to the case
> where no pivoting occurs at all once the pivots are selected during the analysis
> phase.  But aren't there current uses of that phrase to include some minor
> variations on the pivot order determined during the analysis?
>
> Then "restricted pivoting" would be (1) pivot in the analysis phase for reducing
> fill-in and for finding large pivots a priori, followed by (2) minor modifications
> to that pivot order during numerical factorization.
>
> Now, what do you call it when a pivot fails in one frontal matrix, but is
> recovered in, say, its grandparent.  We already call this "delayed pivoting", correct?
> This sounds very much like "restricted pivoting" however.
>
> Thanks,
> Tim
>
> p.s.  Speaking of the naming of pivoting schemes, what do you call a fill-reducing
> ordering followed by a call to a template sparse solver (http://www.netlib.org/templates/),
> when done after 6pm, and when used in the solution of a CFD problem involving argon or
> another one of the noble gases?
>
> answer: "a noble night of the minimum degree order of the templates"
>
> Cleve Ashcraft wrote:
>
>  On Tuesday 12 September 2006 13:33, Tim Davis wrote:
>
>
>  Actually, the fact that there may be row/col swaps hidden inside
> a static scheme fits the word nicely.
>
>
>  Pivoting within a supernode could be called "restricted pivoting".
> _______________________________________________
> Csc mailing list
> Csc-lVPF2Jp+hdv2fBVCVOL8/[email protected]
> http://list.odu.edu/listinfo/csc
>
>
>
>
>