]> https://gitweb.dealii.org/ - dealii.git/log
dealii.git
2 years agoMerge pull request #14447 from tjhei/fix-indent-script
David Wells [Mon, 21 Nov 2022 18:11:13 +0000 (13:11 -0500)]
Merge pull request #14447 from tjhei/fix-indent-script

fix indent script

2 years agoMerge pull request #14446 from bangerth/vtu-5
David Wells [Mon, 21 Nov 2022 17:40:14 +0000 (12:40 -0500)]
Merge pull request #14446 from bangerth/vtu-5

Simplify code

2 years agofix indent script 14447/head
Timo Heister [Mon, 21 Nov 2022 15:44:06 +0000 (10:44 -0500)]
fix indent script

2 years agoMerge pull request #14444 from bangerth/ident
David Wells [Mon, 21 Nov 2022 14:33:26 +0000 (09:33 -0500)]
Merge pull request #14444 from bangerth/ident

Add a bit of documentation.

2 years agoMerge pull request #14429 from masterleinad/cuda_build_tests
Daniel Arndt [Mon, 21 Nov 2022 13:06:02 +0000 (08:06 -0500)]
Merge pull request #14429 from masterleinad/cuda_build_tests

2 years agoConvert yet another function. 14446/head
Wolfgang Bangerth [Sun, 20 Nov 2022 01:23:44 +0000 (18:23 -0700)]
Convert yet another function.

2 years agoSimplify yet more code.
Wolfgang Bangerth [Fri, 18 Nov 2022 23:03:09 +0000 (16:03 -0700)]
Simplify yet more code.

2 years agoSimplify more code.
Wolfgang Bangerth [Fri, 18 Nov 2022 22:58:45 +0000 (15:58 -0700)]
Simplify more code.

2 years agoAdd a bit of documentation. 14444/head
Wolfgang Bangerth [Fri, 18 Nov 2022 16:17:37 +0000 (09:17 -0700)]
Add a bit of documentation.

2 years agoMerge pull request #14424 from masterleinad/kokkos_print_backends_and_archs
Matthias Maier [Sat, 19 Nov 2022 03:27:43 +0000 (21:27 -0600)]
Merge pull request #14424 from masterleinad/kokkos_print_backends_and_archs

Print backends and architectures for Kokkos

2 years agoPrint backends and architectures for Kokkos 14424/head
Daniel Arndt [Tue, 15 Nov 2022 18:15:10 +0000 (12:15 -0600)]
Print backends and architectures for Kokkos

2 years agoBuild CUDA tests in CUDA CI 14429/head
Daniel Arndt [Sat, 19 Nov 2022 00:49:41 +0000 (18:49 -0600)]
Build CUDA tests in CUDA CI

2 years agoSimplify some code.
Wolfgang Bangerth [Fri, 18 Nov 2022 22:51:49 +0000 (15:51 -0700)]
Simplify some code.

2 years agoMerge pull request #14442 from marcfehling/vectorcopy
David Wells [Fri, 18 Nov 2022 04:01:11 +0000 (23:01 -0500)]
Merge pull request #14442 from marcfehling/vectorcopy

Use copy function internally in LA::d::Vector.

2 years agoMerge pull request #14438 from bangerth/petsc
Daniel Arndt [Thu, 17 Nov 2022 23:21:49 +0000 (17:21 -0600)]
Merge pull request #14438 from bangerth/petsc

2 years agoRe-introduce (now deprecated) old constructor. 14438/head
Wolfgang Bangerth [Thu, 17 Nov 2022 20:46:24 +0000 (13:46 -0700)]
Re-introduce (now deprecated) old constructor.

2 years agoMerge pull request #14436 from Rombur/cuda_ice
Timo Heister [Thu, 17 Nov 2022 20:28:15 +0000 (15:28 -0500)]
Merge pull request #14436 from Rombur/cuda_ice

Update CI workflow to CUDA 11.8

2 years agoMerge pull request #14421 from bangerth/vtu-3
Bruno Turcksin [Thu, 17 Nov 2022 16:28:05 +0000 (11:28 -0500)]
Merge pull request #14421 from bangerth/vtu-3

2 years agoInstall boost 14436/head
Bruno Turcksin [Thu, 17 Nov 2022 15:57:14 +0000 (15:57 +0000)]
Install boost

2 years agoAdd a changelog.
Wolfgang Bangerth [Thu, 17 Nov 2022 04:33:33 +0000 (21:33 -0700)]
Add a changelog.

2 years agoSimplify step-40.
Wolfgang Bangerth [Thu, 17 Nov 2022 00:19:50 +0000 (17:19 -0700)]
Simplify step-40.

2 years agoAvoid the use of deprecated constructors.
Wolfgang Bangerth [Thu, 17 Nov 2022 00:19:37 +0000 (17:19 -0700)]
Avoid the use of deprecated constructors.

2 years agoDeprecate the use of PETScWrapper::Solver* constructors that take an MPI communicator.
Wolfgang Bangerth [Thu, 17 Nov 2022 00:13:10 +0000 (17:13 -0700)]
Deprecate the use of PETScWrapper::Solver* constructors that take an MPI communicator.

2 years agoSimplify code.
Wolfgang Bangerth [Wed, 16 Nov 2022 21:58:27 +0000 (14:58 -0700)]
Simplify code.

2 years agoRemove the need to store the MPI communicator in PETScWrappers::SolverBase.
Wolfgang Bangerth [Wed, 16 Nov 2022 21:52:55 +0000 (14:52 -0700)]
Remove the need to store the MPI communicator in PETScWrappers::SolverBase.

2 years agoProvide the ability to query the communicator from PETSc precondiditioners.
Wolfgang Bangerth [Wed, 16 Nov 2022 21:52:19 +0000 (14:52 -0700)]
Provide the ability to query the communicator from PETSc precondiditioners.

2 years agoMerge pull request #14434 from drwells/sparsity-pattern-base-dof-tools
Daniel Arndt [Thu, 17 Nov 2022 15:29:35 +0000 (09:29 -0600)]
Merge pull request #14434 from drwells/sparsity-pattern-base-dof-tools

2 years agoMerge pull request #14439 from bangerth/comment
Timo Heister [Thu, 17 Nov 2022 12:47:41 +0000 (07:47 -0500)]
Merge pull request #14439 from bangerth/comment

Fix a misspelling in a comment.

2 years agoFix a misspelling in a comment. 14439/head
Wolfgang Bangerth [Thu, 17 Nov 2022 04:51:27 +0000 (21:51 -0700)]
Fix a misspelling in a comment.

2 years agoDoFTools: use some ranged-for loops 14434/head
David Wells [Wed, 16 Nov 2022 23:06:22 +0000 (18:06 -0500)]
DoFTools: use some ranged-for loops

2 years agoUpdate CI workflow to CUDA 11.8
Bruno Turcksin [Wed, 16 Nov 2022 20:23:07 +0000 (20:23 +0000)]
Update CI workflow to CUDA 11.8

2 years agoDoFTools: convert some functions to use add_row_entries().
David Wells [Wed, 16 Nov 2022 18:32:52 +0000 (13:32 -0500)]
DoFTools: convert some functions to use add_row_entries().

2 years agoSparsityPattern: add an assertion in input data.
David Wells [Wed, 16 Nov 2022 18:30:31 +0000 (13:30 -0500)]
SparsityPattern: add an assertion in input data.

2 years agoMerge pull request #14427 from marcfehling/reinit
Martin Kronbichler [Wed, 16 Nov 2022 15:38:44 +0000 (16:38 +0100)]
Merge pull request #14427 from marcfehling/reinit

Make TrilinosWrappers::MPI::Vector::reinit accept Utilities::MPI::Partitioner.

2 years agoMerge pull request #14432 from drwells/sparsity-pattern-base-particles
Martin Kronbichler [Wed, 16 Nov 2022 15:37:39 +0000 (16:37 +0100)]
Merge pull request #14432 from drwells/sparsity-pattern-base-particles

Particles: Remove template on sparsity pattern type.

2 years agoMerge pull request #14431 from drwells/sparsity-pattern-base-nonmatching
Daniel Arndt [Wed, 16 Nov 2022 05:02:32 +0000 (23:02 -0600)]
Merge pull request #14431 from drwells/sparsity-pattern-base-nonmatching

2 years agoMerge pull request #14425 from Rombur/belos
Wolfgang Bangerth [Wed, 16 Nov 2022 03:28:37 +0000 (20:28 -0700)]
Merge pull request #14425 from Rombur/belos

Guard Belos code with macros

2 years agoDoFTools: Remove template on sparsity pattern type.
David Wells [Sat, 5 Nov 2022 15:42:24 +0000 (11:42 -0400)]
DoFTools: Remove template on sparsity pattern type.

2 years agoParticles: Remove template on sparsity pattern type. 14432/head
David Wells [Sat, 5 Nov 2022 02:33:46 +0000 (22:33 -0400)]
Particles: Remove template on sparsity pattern type.

2 years agoNonMatching: Remove template on sparsity pattern type. 14431/head
David Wells [Sat, 5 Nov 2022 02:10:56 +0000 (22:10 -0400)]
NonMatching: Remove template on sparsity pattern type.

This is no longer necessary since all of our SparsityPatterns inherit from
SparsityPatternBase. This file was previously the most expensive to compile in
the entire library so cutting down the number of instantiations by a factor of 4
or 6 (with Trilinos) is quite nice.

2 years agoMerge pull request #14415 from tamiko/ninja_workaround
Daniel Arndt [Tue, 15 Nov 2022 21:40:52 +0000 (15:40 -0600)]
Merge pull request #14415 from tamiko/ninja_workaround

2 years agoUse copy function internally in LA::d::Vector. 14442/head
Marc Fehling [Tue, 15 Nov 2022 21:40:19 +0000 (14:40 -0700)]
Use copy function internally in LA::d::Vector.

2 years agoMerge pull request #14423 from tjhei/rename-inst-file
Daniel Arndt [Tue, 15 Nov 2022 20:34:01 +0000 (14:34 -0600)]
Merge pull request #14423 from tjhei/rename-inst-file

2 years agoMake TrilinosWrappers::MPI::Vector::reinit accept Utilities::MPI::Partitioner. 14427/head
Marc Fehling [Tue, 15 Nov 2022 20:13:33 +0000 (13:13 -0700)]
Make TrilinosWrappers::MPI::Vector::reinit accept Utilities::MPI::Partitioner.

2 years agoGuard Belos code with macros 14425/head
Bruno Turcksin [Tue, 15 Nov 2022 19:38:28 +0000 (19:38 +0000)]
Guard Belos code with macros

2 years agoMerge pull request #14422 from tjhei/remove-cu-doc-stuff
Daniel Arndt [Tue, 15 Nov 2022 18:28:30 +0000 (12:28 -0600)]
Merge pull request #14422 from tjhei/remove-cu-doc-stuff

2 years agoMerge pull request #14420 from drwells/de-template-af-local-to-global
Martin Kronbichler [Tue, 15 Nov 2022 17:57:03 +0000 (18:57 +0100)]
Merge pull request #14420 from drwells/de-template-af-local-to-global

De template AffineConstraints::add_entries_local_to_global()

2 years agoTry to use std::make_tuple() to get around compiler warnings. 14421/head
Wolfgang Bangerth [Tue, 15 Nov 2022 16:58:47 +0000 (09:58 -0700)]
Try to use std::make_tuple() to get around compiler warnings.

2 years agorename an inst file 14423/head
Timo Heister [Tue, 15 Nov 2022 15:04:38 +0000 (10:04 -0500)]
rename an inst file

2 years agoremove .cu extension in program2doxygen 14422/head
Timo Heister [Tue, 15 Nov 2022 14:48:34 +0000 (09:48 -0500)]
remove .cu extension in program2doxygen

2 years agoWork around a strange compiler issue.
Wolfgang Bangerth [Tue, 15 Nov 2022 14:47:03 +0000 (07:47 -0700)]
Work around a strange compiler issue.

The MPI build bot complains that using a {...} initializer list would require the use
of a conversion constructor marked as 'explicit':

data_out_base.cc:1156:29: error: converting to ‘std::tuple<unsigned int, unsigned int>’ from initializer list would use explicit constructor ‘constexpr std::tuple<_T1, _T2>::tuple(_U1&&, _U2&&) [with _U1 = unsigned int&; _U2 = unsigned int&; <template-parameter-2-3> = void; _T1 = unsigned int; _T2 = unsigned int]’

That makes no sense, but let's appease the compiler by being explicit.

2 years agoMerge pull request #14377 from masterleinad/move_cuda_to_nvcc_wrapper_kokkos
Timo Heister [Tue, 15 Nov 2022 14:36:47 +0000 (09:36 -0500)]
Merge pull request #14377 from masterleinad/move_cuda_to_nvcc_wrapper_kokkos

Move CUDA compilation to Kokkos with nvcc_wrapper

2 years agoMake a function return its object by value, rather than via reference.
Wolfgang Bangerth [Thu, 10 Nov 2022 22:20:53 +0000 (15:20 -0700)]
Make a function return its object by value, rather than via reference.

2 years agoRename a function and associated objects.
Wolfgang Bangerth [Thu, 10 Nov 2022 21:56:54 +0000 (14:56 -0700)]
Rename a function and associated objects.

2 years agoSimplify the call to a function on a separate task.
Wolfgang Bangerth [Thu, 10 Nov 2022 21:30:45 +0000 (14:30 -0700)]
Simplify the call to a function on a separate task.

2 years agoRename and refactor another function.
Wolfgang Bangerth [Thu, 10 Nov 2022 21:13:27 +0000 (14:13 -0700)]
Rename and refactor another function.

2 years agoSimplify testing whether a cell is a hypercube.
Wolfgang Bangerth [Thu, 10 Nov 2022 20:48:47 +0000 (13:48 -0700)]
Simplify testing whether a cell is a hypercube.

2 years agoDon't explicitly specify template arguments where not necessary.
Wolfgang Bangerth [Thu, 10 Nov 2022 18:57:36 +0000 (11:57 -0700)]
Don't explicitly specify template arguments where not necessary.

2 years agoRename a variable. Unname another.
Wolfgang Bangerth [Thu, 10 Nov 2022 18:56:16 +0000 (11:56 -0700)]
Rename a variable. Unname another.

2 years agoLet a function return by-value, rather than through multiple reference arguments.
Wolfgang Bangerth [Thu, 10 Nov 2022 18:52:53 +0000 (11:52 -0700)]
Let a function return by-value, rather than through multiple reference arguments.

2 years agoAffineConstraints: de-template add_entries_local_to_global() part 2 14420/head
David Wells [Thu, 3 Nov 2022 18:21:18 +0000 (14:21 -0400)]
AffineConstraints: de-template add_entries_local_to_global() part 2

2 years agoAffineConstraints: de-template add_entries_local_to_global() part 1
David Wells [Thu, 3 Nov 2022 17:38:34 +0000 (13:38 -0400)]
AffineConstraints: de-template add_entries_local_to_global() part 1

2 years agoAffineConstraints: de-template set_sparsity_diagonals().
David Wells [Thu, 3 Nov 2022 16:50:51 +0000 (12:50 -0400)]
AffineConstraints: de-template set_sparsity_diagonals().

2 years agoMerge pull request #14418 from drwells/tw-sparsity-pattern-base
Daniel Arndt [Mon, 14 Nov 2022 16:02:28 +0000 (11:02 -0500)]
Merge pull request #14418 from drwells/tw-sparsity-pattern-base

2 years agoMerge pull request #14398 from drwells/remove-bsp-constraint-function
Daniel Arndt [Mon, 14 Nov 2022 15:59:06 +0000 (10:59 -0500)]
Merge pull request #14398 from drwells/remove-bsp-constraint-function

2 years agoMerge pull request #14419 from drwells/bsp-block-sparsity-pattern-base
Martin Kronbichler [Mon, 14 Nov 2022 09:07:17 +0000 (10:07 +0100)]
Merge pull request #14419 from drwells/bsp-block-sparsity-pattern-base

BlockSparsityPatternBase: inherit from SparsityPatternBase.

2 years agoBlockSparsityPatternBase: inherit from SparsityPatternBase. 14419/head
David Wells [Sun, 30 Oct 2022 13:05:50 +0000 (09:05 -0400)]
BlockSparsityPatternBase: inherit from SparsityPatternBase.

2 years agoTrilinosWrappers::SparsityPattern: inherit from SparsityPatternBase. 14418/head
David Wells [Thu, 3 Nov 2022 16:40:14 +0000 (12:40 -0400)]
TrilinosWrappers::SparsityPattern: inherit from SparsityPatternBase.

2 years agoTrilinosWrappers::SparsityPattern: clean up whitespace.
David Wells [Tue, 1 Nov 2022 12:51:53 +0000 (08:51 -0400)]
TrilinosWrappers::SparsityPattern: clean up whitespace.

2 years agoOptimize BlockSparsityPatternBase::add_entries(). 14398/head
David Wells [Sun, 13 Nov 2022 21:06:21 +0000 (16:06 -0500)]
Optimize BlockSparsityPatternBase::add_entries().

We can avoid some expensive parts and looping over data more than once by
utilizing the fact that the DoFs are sorted, which implies that they are also
sorted by block.

2 years agoMerge pull request #14406 from bangerth/vtu-2
David Wells [Sun, 13 Nov 2022 21:12:26 +0000 (16:12 -0500)]
Merge pull request #14406 from bangerth/vtu-2

More refactoring of the VTU output.

2 years agoMerge pull request #14413 from tcojean/fix_findginkgo_cmake
Martin Kronbichler [Sun, 13 Nov 2022 16:23:46 +0000 (17:23 +0100)]
Merge pull request #14413 from tcojean/fix_findginkgo_cmake

Fix FindGINKGO.cmake for Ginkgo 1.5.0

2 years agoMerge pull request #14384 from kronbichler/performance_test
Peter Munch [Sat, 12 Nov 2022 20:40:41 +0000 (21:40 +0100)]
Merge pull request #14384 from kronbichler/performance_test

Add performance test showing compressible Navier-Stokes equations

2 years agoAffineConstraints: remove the block sparsity constraint function.
David Wells [Sun, 6 Nov 2022 13:13:25 +0000 (08:13 -0500)]
AffineConstraints: remove the block sparsity constraint function.

This function predates BlockSparsityPatternBase<T>::add_entries(), which does
the same thing (essentially just updating the individual blocks after computing
offsets). For good measure I added a test which verifies that we get identical
output with constraints and a DoF mask table before and after the switch.

2 years agoMerge pull request #14414 from lpsaavedra/fix_small_typo
Peter Munch [Sat, 12 Nov 2022 06:32:46 +0000 (07:32 +0100)]
Merge pull request #14414 from lpsaavedra/fix_small_typo

Fix small typos in the intro.dox of Step-66

2 years agoadd a changes entry 14415/head
Matthias Maier [Fri, 11 Nov 2022 23:47:10 +0000 (17:47 -0600)]
add a changes entry

2 years agoCMake: add a compile_test_executables target
Matthias Maier [Fri, 11 Nov 2022 23:37:12 +0000 (17:37 -0600)]
CMake: add a compile_test_executables target

This commit adds a top-level compile_test_executables and modifies the
ADD_TEST() macro to let it depend on all executables generated by our
DEAL_II_PICKUP_TESTS() mechanism.

This is mainly intended for user projects and works around the issue
that simultaneous calls to ninja are not reentrant, meaning when
configuring a user project with ninja one is otherwise limited to run
tests in serial.

2 years agoFix small typos in intro.dox 14414/head
Laura Prieto Saavedra [Fri, 11 Nov 2022 19:41:01 +0000 (14:41 -0500)]
Fix small typos in intro.dox

2 years agoMention reason for change. 14413/head
Terry Cojean [Fri, 11 Nov 2022 12:18:19 +0000 (13:18 +0100)]
Mention reason for change.

Co-authored-by: Martin Kronbichler <martin.kronbichler@uni-a.de>
2 years agoFix FindGINKGO.cmake for Ginkgo 1.5.0
Terry Cojean [Fri, 11 Nov 2022 12:04:41 +0000 (13:04 +0100)]
Fix FindGINKGO.cmake for Ginkgo 1.5.0

2 years agoMerge pull request #14396 from drwells/add-sparsity-pattern-base
Martin Kronbichler [Fri, 11 Nov 2022 07:45:18 +0000 (08:45 +0100)]
Merge pull request #14396 from drwells/add-sparsity-pattern-base

(Dynamic)SparsityPattern: inherit from SparsityPatternBase.

2 years agoMerge pull request #14409 from bangerth/x-1
Daniel Arndt [Thu, 10 Nov 2022 23:02:34 +0000 (18:02 -0500)]
Merge pull request #14409 from bangerth/x-1

2 years agoMerge pull request #14410 from bangerth/x-3
David Wells [Thu, 10 Nov 2022 22:55:12 +0000 (17:55 -0500)]
Merge pull request #14410 from bangerth/x-3

Remove comments that don't actually explain anything.

2 years agoUse a range-based for loop. 14406/head
Wolfgang Bangerth [Thu, 10 Nov 2022 21:31:49 +0000 (14:31 -0700)]
Use a range-based for loop.

2 years agoRemove comments that don't actually explain anything. 14410/head
Wolfgang Bangerth [Thu, 10 Nov 2022 18:46:22 +0000 (11:46 -0700)]
Remove comments that don't actually explain anything.

2 years agoDo not mark functions in anonymous namespaces as 'static'. 14409/head
Wolfgang Bangerth [Thu, 10 Nov 2022 18:40:43 +0000 (11:40 -0700)]
Do not mark functions in anonymous namespaces as 'static'.

2 years agoMerge pull request #14407 from bangerth/doc-3
David Wells [Thu, 10 Nov 2022 03:03:41 +0000 (22:03 -0500)]
Merge pull request #14407 from bangerth/doc-3

Add a documentation note.

2 years agoMerge pull request #14408 from bangerth/doc-2
David Wells [Thu, 10 Nov 2022 03:03:12 +0000 (22:03 -0500)]
Merge pull request #14408 from bangerth/doc-2

Fix a mistake in the documentation.

2 years agoRe-factor the use of write_nodes() for VTU output.
Wolfgang Bangerth [Wed, 9 Nov 2022 04:30:58 +0000 (21:30 -0700)]
Re-factor the use of write_nodes() for VTU output.

As part of this, just output the vertex locations as one big array (like all
other data), rather than breaking it every three coordinates onto a
separate line.

2 years agoDo not use VtuStream to output data object; do it directly.
Wolfgang Bangerth [Tue, 8 Nov 2022 22:09:37 +0000 (15:09 -0700)]
Do not use VtuStream to output data object; do it directly.

2 years agoFix a mistake in the documentation. 14408/head
Wolfgang Bangerth [Thu, 10 Nov 2022 00:20:17 +0000 (17:20 -0700)]
Fix a mistake in the documentation.

2 years agoAdd a documentation note. 14407/head
Wolfgang Bangerth [Thu, 10 Nov 2022 00:22:52 +0000 (17:22 -0700)]
Add a documentation note.

2 years agoMerge pull request #14397 from drwells/cleanup-bsp
Martin Kronbichler [Wed, 9 Nov 2022 09:11:41 +0000 (10:11 +0100)]
Merge pull request #14397 from drwells/cleanup-bsp

Cleanup BSP

2 years agoMerge pull request #14404 from drwells/update-cuda-tester-backports
Martin Kronbichler [Wed, 9 Nov 2022 09:10:00 +0000 (10:10 +0100)]
Merge pull request #14404 from drwells/update-cuda-tester-backports

CI: Update backports.

2 years agoMerge pull request #14402 from bangerth/vtu
Martin Kronbichler [Wed, 9 Nov 2022 09:09:39 +0000 (10:09 +0100)]
Merge pull request #14402 from bangerth/vtu

Refactor compressing arrays for VTU output.

2 years agoCI: Update backports. 14404/head
David Wells [Tue, 8 Nov 2022 22:52:48 +0000 (17:52 -0500)]
CI: Update backports.

We may as well use the most recent dependency backports.

2 years agoMerge pull request #14401 from bangerth/doc-3
Daniel Arndt [Tue, 8 Nov 2022 19:32:39 +0000 (14:32 -0500)]
Merge pull request #14401 from bangerth/doc-3

2 years agoRefactor compressing arrays for VTU output. 14402/head
Wolfgang Bangerth [Tue, 8 Nov 2022 17:26:14 +0000 (10:26 -0700)]
Refactor compressing arrays for VTU output.

2 years agoMinor fix to aid understand what a sentence means. 14401/head
Wolfgang Bangerth [Tue, 8 Nov 2022 15:53:48 +0000 (08:53 -0700)]
Minor fix to aid understand what a sentence means.

2 years agoAdd BlockSparsityPatternBase::memory_consumption(). 14397/head
David Wells [Sun, 30 Oct 2022 00:59:46 +0000 (20:59 -0400)]
Add BlockSparsityPatternBase::memory_consumption().

None of the inheriting classes add any class members so this should be here
instead of in BlockSparsityPattern.


In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.