ShyLU  Version of the Day
shylu_symbolic.h
1 
2 //@HEADER
3 // ************************************************************************
4 //
5 // ShyLU: Hybrid preconditioner package
6 // Copyright 2012 Sandia Corporation
7 //
8 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
9 // the U.S. Government retains certain rights in this software.
10 //
11 // Redistribution and use in source and binary forms, with or without
12 // modification, are permitted provided that the following conditions are
13 // met:
14 //
15 // 1. Redistributions of source code must retain the above copyright
16 // notice, this list of conditions and the following disclaimer.
17 //
18 // 2. Redistributions in binary form must reproduce the above copyright
19 // notice, this list of conditions and the following disclaimer in the
20 // documentation and/or other materials provided with the distribution.
21 //
22 // 3. Neither the name of the Corporation nor the names of the
23 // contributors may be used to endorse or promote products derived from
24 // this software without specific prior written permission.
25 //
26 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37 //
38 // Questions? Contact Michael A. Heroux (maherou@sandia.gov)
39 //
40 // ************************************************************************
41 //@HEADER
42 
43 #ifndef SHYLU_SYMBOLIC_H
44 #define SHYLU_SYMBOLIC_H
45 
46 #include "Epetra_CrsMatrix.h"
47 #include "Epetra_Map.h"
48 #include "Epetra_MultiVector.h"
49 #include "Epetra_LinearProblem.h"
50 #include "Isorropia_EpetraProber.hpp"
51 //#include "EpetraExt_Transpose_RowMatrix.h"
52 #include <EpetraExt_Reindex_LinearProblem2.h>
53 #include "Ifpack_Preconditioner.h"
54 
55 // This is NOT just the symbolic structure, needs a better name
56 typedef struct
57 {
58  Teuchos::RCP<Epetra_CrsMatrix> D; // D Matrix
59  //Teuchos::RCP<Epetra_CrsMatrix> DT; // D Transpose Matrix
60  //Teuchos::RCP<EpetraExt::RowMatrix_Transpose> transposer;
61  Teuchos::RCP<Epetra_CrsMatrix> C; // Column separator
62  Teuchos::RCP<Epetra_CrsMatrix> R; // Row separator
63  Teuchos::RCP<Epetra_CrsMatrix> G; // G Matrix (A22 block)
64  Teuchos::RCP<Epetra_LinearProblem> LP; // Local problem to solve D
65  Teuchos::RCP<Epetra_LinearProblem> OrigLP; // Local problem to solve D
66  Teuchos::RCP<EpetraExt::ViewTransform<Epetra_LinearProblem> > ReIdx_LP ;
67  Teuchos::RCP<Epetra_MultiVector> Dlhs;
68  Teuchos::RCP<Epetra_MultiVector> Drhs;
69  Teuchos::RCP<Epetra_MultiVector> Gvec;
70  Teuchos::RCP<Amesos_BaseSolver> Solver; // Local solver for D
71  Teuchos::RCP<Ifpack_Preconditioner> ifSolver; //Local incomplete preconditioner
72  Teuchos::RCP<Epetra_CrsGraph> Sg; // The approximate graph of S
73  // Graph(S) + few diagonals
74  Teuchos::RCP<Isorropia::Epetra::Prober> prober; // Prober for Sbar
76 
77 
78 #endif // SHYLU_SYMBOLIC_H