PetscSolverFeti.h 2.22 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
// ============================================================================
// ==                                                                        ==
// == AMDiS - Adaptive multidimensional simulations                          ==
// ==                                                                        ==
// ==  http://www.amdis-fem.org                                              ==
// ==                                                                        ==
// ============================================================================
//
// Software License for AMDiS
//
// Copyright (c) 2010 Dresden University of Technology 
// All rights reserved.
// Authors: Simon Vey, Thomas Witkowski et al.
//
// This file is part of AMDiS
//
// See also license.opensource.txt in the distribution.



/** \file PetscSolverFeti.h */

#include "parallel/PetscSolver.h"
24
#include "parallel/ParallelTypes.h"
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44

#ifndef AMDIS_PETSC_SOLVER_FETI_H
#define AMDIS_PETSC_SOLVER_FETI_H

namespace AMDiS {

  using namespace std;


#ifdef HAVE_PETSC_DEV
  class PetscSolverFeti : public PetscSolver
  {
  public:
    PetscSolverFeti() 
      : PetscSolver()
    {}

    void fillPetscMatrix(Matrix<DOFMatrix*> *mat, SystemVector *vec);

    void solvePetscMatrix(SystemVector &vec, AdaptInfo *adaptInfo);
45
46
47
48
49
50
51
52
53

    Flag getBoundaryDofRequirement()
    {
      return 
	MeshDistributor::BOUNDARY_SUBOBJ_SORTED |
	MeshDistributor::BOUNDARY_FILL_INFO_SEND_DOFS |
	MeshDistributor::BOUNDARY_FILL_INFO_RECV_DOFS;
    }
  protected:
54
    void updateDofData();
55

56
57
58
59
60
61
62
63
64
    void createPrimals();

    void createDuals();

    void createLagrange();

    void createIndexB();

    void createMatLagrange(int nComponents);
65
66

  protected:
67
    DofIndexSet primals;
68
69

    DofMapping globalPrimalIndex;
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

    int nRankPrimals;

    int nOverallPrimals;

    DofIndexSet duals;
    
    DofMapping globalDualIndex;

    DofIndexToPartitions boundaryDofRanks;

    DofMapping dofFirstLagrange;

    int nRankLagrange;

    int nOverallLagrange;
    
    DofMapping globalIndexB;

    int nRankB;

    int nOverallB;
    
    int rStartB;    

    Mat mat_b_b;

    Mat mat_primal_primal;

    Mat mat_b_primal, mat_primal_b;

    Mat mat_lagrange;

    Vec vec_b;

    Vec vec_primal;
106
107
108
109
110
111
  };
#endif

}

#endif