Branches

debo/mat-taij

View source
debo/mat-taij
Author Commit Message Labels Comments Date
250+ commits behind master.
Debojyoti Ghosh
MatTAIJ: Optimized the backward sweep of MatSOR_SeqTAIJ by saving the application of lower triangular during forward sweep.
Debojyoti Ghosh
MatTAIJ: Implemented MatSOR_SeqTAIJ
Jed Brown
KSP ex60: report time steps and Krylov iterations
Debojyoti Ghosh
MATTAIJ: Fixed two compilation warnings resulting from not calling CHKERRQ(ierr).
Jed Brown
KSP ex60: report final time, update output files
Jed Brown
KSP ex60: add centered-difference discretization of advection
Jed Brown
KSP ex60: make new arbitrary-order Gauss method the default
Jed Brown
KSP ex60: add support for Gauss methods of arbitrary order
Jed Brown
KSP ex60: use PetscFunctionList to factor out 2-stage Gauss method
Jed Brown
KSP ex60: change variable name I to Identity
Jed Brown
KSP ex60: fix __FUNCT__, PETSC_NULL -> NULL, and PetscOptionsBegin
Jed Brown
PCPBJacobi: add PCApplyTranspose (not optimized)
Jed Brown
PCPBJacobi: add PCApply_PBJacobi_N (not unrolled/optimized)
Jed Brown
PCPBJacobi: use VecGetArrayRead for unmodified vector
Jed Brown
PCPBJacobi: fix flop counting
Jed Brown
DTGaussQuadrature: fix quadrature weights
Debojyoti Ghosh
MATTAIJ: Setting block sizes for rectangular S,T also
Debojyoti Ghosh
MatConvert_Basic: Setting block sizes so as to handle rectangular blocks as well.
Debojyoti Ghosh
KSP ex60: Updated the test output for ex60
Debojyoti Ghosh
MATTAIJ: For square S,T, set the blocksize for the TAIJ matrix.
Comments 3
Debojyoti Ghosh
MATTAIJ: Updated the InvertBlockDiagonal for TAIJ to account for T, as well as special cases like S,T being NULL, T being identity.
Debojyoti Ghosh
MAT ex171: Added a test for MATTAIJ where T is specified as an identity matrix; and setting S such that most of its values are not 0!
Debojyoti Ghosh
MATTAIJ: Treating T = identity matrix as a special case since it reduces operations
Debojyoti Ghosh
MATTAIJ and MAT ex 171: Minor corrections in the comments
Debojyoti Ghosh
MATTAIJ: Added if conditions to handle cases where either or both of S and T may be NULL (so that arrays full of zeros are not allocated and operated on)
Debojyoti Ghosh
MAT ex171: Using the matrix "small" for ex171 testing (MATTAIJ functions); testing with nprocs 1,2,3,4 and different square and rectangular sizes for S,T
Debojyoti Ghosh
KSP ex60: Bugfix in the index of local uarr in exact solution function; and removing an unused variable
Debojyoti Ghosh
MATTAIJ: Corrected the comments to be consistent with the mathematical definition of the Kronecker product
Debojyoti Ghosh
MATTAIJ: Implemented MatGetRow (Seq and MPI) for MATTAIJ. Commented out the implementations of MatConvert; MATTAIJ will use MatConvert_Basic (which uses MatGetRow).
Comments 3
Debojyoti Ghosh
MATTAIJ: Passing Vec instead of Vec* in MatMult and MatMultAdd
Debojyoti Ghosh
MATTAIJ: Bugfix in MatMult for non-square I in I \otimes S
Debojyoti Ghosh
MAT_MPIAIJ: In MatMPIAIJGetSeqAIJ, check if pointers are NULL before assigning them. This will allow passing NULL pointers for information that is not needed by the calling function
Debojyoti Ghosh
KSP ex60: Reduced the diffusion coefficient of the test run (runex60) so that the initial sine wave does not decay to round-off errors and updated the test output file
Debojyoti Ghosh
KSP ex60: Updated creation of MATTAIJ to match the new interface, Using the extended TAIJ implementation to create the right-hand-side for solving the stages, and for step completion.
Debojyoti Ghosh
MAT ex171: updated ex171 (MATTAIJ tests) to form TAIJ with S and T; and try non-square S,T; Added tests for MATTAIJ where S or T can be NULL
Debojyoti Ghosh
MATTAIJ: Extended the TAIJ type to handle operations of the form [S \otimes I + T \otimes A] where S and T can be rectangular matrices; If NULL passed as S or T to MatCreateTAIJ, then set them to zero.
Jed Brown
MatTAIJ: fix MatInvertBlockDiagonal_SeqTAIJ, ksp ex60 works with pbjacobi
Jed Brown
KSP ex60: fix for Complex and report errors with %G
Debojyoti Ghosh
KSP ex60: Added a KSP example that demonstrates MATTAIJ by solving a constant-coefficient heat equation (1D) with an implicit Runge-Kutta method.
Debojyoti Ghosh
MAT ex171: Added an example 171 to test the MATTAIJ implementation.
Debojyoti Ghosh
MATTAIJ: Handling MatMult and MatMultAdd in TAIJ through a single function
Debojyoti Ghosh
MATTAIJ: Created and implemented with MatMult, MatMultAdd and InvertBlockDiagonal operations