/usr/include/hypre/HYPRE_LinSysCore.h is in libhypre-dev 2.11.1-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 | /*BHEADER**********************************************************************
* Copyright (c) 2008, Lawrence Livermore National Security, LLC.
* Produced at the Lawrence Livermore National Laboratory.
* This file is part of HYPRE. See file COPYRIGHT for details.
*
* HYPRE is free software; you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License (as published by the Free
* Software Foundation) version 2.1 dated February 1999.
*
* $Revision$
***********************************************************************EHEADER*/
// *************************************************************************
// This is the HYPRE implementation of LinearSystemCore.
// *************************************************************************
#ifndef _HYPRE_LinSysCore_h_
#define _HYPRE_LinSysCore_h_
#define HYPRE_FEI_Version() "FEI/HYPRE 2.7.0R1"
// *************************************************************************
// system libraries used
// -------------------------------------------------------------------------
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <assert.h>
#include <math.h>
#ifdef NOFEI
#undef NOFEI
#endif
// *************************************************************************
// FEI-specific include files
// -------------------------------------------------------------------------
#include "HYPRE_FEI_includes.h"
// *************************************************************************
// local enumerations and defines
// -------------------------------------------------------------------------
enum HYsolverID {HYPCG,HYLSICG,HYGMRES,HYFGMRES,HYCGSTAB,HYCGSTABL,HYTFQMR,
HYBICGS,HYSYMQMR,HYAMG,HYSUPERLU,HYSUPERLUX,HYDSUPERLU,
HYY12M,HYAMGE,HYHYBRID};
enum HYpreconID {HYIDENTITY,HYDIAGONAL,HYPILUT,HYPARASAILS,HYBOOMERAMG,HYML,
HYDDILUT,HYPOLY,HYDDICT,HYSCHWARZ,HYEUCLID,HYBLOCK,HYMLI,
HYUZAWA,HYMLMAXWELL,HYAMS,HYSYSPDE,HYDSLU};
#define HYFEI_HIGHMASK (2147483647-255)
#define HYFEI_SPECIALMASK 255
#define HYFEI_SLIDEREDUCE1 256
#define HYFEI_SLIDEREDUCE2 512
#define HYFEI_SLIDEREDUCE3 1024
#define HYFEI_PRINTMAT 2048
#define HYFEI_PRINTREDMAT 4096
#define HYFEI_PRINTSOL 8192
#define HYFEI_DDILUT 16384
#define HYFEI_SCHURREDUCE1 32768
#define HYFEI_SCHURREDUCE2 65536
#define HYFEI_SCHURREDUCE3 131072
#define HYFEI_PRINTFEINFO 262144
#define HYFEI_AMGDEBUG 524288
#define HYFEI_STOPAFTERPRINT 1048576
#define HYFEI_PRINTPARCSRMAT 2097152
#define HYFEI_IMPOSENOBC 4194304
// *************************************************************************
// substructure definition
// -------------------------------------------------------------------------
typedef struct
{
int *EdgeNodeList_;
int *NodeNumbers_;
int numEdges_;
int numLocalNodes_;
int numNodes_;
double *NodalCoord_;
} HYPRE_FEI_AMSData;
// *************************************************************************
// class definition
// -------------------------------------------------------------------------
class HYPRE_LinSysCore
#ifndef NOFEI
: public LinearSystemCore
#endif
{
public:
HYPRE_LinSysCore(MPI_Comm comm);
virtual ~HYPRE_LinSysCore();
// ----------------------------------------------------------------------
// for creating another one, w/o knowing the run-time type of 'this' one.
// ----------------------------------------------------------------------
#ifndef NOFEI
LinearSystemCore* clone();
#endif
// ----------------------------------------------------------------------
// parameters : for setting generic argc/argv style parameters.
// ----------------------------------------------------------------------
int parameters(int numParams, char** params);
// ----------------------------------------------------------------------
// new functions in FEI 1.5 and above (not implemented here)
// ----------------------------------------------------------------------
int setLookup(Lookup& lookup);
int setConnectivities(GlobalID elemBlock, int numElements,
int numNodesPerElem, const GlobalID* elemIDs,
const int* const* connNodes) ;
int setStiffnessMatrices(GlobalID elemBlock, int numElems,
const GlobalID* elemIDs, const double *const *const *stiff,
int numEqnsPerElem, const int *const * eqnIndices);
int setLoadVectors(GlobalID elemBlock, int numElems,
const GlobalID* elemIDs, const double *const * load,
int numEqnsPerElem, const int *const * eqnIndices);
int setMultCREqns(int multCRSetID, int numCRs, int numNodesPerCR,
int** nodeNumbers, int** eqnNumbers, int* fieldIDs,
int* multiplierEqnNumbers);
int setPenCREqns(int penCRSetID, int numCRs, int numNodesPerCR,
int** nodeNumbers, int** eqnNumbers, int* fieldIDs);
// ----------------------------------------------------------------------
// setGlobalOffsets : provide info for initial creation of
// matrix/vector data, Equation numbers are 1-based, and local sets
// of equation numbers are contiguous.
// ----------------------------------------------------------------------
int setGlobalOffsets(int len, int* nodeOffsets, int* eqnOffsets,
int* blkEqnOffsets);
// ----------------------------------------------------------------------
// setMatrixStructure : provide enough info to allocate the matrix --
// i.e., define the structure.
// ----------------------------------------------------------------------
int setMatrixStructure(int** ptColIndices, int* ptRowLengths,
int** blkColIndices, int* blkRowLengths, int* ptRowsPerBlkRow);
// ----------------------------------------------------------------------
// resetMatrixAndVector : don't destroy the structure of the matrix,
// but set the value 's' throughout the matrix and vectors.
// ----------------------------------------------------------------------
int resetMatrixAndVector(double s);
// ----------------------------------------------------------------------
// reset matrix and vector individually
// ----------------------------------------------------------------------
int resetMatrix(double s);
int resetRHSVector(double s);
// ----------------------------------------------------------------------
// sumIntoSystemMatrix:
// this is the primary assembly function. The coefficients 'values'
// are to be accumumlated into (added to any values already in place)
// global (0-based) equation 'row' of the matrix.
// ----------------------------------------------------------------------
int sumIntoSystemMatrix(int numPtRows, const int* ptRows,
int numPtCols, const int* ptCols, int numBlkRows,
const int* blkRows, int numBlkCols, const int* blkCols,
const double* const* values);
// ----------------------------------------------------------------------
// sumIntoSystemMatrix:
// this is the primary assembly function. The coefficients 'values'
// are to be accumumlated into (added to any values already in place)
// global (1-based) [old - 0-based] equation 'row' of the matrix.
// ----------------------------------------------------------------------
int sumIntoSystemMatrix(int numPtRows, const int* ptRows,
int numPtCols, const int* ptCols,
const double* const* values);
// ----------------------------------------------------------------------
// Point-entry matrix data as for 'sumIntoSystemMatrix', but in this case
// the data should be "put" into the matrix (i.e., overwrite any coefficients
// already present) rather than being "summed" into the matrix.
// ----------------------------------------------------------------------
int putIntoSystemMatrix(int numPtRows, const int* ptRows, int numPtCols,
const int* ptCols, const double* const* values);
// ----------------------------------------------------------------------
// Get the length of a row of the matrix.
// ----------------------------------------------------------------------
int getMatrixRowLength(int row, int& length);
// ----------------------------------------------------------------------
// Obtain the coefficients and indices for a row of the matrix.
// ----------------------------------------------------------------------
int getMatrixRow(int row, double* coefs, int* indices, int len,
int& rowLength);
// ----------------------------------------------------------------------
// sumIntoRHSVector:
// this is the rhs vector equivalent to sumIntoSystemMatrix above.
// ----------------------------------------------------------------------
int sumIntoRHSVector(int num, const double* values, const int* indices);
// ----------------------------------------------------------------------
// For putting coefficients into the rhs vector
// ----------------------------------------------------------------------
int putIntoRHSVector(int num, const double* values, const int* indices);
// ----------------------------------------------------------------------
// For getting coefficients out of the rhs vector
// ----------------------------------------------------------------------
int getFromRHSVector(int num, double* values, const int* indices);
// ----------------------------------------------------------------------
// matrixLoadComplete:
// do any internal synchronization/communication.
// ----------------------------------------------------------------------
int matrixLoadComplete();
// ----------------------------------------------------------------------
// Pass nodal data that probably doesn't mean anything to the FEI
// implementation, but may mean something to the linear solver. Examples:
// geometric coordinates, nullspace data, etc.
// ----------------------------------------------------------------------
int putNodalFieldData(int fieldID, int fieldSize, int* nodeNumbers,
int numNodes, const double* data);
// ----------------------------------------------------------------------
// functions for enforcing boundary conditions.
// ----------------------------------------------------------------------
int enforceEssentialBC(int* globalEqn,double* alpha,double* gamma,int len);
int enforceRemoteEssBCs(int numEqns, int* globalEqns, int** colIndices,
int* colIndLen, double** coefs);
int enforceOtherBC(int* globalEqn, double* alpha, double* beta,
double* gamma, int len);
// ----------------------------------------------------------------------
//functions for getting/setting matrix or vector pointers.
// ----------------------------------------------------------------------
// ----------------------------------------------------------------------
// getMatrixPtr:
// obtain a pointer to the 'A' matrix. This should be considered a
// constant pointer -- i.e., this class remains responsible for the
// matrix (e.g., de-allocation upon destruction).
// ----------------------------------------------------------------------
#ifndef NOFEI
int getMatrixPtr(Data& data);
#endif
// ----------------------------------------------------------------------
// copyInMatrix:
// replaces the internal matrix with a copy of the input argument, scaled
// by the coefficient 'scalar'.
// ----------------------------------------------------------------------
#ifndef NOFEI
int copyInMatrix(double scalar, const Data& data);
#endif
// ----------------------------------------------------------------------
// copyOutMatrix:
// passes out a copy of the internal matrix, scaled by the coefficient
// 'scalar'.
// ----------------------------------------------------------------------
#ifndef NOFEI
int copyOutMatrix(double scalar, Data& data);
#endif
// ----------------------------------------------------------------------
// sumInMatrix:
// accumulate (sum) a copy of the input argument into the internal
// matrix, scaling the input by the coefficient 'scalar'.
// ----------------------------------------------------------------------
#ifndef NOFEI
int sumInMatrix(double scalar, const Data& data);
#endif
// ----------------------------------------------------------------------
// get/setRHSVectorPtr:
// the same semantics apply here as for the matrixPtr functions above.
// ----------------------------------------------------------------------
#ifndef NOFEI
int getRHSVectorPtr(Data& data);
#endif
// ----------------------------------------------------------------------
// copyInRHSVector/copyOutRHSVector/sumInRHSVector:
// the same semantics apply here as for the matrix functions above.
// ----------------------------------------------------------------------
#ifndef NOFEI
int copyInRHSVector(double scalar, const Data& data);
int copyOutRHSVector(double scalar, Data& data);
int sumInRHSVector(double scalar, const Data& data);
#endif
// ----------------------------------------------------------------------
// destroyMatrixData/destroyVectorData:
// Utility function for destroying the matrix (or vector) in Data
// ----------------------------------------------------------------------
#ifndef NOFEI
int destroyMatrixData(Data& data);
int destroyVectorData(Data& data);
#endif
// ----------------------------------------------------------------------
// functions for managing multiple rhs vectors
// ----------------------------------------------------------------------
int setNumRHSVectors(int numRHSs, const int* rhsIDs);
// ----------------------------------------------------------------------
// setRHSID:
// set the 'current' rhs context, assuming multiple rhs vectors.
// ----------------------------------------------------------------------
int setRHSID(int rhsID);
// ----------------------------------------------------------------------
// putInitialGuess:
// function for setting (a subset of) the initial-guess
// solution values (i.e., in the 'x' vector).
// ----------------------------------------------------------------------
int putInitialGuess(const int* eqnNumbers, const double* values,int len);
// ----------------------------------------------------------------------
// function for getting all of the answers ('x' vector).
// ----------------------------------------------------------------------
int getSolution(double* answers, int len);
// ----------------------------------------------------------------------
// function for getting the (single) entry at equation number 'eqnNumber'.
// ----------------------------------------------------------------------
int getSolnEntry(int eqnNumber, double& answer);
// ----------------------------------------------------------------------
// This will be called to request that LinearSystemCore form the residual
// vector r = b - A*x, and pass the coefficients for r back out in the
// 'values' list.
// ----------------------------------------------------------------------
int formResidual(double* values, int len);
double HYPRE_LSC_GetRNorm();
// ----------------------------------------------------------------------
// function for launching the linear solver
// ----------------------------------------------------------------------
int launchSolver(int& solveStatus, int& iterations);
// ----------------------------------------------------------------------
// other functions
// ----------------------------------------------------------------------
int writeSystem(const char *);
// ----------------------------------------------------------------------
// old functions before FEI 1.5 (but still needed here)
// ----------------------------------------------------------------------
int createMatricesAndVectors(int numGlobalEqns, int firstLocalEqn,
int numLocalEqns);
int allocateMatrix(int** colIndices, int* rowLengths);
int sumIntoSystemMatrix(int row, int numValues, const double* values,
const int* scatterIndices);
// ----------------------------------------------------------------------
// HYPRE-specific public functions
// ----------------------------------------------------------------------
void loadConstraintNumbers(int length, int *list);
char *getVersion();
void beginCreateMapFromSoln();
void endCreateMapFromSoln();
void putIntoMappedMatrix(int row, int numValues, const double* values,
const int* scatterIndices);
void getFEDataObject(void **object) { (*object) = feData_; }
int HYPRE_LSC_Matvec(void *x, void *y);
int HYPRE_LSC_Axpby(double, void *, double, void *);
void *HYPRE_LSC_GetRHSVector();
void *HYPRE_LSC_GetSolVector();
void *HYPRE_LSC_GetMatrix();
void *HYPRE_LSC_SetColMap(int, int);
void *HYPRE_LSC_MatMatMult(void *);
// ----------------------------------------------------------------------
// MLI-specific public functions
// ----------------------------------------------------------------------
void FE_initFields(int nFields, int *fieldSizes, int *fieldIDs);
void FE_initElemBlock(int nElems, int nNodesPerElem, int numNodeFields,
int *nodeFieldIDs);
void FE_initElemNodeList(int elemID, int nNodesPerElem, int *nodeIDs);
void FE_initSharedNodes(int nShared, int *sharedIDs, int *sharedPLengs,
int **sharedProcs);
void FE_initComplete();
void FE_loadElemMatrix(int elemID, int nNodes, int *elemNodeList,
int matDim, double **elemMat);
private: //functions
// ----------------------------------------------------------------------
// HYPRE specific private functions
// ----------------------------------------------------------------------
void setupPCGPrecon();
void setupLSICGPrecon();
void setupGMRESPrecon();
void setupFGMRESPrecon();
void setupBiCGSTABPrecon();
void setupBiCGSTABLPrecon();
void setupTFQmrPrecon();
void setupBiCGSPrecon();
void setupSymQMRPrecon();
void setupPreconBoomerAMG();
void setupPreconParaSails();
void setupPreconDDICT();
void setupPreconDDILUT();
void setupPreconPILUT();
void setupPreconPoly();
void setupPreconSchwarz();
void setupPreconML();
void setupPreconMLMaxwell();
void setupPreconAMS();
void setupPreconBlock();
void setupPreconEuclid();
void setupPreconSysPDE();
void solveUsingBoomeramg(int&);
double solveUsingSuperLU(int&);
double solveUsingSuperLUX(int&);
double solveUsingDSuperLU(int&);
void solveUsingY12M(int&);
void solveUsingAMGe(int&);
void buildSlideReducedSystem();
void buildSlideReducedSystem2();
double buildSlideReducedSoln();
double buildSlideReducedSoln2();
void buildSchurReducedSystem();
void buildSchurReducedSystem2();
void buildSlideReducedSystemPartA(int*,int*,int,int,int*,int*);
void buildSlideReducedSystemPartB(int*,int*,int,int,int*,int*,
HYPRE_ParCSRMatrix *);
void buildSlideReducedSystemPartC(int*,int*,int,int,int*,int*,
HYPRE_ParCSRMatrix);
void buildSchurReducedRHS();
void buildSchurInitialGuess();
double buildSchurReducedSoln();
void computeAConjProjection(HYPRE_ParCSRMatrix A_csr, HYPRE_ParVector x,
HYPRE_ParVector b);
void computeMinResProjection(HYPRE_ParCSRMatrix A_csr, HYPRE_ParVector x,
HYPRE_ParVector b);
void addToAConjProjectionSpace(HYPRE_IJVector x, HYPRE_IJVector b);
void addToMinResProjectionSpace(HYPRE_IJVector x, HYPRE_IJVector b);
int HYPRE_Schur_Search(int,int,int*,int*,int,int);
void HYPRE_LSI_BuildNodalCoordinates();
// ----------------------------------------------------------------------
// private functions for selecting solver/preconditioner
// ----------------------------------------------------------------------
void selectSolver(char* name);
void selectPreconditioner(char* name);
private: //variables
// ----------------------------------------------------------------------
// parallel communication information and output levels
// ----------------------------------------------------------------------
MPI_Comm comm_; // MPI communicator
int numProcs_; // number of processors
int mypid_; // my processor ID
int HYOutputLevel_; // control print information
int memOptimizerFlag_; // turn on memory optimizer
// ----------------------------------------------------------------------
// for storing information about how to load matrix directly (bypass FEI)
// ----------------------------------------------------------------------
int mapFromSolnFlag_;
int mapFromSolnLeng_;
int mapFromSolnLengMax_;
int *mapFromSolnList_;
int *mapFromSolnList2_;
// ----------------------------------------------------------------------
// matrix and vectors
// ----------------------------------------------------------------------
HYPRE_IJMatrix HYA_; // the system matrix
HYPRE_IJMatrix HYnormalA_; // normalized system matrix
HYPRE_IJVector HYb_; // the current RHS
HYPRE_IJVector HYnormalB_; // normalized system rhs
HYPRE_IJVector *HYbs_; // an array of RHSs
HYPRE_IJVector HYx_; // the solution vector
HYPRE_IJVector HYr_; // temporary vector for residual
HYPRE_IJVector *HYpxs_; // an array of previous solutions
HYPRE_IJVector *HYpbs_; // an array of previous rhs
int numGlobalRows_;
int localStartRow_;
int localEndRow_;
int localStartCol_;
int localEndCol_;
int *rowLengths_;
int **colIndices_;
double **colValues_;
double truncThresh_;
double rnorm_;
// ----------------------------------------------------------------------
// matrix and vectors for reduction
// ----------------------------------------------------------------------
HYPRE_IJMatrix reducedA_; // matrix for reduction
HYPRE_IJVector reducedB_; // RHS vector for reduction
HYPRE_IJVector reducedX_; // solution vector for reduction
HYPRE_IJVector reducedR_; // temporary vector for reduction
HYPRE_IJMatrix HYA21_; // (2,1) block in reduction
HYPRE_IJMatrix HYA12_; // (1,2) block in reduction
int A21NRows_; // number of rows in (2,1) block
int A21NCols_; // number of cols in (2,1) block
int reducedAStartRow_; // Nrows in reduced system
HYPRE_IJMatrix HYinvA22_; // inv(A22) in slide reduction
// ----------------------------------------------------------------------
// pointers to current matrix and vectors for solver
// ----------------------------------------------------------------------
HYPRE_IJMatrix currA_;
HYPRE_IJVector currB_;
HYPRE_IJVector currX_;
HYPRE_IJVector currR_;
int currentRHS_;
int *rhsIDs_;
int numRHSs_;
int nStored_;
int *storedIndices_;
int *auxStoredIndices_;
int mRHSFlag_;
int mRHSNumGEqns_;
int *mRHSGEqnIDs_;
int *mRHSNEntries_;
int *mRHSBCType_;
int **mRHSRowInds_;
double **mRHSRowVals_;
// ----------------------------------------------------------------------
// flags for matrix assembly, various reductions, and projections
// ----------------------------------------------------------------------
int matrixVectorsCreated_;
int systemAssembled_;
int slideReduction_;
double slideReductionMinNorm_;
int slideReductionScaleMatrix_;
int schurReduction_;
int schurReductionCreated_;
int projectionScheme_;
int projectSize_;
int projectCurrSize_;
double **projectionMatrix_;
int normalEqnFlag_;
void *slideObj_;
// ----------------------------------------------------------------------
// variables for slide and Schur reduction
// ----------------------------------------------------------------------
int *selectedList_;
int *selectedListAux_;
int nConstraints_;
int *constrList_;
int matrixPartition_;
// ----------------------------------------------------------------------
// variables for the selected solver and preconditioner
// ----------------------------------------------------------------------
char *HYSolverName_;
HYPRE_Solver HYSolver_;
HYsolverID HYSolverID_;
int gmresDim_;
int fgmresUpdateTol_;
int maxIterations_;
double tolerance_;
int normAbsRel_;
int pcgRecomputeRes_;
char *HYPreconName_;
HYPRE_Solver HYPrecon_;
HYpreconID HYPreconID_;
int HYPreconReuse_;
int HYPreconSetup_;
// ----------------------------------------------------------------------
// preconditioner specific variables
// ----------------------------------------------------------------------
int amgMaxLevels_;
int amgCoarsenType_;
int amgMaxIter_;
int amgMeasureType_;
int amgNumSweeps_[4];
int amgRelaxType_[4];
int amgGridRlxType_;
double amgRelaxWeight_[25];
double amgRelaxOmega_[25];
double amgStrongThreshold_;
int amgSystemSize_;
int amgSmoothType_;
int amgSmoothNumLevels_;
int amgSmoothNumSweeps_;
int amgCGSmoothNumSweeps_;
double amgSchwarzRelaxWt_;
int amgSchwarzVariant_;
int amgSchwarzOverlap_;
int amgSchwarzDomainType_;
int amgUseGSMG_;
int amgGSMGNSamples_;
int amgAggLevels_;
int amgInterpType_;
int amgPmax_;
int pilutFillin_;
double pilutDropTol_;
int pilutMaxNnzPerRow_;
int parasailsSym_;
double parasailsThreshold_;
int parasailsNlevels_;
double parasailsFilter_;
double parasailsLoadbal_;
int parasailsReuse_;
int mlMethod_;
int mlNumPreSweeps_;
int mlNumPostSweeps_;
int mlPresmootherType_;
int mlPostsmootherType_;
double mlRelaxWeight_;
double mlStrongThreshold_;
int mlCoarseSolver_;
int mlCoarsenScheme_;
int mlNumPDEs_;
int superluOrdering_;
char superluScale_[1];
double ddilutFillin_;
double ddilutDropTol_;
int ddilutOverlap_;
int ddilutReorder_;
double ddictFillin_;
double ddictDropTol_;
double schwarzFillin_;
int schwarzNblocks_;
int schwarzBlksize_;
int polyOrder_;
int euclidargc_;
char **euclidargv_;
HYPRE_IJVector amsX_;
HYPRE_IJVector amsY_;
HYPRE_IJVector amsZ_;
int localStartRowAMSV_;
int localEndRowAMSV_;
HYPRE_IJMatrix amsG_;
HYPRE_IJMatrix amsD0_;
HYPRE_IJMatrix amsD1_;
int localStartRowAMSG_;
int localEndRowAMSG_;
int localStartColAMSG_;
int localEndColAMSG_;
HYPRE_ParCSRMatrix amsBetaPoisson_;
int amsNumPDEs_;
int amsMaxIter_;
double amsTol_;
int amsCycleType_;
int amsRelaxType_;
int amsRelaxTimes_;
double amsRelaxWt_;
double amsRelaxOmega_;
int amsPrintLevel_;
int amsAlphaCoarsenType_;
int amsAlphaAggLevels_;
int amsAlphaRelaxType_;
double amsAlphaStrengthThresh_;
int amsAlphaInterpType_;
int amsAlphaPmax_;
int amsBetaCoarsenType_;
int amsBetaAggLevels_;
int amsBetaRelaxType_;
double amsBetaStrengthThresh_;
int amsBetaInterpType_;
int amsBetaPmax_;
int sysPDEMethod_;
int sysPDEFormat_;
double sysPDETol_;
int sysPDEMaxIter_;
int sysPDENumPre_;
int sysPDENumPost_;
int sysPDENVars_;
// ----------------------------------------------------------------------
// FEI and MLI variables
// ----------------------------------------------------------------------
void *feData_;
int haveFEData_;
Lookup *lookup_;
int haveLookup_;
int MLI_NumNodes_;
int MLI_FieldSize_;
int *MLI_EqnNumbers_;
double *MLI_NodalCoord_;
int MLI_Hybrid_NSIncr_;
int MLI_Hybrid_GSA_;
int MLI_Hybrid_MaxIter_;
double MLI_Hybrid_ConvRate_;
int MLI_Hybrid_NTrials_;
HYPRE_FEI_AMSData AMSData_;
int FEI_mixedDiagFlag_;
double *FEI_mixedDiag_;
// ----------------------------------------------------------------------
// ML Maxwell variables
// ----------------------------------------------------------------------
HYPRE_ParCSRMatrix maxwellANN_; // Maxwell nodal matrix
HYPRE_ParCSRMatrix maxwellGEN_; // Maxwell gradient matrix
// ----------------------------------------------------------------------
// temporary functions for testing purposes
// ----------------------------------------------------------------------
friend void fei_hypre_test(int argc, char *argv[]);
friend void fei_hypre_domaindecomposition(int argc, char *argv[]);
};
#endif
|