MueLu  Version of the Day
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
MueLu_RefMaxwell_decl.hpp
Go to the documentation of this file.
1 // @HEADER
2 // *****************************************************************************
3 // MueLu: A package for multigrid based preconditioning
4 //
5 // Copyright 2012 NTESS and the MueLu contributors.
6 // SPDX-License-Identifier: BSD-3-Clause
7 // *****************************************************************************
8 // @HEADER
9 
10 #ifndef MUELU_REFMAXWELL_DECL_HPP
11 #define MUELU_REFMAXWELL_DECL_HPP
12 
13 #include "MueLu_ConfigDefs.hpp"
14 #include "MueLu_BaseClass.hpp"
15 
16 #include "MueLu_CoalesceDropFactory_fwd.hpp"
17 #include "MueLu_CoarseMapFactory_fwd.hpp"
18 #include "MueLu_CoordinatesTransferFactory_fwd.hpp"
19 #include "MueLu_TentativePFactory_fwd.hpp"
20 #include "MueLu_SaPFactory_fwd.hpp"
21 #include "MueLu_UncoupledAggregationFactory_fwd.hpp"
22 #include "MueLu_AggregationExportFactory_fwd.hpp"
23 #include "MueLu_Utilities_fwd.hpp"
24 
25 #include "MueLu_MasterList.hpp"
26 #include "MueLu_FactoryManager.hpp"
27 #include "MueLu_SmootherBase_fwd.hpp"
28 #include "MueLu_SmootherPrototype_fwd.hpp"
29 #include "MueLu_PerfUtils_fwd.hpp"
30 #include "MueLu_RAPFactory_fwd.hpp"
31 
32 #include "MueLu_CoalesceDropFactory_kokkos_fwd.hpp"
33 #include "MueLu_TentativePFactory_kokkos_fwd.hpp"
34 
35 #include "MueLu_ZoltanInterface_fwd.hpp"
36 #include "MueLu_Zoltan2Interface_fwd.hpp"
37 #include "MueLu_RepartitionHeuristicFactory_fwd.hpp"
38 #include "MueLu_RepartitionFactory_fwd.hpp"
39 #include "MueLu_RebalanceAcFactory_fwd.hpp"
40 #include "MueLu_RebalanceTransferFactory_fwd.hpp"
41 
42 #include "MueLu_SmootherFactory_fwd.hpp"
43 #include "MueLu_TrilinosSmoother_fwd.hpp"
44 #include "MueLu_Hierarchy_fwd.hpp"
45 
46 #include "Xpetra_Operator.hpp"
47 #include "Xpetra_Map_fwd.hpp"
48 #include "Xpetra_Matrix_fwd.hpp"
49 #include "Xpetra_MapFactory_fwd.hpp"
50 #include "Xpetra_MatrixFactory_fwd.hpp"
51 #include "Xpetra_MultiVectorFactory_fwd.hpp"
52 #include "Xpetra_VectorFactory_fwd.hpp"
53 #include "Xpetra_CrsMatrixWrap_fwd.hpp"
54 
55 namespace MueLu {
56 
216 template <class Scalar,
217  class LocalOrdinal,
218  class GlobalOrdinal,
219  class Node>
220 class RefMaxwell : public VerboseObject, public Xpetra::Operator<Scalar, LocalOrdinal, GlobalOrdinal, Node> {
221 #undef MUELU_REFMAXWELL_SHORT
222 #include "MueLu_UseShortNames.hpp"
223 
224  public:
225  typedef typename Teuchos::ScalarTraits<Scalar>::magnitudeType magnitudeType;
226  typedef typename Teuchos::ScalarTraits<Scalar>::coordinateType coordinateType;
227  typedef typename Xpetra::MultiVector<coordinateType, LO, GO, NO> RealValuedMultiVector;
228 
230  RefMaxwell()
231  : HierarchyCoarse11_(Teuchos::null)
232  , Hierarchy22_(Teuchos::null)
233  , disable_addon_(MasterList::getDefault<bool>("refmaxwell: disable addon"))
234  , mode_(MasterList::getDefault<std::string>("refmaxwell: mode")) {
235  }
236 
238  RefMaxwell(Teuchos::RCP<Hierarchy> HH, Teuchos::RCP<Hierarchy> H22)
239  : HierarchyCoarse11_(HH)
240  , Hierarchy22_(H22)
241  , disable_addon_(MasterList::getDefault<bool>("refmaxwell: disable addon"))
242  , mode_(MasterList::getDefault<std::string>("refmaxwell: mode")) {
243  }
244 
245  RefMaxwell(const Teuchos::RCP<Matrix> &SM_Matrix,
246  const Teuchos::RCP<Matrix> &Dk_1,
247  const Teuchos::RCP<Matrix> &Dk_2,
248  const Teuchos::RCP<Matrix> &D0,
249  const Teuchos::RCP<Matrix> &M1_beta,
250  const Teuchos::RCP<Matrix> &M1_alpha,
251  const Teuchos::RCP<Matrix> &Mk_one,
252  const Teuchos::RCP<Matrix> &Mk_1_one,
253  const Teuchos::RCP<Matrix> &invMk_1_invBeta,
254  const Teuchos::RCP<Matrix> &invMk_2_invAlpha,
255  const Teuchos::RCP<MultiVector> &Nullspace11,
256  const Teuchos::RCP<MultiVector> &Nullspace22,
257  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
258  Teuchos::ParameterList &List,
259  bool ComputePrec = true) {
260  int spaceNumber = List.get<int>("refmaxwell: space number", 1);
261  initialize(spaceNumber,
262  Dk_1, Dk_2, D0,
263  M1_beta, M1_alpha,
264  Mk_one, Mk_1_one,
265  invMk_1_invBeta, invMk_2_invAlpha,
266  Nullspace11, Nullspace22, NodalCoords,
267  Teuchos::null, Teuchos::null,
268  List);
269  resetMatrix(SM_Matrix, ComputePrec);
270  }
271 
272  RefMaxwell(const Teuchos::RCP<Matrix> &SM_Matrix,
273  const Teuchos::RCP<Matrix> &Dk_1,
274  const Teuchos::RCP<Matrix> &Dk_2,
275  const Teuchos::RCP<Matrix> &D0,
276  const Teuchos::RCP<Matrix> &M1_beta,
277  const Teuchos::RCP<Matrix> &M1_alpha,
278  const Teuchos::RCP<Matrix> &Mk_one,
279  const Teuchos::RCP<Matrix> &Mk_1_one,
280  const Teuchos::RCP<Matrix> &invMk_1_invBeta,
281  const Teuchos::RCP<Matrix> &invMk_2_invAlpha,
282  const Teuchos::RCP<MultiVector> &Nullspace11,
283  const Teuchos::RCP<MultiVector> &Nullspace22,
284  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
285  const Teuchos::RCP<MultiVector> &Material_beta,
286  const Teuchos::RCP<MultiVector> &Material_alpha,
287  Teuchos::ParameterList &List,
288  bool ComputePrec = true) {
289  int spaceNumber = List.get<int>("refmaxwell: space number", 1);
290  initialize(spaceNumber,
291  Dk_1, Dk_2, D0,
292  M1_beta, M1_alpha,
293  Mk_one, Mk_1_one,
294  invMk_1_invBeta, invMk_2_invAlpha,
295  Nullspace11, Nullspace22, NodalCoords,
296  Material_beta, Material_alpha,
297  List);
298  resetMatrix(SM_Matrix, ComputePrec);
299  }
300 
313  RefMaxwell(const Teuchos::RCP<Matrix> &SM_Matrix,
314  const Teuchos::RCP<Matrix> &D0_Matrix,
315  const Teuchos::RCP<Matrix> &Ms_Matrix,
316  const Teuchos::RCP<Matrix> &M0inv_Matrix,
317  const Teuchos::RCP<Matrix> &M1_Matrix,
318  const Teuchos::RCP<MultiVector> &Nullspace11,
319  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
320  Teuchos::ParameterList &List,
321  bool ComputePrec = true) {
322  initialize(D0_Matrix, Ms_Matrix, M0inv_Matrix, M1_Matrix, Nullspace11, NodalCoords, Teuchos::null, List);
323  resetMatrix(SM_Matrix, ComputePrec);
324  }
325 
326  RefMaxwell(const Teuchos::RCP<Matrix> &SM_Matrix,
327  const Teuchos::RCP<Matrix> &D0_Matrix,
328  const Teuchos::RCP<Matrix> &Ms_Matrix,
329  const Teuchos::RCP<Matrix> &M0inv_Matrix,
330  const Teuchos::RCP<Matrix> &M1_Matrix,
331  const Teuchos::RCP<MultiVector> &Nullspace11,
332  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
333  const Teuchos::RCP<MultiVector> &Material,
334  Teuchos::ParameterList &List,
335  bool ComputePrec = true) {
336  initialize(D0_Matrix, Ms_Matrix, M0inv_Matrix, M1_Matrix, Nullspace11, NodalCoords, Material, List);
337  resetMatrix(SM_Matrix, ComputePrec);
338  }
339 
351  RefMaxwell(const Teuchos::RCP<Matrix> &SM_Matrix,
352  const Teuchos::RCP<Matrix> &D0_Matrix,
353  const Teuchos::RCP<Matrix> &M0inv_Matrix,
354  const Teuchos::RCP<Matrix> &M1_Matrix,
355  const Teuchos::RCP<MultiVector> &Nullspace11,
356  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
357  Teuchos::ParameterList &List,
358  bool ComputePrec = true) {
359  initialize(D0_Matrix, M1_Matrix, M0inv_Matrix, M1_Matrix, Nullspace11, NodalCoords, Teuchos::null, List);
360  resetMatrix(SM_Matrix, ComputePrec);
361  }
362 
372  RefMaxwell(const Teuchos::RCP<Matrix> &D0_Matrix,
373  const Teuchos::RCP<Matrix> &M0inv_Matrix,
374  const Teuchos::RCP<Matrix> &M1_Matrix,
375  const Teuchos::RCP<MultiVector> &Nullspace11,
376  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
377  Teuchos::ParameterList &List)
378  : SM_Matrix_(Teuchos::null) {
379  initialize(D0_Matrix, M1_Matrix, M0inv_Matrix, M1_Matrix, Nullspace11, NodalCoords, Teuchos::null, List);
380  }
381 
392  RefMaxwell(const Teuchos::RCP<Matrix> &SM_Matrix,
393  const Teuchos::RCP<Matrix> &D0_Matrix,
394  const Teuchos::RCP<Matrix> &M1_Matrix,
395  const Teuchos::RCP<MultiVector> &Nullspace11,
396  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
397  Teuchos::ParameterList &List,
398  bool ComputePrec) {
399  initialize(D0_Matrix, M1_Matrix, Teuchos::null, M1_Matrix, Nullspace11, NodalCoords, Teuchos::null, List);
400  resetMatrix(SM_Matrix, ComputePrec);
401  }
402 
411  RefMaxwell(const Teuchos::RCP<Matrix> &D0_Matrix,
412  const Teuchos::RCP<Matrix> &M1_Matrix,
413  const Teuchos::RCP<MultiVector> &Nullspace11,
414  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
415  Teuchos::ParameterList &List)
416  : SM_Matrix_(Teuchos::null) {
417  initialize(D0_Matrix, M1_Matrix, Teuchos::null, M1_Matrix, Nullspace11, NodalCoords, Teuchos::null, List);
418  }
419 
426  RefMaxwell(const Teuchos::RCP<Matrix> &SM_Matrix,
427  Teuchos::ParameterList &List,
428  bool ComputePrec = true);
429 
431  virtual ~RefMaxwell() {}
432 
434  const Teuchos::RCP<const Map> getDomainMap() const;
435 
437  const Teuchos::RCP<const Map> getRangeMap() const;
438 
440  const Teuchos::RCP<Matrix> &getJacobian() const {
441  return SM_Matrix_;
442  }
443 
445  void setParameters(Teuchos::ParameterList &list);
446 
448  void compute(bool reuse = false);
449 
451  void resetMatrix(Teuchos::RCP<Matrix> SM_Matrix_new, bool ComputePrec = true);
452 
456  void apply(const MultiVector &X, MultiVector &Y,
457  Teuchos::ETransp mode = Teuchos::NO_TRANS,
458  Scalar alpha = Teuchos::ScalarTraits<Scalar>::one(),
459  Scalar beta = Teuchos::ScalarTraits<Scalar>::zero()) const;
460 
462  bool hasTransposeApply() const;
463 
464  void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel = Teuchos::VERB_HIGH) const;
465 
467  void residual(const MultiVector &X,
468  const MultiVector &B,
469  MultiVector &R) const {
470  using STS = Teuchos::ScalarTraits<Scalar>;
471  R.update(STS::one(), B, STS::zero());
472  this->apply(X, R, Teuchos::NO_TRANS, -STS::one(), STS::one());
473  }
474 
475  private:
476  Teuchos::RCP<Teuchos::ParameterList> getValidParamterList();
477 
490  void initialize(const Teuchos::RCP<Matrix> &D0_Matrix,
491  const Teuchos::RCP<Matrix> &Ms_Matrix,
492  const Teuchos::RCP<Matrix> &M0inv_Matrix,
493  const Teuchos::RCP<Matrix> &M1_Matrix,
494  const Teuchos::RCP<MultiVector> &Nullspace11,
495  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
496  const Teuchos::RCP<MultiVector> &Material,
497  Teuchos::ParameterList &List);
498 
516  void initialize(const int k,
517  const Teuchos::RCP<Matrix> &Dk_1,
518  const Teuchos::RCP<Matrix> &Dk_2,
519  const Teuchos::RCP<Matrix> &D0,
520  const Teuchos::RCP<Matrix> &M1_beta,
521  const Teuchos::RCP<Matrix> &M1_alpha,
522  const Teuchos::RCP<Matrix> &Mk_one,
523  const Teuchos::RCP<Matrix> &Mk_1_one,
524  const Teuchos::RCP<Matrix> &invMk_1_invBeta,
525  const Teuchos::RCP<Matrix> &invMk_2_invAlpha,
526  const Teuchos::RCP<MultiVector> &Nullspace11,
527  const Teuchos::RCP<MultiVector> &Nullspace22,
528  const Teuchos::RCP<RealValuedMultiVector> &NodalCoords,
529  const Teuchos::RCP<MultiVector> &Material_beta,
530  const Teuchos::RCP<MultiVector> &Material_alpha,
531  Teuchos::ParameterList &List);
532 
534  void determineSubHierarchyCommSizes(bool &doRebalancing, int &rebalanceStriding, int &numProcsCoarseA11, int &numProcsA22);
535 
537  void buildCoarse11Matrix();
538 
540  void rebalanceCoarse11Matrix(const int rebalanceStriding, const int numProcsCoarseA11);
541 
543  void build22Matrix(const bool reuse, const bool doRebalancing, const int rebalanceStriding, const int numProcsA22);
544 
545  public: // due to Cuda build errors otherwise
547  RCP<MultiVector> buildNullspace(const int spaceNumber, const Kokkos::View<bool *, typename Node::device_type> &bcs, const bool applyBCs);
548 
550  Teuchos::RCP<Matrix> buildProjection(const int spaceNumber, const RCP<MultiVector> &EdgeNullspace) const;
551 
558  void buildNodalProlongator(const Teuchos::RCP<Matrix> &A_nodal,
559  Teuchos::RCP<Matrix> &P_nodal,
560  Teuchos::RCP<MultiVector> &Nullspace_nodal,
561  Teuchos::RCP<RealValuedMultiVector> &Coords_nodal) const;
562 
568  RCP<Matrix> buildVectorNodalProlongator(const Teuchos::RCP<Matrix> &P_nodal) const;
569 
579  void buildProlongator(const int spaceNumber,
580  const Teuchos::RCP<Matrix> &A_nodal_Matrix,
581  const RCP<MultiVector> &EdgeNullspace,
582  Teuchos::RCP<Matrix> &edgeProlongator,
583  Teuchos::RCP<MultiVector> &coarseEdgeNullspace,
584  Teuchos::RCP<RealValuedMultiVector> &coarseNodalCoords) const;
585 
590  RCP<Matrix> buildAddon(const int spaceNumber);
591 
592  private:
594  void setupSubSolve(Teuchos::RCP<Hierarchy> &hierarchy,
595  Teuchos::RCP<Operator> &thyraPrecOp,
596  const Teuchos::RCP<Matrix> &A,
597  const Teuchos::RCP<MultiVector> &Nullspace,
598  const Teuchos::RCP<RealValuedMultiVector> &Coords,
599  const Teuchos::RCP<MultiVector> &Material,
600  Teuchos::ParameterList &params,
601  std::string &label,
602  const bool reuse,
603  const bool isSingular = false);
604 
606  void setFineLevelSmoother11();
607 
609  void applyInverseAdditive(const MultiVector &RHS, MultiVector &X) const;
610 
612  void solveH(const MultiVector &RHS, MultiVector &X) const;
613 
615  void solve22(const MultiVector &RHS, MultiVector &X) const;
616 
618  void allocateMemory(int numVectors) const;
619 
621  void dump(const RCP<Matrix> &A, std::string name) const;
622 
624  void dump(const RCP<MultiVector> &X, std::string name) const;
625 
627  void dumpCoords(const RCP<RealValuedMultiVector> &X, std::string name) const;
628 
630  void dump(const Teuchos::ArrayRCP<bool> &v, std::string name) const;
631 
633  void dump(const Kokkos::View<bool *, typename Node::device_type> &v, std::string name) const;
634 
636  Teuchos::RCP<Teuchos::TimeMonitor> getTimer(std::string name, RCP<const Teuchos::Comm<int> > comm = Teuchos::null) const;
637 
639  Teuchos::RCP<Hierarchy> HierarchyCoarse11_, Hierarchy22_;
640  Teuchos::RCP<SmootherBase> PreSmoother11_, PostSmoother11_;
641  Teuchos::RCP<SmootherPrototype> PreSmootherData11_, PostSmootherData11_;
642  RCP<Operator> thyraPrecOpH_, thyraPrecOp22_;
644  int spaceNumber_;
646  std::string solverName_;
648  size_t dim_;
650  Teuchos::RCP<Matrix> SM_Matrix_;
652  Teuchos::RCP<Matrix> Dk_1_, Dk_1_T_;
654  Teuchos::RCP<Matrix> Dk_2_;
656  Teuchos::RCP<Matrix> D0_;
658  Teuchos::RCP<Matrix> invMk_1_invBeta_, invMk_2_invAlpha_;
660  Teuchos::RCP<Matrix> Mk_one_, Mk_1_one_;
662  Teuchos::RCP<Matrix> M1_beta_, M1_alpha_;
664  Teuchos::RCP<MultiVector> Material_beta_, Material_alpha_;
666  Teuchos::RCP<Matrix> P11_, R11_;
668  Teuchos::RCP<Matrix> P22_, R22_;
670  Teuchos::RCP<Matrix> coarseA11_, A22_, coarseA22_;
672  Teuchos::RCP<Matrix> Addon11_;
674  Teuchos::RCP<Matrix> Addon22_;
675  Teuchos::RCP<const Map> DorigDomainMap_;
676  Teuchos::RCP<const Import> DorigImporter_;
678  Kokkos::View<bool *, typename Node::device_type::memory_space> BCrows11_, BCcols22_, BCdomain22_;
679  int globalNumberBoundaryUnknowns11_, globalNumberBoundaryUnknowns22_;
681  Teuchos::RCP<MultiVector> Nullspace11_, Nullspace22_;
683  Teuchos::RCP<RealValuedMultiVector> NodalCoords_, CoordsCoarse11_, Coords22_;
685  Teuchos::RCP<MultiVector> NullspaceCoarse11_;
687  Teuchos::RCP<MultiVector> CoarseNullspace22_;
689  Teuchos::RCP<const Import> ImporterCoarse11_, Importer22_;
690  bool Dk_1_T_R11_colMapsMatch_, asyncTransfers_;
691  bool onlyBoundary11_, onlyBoundary22_;
693  Teuchos::ParameterList parameterList_, precList11_, precList22_;
694  Teuchos::RCP<Teuchos::ParameterList> coarseA11_AP_reuse_data_, coarseA11_RAP_reuse_data_;
695  Teuchos::RCP<Teuchos::ParameterList> A22_AP_reuse_data_, A22_RAP_reuse_data_;
697  bool disable_addon_, disable_addon_22_, dump_matrices_, useKokkos_, use_as_preconditioner_, implicitTranspose_, fuseProlongationAndUpdate_, syncTimers_, enable_reuse_, skipFirst11Level_, skipFirst22Level_;
698  bool applyBCsToAnodal_, applyBCsToCoarse11_, applyBCsTo22_;
699  int numItersCoarse11_, numIters22_;
700  std::string mode_;
701 
703  mutable Teuchos::RCP<MultiVector> P11res_, P11x_, P11resSubComm_, P11xSubComm_, DresIntermediate_, Dres_, DxIntermediate_, Dx_, DresSubComm_, DxSubComm_, residual_, P11resTmp_, DresTmp_, DTR11Tmp_, P11x_colmap_, Dx_colmap_;
704 };
705 
706 } // namespace MueLu
707 
708 #define MUELU_REFMAXWELL_SHORT
709 #endif // MUELU_REFMAXWELL_DECL_HPP
MueLu_SmootherFactory_fwd.hpp
MueLu::RefMaxwell::getJacobian
const Teuchos::RCP< Matrix > & getJacobian() const
Returns Jacobian matrix SM.
Definition: MueLu_RefMaxwell_decl.hpp:440
MueLu::RefMaxwell::Importer22_
Teuchos::RCP< const Import > Importer22_
Definition: MueLu_RefMaxwell_decl.hpp:689
MueLu::RefMaxwell::dim_
size_t dim_
The spatial dimension.
Definition: MueLu_RefMaxwell_decl.hpp:648
MueLu::RefMaxwell::buildAddon
RCP< Matrix > buildAddon(const int spaceNumber)
Definition: MueLu_RefMaxwell_def.hpp:701
MueLu::RefMaxwell::PreSmootherData11_
Teuchos::RCP< SmootherPrototype > PreSmootherData11_
Definition: MueLu_RefMaxwell_decl.hpp:641
MueLu_RebalanceTransferFactory_fwd.hpp
MueLu_TrilinosSmoother_fwd.hpp
MueLu::RefMaxwell::R22_
Teuchos::RCP< Matrix > R22_
Definition: MueLu_RefMaxwell_decl.hpp:668
LocalOrdinal
MueLu::DefaultLocalOrdinal LocalOrdinal
Definition: MueLu_UseDefaultTypes.hpp:13
MueLu::RefMaxwell::BCdomain22_
Kokkos::View< bool *, typename Node::device_type::memory_space > BCdomain22_
Definition: MueLu_RefMaxwell_decl.hpp:678
MueLu::RefMaxwell::solve22
void solve22(const MultiVector &RHS, MultiVector &X) const
apply solve to 2-2 block only
Definition: MueLu_RefMaxwell_def.hpp:2385
MueLu::RefMaxwell::PreSmoother11_
Teuchos::RCP< SmootherBase > PreSmoother11_
Definition: MueLu_RefMaxwell_decl.hpp:640
MueLu::RefMaxwell::A22_RAP_reuse_data_
Teuchos::RCP< Teuchos::ParameterList > A22_RAP_reuse_data_
Definition: MueLu_RefMaxwell_decl.hpp:695
MueLu::RefMaxwell::dumpCoords
void dumpCoords(const RCP< RealValuedMultiVector > &X, std::string name) const
dump out real-valued multivector
Definition: MueLu_RefMaxwell_def.hpp:1350
MueLu::RefMaxwell::getRangeMap
const Teuchos::RCP< const Map > getRangeMap() const
Returns the Xpetra::Map object associated with the range of this operator.
Definition: MueLu_RefMaxwell_def.hpp:91
MueLu_Zoltan2Interface_fwd.hpp
MueLu_TentativePFactory_kokkos_fwd.hpp
MueLu::RefMaxwell::RealValuedMultiVector
Xpetra::MultiVector< coordinateType, LO, GO, NO > RealValuedMultiVector
Definition: MueLu_RefMaxwell_decl.hpp:227
MueLu::RefMaxwell::Material_alpha_
Teuchos::RCP< MultiVector > Material_alpha_
Definition: MueLu_RefMaxwell_decl.hpp:664
MueLu_CoordinatesTransferFactory_fwd.hpp
Xpetra_MapFactory_fwd.hpp
MueLu::RefMaxwell::BCrows11_
Kokkos::View< bool *, typename Node::device_type::memory_space > BCrows11_
Vectors for BCs.
Definition: MueLu_RefMaxwell_decl.hpp:678
MueLu::RefMaxwell::CoarseNullspace22_
Teuchos::RCP< MultiVector > CoarseNullspace22_
Nullspace for coarse (2,2) problem.
Definition: MueLu_RefMaxwell_decl.hpp:687
Teuchos::ParameterList::get
T & get(const std::string &name, T def_value)
MueLu::RefMaxwell::setFineLevelSmoother11
void setFineLevelSmoother11()
Set the fine level smoother.
Definition: MueLu_RefMaxwell_def.hpp:1185
MueLu::RefMaxwell::SM_Matrix_
Teuchos::RCP< Matrix > SM_Matrix_
The system that is getting preconditioned.
Definition: MueLu_RefMaxwell_decl.hpp:650
MueLu_SaPFactory_fwd.hpp
Xpetra::Operator::MultiVector
Xpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > MultiVector
MueLu::RefMaxwell::Mk_one_
Teuchos::RCP< Matrix > Mk_one_
mass matrices with unit weight on k-th and (k-1)-th spaces
Definition: MueLu_RefMaxwell_decl.hpp:660
MueLu_TentativePFactory_fwd.hpp
MueLu_RebalanceAcFactory_fwd.hpp
MueLu_BaseClass.hpp
MueLu::RefMaxwell::solverName_
std::string solverName_
The name of the solver.
Definition: MueLu_RefMaxwell_decl.hpp:646
MueLu::RefMaxwell::buildProjection
Teuchos::RCP< Matrix > buildProjection(const int spaceNumber, const RCP< MultiVector > &EdgeNullspace) const
Builds a projection from a vector values space into a vector valued nodal space.
Definition: MueLu_RefMaxwell_def.hpp:1554
MueLu::RefMaxwell::DxIntermediate_
Teuchos::RCP< MultiVector > DxIntermediate_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::buildNodalProlongator
void buildNodalProlongator(const Teuchos::RCP< Matrix > &A_nodal, Teuchos::RCP< Matrix > &P_nodal, Teuchos::RCP< MultiVector > &Nullspace_nodal, Teuchos::RCP< RealValuedMultiVector > &Coords_nodal) const
Definition: MueLu_RefMaxwell_def.hpp:1679
MueLu::RefMaxwell::A22_
Teuchos::RCP< Matrix > A22_
Definition: MueLu_RefMaxwell_decl.hpp:670
Xpetra_VectorFactory_fwd.hpp
MueLu_Utilities_fwd.hpp
MueLu_UseShortNames.hpp
MueLu::RefMaxwell::determineSubHierarchyCommSizes
void determineSubHierarchyCommSizes(bool &doRebalancing, int &rebalanceStriding, int &numProcsCoarseA11, int &numProcsA22)
Determine how large the sub-communicators for the two hierarchies should be.
Definition: MueLu_RefMaxwell_def.hpp:606
MueLu::RefMaxwell::PostSmoother11_
Teuchos::RCP< SmootherBase > PostSmoother11_
Definition: MueLu_RefMaxwell_decl.hpp:640
Teuchos::EVerbosityLevel
EVerbosityLevel
MueLu::RefMaxwell::solveH
void solveH(const MultiVector &RHS, MultiVector &X) const
apply solve to 1-1 block only
Definition: MueLu_RefMaxwell_def.hpp:2354
MueLu::RefMaxwell::Addon11_
Teuchos::RCP< Matrix > Addon11_
the addon for the 11 block
Definition: MueLu_RefMaxwell_decl.hpp:672
MueLu::RefMaxwell::Material_beta_
Teuchos::RCP< MultiVector > Material_beta_
material for first space
Definition: MueLu_RefMaxwell_decl.hpp:664
MueLu::RefMaxwell::Mk_1_one_
Teuchos::RCP< Matrix > Mk_1_one_
Definition: MueLu_RefMaxwell_decl.hpp:660
MueLu::RefMaxwell::setParameters
void setParameters(Teuchos::ParameterList &list)
Set parameters.
Definition: MueLu_RefMaxwell_def.hpp:193
MueLu_ConfigDefs.hpp
Node
MueLu::DefaultNode Node
Definition: MueLu_UseDefaultTypes.hpp:15
MueLu::RefMaxwell::invMk_2_invAlpha_
Teuchos::RCP< Matrix > invMk_2_invAlpha_
Definition: MueLu_RefMaxwell_decl.hpp:658
MueLu::RefMaxwell::magnitudeType
Teuchos::ScalarTraits< Scalar >::magnitudeType magnitudeType
Definition: MueLu_RefMaxwell_decl.hpp:225
MueLu::RefMaxwell::D0_
Teuchos::RCP< Matrix > D0_
D_0 matrix.
Definition: MueLu_RefMaxwell_decl.hpp:656
MueLu::RefMaxwell::P11res_
Teuchos::RCP< MultiVector > P11res_
Temporary memory.
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::buildNullspace
RCP< MultiVector > buildNullspace(const int spaceNumber, const Kokkos::View< bool *, typename Node::device_type > &bcs, const bool applyBCs)
Builds a nullspace.
Definition: MueLu_RefMaxwell_def.hpp:1398
MueLu::VerboseObject
Verbose class for MueLu classes.
Definition: MueLu_VerboseObject.hpp:43
MueLu::MasterList
Static class that holds the complete list of valid MueLu parameters.
Definition: MueLu_MasterList.hpp:37
MueLu_UncoupledAggregationFactory_fwd.hpp
MueLu_MasterList.hpp
MueLu::RefMaxwell::P11x_colmap_
Teuchos::RCP< MultiVector > P11x_colmap_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell
Preconditioner (wrapped as a Xpetra::Operator) for Maxwell&#39;s equations in curl-curl form...
Definition: MueLu_RefMaxwell_decl.hpp:220
MueLu::RefMaxwell::ImporterCoarse11_
Teuchos::RCP< const Import > ImporterCoarse11_
Importer to coarse (1,1) hierarchy.
Definition: MueLu_RefMaxwell_decl.hpp:689
MueLu::RefMaxwell::RefMaxwell
RefMaxwell(const Teuchos::RCP< Matrix > &SM_Matrix, const Teuchos::RCP< Matrix > &D0_Matrix, const Teuchos::RCP< Matrix > &M1_Matrix, const Teuchos::RCP< MultiVector > &Nullspace11, const Teuchos::RCP< RealValuedMultiVector > &NodalCoords, Teuchos::ParameterList &List, bool ComputePrec)
Definition: MueLu_RefMaxwell_decl.hpp:392
MueLu::RefMaxwell::P11resSubComm_
Teuchos::RCP< MultiVector > P11resSubComm_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::RefMaxwell
RefMaxwell(const Teuchos::RCP< Matrix > &SM_Matrix, const Teuchos::RCP< Matrix > &D0_Matrix, const Teuchos::RCP< Matrix > &M0inv_Matrix, const Teuchos::RCP< Matrix > &M1_Matrix, const Teuchos::RCP< MultiVector > &Nullspace11, const Teuchos::RCP< RealValuedMultiVector > &NodalCoords, Teuchos::ParameterList &List, bool ComputePrec=true)
Definition: MueLu_RefMaxwell_decl.hpp:351
MueLu::RefMaxwell::allocateMemory
void allocateMemory(int numVectors) const
allocate multivectors for solve
Definition: MueLu_RefMaxwell_def.hpp:1258
MueLu::RefMaxwell::coarseA22_
Teuchos::RCP< Matrix > coarseA22_
Definition: MueLu_RefMaxwell_decl.hpp:670
MueLu::RefMaxwell::residual_
Teuchos::RCP< MultiVector > residual_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::DxSubComm_
Teuchos::RCP< MultiVector > DxSubComm_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::initialize
void initialize(const Teuchos::RCP< Matrix > &D0_Matrix, const Teuchos::RCP< Matrix > &Ms_Matrix, const Teuchos::RCP< Matrix > &M0inv_Matrix, const Teuchos::RCP< Matrix > &M1_Matrix, const Teuchos::RCP< MultiVector > &Nullspace11, const Teuchos::RCP< RealValuedMultiVector > &NodalCoords, const Teuchos::RCP< MultiVector > &Material, Teuchos::ParameterList &List)
Definition: MueLu_RefMaxwell_def.hpp:2572
MueLu::RefMaxwell::applyInverseAdditive
void applyInverseAdditive(const MultiVector &RHS, MultiVector &X) const
apply additive algorithm for 2x2 solve
Definition: MueLu_RefMaxwell_def.hpp:2081
MueLu::RefMaxwell::P11resTmp_
Teuchos::RCP< MultiVector > P11resTmp_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::Dk_2_
Teuchos::RCP< Matrix > Dk_2_
D_{k-2} matrix.
Definition: MueLu_RefMaxwell_decl.hpp:654
MueLu::RefMaxwell::invMk_1_invBeta_
Teuchos::RCP< Matrix > invMk_1_invBeta_
inverse of mass matrices on (k-1)-th and (k-2)-th space with weights 1/beta and 1/alpha respectively ...
Definition: MueLu_RefMaxwell_decl.hpp:658
Xpetra_CrsMatrixWrap_fwd.hpp
Teuchos::VERB_HIGH
MueLu::RefMaxwell::BCcols22_
Kokkos::View< bool *, typename Node::device_type::memory_space > BCcols22_
Definition: MueLu_RefMaxwell_decl.hpp:678
Scalar
MueLu::DefaultScalar Scalar
Definition: MueLu_UseDefaultTypes.hpp:12
MueLu::RefMaxwell::resetMatrix
void resetMatrix(Teuchos::RCP< Matrix > SM_Matrix_new, bool ComputePrec=true)
Reset system matrix.
Definition: MueLu_RefMaxwell_def.hpp:2072
MueLu_AggregationExportFactory_fwd.hpp
GlobalOrdinal
MueLu::DefaultGlobalOrdinal GlobalOrdinal
Definition: MueLu_UseDefaultTypes.hpp:14
MueLu::RefMaxwell::rebalanceCoarse11Matrix
void rebalanceCoarse11Matrix(const int rebalanceStriding, const int numProcsCoarseA11)
rebalance the coarse A11 matrix, as well as P11, CoordsCoarse11 and Addon11
Definition: MueLu_RefMaxwell_def.hpp:874
MueLu::RefMaxwell::precList11_
Teuchos::ParameterList precList11_
Definition: MueLu_RefMaxwell_decl.hpp:693
MueLu_CoalesceDropFactory_kokkos_fwd.hpp
Teuchos::basic_FancyOStream
MueLu::RefMaxwell::DorigImporter_
Teuchos::RCP< const Import > DorigImporter_
Definition: MueLu_RefMaxwell_decl.hpp:676
MueLu::RefMaxwell::P11_
Teuchos::RCP< Matrix > P11_
special prolongator for 11 block and its transpose
Definition: MueLu_RefMaxwell_decl.hpp:666
MueLu_RepartitionFactory_fwd.hpp
Xpetra_MatrixFactory_fwd.hpp
MueLu::RefMaxwell::Hierarchy22_
Teuchos::RCP< Hierarchy > Hierarchy22_
Definition: MueLu_RefMaxwell_decl.hpp:639
MueLu::RefMaxwell::buildVectorNodalProlongator
RCP< Matrix > buildVectorNodalProlongator(const Teuchos::RCP< Matrix > &P_nodal) const
Definition: MueLu_RefMaxwell_def.hpp:1795
MueLu::RefMaxwell::Dres_
Teuchos::RCP< MultiVector > Dres_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::dump
void dump(const RCP< Matrix > &A, std::string name) const
dump out matrix
Definition: MueLu_RefMaxwell_def.hpp:1334
MueLu::RefMaxwell::buildCoarse11Matrix
void buildCoarse11Matrix()
Compute coarseA11 = P11^{T}*SM*P11 + addon efficiently.
Definition: MueLu_RefMaxwell_def.hpp:787
MueLu::RefMaxwell::spaceNumber_
int spaceNumber_
The number of the space in the deRham complex.
Definition: MueLu_RefMaxwell_decl.hpp:644
MueLu_PerfUtils_fwd.hpp
MueLu::RefMaxwell::coarseA11_RAP_reuse_data_
Teuchos::RCP< Teuchos::ParameterList > coarseA11_RAP_reuse_data_
Definition: MueLu_RefMaxwell_decl.hpp:694
MueLu_CoalesceDropFactory_fwd.hpp
MueLu::RefMaxwell::HierarchyCoarse11_
Teuchos::RCP< Hierarchy > HierarchyCoarse11_
Two hierarchies: one for the coarse (1,1)-block, another for the (2,2)-block.
Definition: MueLu_RefMaxwell_decl.hpp:639
MueLu::RefMaxwell::DTR11Tmp_
Teuchos::RCP< MultiVector > DTR11Tmp_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::RefMaxwell
RefMaxwell(const Teuchos::RCP< Matrix > &SM_Matrix, const Teuchos::RCP< Matrix > &Dk_1, const Teuchos::RCP< Matrix > &Dk_2, const Teuchos::RCP< Matrix > &D0, const Teuchos::RCP< Matrix > &M1_beta, const Teuchos::RCP< Matrix > &M1_alpha, const Teuchos::RCP< Matrix > &Mk_one, const Teuchos::RCP< Matrix > &Mk_1_one, const Teuchos::RCP< Matrix > &invMk_1_invBeta, const Teuchos::RCP< Matrix > &invMk_2_invAlpha, const Teuchos::RCP< MultiVector > &Nullspace11, const Teuchos::RCP< MultiVector > &Nullspace22, const Teuchos::RCP< RealValuedMultiVector > &NodalCoords, Teuchos::ParameterList &List, bool ComputePrec=true)
Definition: MueLu_RefMaxwell_decl.hpp:245
MueLu::RefMaxwell::Dk_1_T_
Teuchos::RCP< Matrix > Dk_1_T_
Definition: MueLu_RefMaxwell_decl.hpp:652
MueLu::RefMaxwell::thyraPrecOp22_
RCP< Operator > thyraPrecOp22_
Definition: MueLu_RefMaxwell_decl.hpp:642
MueLu::RefMaxwell::build22Matrix
void build22Matrix(const bool reuse, const bool doRebalancing, const int rebalanceStriding, const int numProcsA22)
Setup A22 = D0^T SM D0 and rebalance it, as well as D0 and Coords_.
Definition: MueLu_RefMaxwell_def.hpp:993
MueLu::RefMaxwell::CoordsCoarse11_
Teuchos::RCP< RealValuedMultiVector > CoordsCoarse11_
Definition: MueLu_RefMaxwell_decl.hpp:683
MueLu::RefMaxwell::NullspaceCoarse11_
Teuchos::RCP< MultiVector > NullspaceCoarse11_
Nullspace for coarse (1,1) problem.
Definition: MueLu_RefMaxwell_decl.hpp:685
MueLu_RepartitionHeuristicFactory_fwd.hpp
MueLu::RefMaxwell::M1_alpha_
Teuchos::RCP< Matrix > M1_alpha_
Definition: MueLu_RefMaxwell_decl.hpp:662
MueLu::RefMaxwell::getTimer
Teuchos::RCP< Teuchos::TimeMonitor > getTimer(std::string name, RCP< const Teuchos::Comm< int > > comm=Teuchos::null) const
get a (synced) timer
Definition: MueLu_RefMaxwell_def.hpp:1382
MueLu::RefMaxwell::compute
void compute(bool reuse=false)
Setup the preconditioner.
Definition: MueLu_RefMaxwell_def.hpp:287
MueLu_RAPFactory_fwd.hpp
MueLu::RefMaxwell::P11xSubComm_
Teuchos::RCP< MultiVector > P11xSubComm_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu_CoarseMapFactory_fwd.hpp
MueLu::RefMaxwell::DorigDomainMap_
Teuchos::RCP< const Map > DorigDomainMap_
Definition: MueLu_RefMaxwell_decl.hpp:675
MueLu::RefMaxwell::getValidParamterList
Teuchos::RCP< Teuchos::ParameterList > getValidParamterList()
Definition: MueLu_RefMaxwell_def.hpp:98
MueLu::RefMaxwell::describe
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::VERB_HIGH) const
Definition: MueLu_RefMaxwell_def.hpp:2803
MueLu::RefMaxwell::disable_addon_
bool disable_addon_
Some options.
Definition: MueLu_RefMaxwell_decl.hpp:697
MueLu::RefMaxwell::RefMaxwell
RefMaxwell(const Teuchos::RCP< Matrix > &D0_Matrix, const Teuchos::RCP< Matrix > &M1_Matrix, const Teuchos::RCP< MultiVector > &Nullspace11, const Teuchos::RCP< RealValuedMultiVector > &NodalCoords, Teuchos::ParameterList &List)
Definition: MueLu_RefMaxwell_decl.hpp:411
MueLu::RefMaxwell::RefMaxwell
RefMaxwell(const Teuchos::RCP< Matrix > &SM_Matrix, const Teuchos::RCP< Matrix > &D0_Matrix, const Teuchos::RCP< Matrix > &Ms_Matrix, const Teuchos::RCP< Matrix > &M0inv_Matrix, const Teuchos::RCP< Matrix > &M1_Matrix, const Teuchos::RCP< MultiVector > &Nullspace11, const Teuchos::RCP< RealValuedMultiVector > &NodalCoords, Teuchos::ParameterList &List, bool ComputePrec=true)
Definition: MueLu_RefMaxwell_decl.hpp:313
MueLu::RefMaxwell::~RefMaxwell
virtual ~RefMaxwell()
Destructor.
Definition: MueLu_RefMaxwell_decl.hpp:431
MueLu::RefMaxwell::M1_beta_
Teuchos::RCP< Matrix > M1_beta_
mass matrices on first space with weights beta and alpha respectively
Definition: MueLu_RefMaxwell_decl.hpp:662
MueLu::RefMaxwell::coarseA11_
Teuchos::RCP< Matrix > coarseA11_
coarse 11, 22 and coarse 22 blocks
Definition: MueLu_RefMaxwell_decl.hpp:670
MueLu::RefMaxwell::parameterList_
Teuchos::ParameterList parameterList_
Parameter lists.
Definition: MueLu_RefMaxwell_decl.hpp:693
MueLu::RefMaxwell::P22_
Teuchos::RCP< Matrix > P22_
special prolongator for 22 block and its transpose
Definition: MueLu_RefMaxwell_decl.hpp:668
MueLu::RefMaxwell::Coords22_
Teuchos::RCP< RealValuedMultiVector > Coords22_
Definition: MueLu_RefMaxwell_decl.hpp:683
MueLu_Hierarchy_fwd.hpp
MueLu::RefMaxwell::Nullspace11_
Teuchos::RCP< MultiVector > Nullspace11_
Nullspace for (1.1) block.
Definition: MueLu_RefMaxwell_decl.hpp:681
MueLu::RefMaxwell::Dx_colmap_
Teuchos::RCP< MultiVector > Dx_colmap_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::Nullspace22_
Teuchos::RCP< MultiVector > Nullspace22_
Definition: MueLu_RefMaxwell_decl.hpp:681
MueLu::RefMaxwell::DresIntermediate_
Teuchos::RCP< MultiVector > DresIntermediate_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::P11x_
Teuchos::RCP< MultiVector > P11x_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::R11_
Teuchos::RCP< Matrix > R11_
Definition: MueLu_RefMaxwell_decl.hpp:666
MueLu::RefMaxwell::apply
void apply(const MultiVector &X, MultiVector &Y, Teuchos::ETransp mode=Teuchos::NO_TRANS, Scalar alpha=Teuchos::ScalarTraits< Scalar >::one(), Scalar beta=Teuchos::ScalarTraits< Scalar >::zero()) const
Definition: MueLu_RefMaxwell_def.hpp:2419
MueLu::RefMaxwell::RefMaxwell
RefMaxwell(const Teuchos::RCP< Matrix > &D0_Matrix, const Teuchos::RCP< Matrix > &M0inv_Matrix, const Teuchos::RCP< Matrix > &M1_Matrix, const Teuchos::RCP< MultiVector > &Nullspace11, const Teuchos::RCP< RealValuedMultiVector > &NodalCoords, Teuchos::ParameterList &List)
Definition: MueLu_RefMaxwell_decl.hpp:372
MueLu::RefMaxwell::setupSubSolve
void setupSubSolve(Teuchos::RCP< Hierarchy > &hierarchy, Teuchos::RCP< Operator > &thyraPrecOp, const Teuchos::RCP< Matrix > &A, const Teuchos::RCP< MultiVector > &Nullspace, const Teuchos::RCP< RealValuedMultiVector > &Coords, const Teuchos::RCP< MultiVector > &Material, Teuchos::ParameterList &params, std::string &label, const bool reuse, const bool isSingular=false)
Setup a subsolve.
Definition: MueLu_RefMaxwell_def.hpp:1995
MueLu::RefMaxwell::getDomainMap
const Teuchos::RCP< const Map > getDomainMap() const
Returns the Xpetra::Map object associated with the domain of this operator.
Definition: MueLu_RefMaxwell_def.hpp:86
MueLu::RefMaxwell::residual
void residual(const MultiVector &X, const MultiVector &B, MultiVector &R) const
Compute a residual R = B - (*this) * X.
Definition: MueLu_RefMaxwell_decl.hpp:467
MueLu_SmootherBase_fwd.hpp
MueLu::RefMaxwell::DresSubComm_
Teuchos::RCP< MultiVector > DresSubComm_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::RefMaxwell
RefMaxwell(const Teuchos::RCP< Matrix > &SM_Matrix, const Teuchos::RCP< Matrix > &D0_Matrix, const Teuchos::RCP< Matrix > &Ms_Matrix, const Teuchos::RCP< Matrix > &M0inv_Matrix, const Teuchos::RCP< Matrix > &M1_Matrix, const Teuchos::RCP< MultiVector > &Nullspace11, const Teuchos::RCP< RealValuedMultiVector > &NodalCoords, const Teuchos::RCP< MultiVector > &Material, Teuchos::ParameterList &List, bool ComputePrec=true)
Definition: MueLu_RefMaxwell_decl.hpp:326
MueLu::RefMaxwell::coarseA11_AP_reuse_data_
Teuchos::RCP< Teuchos::ParameterList > coarseA11_AP_reuse_data_
Definition: MueLu_RefMaxwell_decl.hpp:694
MueLu::RefMaxwell::DresTmp_
Teuchos::RCP< MultiVector > DresTmp_
Definition: MueLu_RefMaxwell_decl.hpp:703
MueLu::RefMaxwell::Addon22_
Teuchos::RCP< Matrix > Addon22_
the addon for the 22 block
Definition: MueLu_RefMaxwell_decl.hpp:674
MueLu_SmootherPrototype_fwd.hpp
MueLu::RefMaxwell::RefMaxwell
RefMaxwell(const Teuchos::RCP< Matrix > &SM_Matrix, const Teuchos::RCP< Matrix > &Dk_1, const Teuchos::RCP< Matrix > &Dk_2, const Teuchos::RCP< Matrix > &D0, const Teuchos::RCP< Matrix > &M1_beta, const Teuchos::RCP< Matrix > &M1_alpha, const Teuchos::RCP< Matrix > &Mk_one, const Teuchos::RCP< Matrix > &Mk_1_one, const Teuchos::RCP< Matrix > &invMk_1_invBeta, const Teuchos::RCP< Matrix > &invMk_2_invAlpha, const Teuchos::RCP< MultiVector > &Nullspace11, const Teuchos::RCP< MultiVector > &Nullspace22, const Teuchos::RCP< RealValuedMultiVector > &NodalCoords, const Teuchos::RCP< MultiVector > &Material_beta, const Teuchos::RCP< MultiVector > &Material_alpha, Teuchos::ParameterList &List, bool ComputePrec=true)
Definition: MueLu_RefMaxwell_decl.hpp:272
MueLu::RefMaxwell::thyraPrecOpH_
RCP< Operator > thyraPrecOpH_
Definition: MueLu_RefMaxwell_decl.hpp:642
MueLu_ZoltanInterface_fwd.hpp
MueLu::RefMaxwell::coordinateType
Teuchos::ScalarTraits< Scalar >::coordinateType coordinateType
Definition: MueLu_RefMaxwell_decl.hpp:226
MueLu::RefMaxwell::RefMaxwell
RefMaxwell(Teuchos::RCP< Hierarchy > HH, Teuchos::RCP< Hierarchy > H22)
Constructor with Hierarchies.
Definition: MueLu_RefMaxwell_decl.hpp:238
MueLu::RefMaxwell::buildProlongator
void buildProlongator(const int spaceNumber, const Teuchos::RCP< Matrix > &A_nodal_Matrix, const RCP< MultiVector > &EdgeNullspace, Teuchos::RCP< Matrix > &edgeProlongator, Teuchos::RCP< MultiVector > &coarseEdgeNullspace, Teuchos::RCP< RealValuedMultiVector > &coarseNodalCoords) const
Definition: MueLu_RefMaxwell_def.hpp:1861
Xpetra_Map_fwd.hpp
MueLu::RefMaxwell::A22_AP_reuse_data_
Teuchos::RCP< Teuchos::ParameterList > A22_AP_reuse_data_
Definition: MueLu_RefMaxwell_decl.hpp:695
MueLu::RefMaxwell::PostSmootherData11_
Teuchos::RCP< SmootherPrototype > PostSmootherData11_
Definition: MueLu_RefMaxwell_decl.hpp:641
MueLu::RefMaxwell::Dx_
Teuchos::RCP< MultiVector > Dx_
Definition: MueLu_RefMaxwell_decl.hpp:703
Xpetra_Operator.hpp
Xpetra_MultiVectorFactory_fwd.hpp
MueLu::RefMaxwell::NodalCoords_
Teuchos::RCP< RealValuedMultiVector > NodalCoords_
Coordinates.
Definition: MueLu_RefMaxwell_decl.hpp:683
Xpetra_Matrix_fwd.hpp
MueLu::RefMaxwell::precList22_
Teuchos::ParameterList precList22_
Definition: MueLu_RefMaxwell_decl.hpp:693
MueLu::RefMaxwell::Dk_1_
Teuchos::RCP< Matrix > Dk_1_
D_{k-1} matrix and its transpose.
Definition: MueLu_RefMaxwell_decl.hpp:652
MueLu::RefMaxwell::hasTransposeApply
bool hasTransposeApply() const
Indicates whether this operator supports applying the adjoint operator.
Definition: MueLu_RefMaxwell_def.hpp:2481