merge

Hg-commit: b723e36a33b01b236644b171b2ed49e502be881c

remove execute permissions from sourcefiles

Hg-commit: 51e6d6c8360e0f44051ea90e634ee2c0b115205d

add support MatMult() for Hermitian matrix in mpisbaij format
add MatMultHermitianTranspose() and MatMultHermitianTransposeAdd()

Hg-commit: bd3ee2c07764767e1715e99fb8df83bd35b221a3

Added SeqBAIJ LU factorization routines with Natural Ordering for bs=2,3,4 (revised new data structure)

Hg-commit: a6e02cc316d3c462357b99f997079b49ffd57bdd

changed PetscMap to PetscLayout

Hg-commit: b8f4632bd00d1563dcf8108fd4199169f8fd320b

changed PetscMap to be pointer to _p_PetscMap
changed PetscMapInitialize() to PetscMapCreate() and also have it malloc the space

Hg-commit: 2d1edc265bc1b52b0bfe730b4711ac37a11744fd

Finished ILU Factorization routines with revised new data structure

Hg-commit: e99636f90e9814b7b9a9e351e7ac0941d7fd93b9

add icc(0) with matorderings under new data structure

Hg-commit: 5b99ce0799b371a259538425a489fed8c02ef64a

add comments
cleanup PetscFreeSpaceContiguous_newdatastruct()

Hg-commit: eb40bb9b7a50536f0c3db779c3d4b00117ea3424

bugfix for MatLUFactorSymbolic_SeqAIJ_newdatastruct()

Hg-commit: f200fccfeabc1a15087bc06a9bc6cbc3abcd7109

mv PetscFreeSpaceContiguous_newdatastruct() from aij files into freespace.*

Hg-commit: 3af95641839072d4a1c681049c89fefe2ddfe79c

commit after merge of patch of 3.0.0 for complex numbrers == 0.0

Hg-commit: da7a9994d64b9371e60d96135ce6b091627eca3f

have to test complex numbers with == 0.0 not !

Hg-commit: 0ed12323236f012a7708363132f4e32693d2f1b9

commit after merge of fixes for minlambda in output files from 3.0.0

Hg-commit: a33729fb1502ca9ab82abe42ce283db7faa76ff2

fixed bug of not setting cnt to zero at start of routine

Hg-commit: 5bb239b6a71ddc130d69f4fd874a599d2e44c61a

pull bug fix for The situation is the following: I'm calling MatSetValues with the flag ADD_VALUES and with matrix option MAT_IGNORE_ZERO_ENTRIES. I'm inserting multiple values at once, multiple colu

pull bug fix for The situation is the following: I'm calling MatSetValues with the flag ADD_VALUES and with matrix option MAT_IGNORE_ZERO_ENTRIES. I'm inserting multiple values at once, multiple columns and rows, so I provide a rank-2 matrix of values. As I'm calling this from fortran I'm also using MAT_COLUMN_ORIENTED. Now for provided rows that are not owned by the process, it jumps to mpiaij.c:394 (line numbers as in petsc-dev). On line 399, it checks for zero entries, but only checks the very first entry of the (non-owned) row. If however other entries of that same row are nonzero, the entire row is still dropped. Note that this is independent of row_oriented/column_oriented as line 396 does exactly the same.

If I don't set the option MAT_IGNORE_ZERO_ENTRIES the problem disappears. In that case however we would either have to preallocate substantially more nonzeros, or complicate the matrix assembly in our code by taking out the zero entries ourselves and call MatSetValues for each entry seperately.

Hg-commit: f09f0ab443e582d40fb2b646b2461bcc61f6e436

show more ...

From: s.kramer@imperial.ac.uk
Subject: non-local values being dropped in MatSetValues
Date: July 17, 2009 7:44:09 AM CDT
To: petsc-users@mcs.anl.gov
Reply-To: petsc-users@mcs.anl.gov

The si

From: s.kramer@imperial.ac.uk
Subject: non-local values being dropped in MatSetValues
Date: July 17, 2009 7:44:09 AM CDT
To: petsc-users@mcs.anl.gov
Reply-To: petsc-users@mcs.anl.gov

The situation is the following: I'm calling MatSetValues with the flag ADD_VALUES and with matrix option MAT_IGNORE_ZERO_ENTRIES. I'm inserting multiple values at once, multiple columns and rows, so I provide a rank-2 matrix of values. As I'm calling this from fortran I'm also using MAT_COLUMN_ORIENTED. Now for provided rows that are not owned by the process, it jumps to mpiaij.c:394 (line numbers as in petsc-dev). On line 399, it checks for zero entries, but only checks the very first entry of the (non-owned) row. If however other entries of that same row are nonzero, the entire row is still dropped. Note that this is independent of row_oriented/column_oriented as line 396 does exactly the same.

If I don't set the option MAT_IGNORE_ZERO_ENTRIES the problem disappears. In that case however we would either have to preallocate substantially more nonzeros, or complicate the matrix assembly in our code by taking out the zero entries ourselves and call MatSetValues for each entry seperately.

Hg-commit: d28deff5e1b8c4e99cf2d919d60e6073ac99448f

show more ...

more information in manual page for MatGetColumnNorms()

Hg-commit: 0bc9bf8bb334ee80b6d4a99576c4c3629a96cb4b

added MatGetColumnNorms() as requested

Hg-commit: 5b6f66184e67c7b90ce8f16c13397f5857ca7adc

Use MatGetVecs instead of VecCreateMPI in the matrix interface

Trivial in:
* MatNullSpaceTest
* MatMatSolve_Basic

Slightly deeper:
* MatComputeExplicitOperator
* Should now work with non-square m

Use MatGetVecs instead of VecCreateMPI in the matrix interface

Trivial in:
* MatNullSpaceTest
* MatMatSolve_Basic

Slightly deeper:
* MatComputeExplicitOperator
* Should now work with non-square matrices (but not tested).
* A little faster by using VEC_IGNORE_OFF_PROC_ENTRIES when assemblying columns of the identity.

Hg-commit: fc9228340b61fbb40680760a26a15b33ee58845a

show more ...

changed PetscOptionsName() and PetscOptionsHasName() to return true if option is set, regardless of the value it is set to, include 0 and false
changed most uses of these routines to the Truth vari

changed PetscOptionsName() and PetscOptionsHasName() to return true if option is set, regardless of the value it is set to, include 0 and false
changed most uses of these routines to the Truth varient

Hg-commit: fd6366abd7eb1f3ec925faa0ecd078d5cc025eb2

show more ...

removed -I$PETSC_DIR from make system and changed the corresponding
#include "src/dsdsdsds" to "../src/sdsdsds" in the source code

Hg-commit: 6e3c4cc4b4bfcb366bae9b71afcb5009c1f0429d

From: dalcinl@gmail.com
Subject: about viewer formats
Date: November 3, 2008 7:50:10 AM CST
To: petsc-dev@mcs.anl.gov
Reply-To: petsc-dev@mcs.anl.gov
Barry, I saw you commit removing PETSC_V

From: dalcinl@gmail.com
Subject: about viewer formats
Date: November 3, 2008 7:50:10 AM CST
To: petsc-dev@mcs.anl.gov
Reply-To: petsc-dev@mcs.anl.gov
Barry, I saw you commit removing PETSC_VIEWER_BINARY_NATIVE. This
makes me think that PETSC_VIEWER_ASCII_DEFAULT and
PETSC_VIEWER_BINARY_DEFAULT (which is NOT used in sources) are a bit
redundat, Why not to remove both and add a new PETSC_VIEWER_DEFAULT?
Does it make sense?

Hg-commit: 9957a5666f6211584952873f183a5f14ec823bbe

show more ...

added MatCompositeSetType() MatComposite can now do multiply in addition to add

Hg-commit: 4c5e06c87351afd0630cd1b00305dc0ca4344355