deal.II version GIT relicensing-2206-gaa53ff9447 2024-12-02 09:10:00+00:00
\(\newcommand{\dealvcentcolon}{\mathrel{\mathop{:}}}\) \(\newcommand{\dealcoloneq}{\dealvcentcolon\mathrel{\mkern-1.2mu}=}\) \(\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}\) \(\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}\)
Loading...
Searching...
No Matches
trilinos_precondition.cc
Go to the documentation of this file.
1// ------------------------------------------------------------------------
2//
3// SPDX-License-Identifier: LGPL-2.1-or-later
4// Copyright (C) 2008 - 2023 by the deal.II authors
5//
6// This file is part of the deal.II library.
7//
8// Part of the source code is dual licensed under Apache-2.0 WITH
9// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
10// governing the source code and code contributions can be found in
11// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
12//
13// ------------------------------------------------------------------------
14
15#include <deal.II/lac/trilinos_precondition.h>
16
17#ifdef DEAL_II_WITH_TRILINOS
18
19# include <deal.II/lac/sparse_matrix.h>
20# include <deal.II/lac/trilinos_sparse_matrix.h>
21# include <deal.II/lac/vector.h>
22
23# include <Epetra_MultiVector.h>
24# include <Ifpack.h>
25# include <Ifpack_Chebyshev.h>
26# include <Teuchos_ParameterList.hpp>
27# include <Teuchos_RCP.hpp>
28
29DEAL_II_NAMESPACE_OPEN
30
31namespace TrilinosWrappers
32{
33 PreconditionBase::PreconditionBase()
34 : communicator(MPI_COMM_SELF)
35 {}
36
37
38 void
39 PreconditionBase::clear()
40 {
41 preconditioner.reset();
42 communicator = MPI_COMM_SELF;
43 }
44
45
46 MPI_Comm
47 PreconditionBase::get_mpi_communicator() const
48 {
49 return communicator.Comm();
50 }
51
52
53 Epetra_Operator &
54 PreconditionBase::trilinos_operator() const
55 {
56 AssertThrow(!preconditioner.is_null(),
57 ExcMessage("Trying to dereference a null pointer."));
58 return (*preconditioner);
59 }
60
61
62 IndexSet
63 PreconditionBase::locally_owned_domain_indices() const
64 {
65 return IndexSet(preconditioner->OperatorDomainMap());
66 }
67
68
69 IndexSet
70 PreconditionBase::locally_owned_range_indices() const
71 {
72 return IndexSet(preconditioner->OperatorRangeMap());
73 }
74
75 /* -------------------------- PreconditionJacobi -------------------------- */
76
77 PreconditionJacobi::AdditionalData::AdditionalData(
78 const double omega,
79 const double min_diagonal,
80 const unsigned int n_sweeps)
81 : omega(omega)
82 , min_diagonal(min_diagonal)
83 , n_sweeps(n_sweeps)
84 {}
85
86
87
88 void
89 PreconditionJacobi::initialize(const SparseMatrix &matrix,
90 const AdditionalData &additional_data)
91 {
92 // release memory before reallocation
93 preconditioner.reset();
94 preconditioner.reset(
95 Ifpack().Create("point relaxation",
96 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
97 0));
98
99 Ifpack_Preconditioner *ifpack =
100 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
101 Assert(ifpack != nullptr,
102 ExcMessage("Trilinos could not create this "
103 "preconditioner"));
104
105 int ierr;
106
107 Teuchos::ParameterList parameter_list;
108 parameter_list.set("relaxation: sweeps",
109 static_cast<int>(additional_data.n_sweeps));
110 parameter_list.set("relaxation: type", "Jacobi");
111 parameter_list.set("relaxation: damping factor", additional_data.omega);
112 parameter_list.set("relaxation: min diagonal value",
113 additional_data.min_diagonal);
114
115 ierr = ifpack->SetParameters(parameter_list);
116 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
117
118 ierr = ifpack->Initialize();
119 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
120
121 ierr = ifpack->Compute();
122 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
123 }
124
125
126
127 /* -------------------------- PreconditionSSOR -------------------------- */
128
129 PreconditionSSOR::AdditionalData::AdditionalData(const double omega,
130 const double min_diagonal,
131 const unsigned int overlap,
132 const unsigned int n_sweeps)
133 : omega(omega)
134 , min_diagonal(min_diagonal)
135 , overlap(overlap)
136 , n_sweeps(n_sweeps)
137 {}
138
139
140
141 void
142 PreconditionSSOR::initialize(const SparseMatrix &matrix,
143 const AdditionalData &additional_data)
144 {
145 preconditioner.reset();
146 preconditioner.reset(
147 Ifpack().Create("point relaxation",
148 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
149 additional_data.overlap));
150
151 Ifpack_Preconditioner *ifpack =
152 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
153 Assert(ifpack != nullptr,
154 ExcMessage("Trilinos could not create this "
155 "preconditioner"));
156
157 int ierr;
158
159 Teuchos::ParameterList parameter_list;
160 parameter_list.set("relaxation: sweeps",
161 static_cast<int>(additional_data.n_sweeps));
162 parameter_list.set("relaxation: type", "symmetric Gauss-Seidel");
163 parameter_list.set("relaxation: damping factor", additional_data.omega);
164 parameter_list.set("relaxation: min diagonal value",
165 additional_data.min_diagonal);
166 parameter_list.set("schwarz: combine mode", "Add");
167
168 ierr = ifpack->SetParameters(parameter_list);
169 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
170
171 ierr = ifpack->Initialize();
172 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
173
174 ierr = ifpack->Compute();
175 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
176 }
177
178
179
180 /* -------------------------- PreconditionSOR -------------------------- */
181
182 PreconditionSOR::AdditionalData::AdditionalData(const double omega,
183 const double min_diagonal,
184 const unsigned int overlap,
185 const unsigned int n_sweeps)
186 : omega(omega)
187 , min_diagonal(min_diagonal)
188 , overlap(overlap)
189 , n_sweeps(n_sweeps)
190 {}
191
192
193
194 void
195 PreconditionSOR::initialize(const SparseMatrix &matrix,
196 const AdditionalData &additional_data)
197 {
198 preconditioner.reset();
199 preconditioner.reset(
200 Ifpack().Create("point relaxation",
201 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
202 additional_data.overlap));
203
204 Ifpack_Preconditioner *ifpack =
205 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
206 Assert(ifpack != nullptr,
207 ExcMessage("Trilinos could not create this "
208 "preconditioner"));
209
210 int ierr;
211
212 Teuchos::ParameterList parameter_list;
213 parameter_list.set("relaxation: sweeps",
214 static_cast<int>(additional_data.n_sweeps));
215 parameter_list.set("relaxation: type", "Gauss-Seidel");
216 parameter_list.set("relaxation: damping factor", additional_data.omega);
217 parameter_list.set("relaxation: min diagonal value",
218 additional_data.min_diagonal);
219 parameter_list.set("schwarz: combine mode", "Add");
220
221 ierr = ifpack->SetParameters(parameter_list);
222 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
223
224 ierr = ifpack->Initialize();
225 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
226
227 ierr = ifpack->Compute();
228 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
229 }
230
231
232
233 /* ----------------------- PreconditionBlockJacobi ---------------------- */
234
235 PreconditionBlockJacobi::AdditionalData::AdditionalData(
236 const unsigned int block_size,
237 const std::string &block_creation_type,
238 const double omega,
239 const double min_diagonal,
240 const unsigned int n_sweeps)
241 : block_size(block_size)
242 , block_creation_type(block_creation_type)
243 , omega(omega)
244 , min_diagonal(min_diagonal)
245 , n_sweeps(n_sweeps)
246 {}
247
248
249
250 void
251 PreconditionBlockJacobi::initialize(const SparseMatrix &matrix,
252 const AdditionalData &additional_data)
253 {
254 // release memory before reallocation
255 preconditioner.reset();
256
257 // Block relaxation setup fails if we have no locally owned rows. As a
258 // work-around we just pretend to use point relaxation on those processors:
259 preconditioner.reset(
260 Ifpack().Create((matrix.trilinos_matrix().NumMyRows() == 0) ?
261 "point relaxation" :
262 "block relaxation",
263 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
264 0));
265
266 Ifpack_Preconditioner *ifpack =
267 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
268 Assert(ifpack != nullptr,
269 ExcMessage("Trilinos could not create this "
270 "preconditioner"));
271
272 int ierr;
273
274 Teuchos::ParameterList parameter_list;
275 parameter_list.set("relaxation: sweeps",
276 static_cast<int>(additional_data.n_sweeps));
277 parameter_list.set("relaxation: type", "Jacobi");
278 parameter_list.set("relaxation: damping factor", additional_data.omega);
279 parameter_list.set("relaxation: min diagonal value",
280 additional_data.min_diagonal);
281 parameter_list.set("partitioner: type",
282 additional_data.block_creation_type);
283 int n_local_parts =
284 (matrix.trilinos_matrix().NumMyRows() + additional_data.block_size - 1) /
285 additional_data.block_size;
286 parameter_list.set("partitioner: local parts", n_local_parts);
287
288 ierr = ifpack->SetParameters(parameter_list);
289 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
290
291 ierr = ifpack->Initialize();
292 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
293
294 ierr = ifpack->Compute();
295 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
296 }
297
298
299
300 /* ----------------------- PreconditionBlockSSOR ------------------------ */
301
302 PreconditionBlockSSOR::AdditionalData::AdditionalData(
303 const unsigned int block_size,
304 const std::string &block_creation_type,
305 const double omega,
306 const double min_diagonal,
307 const unsigned int overlap,
308 const unsigned int n_sweeps)
309 : block_size(block_size)
310 , block_creation_type(block_creation_type)
311 , omega(omega)
312 , min_diagonal(min_diagonal)
313 , overlap(overlap)
314 , n_sweeps(n_sweeps)
315 {}
316
317
318
319 void
320 PreconditionBlockSSOR::initialize(const SparseMatrix &matrix,
321 const AdditionalData &additional_data)
322 {
323 preconditioner.reset();
324
325 // Block relaxation setup fails if we have no locally owned rows. As a
326 // work-around we just pretend to use point relaxation on those processors:
327 preconditioner.reset(
328 Ifpack().Create((matrix.trilinos_matrix().NumMyRows() == 0) ?
329 "point relaxation" :
330 "block relaxation",
331 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
332 additional_data.overlap));
333
334 Ifpack_Preconditioner *ifpack =
335 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
336 Assert(ifpack != nullptr,
337 ExcMessage("Trilinos could not create this "
338 "preconditioner"));
339
340 int ierr;
341
342 Teuchos::ParameterList parameter_list;
343 parameter_list.set("relaxation: sweeps",
344 static_cast<int>(additional_data.n_sweeps));
345 parameter_list.set("relaxation: type", "symmetric Gauss-Seidel");
346 parameter_list.set("relaxation: damping factor", additional_data.omega);
347 parameter_list.set("relaxation: min diagonal value",
348 additional_data.min_diagonal);
349 parameter_list.set("schwarz: combine mode", "Add");
350 parameter_list.set("partitioner: type",
351 additional_data.block_creation_type);
352 int n_local_parts =
353 (matrix.trilinos_matrix().NumMyRows() + additional_data.block_size - 1) /
354 additional_data.block_size;
355 parameter_list.set("partitioner: local parts", n_local_parts);
356
357 ierr = ifpack->SetParameters(parameter_list);
358 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
359
360 ierr = ifpack->Initialize();
361 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
362
363 ierr = ifpack->Compute();
364 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
365 }
366
367
368
369 /* ------------------------ PreconditionBlockSOR ------------------------ */
370
371 PreconditionBlockSOR::AdditionalData::AdditionalData(
372 const unsigned int block_size,
373 const std::string &block_creation_type,
374 const double omega,
375 const double min_diagonal,
376 const unsigned int overlap,
377 const unsigned int n_sweeps)
378 : block_size(block_size)
379 , block_creation_type(block_creation_type)
380 , omega(omega)
381 , min_diagonal(min_diagonal)
382 , overlap(overlap)
383 , n_sweeps(n_sweeps)
384 {}
385
386
387
388 void
389 PreconditionBlockSOR::initialize(const SparseMatrix &matrix,
390 const AdditionalData &additional_data)
391 {
392 preconditioner.reset();
393
394 // Block relaxation setup fails if we have no locally owned rows. As a
395 // work-around we just pretend to use point relaxation on those processors:
396 preconditioner.reset(
397 Ifpack().Create((matrix.trilinos_matrix().NumMyRows() == 0) ?
398 "point relaxation" :
399 "block relaxation",
400 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
401 additional_data.overlap));
402
403 Ifpack_Preconditioner *ifpack =
404 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
405 Assert(ifpack != nullptr,
406 ExcMessage("Trilinos could not create this "
407 "preconditioner"));
408
409 int ierr;
410
411 Teuchos::ParameterList parameter_list;
412 parameter_list.set("relaxation: sweeps",
413 static_cast<int>(additional_data.n_sweeps));
414 parameter_list.set("relaxation: type", "Gauss-Seidel");
415 parameter_list.set("relaxation: damping factor", additional_data.omega);
416 parameter_list.set("relaxation: min diagonal value",
417 additional_data.min_diagonal);
418 parameter_list.set("schwarz: combine mode", "Add");
419 parameter_list.set("partitioner: type",
420 additional_data.block_creation_type);
421 int n_local_parts =
422 (matrix.trilinos_matrix().NumMyRows() + additional_data.block_size - 1) /
423 additional_data.block_size;
424 parameter_list.set("partitioner: local parts", n_local_parts);
425
426 ierr = ifpack->SetParameters(parameter_list);
427 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
428
429 ierr = ifpack->Initialize();
430 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
431
432 ierr = ifpack->Compute();
433 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
434 }
435
436
437
438 /* -------------------------- PreconditionIC -------------------------- */
439
440 PreconditionIC::AdditionalData::AdditionalData(const unsigned int ic_fill,
441 const double ic_atol,
442 const double ic_rtol,
443 const unsigned int overlap)
444 : ic_fill(ic_fill)
445 , ic_atol(ic_atol)
446 , ic_rtol(ic_rtol)
447 , overlap(overlap)
448 {}
449
450
451
452 void
453 PreconditionIC::initialize(const SparseMatrix &matrix,
454 const AdditionalData &additional_data)
455 {
456 preconditioner.reset();
457 preconditioner.reset(
458 Ifpack().Create("IC",
459 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
460 additional_data.overlap));
461
462 Ifpack_Preconditioner *ifpack =
463 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
464 Assert(ifpack != nullptr,
465 ExcMessage("Trilinos could not create this "
466 "preconditioner"));
467
468 int ierr;
469
470 Teuchos::ParameterList parameter_list;
471 parameter_list.set("fact: level-of-fill", additional_data.ic_fill);
472 parameter_list.set("fact: absolute threshold", additional_data.ic_atol);
473 parameter_list.set("fact: relative threshold", additional_data.ic_rtol);
474 parameter_list.set("schwarz: combine mode", "Add");
475
476 ierr = ifpack->SetParameters(parameter_list);
477 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
478
479 ierr = ifpack->Initialize();
480 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
481
482 ierr = ifpack->Compute();
483 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
484 }
485
486
487
488 /* -------------------------- PreconditionILU -------------------------- */
489
490 PreconditionILU::AdditionalData::AdditionalData(const unsigned int ilu_fill,
491 const double ilu_atol,
492 const double ilu_rtol,
493 const unsigned int overlap)
494 : ilu_fill(ilu_fill)
495 , ilu_atol(ilu_atol)
496 , ilu_rtol(ilu_rtol)
497 , overlap(overlap)
498 {}
499
500
501
502 void
503 PreconditionILU::initialize(const SparseMatrix &matrix,
504 const AdditionalData &additional_data)
505 {
506 preconditioner.reset();
507 preconditioner.reset(
508 Ifpack().Create("ILU",
509 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
510 additional_data.overlap));
511
512 Ifpack_Preconditioner *ifpack =
513 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
514 Assert(ifpack != nullptr,
515 ExcMessage("Trilinos could not create this "
516 "preconditioner"));
517
518 int ierr;
519
520 Teuchos::ParameterList parameter_list;
521 parameter_list.set("fact: level-of-fill",
522 static_cast<int>(additional_data.ilu_fill));
523 parameter_list.set("fact: absolute threshold", additional_data.ilu_atol);
524 parameter_list.set("fact: relative threshold", additional_data.ilu_rtol);
525 parameter_list.set("schwarz: combine mode", "Add");
526
527 ierr = ifpack->SetParameters(parameter_list);
528 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
529
530 ierr = ifpack->Initialize();
531 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
532
533 ierr = ifpack->Compute();
534 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
535 }
536
537
538
539 /* -------------------------- PreconditionILUT -------------------------- */
540
541 PreconditionILUT::AdditionalData::AdditionalData(const double ilut_drop,
542 const unsigned int ilut_fill,
543 const double ilut_atol,
544 const double ilut_rtol,
545 const unsigned int overlap)
546 : ilut_drop(ilut_drop)
547 , ilut_fill(ilut_fill)
548 , ilut_atol(ilut_atol)
549 , ilut_rtol(ilut_rtol)
550 , overlap(overlap)
551 {}
552
553
554
555 void
556 PreconditionILUT::initialize(const SparseMatrix &matrix,
557 const AdditionalData &additional_data)
558 {
559 preconditioner.reset();
560 preconditioner.reset(
561 Ifpack().Create("ILUT",
562 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
563 additional_data.overlap));
564
565 Ifpack_Preconditioner *ifpack =
566 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
567 Assert(ifpack != nullptr,
568 ExcMessage("Trilinos could not create this "
569 "preconditioner"));
570
571 int ierr;
572
573 Teuchos::ParameterList parameter_list;
574 parameter_list.set("fact: drop value", additional_data.ilut_drop);
575 parameter_list.set("fact: level-of-fill",
576 static_cast<int>(additional_data.ilut_fill));
577 parameter_list.set("fact: absolute threshold", additional_data.ilut_atol);
578 parameter_list.set("fact: relative threshold", additional_data.ilut_rtol);
579 parameter_list.set("schwarz: combine mode", "Add");
580
581 ierr = ifpack->SetParameters(parameter_list);
582 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
583
584 ierr = ifpack->Initialize();
585 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
586
587 ierr = ifpack->Compute();
588 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
589 }
590
591
592
593 /* ---------------------- PreconditionBlockDirect --------------------- */
594
595 PreconditionBlockwiseDirect::AdditionalData::AdditionalData(
596 const unsigned int overlap)
597 : overlap(overlap)
598 {}
599
600
601
602 void
603 PreconditionBlockwiseDirect::initialize(const SparseMatrix &matrix,
604 const AdditionalData &additional_data)
605 {
606 preconditioner.reset();
607 preconditioner.reset(
608 Ifpack().Create("Amesos",
609 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
610 additional_data.overlap));
611
612 Ifpack_Preconditioner *ifpack =
613 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
614 Assert(ifpack != nullptr,
615 ExcMessage("Trilinos could not create this "
616 "preconditioner"));
617
618 int ierr;
619
620 Teuchos::ParameterList parameter_list;
621 parameter_list.set("schwarz: combine mode", "Add");
622
623 ierr = ifpack->SetParameters(parameter_list);
624 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
625
626 ierr = ifpack->Initialize();
627 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
628
629 ierr = ifpack->Compute();
630 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
631 }
632
633
634
635 /* ---------------------- PreconditionBlockDirect --------------------- */
636
637 PreconditionChebyshev::AdditionalData::AdditionalData(
638 const unsigned int degree,
639 const double max_eigenvalue,
640 const double eigenvalue_ratio,
641 const double min_eigenvalue,
642 const double min_diagonal,
643 const bool nonzero_starting)
644 : degree(degree)
645 , max_eigenvalue(max_eigenvalue)
646 , eigenvalue_ratio(eigenvalue_ratio)
647 , min_eigenvalue(min_eigenvalue)
648 , min_diagonal(min_diagonal)
649 , nonzero_starting(nonzero_starting)
650 {}
651
652
653
654 void
655 PreconditionChebyshev::initialize(const SparseMatrix &matrix,
656 const AdditionalData &additional_data)
657 {
658 preconditioner =
659 Teuchos::rcp(new Ifpack_Chebyshev(&matrix.trilinos_matrix()));
660
661 Ifpack_Chebyshev *ifpack =
662 dynamic_cast<Ifpack_Chebyshev *>(preconditioner.get());
663 Assert(ifpack != nullptr,
664 ExcMessage("Trilinos could not create this "
665 "preconditioner"));
666
667 int ierr;
668
669 Teuchos::ParameterList parameter_list;
670 parameter_list.set("chebyshev: ratio eigenvalue",
671 additional_data.eigenvalue_ratio);
672 parameter_list.set("chebyshev: min eigenvalue",
673 additional_data.min_eigenvalue);
674 parameter_list.set("chebyshev: max eigenvalue",
675 additional_data.max_eigenvalue);
676 parameter_list.set("chebyshev: degree",
677 static_cast<int>(additional_data.degree));
678 parameter_list.set("chebyshev: min diagonal value",
679 additional_data.min_diagonal);
680 parameter_list.set("chebyshev: zero starting solution",
681 !additional_data.nonzero_starting);
682
683 ierr = ifpack->SetParameters(parameter_list);
684 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
685
686 ierr = ifpack->Initialize();
687 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
688
689 ierr = ifpack->Compute();
690 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
691 }
692
693
694
695 /* -------------------------- PreconditionIdentity --------------------- */
696
697 void
698 PreconditionIdentity::initialize(const SparseMatrix &matrix,
699 const AdditionalData &)
700 {
701 // What follows just configures a dummy preconditioner that
702 // sets up the domain and range maps, as well as the communicator.
703 // It is never used as the vmult, Tvmult operations are
704 // given a custom definition.
705 // Note: This is only required in order to wrap this
706 // preconditioner in a LinearOperator without an exemplar
707 // matrix.
708
709 // From PreconditionJacobi:
710 // release memory before reallocation
711 preconditioner.reset();
712 preconditioner.reset(
713 Ifpack().Create("point relaxation",
714 const_cast<Epetra_CrsMatrix *>(&matrix.trilinos_matrix()),
715 0));
716
717 Ifpack_Preconditioner *ifpack =
718 dynamic_cast<Ifpack_Preconditioner *>(preconditioner.get());
719 Assert(ifpack != nullptr,
720 ExcMessage("Trilinos could not create this "
721 "preconditioner"));
722
723 int ierr;
724
725 Teuchos::ParameterList parameter_list;
726 parameter_list.set("relaxation: sweeps", 1);
727 parameter_list.set("relaxation: type", "Jacobi");
728 parameter_list.set("relaxation: damping factor", 1.0);
729 parameter_list.set("relaxation: min diagonal value", 0.0);
730
731 ierr = ifpack->SetParameters(parameter_list);
732 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
733
734 ierr = ifpack->Initialize();
735 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
736
737 ierr = ifpack->Compute();
738 AssertThrow(ierr == 0, ExcTrilinosError(ierr));
739 }
740
741 void
742 PreconditionIdentity::vmult(MPI::Vector &dst, const MPI::Vector &src) const
743 {
744 dst = src;
745 }
746
747 void
748 PreconditionIdentity::Tvmult(MPI::Vector &dst, const MPI::Vector &src) const
749 {
750 dst = src;
751 }
752
753 void
754 PreconditionIdentity::vmult(::Vector<double> &dst,
755 const ::Vector<double> &src) const
756 {
757 dst = src;
758 }
759
760 void
761 PreconditionIdentity::Tvmult(::Vector<double> &dst,
762 const ::Vector<double> &src) const
763 {
764 dst = src;
765 }
766
767# ifndef DOXYGEN
768 void
769 PreconditionIdentity::vmult(
770 LinearAlgebra::distributed::Vector<double> &dst,
771 const LinearAlgebra::distributed::Vector<double> &src) const
772 {
773 dst = src;
774 }
775
776 void
777 PreconditionIdentity::Tvmult(
778 LinearAlgebra::distributed::Vector<double> &dst,
779 const LinearAlgebra::distributed::Vector<double> &src) const
780 {
781 dst = src;
782 }
783# endif // DOXYGEN
784} // namespace TrilinosWrappers
785
786DEAL_II_NAMESPACE_CLOSE
787
788#endif // DEAL_II_WITH_TRILINOS
Epetra_Operator
MPI_Comm
TrilinosWrappers::PreconditionBase::trilinos_operator
Epetra_Operator & trilinos_operator() const
Definition trilinos_precondition.cc:54
TrilinosWrappers::PreconditionBase::preconditioner
Teuchos::RCP< Epetra_Operator > preconditioner
Definition trilinos_precondition.h:222
TrilinosWrappers::PreconditionBlockJacobi::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:251
TrilinosWrappers::PreconditionBlockSOR::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:389
TrilinosWrappers::PreconditionBlockSSOR::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:320
TrilinosWrappers::PreconditionBlockwiseDirect::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:603
TrilinosWrappers::PreconditionChebyshev::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:655
TrilinosWrappers::PreconditionIC::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:453
TrilinosWrappers::PreconditionILUT::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:556
TrilinosWrappers::PreconditionILU::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:503
TrilinosWrappers::PreconditionIdentity::vmult
void vmult(MPI::Vector &dst, const MPI::Vector &src) const override
Definition trilinos_precondition.cc:742
TrilinosWrappers::PreconditionIdentity::Tvmult
void Tvmult(MPI::Vector &dst, const MPI::Vector &src) const override
Definition trilinos_precondition.cc:748
TrilinosWrappers::PreconditionIdentity::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:698
TrilinosWrappers::PreconditionJacobi::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:89
TrilinosWrappers::PreconditionSOR::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:195
TrilinosWrappers::PreconditionSSOR::initialize
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
Definition trilinos_precondition.cc:142
DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_OPEN
Definition config.h:498
DEAL_II_NAMESPACE_CLOSE
#define DEAL_II_NAMESPACE_CLOSE
Definition config.h:499
Assert
#define Assert(cond, exc)
Definition exceptions.h:1659
LACExceptions::ExcTrilinosError
static ::ExceptionBase & ExcTrilinosError(int arg1)
StandardExceptions::ExcMessage
static ::ExceptionBase & ExcMessage(std::string arg1)
AssertThrow
#define AssertThrow(cond, exc)
Definition exceptions.h:1780
sparse_matrix.h
TrilinosWrappers::PreconditionBlockJacobi::AdditionalData::AdditionalData
AdditionalData(const unsigned int block_size=1, const std::string &block_creation_type="linear", const double omega=1, const double min_diagonal=0, const unsigned int n_sweeps=1)
Definition trilinos_precondition.cc:235
TrilinosWrappers::PreconditionBlockSOR::AdditionalData::AdditionalData
AdditionalData(const unsigned int block_size=1, const std::string &block_creation_type="linear", const double omega=1, const double min_diagonal=0, const unsigned int overlap=0, const unsigned int n_sweeps=1)
Definition trilinos_precondition.cc:371
TrilinosWrappers::PreconditionBlockSSOR::AdditionalData::AdditionalData
AdditionalData(const unsigned int block_size=1, const std::string &block_creation_type="linear", const double omega=1, const double min_diagonal=0, const unsigned int overlap=0, const unsigned int n_sweeps=1)
Definition trilinos_precondition.cc:302
TrilinosWrappers::PreconditionChebyshev::AdditionalData::AdditionalData
AdditionalData(const unsigned int degree=1, const double max_eigenvalue=10., const double eigenvalue_ratio=30., const double min_eigenvalue=1., const double min_diagonal=1e-12, const bool nonzero_starting=false)
Definition trilinos_precondition.cc:637
TrilinosWrappers::PreconditionIC::AdditionalData::AdditionalData
AdditionalData(const unsigned int ic_fill=0, const double ic_atol=0., const double ic_rtol=1., const unsigned int overlap=0)
Definition trilinos_precondition.cc:440
TrilinosWrappers::PreconditionILUT::AdditionalData::AdditionalData
AdditionalData(const double ilut_drop=0., const unsigned int ilut_fill=0, const double ilut_atol=0., const double ilut_rtol=1., const unsigned int overlap=0)
Definition trilinos_precondition.cc:541
TrilinosWrappers::PreconditionILU::AdditionalData::AdditionalData
AdditionalData(const unsigned int ilu_fill=0, const double ilu_atol=0., const double ilu_rtol=1., const unsigned int overlap=0)
Definition trilinos_precondition.cc:490
TrilinosWrappers::PreconditionJacobi::AdditionalData::AdditionalData
AdditionalData(const double omega=1, const double min_diagonal=0, const unsigned int n_sweeps=1)
Definition trilinos_precondition.cc:77
TrilinosWrappers::PreconditionSOR::AdditionalData::AdditionalData
AdditionalData(const double omega=1, const double min_diagonal=0, const unsigned int overlap=0, const unsigned int n_sweeps=1)
Definition trilinos_precondition.cc:182
TrilinosWrappers::PreconditionSSOR::AdditionalData::AdditionalData
AdditionalData(const double omega=1, const double min_diagonal=0, const unsigned int overlap=0, const unsigned int n_sweeps=1)
Definition trilinos_precondition.cc:129
trilinos_precondition.h
trilinos_sparse_matrix.h
Generated by doxygen 1.9.8