mform/m_mform_constraint_boundary.F90
| Line | Branch | Exec | Source |
|---|---|---|---|
| 1 | !> @brief Module containing implementations of various local m-form constraints, such as: | ||
| 2 | !> - MFormDirichlet, which imposes a Dirichlet boundary condition on an m-form | ||
| 3 | !> - MFormNeumann, which imposes a Neumann boundary condition on a 0-form | ||
| 4 | !> - MFormUniqueness constraint, which imposes a uniqueness constraint on an m-form (e.g., to fix the constant mode in a purely periodic Poisson problem) | ||
| 5 | module m_mform_constraint_boundary | ||
| 6 | use m_bspline, only: BSplineFun | ||
| 7 | use m_common, only: wp | ||
| 8 | use m_mform_basis, only: MFormSpace | ||
| 9 | use m_sparsemat, only: SparseMat | ||
| 10 | use m_tensorprod_indices, only: TensorProdIndices, TensorProdIndexList | ||
| 11 | use m_mform_constraint_abstract, only: MFormConstraintLocal | ||
| 12 | #include "petsc.fi" | ||
| 13 | #include "timer.fi" | ||
| 14 | |||
| 15 | implicit none | ||
| 16 | |||
| 17 | private | ||
| 18 | public :: MFormDirichlet, MFormNeumann, MFormUniqueness | ||
| 19 | |||
| 20 | !> The threshold under which extraction/projection coefficients are considered zero | ||
| 21 | real(wp), parameter :: COEFF_THR = 1e-14_wp | ||
| 22 | |||
| 23 | ! ! > A wrapper for a vector of real numbers | ||
| 24 | ! type VectorWrapper | ||
| 25 | ! !> The data of the vector | ||
| 26 | ! real(wp), allocatable :: data(:) | ||
| 27 | ! end type VectorWrapper | ||
| 28 | |||
| 29 | !> An object that represents a Dirichlet constraint on an m-form space | ||
| 30 | type, extends(MFormConstraintLocal) :: MFormDirichlet | ||
| 31 | |||
| 32 | contains | ||
| 33 | procedure :: init => init_mform_dirichlet | ||
| 34 | procedure :: destroy_derived => destroy_mform_dirichlet | ||
| 35 | end type | ||
| 36 | |||
| 37 | !> An object that represents a Neumann constraint on an m-form space | ||
| 38 | type, extends(MFormConstraintLocal) :: MFormNeumann | ||
| 39 | !> Whether to fix the constant mode in e.g. a Poisson problem | ||
| 40 | logical :: fix_constant | ||
| 41 | contains | ||
| 42 | procedure :: init => init_mform_neumann | ||
| 43 | procedure :: destroy_derived => destroy_mform_neumann | ||
| 44 | end type | ||
| 45 | |||
| 46 | !> An object that represents a uniqueness constraint on an m-form space | ||
| 47 | type, extends(MFormConstraintLocal) :: MFormUniqueness | ||
| 48 | contains | ||
| 49 | procedure :: init => init_mform_uniqueness | ||
| 50 | procedure :: destroy_derived => destroy_mform_uniqueness | ||
| 51 | end type | ||
| 52 | |||
| 53 | ! The constructors are used for derived-type specific parameters, whereas the deferred init method actually initializes the object | ||
| 54 | interface MFormDirichlet | ||
| 55 | module procedure construct_mform_dirichlet | ||
| 56 | end interface | ||
| 57 | |||
| 58 | interface MFormNeumann | ||
| 59 | module procedure construct_mform_neumann | ||
| 60 | end interface | ||
| 61 | |||
| 62 | interface MFormUniqueness | ||
| 63 | module procedure construct_mform_uniqueness | ||
| 64 | end interface | ||
| 65 | |||
| 66 | contains | ||
| 67 | |||
| 68 | !> @brief Construct a Dirichlet constraint on an m-form | ||
| 69 | !> | ||
| 70 | !> @param[in] _(optional)_ bdy_function The function to which the Dirichlet condition is applied (for scalar m-forms) | ||
| 71 | !> @param[in] _(optional)_ bdy_function_x,bdy_function_y,bdy_function_z The functions to which the Dirichlet condition is applied | ||
| 72 | !> (for vector m-forms) | ||
| 73 | !> | ||
| 74 | !> @return The constructed MFormDirichlet object | ||
| 75 | 124 | pure type(MFormDirichlet) function construct_mform_dirichlet(bdy_function, bdy_function_x, bdy_function_y, bdy_function_z) & | |
| 76 | result(this) | ||
| 77 | |||
| 78 | use m_common, only: user_function_3d_interface | ||
| 79 | implicit none | ||
| 80 | |||
| 81 | procedure(user_function_3d_interface), optional :: bdy_function, bdy_function_x, bdy_function_y, bdy_function_z | ||
| 82 | |||
| 83 | ! TODO inhomogeneous Dirichlet condition | ||
| 84 | this%is_homogeneous = .not. (present(bdy_function_x) .or. present(bdy_function_y) .or. present(bdy_function_z) .or. & | ||
| 85 |
2/4✓ Branch 2 → 3 taken 124 times.
✗ Branch 2 → 4 not taken.
✗ Branch 3 → 4 not taken.
✓ Branch 3 → 5 taken 124 times.
|
124 | & present(bdy_function)) |
| 86 |
1/2✗ Branch 5 → 6 not taken.
✓ Branch 5 → 7 taken 124 times.
|
124 | if (.not. this%is_homogeneous) then |
| 87 | ! TODO inhomogeneous Dirichlet condition requires some kind of L2 boundary projection | ||
| 88 | ✗ | error stop "MFormDirichlet::init_mform_dirichlet: inhomogeneous Dirichlet condition is not implemented yet" | |
| 89 | ! allocate(this%inhomogeneous_values(nr_blocks)) | ||
| 90 | end if | ||
| 91 | 124 | end function construct_mform_dirichlet | |
| 92 | |||
| 93 | !> @brief Construct a Neumann constraint on an m-form | ||
| 94 | !> | ||
| 95 | !> @param[in] _(optional)_ fix_constant Whether to fix the constant mode (default: .false.) | ||
| 96 | !> | ||
| 97 | !> @return The constructed MFormNeumann object | ||
| 98 | 105 | pure type(MFormNeumann) function construct_mform_neumann(fix_constant) result(this) | |
| 99 | implicit none | ||
| 100 | |||
| 101 | logical, intent(in), optional :: fix_constant | ||
| 102 | |||
| 103 |
1/2✗ Branch 2 → 3 not taken.
✓ Branch 2 → 4 taken 105 times.
|
105 | if (present(fix_constant)) then |
| 104 | ✗ | this%fix_constant = fix_constant | |
| 105 | else | ||
| 106 | this%fix_constant = .false. | ||
| 107 | end if | ||
| 108 | 105 | end function construct_mform_neumann | |
| 109 | |||
| 110 | !> @brief Construct a uniqueness constraint on an m-form | ||
| 111 | !> @return The constructed MFormUniqueness object | ||
| 112 | 20 | pure type(MFormUniqueness) function construct_mform_uniqueness() result(this) | |
| 113 | implicit none | ||
| 114 | |||
| 115 | 20 | end function construct_mform_uniqueness | |
| 116 | |||
| 117 | !> @brief Initialize a Dirichlet constraint on an m-form space | ||
| 118 | !> | ||
| 119 | !> @param[inout] this The MFormDirichlet object to initialize | ||
| 120 | !> @param[in] space The MFormSpace to which the Dirichlet constraint applies | ||
| 121 | !> @param[in] _(optional)_ coord_transform The coordinate transform associated with the m-form space | ||
| 122 | 176 | subroutine init_mform_dirichlet(this, space, coord_transform) | |
| 123 | use m_mform_basis, only: MFormSpace | ||
| 124 | use m_coord_transform_abstract, only: CoordTransformAbstract | ||
| 125 | implicit none | ||
| 126 | |||
| 127 | class(MFormDirichlet), intent(inout) :: this | ||
| 128 | type(MFormSpace), intent(in) :: space | ||
| 129 | |||
| 130 | class(CoordTransformAbstract), optional, intent(in) :: coord_transform | ||
| 131 | logical :: coord_has_polar_xy_singularity_ | ||
| 132 | integer :: s1, s2, nr_blocks | ||
| 133 | |||
| 134 | TIMER_DECL(t_total) | ||
| 135 | TIMER_INIT_START(t_total, "MFormDirichlet::init") | ||
| 136 | |||
| 137 | 176 | call this%destroy() | |
| 138 | |||
| 139 |
5/8✓ Branch 3 → 4 taken 176 times.
✗ Branch 3 → 7 not taken.
✓ Branch 4 → 5 taken 176 times.
✗ Branch 4 → 6 not taken.
✓ Branch 7 → 8 taken 176 times.
✗ Branch 7 → 10 not taken.
✓ Branch 8 → 9 taken 42 times.
✓ Branch 8 → 10 taken 134 times.
|
176 | this%space = space |
| 140 |
2/2✓ Branch 11 → 12 taken 528 times.
✓ Branch 11 → 13 taken 176 times.
|
704 | this%is_mpi_shared = .false. |
| 141 | |||
| 142 | coord_has_polar_xy_singularity_ = .false. | ||
| 143 |
3/4✓ Branch 13 → 14 taken 86 times.
✓ Branch 13 → 16 taken 90 times.
✓ Branch 14 → 15 taken 86 times.
✗ Branch 14 → 16 not taken.
|
176 | if (present(coord_transform)) coord_has_polar_xy_singularity_ = coord_transform%has_polar_xy_singularity |
| 144 | |||
| 145 | 292 | select case (space%m) | |
| 146 | case (0) | ||
| 147 | nr_blocks = 1 | ||
| 148 |
4/6✗ Branch 17 → 18 not taken.
✓ Branch 17 → 19 taken 116 times.
✗ Branch 19 → 20 not taken.
✓ Branch 19 → 21 taken 116 times.
✓ Branch 22 → 23 taken 116 times.
✓ Branch 22 → 24 taken 116 times.
|
232 | allocate (this%index_list(nr_blocks)) |
| 149 | call init_tp_boundary_indices(this%index_list(1), space%tp_spaces(1), & | ||
| 150 | 116 | include_x_left=.not. coord_has_polar_xy_singularity_) | |
| 151 | case (1) | ||
| 152 | nr_blocks = 3 | ||
| 153 | ! Boundary condition is applied to the two tangential components only | ||
| 154 |
4/6✗ Branch 27 → 28 not taken.
✓ Branch 27 → 29 taken 40 times.
✗ Branch 29 → 30 not taken.
✓ Branch 29 → 31 taken 40 times.
✓ Branch 32 → 33 taken 120 times.
✓ Branch 32 → 34 taken 40 times.
|
160 | allocate (this%index_list(nr_blocks)) |
| 155 |
2/2✓ Branch 35 → 36 taken 40 times.
✓ Branch 35 → 38 taken 120 times.
|
160 | do s1 = 1, 3 |
| 156 | call init_tp_boundary_indices(this%index_list(s1), space%tp_spaces(s1), & | ||
| 157 | include_x_left=.not. coord_has_polar_xy_singularity_ .and. s1 /= 1, include_x_right=s1 /= 1, & | ||
| 158 | include_y_left=s1 /= 2, include_y_right=s1 /= 2, & | ||
| 159 | 160 | include_z_left=s1 /= 3, include_z_right=s1 /= 3) | |
| 160 | end do | ||
| 161 | case (2) | ||
| 162 | nr_blocks = 3 | ||
| 163 | ! Boundary condition is applied to the normal component only | ||
| 164 |
4/6✗ Branch 40 → 41 not taken.
✓ Branch 40 → 42 taken 12 times.
✗ Branch 42 → 43 not taken.
✓ Branch 42 → 44 taken 12 times.
✓ Branch 45 → 46 taken 36 times.
✓ Branch 45 → 47 taken 12 times.
|
48 | allocate (this%index_list(nr_blocks)) |
| 165 |
2/2✓ Branch 48 → 49 taken 12 times.
✓ Branch 48 → 50 taken 36 times.
|
48 | do s1 = 1, 3 |
| 166 | call init_tp_boundary_indices(this%index_list(s1), space%tp_spaces(s1), & | ||
| 167 | include_x_left=.not. coord_has_polar_xy_singularity_ .and. s1 == 1, include_x_right=s1 == 1, & | ||
| 168 | include_y_left=s1 == 2, include_y_right=s1 == 2, & | ||
| 169 | 48 | include_z_left=s1 == 3, include_z_right=s1 == 3) | |
| 170 | end do | ||
| 171 | case (3) | ||
| 172 | nr_blocks = 1 | ||
| 173 | ! No boundary condition | ||
| 174 |
4/6✗ Branch 52 → 53 not taken.
✓ Branch 52 → 54 taken 8 times.
✗ Branch 54 → 55 not taken.
✓ Branch 54 → 56 taken 8 times.
✓ Branch 57 → 58 taken 8 times.
✓ Branch 57 → 59 taken 8 times.
|
16 | allocate (this%index_list(nr_blocks)) |
| 175 | call init_tp_boundary_indices(this%index_list(1), space%tp_spaces(1), & | ||
| 176 | include_x_left=.false., include_x_right=.false., & | ||
| 177 | include_y_left=.false., include_y_right=.false., & | ||
| 178 |
4/5✓ Branch 16 → 17 taken 116 times.
✓ Branch 16 → 27 taken 40 times.
✓ Branch 16 → 40 taken 12 times.
✓ Branch 16 → 52 taken 8 times.
✗ Branch 16 → 60 not taken.
|
184 | include_z_left=.false., include_z_right=.false.) |
| 179 | end select | ||
| 180 | |||
| 181 | ! if (nr_blocks == 1 .and. (present(bdy_function_x) .or. present(bdy_function_y) .or. present(bdy_function_z))) then | ||
| 182 | ! error stop "MFormDirichlet::init_mform_dirichlet: inhomogeneous Dirichlet condition for scalar m-form requires the"// & | ||
| 183 | ! & " optional bdy_function argument" | ||
| 184 | ! end if | ||
| 185 | ! if (nr_blocks == 3 .and. present(bdy_function)) then | ||
| 186 | ! error stop "MFormDirichlet::init_mform_dirichlet: inhomogeneous Dirichlet condition for vector m-form requires the"// & | ||
| 187 | ! & " optional bdy_function_x, bdy_function_y and/or bdy_function_z arguments" | ||
| 188 | ! end if | ||
| 189 | |||
| 190 |
10/16✓ Branch 60 → 61 taken 176 times.
✗ Branch 60 → 62 not taken.
✓ Branch 62 → 63 taken 176 times.
✗ Branch 62 → 64 not taken.
✓ Branch 64 → 65 taken 176 times.
✗ Branch 64 → 66 not taken.
✗ Branch 66 → 67 not taken.
✓ Branch 66 → 68 taken 176 times.
✗ Branch 68 → 69 not taken.
✓ Branch 68 → 70 taken 176 times.
✗ Branch 70 → 71 not taken.
✓ Branch 70 → 72 taken 176 times.
✓ Branch 73 → 74 taken 280 times.
✓ Branch 73 → 78 taken 176 times.
✓ Branch 75 → 76 taken 592 times.
✓ Branch 75 → 77 taken 280 times.
|
1576 | allocate (this%projector(nr_blocks, nr_blocks)) |
| 191 |
5/8✗ Branch 78 → 79 not taken.
✓ Branch 78 → 80 taken 176 times.
✗ Branch 80 → 81 not taken.
✓ Branch 80 → 82 taken 176 times.
✗ Branch 82 → 83 not taken.
✓ Branch 82 → 84 taken 176 times.
✓ Branch 85 → 86 taken 280 times.
✓ Branch 85 → 87 taken 176 times.
|
456 | allocate (this%extractor(nr_blocks)) |
| 192 | |||
| 193 |
2/2✓ Branch 88 → 89 taken 280 times.
✓ Branch 88 → 98 taken 176 times.
|
456 | do s1 = 1, nr_blocks |
| 194 |
2/2✓ Branch 90 → 91 taken 592 times.
✓ Branch 90 → 94 taken 280 times.
|
872 | do s2 = 1, nr_blocks |
| 195 | |||
| 196 | ! The projector is simply the zero-matrix | ||
| 197 | 592 | call this%projector(s1, s2)%init(nr_rows=size(this%index_list(s1)%lin), nr_cols=size(this%index_list(s2)%lin), nr_nonzero=0) | |
| 198 | 872 | call this%projector(s1, s2)%finalize() | |
| 199 | end do | ||
| 200 | |||
| 201 | 280 | call this%extractor(s1)%init(nr_rows=size(this%index_list(s1)%lin), nr_cols=0, nr_nonzero=0) | |
| 202 | 456 | call this%extractor(s1)%finalize() | |
| 203 | end do | ||
| 204 | |||
| 205 | 176 | this%is_symmetric = .true. | |
| 206 | |||
| 207 | TIMER_STOP(t_total) | ||
| 208 | 176 | end subroutine init_mform_dirichlet | |
| 209 | |||
| 210 | !> @brief Initialize a Neumann constraint on an m-form space | ||
| 211 | !> | ||
| 212 | !> @param[inout] this The MFormNeumann object to initialize | ||
| 213 | !> @param[in] space The MFormSpace to which the Neumann constraint applies | ||
| 214 | !> @param[in] _(optional)_ coord_transform The coordinate transform associated with the m-form space | ||
| 215 | 105 | subroutine init_mform_neumann(this, space, coord_transform) | |
| 216 | use m_mform_basis, only: MFormSpace | ||
| 217 | use m_tensorprod_basis, only: size | ||
| 218 | use m_tensorprod, only: TensorProdIndices, intersect | ||
| 219 | use m_coord_transform_abstract, only: CoordTransformAbstract | ||
| 220 | implicit none | ||
| 221 | |||
| 222 | class(MFormNeumann), intent(inout) :: this | ||
| 223 | type(MFormSpace), intent(in) :: space | ||
| 224 | class(CoordTransformAbstract), optional, intent(in) :: coord_transform | ||
| 225 | |||
| 226 | logical :: coord_has_polar_xy_singularity_ | ||
| 227 | |||
| 228 | integer :: nr_blocks, nr_bsplines_x, nr_bsplines_y, nr_bsplines_z, list_size | ||
| 229 | logical :: impose_x_left, impose_x_right | ||
| 230 | integer :: row, col, j, k, rep, nr_per_jk | ||
| 231 | |||
| 232 | type(TensorProdIndices) :: inds, inds_tmp | ||
| 233 | |||
| 234 | TIMER_DECL(t_total) | ||
| 235 | TIMER_INIT_START(t_total, "MFormNeumann::init") | ||
| 236 | |||
| 237 | 105 | call this%destroy() | |
| 238 | |||
| 239 |
4/8✓ Branch 3 → 4 taken 105 times.
✗ Branch 3 → 7 not taken.
✓ Branch 4 → 5 taken 105 times.
✗ Branch 4 → 6 not taken.
✓ Branch 7 → 8 taken 105 times.
✗ Branch 7 → 10 not taken.
✗ Branch 8 → 9 not taken.
✓ Branch 8 → 10 taken 105 times.
|
105 | this%space = space |
| 240 |
2/2✓ Branch 11 → 12 taken 315 times.
✓ Branch 11 → 13 taken 105 times.
|
420 | this%is_mpi_shared = .false. |
| 241 | |||
| 242 |
1/2✗ Branch 13 → 14 not taken.
✓ Branch 13 → 15 taken 105 times.
|
105 | if (this%space%m /= 0) then |
| 243 | ✗ | error stop "MFormNeumann::init_mform_neumann: Neumann constraint is only defined for 0-forms" | |
| 244 | end if | ||
| 245 | nr_blocks = 1 | ||
| 246 | |||
| 247 | ! For simplicity we don't consider domains with corners | ||
| 248 |
1/2✗ Branch 15 → 16 not taken.
✓ Branch 15 → 17 taken 105 times.
|
105 | if (.not. this%space%tp_spaces(1)%spaces(2)%is_periodic) then |
| 249 | ✗ | error stop "MFormNeumann::init_mform_neumann: Neumann constraint is only compatible with periodicity in the y-direction" | |
| 250 | end if | ||
| 251 | |||
| 252 |
1/4✗ Branch 17 → 18 not taken.
✓ Branch 17 → 20 taken 105 times.
✗ Branch 18 → 19 not taken.
✗ Branch 18 → 20 not taken.
|
105 | if (.not. (this%space%dimensionality == 2 .or. .not. this%space%tp_spaces(1)%spaces(3)%is_periodic)) then |
| 253 | error stop "MFormNeumann::init_mform_neumann: 3D Neumann constraint is only compatible with periodicity in the"// & | ||
| 254 | ✗ | & " z-direction" | |
| 255 | end if | ||
| 256 | |||
| 257 | 105 | this%is_homogeneous = .true. | |
| 258 | 105 | this%is_symmetric = .true. | |
| 259 | |||
| 260 | coord_has_polar_xy_singularity_ = .false. | ||
| 261 |
2/4✓ Branch 20 → 21 taken 105 times.
✗ Branch 20 → 23 not taken.
✓ Branch 21 → 22 taken 105 times.
✗ Branch 21 → 23 not taken.
|
105 | if (present(coord_transform)) coord_has_polar_xy_singularity_ = coord_transform%has_polar_xy_singularity |
| 262 | |||
| 263 |
1/2✓ Branch 23 → 24 taken 105 times.
✗ Branch 23 → 25 not taken.
|
105 | if (this%space%tp_spaces(1)%spaces(1)%is_periodic) then |
| 264 | impose_x_left = .false. | ||
| 265 | impose_x_right = .false. | ||
| 266 | else | ||
| 267 | 105 | impose_x_left = .not. coord_has_polar_xy_singularity_ | |
| 268 | impose_x_right = .true. | ||
| 269 | end if | ||
| 270 | |||
| 271 |
4/6✗ Branch 25 → 26 not taken.
✓ Branch 25 → 27 taken 105 times.
✗ Branch 27 → 28 not taken.
✓ Branch 27 → 29 taken 105 times.
✓ Branch 30 → 31 taken 105 times.
✓ Branch 30 → 32 taken 105 times.
|
210 | allocate (this%index_list(nr_blocks)) |
| 272 |
6/8✗ Branch 32 → 33 not taken.
✓ Branch 32 → 34 taken 105 times.
✗ Branch 34 → 35 not taken.
✓ Branch 34 → 36 taken 105 times.
✓ Branch 37 → 38 taken 105 times.
✓ Branch 37 → 41 taken 105 times.
✓ Branch 39 → 37 taken 105 times.
✓ Branch 39 → 40 taken 105 times.
|
315 | allocate (this%projector(nr_blocks, nr_blocks)) |
| 273 |
4/6✗ Branch 41 → 42 not taken.
✓ Branch 41 → 43 taken 105 times.
✗ Branch 43 → 44 not taken.
✓ Branch 43 → 45 taken 105 times.
✓ Branch 46 → 47 taken 105 times.
✓ Branch 46 → 48 taken 105 times.
|
210 | allocate (this%extractor(nr_blocks)) |
| 274 | |||
| 275 |
1/2✗ Branch 48 → 49 not taken.
✓ Branch 48 → 56 taken 105 times.
|
105 | if (.not. impose_x_left .and. .not. impose_x_right) then |
| 276 | |||
| 277 | ! No Neumann constraint to impose | ||
| 278 | ✗ | call inds%empty() | |
| 279 | ✗ | call this%index_list(1)%init(inds, this%space%tp_spaces(1)%rank_resp_bspline, this%space%tp_spaces(1)%rank_l0) | |
| 280 | |||
| 281 | ✗ | call this%extractor(1)%init(nr_rows=0, nr_cols=0, nr_nonzero=0) | |
| 282 | ✗ | call this%extractor(1)%finalize() | |
| 283 | |||
| 284 | ✗ | call this%projector(1, 1)%init(nr_rows=0, nr_cols=0, nr_nonzero=0) | |
| 285 | ✗ | call this%projector(1, 1)%finalize() | |
| 286 | else | ||
| 287 | 105 | nr_bsplines_x = size(this%space%tp_spaces(1), 1) | |
| 288 | 105 | nr_bsplines_y = size(this%space%tp_spaces(1), 2) | |
| 289 | 105 | nr_bsplines_z = size(this%space%tp_spaces(1), 3) | |
| 290 | |||
| 291 | call inds_tmp%init(merge(2, 0, impose_x_left), nr_bsplines_x - 3, 0, nr_bsplines_y - 1, 0, nr_bsplines_z - 1, & | ||
| 292 |
1/2✓ Branch 56 → 57 taken 105 times.
✗ Branch 56 → 58 not taken.
|
210 | this%space%tp_spaces(1)%spaces) |
| 293 | 105 | inds = intersect(inds_tmp, this%space%tp_spaces(1)%rank_resp_bspline) | |
| 294 | 105 | inds%complement = .true. | |
| 295 | |||
| 296 | 105 | call this%index_list(1)%init(inds, this%space%tp_spaces(1)%rank_resp_bspline, this%space%tp_spaces(1)%rank_l0) | |
| 297 | |||
| 298 | 105 | list_size = size(this%index_list(1)%lin) | |
| 299 | 105 | call this%extractor(1)%init(nr_rows=list_size, nr_cols=list_size / 2, nr_nonzero=list_size) | |
| 300 | |||
| 301 | nr_per_jk = 1 + merge(1, 0, impose_x_left) | ||
| 302 | |||
| 303 | 105 | col = 0 | |
| 304 | ! The extractor is a sparse matrix where the i-th column has the (2i - 1)-th and 2i-th row equal to 1 | ||
| 305 |
2/2✓ Branch 62 → 63 taken 3002 times.
✓ Branch 62 → 66 taken 105 times.
|
3107 | do col = 0, list_size / 2 - 1 |
| 306 | 3002 | call this%extractor(1)%insert(2 * col, col, 1.0_wp) | |
| 307 | 3107 | call this%extractor(1)%insert(2 * col + 1, col, 1.0_wp) | |
| 308 | end do | ||
| 309 | 105 | call this%extractor(1)%finalize() | |
| 310 | |||
| 311 | 105 | call this%projector(1, 1)%init(nr_rows=list_size, nr_cols=list_size, nr_nonzero=list_size * 2) | |
| 312 | ! The projector is block diagonal with 2x2 blocks of the form [1 1; 1 1] / sqrt(2) | ||
| 313 |
3/4✗ Branch 69 → 70 not taken.
✓ Branch 69 → 71 taken 105 times.
✓ Branch 75 → 71 taken 2897 times.
✓ Branch 75 → 76 taken 105 times.
|
3107 | do row = 0, list_size - 1, 2 |
| 314 | 3002 | call this%projector(1, 1)%insert(row, row, 1.0_wp / sqrt(2.0_wp)) | |
| 315 | 3002 | call this%projector(1, 1)%insert(row, row + 1, 1.0_wp / sqrt(2.0_wp)) | |
| 316 | 3002 | call this%projector(1, 1)%insert(row + 1, row, 1.0_wp / sqrt(2.0_wp)) | |
| 317 | 3107 | call this%projector(1, 1)%insert(row + 1, row + 1, 1.0_wp / sqrt(2.0_wp)) | |
| 318 | end do | ||
| 319 | |||
| 320 | 105 | call this%projector(1, 1)%finalize() | |
| 321 | end if | ||
| 322 | |||
| 323 | TIMER_STOP(t_total) | ||
| 324 | 105 | end subroutine init_mform_neumann | |
| 325 | |||
| 326 | !> @brief Initialize a uniqueness constraint on an m-form space by fixing one global dof to zero. | ||
| 327 | !> | ||
| 328 | !> The constrained dof is the coefficient at (i,j,k)=(0,0,0), applied to only one m-form component. | ||
| 329 | !> | ||
| 330 | !> @param[inout] this The MFormUniqueness object to initialize | ||
| 331 | !> @param[in] space The MFormSpace to which the uniqueness constraint applies | ||
| 332 | !> @param[in] _(optional)_ coord_transform The coordinate transform associated with the m-form space (unused) | ||
| 333 | 20 | subroutine init_mform_uniqueness(this, space, coord_transform) | |
| 334 | use m_mform_basis, only: MFormSpace | ||
| 335 | use m_coord_transform_abstract, only: CoordTransformAbstract | ||
| 336 | use m_tensorprod_basis, only: size | ||
| 337 | implicit none | ||
| 338 | |||
| 339 | class(MFormUniqueness), intent(inout) :: this | ||
| 340 | type(MFormSpace), intent(in) :: space | ||
| 341 | class(CoordTransformAbstract), optional, intent(in) :: coord_transform | ||
| 342 | |||
| 343 | integer :: nr_blocks, blk_row, blk_col, constrained_blk | ||
| 344 | logical :: owns_constrained_dof | ||
| 345 | logical :: apply_uniqueness | ||
| 346 | logical :: has_y, has_z | ||
| 347 | type(TensorProdIndices) :: inds | ||
| 348 | |||
| 349 | TIMER_DECL(t_total) | ||
| 350 | TIMER_INIT_START(t_total, "MFormUniqueness::init") | ||
| 351 | |||
| 352 | 20 | call this%destroy() | |
| 353 | |||
| 354 |
4/8✓ Branch 3 → 4 taken 20 times.
✗ Branch 3 → 7 not taken.
✓ Branch 4 → 5 taken 20 times.
✗ Branch 4 → 6 not taken.
✓ Branch 7 → 8 taken 20 times.
✗ Branch 7 → 10 not taken.
✗ Branch 8 → 9 not taken.
✓ Branch 8 → 10 taken 20 times.
|
20 | this%space = space |
| 355 | 20 | this%is_homogeneous = .true. | |
| 356 | 20 | this%is_symmetric = .true. | |
| 357 |
2/2✓ Branch 11 → 12 taken 60 times.
✓ Branch 11 → 13 taken 20 times.
|
80 | this%is_mpi_shared = .false. |
| 358 | |||
| 359 | 20 | nr_blocks = size(space%tp_spaces) | |
| 360 | constrained_blk = 1 | ||
| 361 | |||
| 362 | 20 | has_y = size(space%tp_spaces(1), 2) > 1 | |
| 363 | 20 | has_z = size(space%tp_spaces(1), 3) > 1 | |
| 364 | apply_uniqueness = space%tp_spaces(1)%spaces(1)%is_periodic .and. & | ||
| 365 | (.not. has_y .or. space%tp_spaces(1)%spaces(2)%is_periodic) .and. & | ||
| 366 |
3/10✓ Branch 13 → 14 taken 20 times.
✗ Branch 13 → 19 not taken.
✗ Branch 14 → 15 not taken.
✓ Branch 14 → 16 taken 20 times.
✗ Branch 15 → 16 not taken.
✗ Branch 15 → 19 not taken.
✗ Branch 16 → 17 not taken.
✓ Branch 16 → 18 taken 20 times.
✗ Branch 17 → 18 not taken.
✗ Branch 17 → 19 not taken.
|
20 | (.not. has_z .or. space%tp_spaces(1)%spaces(3)%is_periodic) |
| 367 | |||
| 368 |
9/16✗ Branch 19 → 20 not taken.
✓ Branch 19 → 21 taken 20 times.
✓ Branch 21 → 20 taken 20 times.
✗ Branch 21 → 22 not taken.
✓ Branch 22 → 23 taken 20 times.
✗ Branch 22 → 24 not taken.
✓ Branch 24 → 25 taken 20 times.
✗ Branch 24 → 26 not taken.
✗ Branch 26 → 27 not taken.
✓ Branch 26 → 28 taken 20 times.
✗ Branch 28 → 29 not taken.
✓ Branch 28 → 30 taken 20 times.
✗ Branch 30 → 31 not taken.
✓ Branch 30 → 32 taken 20 times.
✓ Branch 33 → 34 taken 20 times.
✓ Branch 33 → 35 taken 20 times.
|
80 | allocate (this%index_list(nr_blocks)) |
| 369 |
13/22✗ Branch 35 → 36 not taken.
✓ Branch 35 → 37 taken 20 times.
✓ Branch 37 → 36 taken 20 times.
✗ Branch 37 → 38 not taken.
✗ Branch 38 → 39 not taken.
✓ Branch 38 → 40 taken 20 times.
✓ Branch 40 → 39 taken 20 times.
✗ Branch 40 → 41 not taken.
✓ Branch 41 → 42 taken 20 times.
✗ Branch 41 → 43 not taken.
✓ Branch 43 → 44 taken 20 times.
✗ Branch 43 → 45 not taken.
✗ Branch 45 → 46 not taken.
✓ Branch 45 → 47 taken 20 times.
✗ Branch 47 → 48 not taken.
✓ Branch 47 → 49 taken 20 times.
✗ Branch 49 → 50 not taken.
✓ Branch 49 → 51 taken 20 times.
✓ Branch 52 → 53 taken 20 times.
✓ Branch 52 → 57 taken 20 times.
✓ Branch 54 → 55 taken 20 times.
✓ Branch 54 → 56 taken 20 times.
|
120 | allocate (this%projector(nr_blocks, nr_blocks)) |
| 370 |
9/16✗ Branch 57 → 58 not taken.
✓ Branch 57 → 59 taken 20 times.
✓ Branch 59 → 58 taken 20 times.
✗ Branch 59 → 60 not taken.
✓ Branch 60 → 61 taken 20 times.
✗ Branch 60 → 62 not taken.
✓ Branch 62 → 63 taken 20 times.
✗ Branch 62 → 64 not taken.
✗ Branch 64 → 65 not taken.
✓ Branch 64 → 66 taken 20 times.
✗ Branch 66 → 67 not taken.
✓ Branch 66 → 68 taken 20 times.
✗ Branch 68 → 69 not taken.
✓ Branch 68 → 70 taken 20 times.
✓ Branch 71 → 72 taken 20 times.
✓ Branch 71 → 73 taken 20 times.
|
80 | allocate (this%extractor(nr_blocks)) |
| 371 | |||
| 372 |
2/2✓ Branch 74 → 75 taken 20 times.
✓ Branch 74 → 98 taken 20 times.
|
40 | do blk_row = 1, nr_blocks |
| 373 |
1/2✓ Branch 75 → 76 taken 20 times.
✗ Branch 75 → 86 not taken.
|
20 | if (apply_uniqueness .and. blk_row == constrained_blk) then |
| 374 | owns_constrained_dof = all([ & | ||
| 375 | space%tp_spaces(blk_row)%rank_resp_bspline%i0 <= 0, 0 <= space%tp_spaces(blk_row)%rank_resp_bspline%i1, & | ||
| 376 | space%tp_spaces(blk_row)%rank_resp_bspline%j0 <= 0, 0 <= space%tp_spaces(blk_row)%rank_resp_bspline%j1, & | ||
| 377 |
6/12✓ Branch 76 → 77 taken 20 times.
✗ Branch 76 → 84 not taken.
✓ Branch 77 → 78 taken 20 times.
✗ Branch 77 → 84 not taken.
✓ Branch 78 → 79 taken 20 times.
✗ Branch 78 → 84 not taken.
✓ Branch 79 → 80 taken 20 times.
✗ Branch 79 → 84 not taken.
✓ Branch 80 → 81 taken 20 times.
✗ Branch 80 → 84 not taken.
✓ Branch 81 → 82 taken 20 times.
✗ Branch 81 → 84 not taken.
|
20 | space%tp_spaces(blk_row)%rank_resp_bspline%k0 <= 0, 0 <= space%tp_spaces(blk_row)%rank_resp_bspline%k1 ]) |
| 378 | |||
| 379 | if (owns_constrained_dof) then | ||
| 380 | 20 | call inds%init(0, 0, 0, 0, 0, 0, space%tp_spaces(blk_row)%spaces) | |
| 381 | else | ||
| 382 | ✗ | call inds%empty() | |
| 383 | end if | ||
| 384 | else | ||
| 385 | ✗ | call inds%empty() | |
| 386 | end if | ||
| 387 | |||
| 388 | 20 | call this%index_list(blk_row)%init(inds, space%tp_spaces(blk_row)%rank_resp_bspline, space%tp_spaces(blk_row)%rank_l0) | |
| 389 | |||
| 390 | 20 | call this%extractor(blk_row)%init(nr_rows=size(this%index_list(blk_row)%lin), nr_cols=0, nr_nonzero=0) | |
| 391 | 20 | call this%extractor(blk_row)%finalize() | |
| 392 | |||
| 393 |
2/2✓ Branch 92 → 93 taken 20 times.
✓ Branch 92 → 96 taken 20 times.
|
60 | do blk_col = 1, nr_blocks |
| 394 | call this%projector(blk_row, blk_col)%init(nr_rows=size(this%index_list(blk_row)%lin), & | ||
| 395 | 20 | nr_cols=size(this%index_list(blk_col)%lin), nr_nonzero=0) | |
| 396 | 40 | call this%projector(blk_row, blk_col)%finalize() | |
| 397 | end do | ||
| 398 | end do | ||
| 399 | |||
| 400 | TIMER_STOP(t_total) | ||
| 401 | 20 | end subroutine init_mform_uniqueness | |
| 402 | |||
| 403 | !> @brief Initialize my interior indices (i.e., all indices on my rank minus the boundary) of the B-spline tensor product space | ||
| 404 | !> | ||
| 405 | !> @param[in] tp_bspline The B-spline tensor product space | ||
| 406 | !> @param[in] include_x_left Whether to have a boundary on the left-hand side in the x-direction (default: .true.) | ||
| 407 | !> @param[in] include_x_right Whether to have a boundary on the right-hand side in the x-direction (default: .true.) | ||
| 408 | !> @param[in] include_y_left Whether to have a boundary on the left-hand side in the y-direction (default: .true.) | ||
| 409 | !> @param[in] include_y_right Whether to have a boundary on the right-hand side in the y-direction (default: .true.) | ||
| 410 | !> @param[in] include_z_left Whether to have a boundary on the left-hand side in the z-direction (default: .true.) | ||
| 411 | !> @param[in] include_z_right Whether to have a boundary on the right-hand side in the z-direction (default: .true.) | ||
| 412 | 280 | pure type(TensorProdIndices) function my_interior_index_range(tp_bspline, include_x_left, include_x_right, & | |
| 413 | & include_y_left, include_y_right, include_z_left, include_z_right) result(inds) | ||
| 414 | |||
| 415 | use m_tensorprod, only: TensorProdSpace, size, TensorProdIndices, intersect | ||
| 416 | implicit none | ||
| 417 | |||
| 418 | type(TensorProdSpace), intent(in) :: tp_bspline | ||
| 419 | logical, optional, intent(in) :: include_x_left, include_x_right, & | ||
| 420 | include_y_left, include_y_right, & | ||
| 421 | include_z_left, include_z_right | ||
| 422 | |||
| 423 | logical :: include_x_left_, include_x_right_, & | ||
| 424 | include_y_left_, include_y_right_, & | ||
| 425 | include_z_left_, include_z_right_ | ||
| 426 | integer :: i0, i1, j0, j1, k0, k1 | ||
| 427 | type(TensorProdIndices) :: inds_tmp | ||
| 428 | |||
| 429 | ! By defualt all boundaries are included (unless is_periodic) | ||
| 430 | include_x_left_ = .true. | ||
| 431 |
1/2✓ Branch 2 → 3 taken 280 times.
✗ Branch 2 → 4 not taken.
|
280 | if (present(include_x_left)) include_x_left_ = include_x_left |
| 432 | include_x_right_ = .true. | ||
| 433 |
2/2✓ Branch 4 → 5 taken 164 times.
✓ Branch 4 → 6 taken 116 times.
|
280 | if (present(include_x_right)) include_x_right_ = include_x_right |
| 434 | include_y_left_ = .true. | ||
| 435 |
2/2✓ Branch 6 → 7 taken 164 times.
✓ Branch 6 → 8 taken 116 times.
|
280 | if (present(include_y_left)) include_y_left_ = include_y_left |
| 436 | include_y_right_ = .true. | ||
| 437 |
2/2✓ Branch 8 → 9 taken 164 times.
✓ Branch 8 → 10 taken 116 times.
|
280 | if (present(include_y_right)) include_y_right_ = include_y_right |
| 438 | include_z_left_ = .true. | ||
| 439 |
2/2✓ Branch 10 → 11 taken 164 times.
✓ Branch 10 → 12 taken 116 times.
|
280 | if (present(include_z_left)) include_z_left_ = include_z_left |
| 440 | include_z_right_ = .true. | ||
| 441 |
2/2✓ Branch 12 → 13 taken 164 times.
✓ Branch 12 → 14 taken 116 times.
|
280 | if (present(include_z_right)) include_z_right_ = include_z_right |
| 442 | |||
| 443 | ! Handle lower dim problem (recognized by less than or equal to 1 basis function in that direction): | ||
| 444 |
2/2✓ Branch 14 → 15 taken 112 times.
✓ Branch 14 → 17 taken 168 times.
|
280 | if (size(tp_bspline, 3) <= 1) then |
| 445 | include_z_left_ = .false. | ||
| 446 | include_z_right_ = .false. | ||
| 447 |
2/2✓ Branch 15 → 16 taken 10 times.
✓ Branch 15 → 17 taken 102 times.
|
112 | if (size(tp_bspline, 2) <= 1) then |
| 448 | include_y_left_ = .false. | ||
| 449 | include_y_right_ = .false. | ||
| 450 | end if | ||
| 451 | end if | ||
| 452 | |||
| 453 |
4/4✓ Branch 17 → 18 taken 86 times.
✓ Branch 17 → 19 taken 194 times.
✓ Branch 19 → 18 taken 128 times.
✓ Branch 19 → 20 taken 66 times.
|
280 | i0 = merge(0, 1, tp_bspline%spaces(1)%is_periodic .or. .not. include_x_left_) |
| 454 |
6/6✓ Branch 20 → 21 taken 194 times.
✓ Branch 20 → 22 taken 86 times.
✓ Branch 21 → 22 taken 42 times.
✓ Branch 21 → 23 taken 152 times.
✓ Branch 23 → 24 taken 152 times.
✓ Branch 23 → 25 taken 128 times.
|
408 | i1 = size(tp_bspline, 1) - 1 - merge(0, 1, tp_bspline%spaces(1)%is_periodic .or. .not. include_x_right_) |
| 455 |
4/4✓ Branch 25 → 26 taken 131 times.
✓ Branch 25 → 27 taken 149 times.
✓ Branch 26 → 27 taken 43 times.
✓ Branch 26 → 28 taken 88 times.
|
280 | j0 = merge(0, 1, tp_bspline%spaces(2)%is_periodic .or. .not. include_y_left_) |
| 456 |
6/6✓ Branch 28 → 29 taken 131 times.
✓ Branch 28 → 30 taken 149 times.
✓ Branch 29 → 30 taken 43 times.
✓ Branch 29 → 31 taken 88 times.
✓ Branch 31 → 32 taken 88 times.
✓ Branch 31 → 33 taken 192 times.
|
472 | j1 = size(tp_bspline, 2) - 1 - merge(0, 1, tp_bspline%spaces(2)%is_periodic .or. .not. include_y_right_) |
| 457 |
4/4✓ Branch 33 → 34 taken 155 times.
✓ Branch 33 → 35 taken 125 times.
✓ Branch 34 → 35 taken 115 times.
✓ Branch 34 → 36 taken 40 times.
|
280 | k0 = merge(0, 1, tp_bspline%spaces(3)%is_periodic .or. .not. include_z_left_) |
| 458 |
6/6✓ Branch 36 → 37 taken 155 times.
✓ Branch 36 → 38 taken 125 times.
✓ Branch 37 → 38 taken 115 times.
✓ Branch 37 → 39 taken 40 times.
✓ Branch 39 → 40 taken 40 times.
✓ Branch 39 → 41 taken 240 times.
|
520 | k1 = size(tp_bspline, 3) - 1 - merge(0, 1, tp_bspline%spaces(3)%is_periodic .or. .not. include_z_right_) |
| 459 | |||
| 460 | 280 | call inds_tmp%init(i0, i1, j0, j1, k0, k1, tp_bspline%spaces) | |
| 461 | 280 | inds = intersect(inds_tmp, tp_bspline%rank_resp_bspline) | |
| 462 | 280 | end function my_interior_index_range | |
| 463 | |||
| 464 | !> @brief Initialize my boundary indices (i.e., the indices of the boundary on the current rank) of the B-spline tensor product | ||
| 465 | !> space | ||
| 466 | !> | ||
| 467 | !> @param[inout] boundary The TensorProdIndexList to initialize | ||
| 468 | !> @param[in] tp_bspline The B-spline tensor product space | ||
| 469 | !> @param[in] include_x_left Whether to have a boundary on the left-hand side in the x-direction | ||
| 470 | !> @param[in] include_x_right Whether to have a boundary on the right-hand side in the x-direction | ||
| 471 | !> @param[in] include_y_left Whether to have a boundary on the left-hand side in the y-direction | ||
| 472 | !> @param[in] include_y_right Whether to have a boundary on the right-hand side in the y-direction | ||
| 473 | !> @param[in] include_z_left Whether to have a boundary on the left-hand side in the z-direction | ||
| 474 | !> @param[in] include_z_right Whether to have a boundary on the right-hand side in the z-direction | ||
| 475 | 280 | subroutine init_tp_boundary_indices(boundary, tp_bspline, include_x_left, include_x_right, & | |
| 476 | include_y_left, include_y_right, include_z_left, include_z_right) | ||
| 477 | |||
| 478 | use m_tensorprod, only: TensorProdSpace, TensorProdIndices | ||
| 479 | implicit none | ||
| 480 | |||
| 481 | type(TensorProdIndexList), intent(inout) :: boundary | ||
| 482 | type(TensorProdSpace), intent(in) :: tp_bspline | ||
| 483 | logical, optional, intent(in) :: include_x_left, include_x_right, & | ||
| 484 | include_y_left, include_y_right, & | ||
| 485 | include_z_left, include_z_right | ||
| 486 | |||
| 487 | ! Start by determining the interior indices | ||
| 488 | type(TensorProdIndices) :: inds | ||
| 489 | |||
| 490 | inds = my_interior_index_range(tp_bspline, & | ||
| 491 | include_x_left=include_x_left, include_x_right=include_x_right, & | ||
| 492 | include_y_left=include_y_left, include_y_right=include_y_right, & | ||
| 493 | 280 | include_z_left=include_z_left, include_z_right=include_z_right) | |
| 494 | |||
| 495 | 280 | inds%complement = .true. ! NOTE: complement is relative to local B-spline indices of this rank | |
| 496 | |||
| 497 | 280 | call boundary%init(inds, tp_bspline%rank_resp_bspline, tp_bspline%rank_l0) | |
| 498 | 280 | end subroutine init_tp_boundary_indices | |
| 499 | |||
| 500 | !> @brief Destroy the MFormDirichlet object | ||
| 501 | !> | ||
| 502 | !> @param[inout] this The MFormDirichlet object to destroy | ||
| 503 | 308 | subroutine destroy_mform_dirichlet(this) | |
| 504 | implicit none | ||
| 505 | class(MFormDirichlet), intent(inout) :: this | ||
| 506 | |||
| 507 | 308 | end subroutine destroy_mform_dirichlet | |
| 508 | |||
| 509 | !> @brief Destroy the MFormNeumann object | ||
| 510 | !> | ||
| 511 | !> @param[inout] this The MFormNeumann object to destroy | ||
| 512 | 210 | subroutine destroy_mform_neumann(this) | |
| 513 | implicit none | ||
| 514 | class(MFormNeumann), intent(inout) :: this | ||
| 515 | |||
| 516 | 210 | end subroutine destroy_mform_neumann | |
| 517 | |||
| 518 | !> @brief Destroy the MFormUniqueness object | ||
| 519 | !> | ||
| 520 | !> @param[inout] this The MFormUniqueness object to destroy | ||
| 521 | 40 | subroutine destroy_mform_uniqueness(this) | |
| 522 | implicit none | ||
| 523 | class(MFormUniqueness), intent(inout) :: this | ||
| 524 | |||
| 525 | 40 | end subroutine destroy_mform_uniqueness | |
| 526 |
51/264__m_mform_constraint_boundary_MOD___copy_m_mform_constraint_boundary_Mformdirichlet:
✓ Branch 2 → 3 taken 132 times.
✗ Branch 2 → 51 not taken.
✗ Branch 3 → 4 not taken.
✓ Branch 3 → 5 taken 132 times.
✗ Branch 6 → 7 not taken.
✓ Branch 6 → 22 taken 132 times.
✗ Branch 8 → 9 not taken.
✗ Branch 8 → 23 not taken.
✗ Branch 9 → 10 not taken.
✗ Branch 9 → 11 not taken.
✗ Branch 12 → 13 not taken.
✗ Branch 12 → 14 not taken.
✗ Branch 15 → 16 not taken.
✗ Branch 15 → 17 not taken.
✗ Branch 18 → 19 not taken.
✗ Branch 18 → 20 not taken.
✗ Branch 23 → 24 not taken.
✓ Branch 23 → 36 taken 132 times.
✗ Branch 25 → 26 not taken.
✗ Branch 25 → 37 not taken.
✗ Branch 26 → 27 not taken.
✗ Branch 26 → 28 not taken.
✗ Branch 29 → 30 not taken.
✗ Branch 29 → 31 not taken.
✗ Branch 32 → 33 not taken.
✗ Branch 32 → 34 not taken.
✗ Branch 37 → 38 not taken.
✓ Branch 37 → 50 taken 132 times.
✗ Branch 39 → 40 not taken.
✗ Branch 39 → 51 not taken.
✗ Branch 40 → 41 not taken.
✗ Branch 40 → 42 not taken.
✗ Branch 43 → 44 not taken.
✗ Branch 43 → 45 not taken.
✗ Branch 46 → 47 not taken.
✗ Branch 46 → 48 not taken.
__m_mform_constraint_boundary_MOD___copy_m_mform_constraint_boundary_Mformneumann:
✓ Branch 2 → 3 taken 105 times.
✗ Branch 2 → 51 not taken.
✗ Branch 3 → 4 not taken.
✓ Branch 3 → 5 taken 105 times.
✗ Branch 6 → 7 not taken.
✓ Branch 6 → 22 taken 105 times.
✗ Branch 8 → 9 not taken.
✗ Branch 8 → 23 not taken.
✗ Branch 9 → 10 not taken.
✗ Branch 9 → 11 not taken.
✗ Branch 12 → 13 not taken.
✗ Branch 12 → 14 not taken.
✗ Branch 15 → 16 not taken.
✗ Branch 15 → 17 not taken.
✗ Branch 18 → 19 not taken.
✗ Branch 18 → 20 not taken.
✗ Branch 23 → 24 not taken.
✓ Branch 23 → 36 taken 105 times.
✗ Branch 25 → 26 not taken.
✗ Branch 25 → 37 not taken.
✗ Branch 26 → 27 not taken.
✗ Branch 26 → 28 not taken.
✗ Branch 29 → 30 not taken.
✗ Branch 29 → 31 not taken.
✗ Branch 32 → 33 not taken.
✗ Branch 32 → 34 not taken.
✗ Branch 37 → 38 not taken.
✓ Branch 37 → 50 taken 105 times.
✗ Branch 39 → 40 not taken.
✗ Branch 39 → 51 not taken.
✗ Branch 40 → 41 not taken.
✗ Branch 40 → 42 not taken.
✗ Branch 43 → 44 not taken.
✗ Branch 43 → 45 not taken.
✗ Branch 46 → 47 not taken.
✗ Branch 46 → 48 not taken.
__m_mform_constraint_boundary_MOD___copy_m_mform_constraint_boundary_Mformuniqueness:
✓ Branch 2 → 3 taken 20 times.
✗ Branch 2 → 51 not taken.
✗ Branch 3 → 4 not taken.
✓ Branch 3 → 5 taken 20 times.
✗ Branch 6 → 7 not taken.
✓ Branch 6 → 22 taken 20 times.
✗ Branch 8 → 9 not taken.
✗ Branch 8 → 23 not taken.
✗ Branch 9 → 10 not taken.
✗ Branch 9 → 11 not taken.
✗ Branch 12 → 13 not taken.
✗ Branch 12 → 14 not taken.
✗ Branch 15 → 16 not taken.
✗ Branch 15 → 17 not taken.
✗ Branch 18 → 19 not taken.
✗ Branch 18 → 20 not taken.
✗ Branch 23 → 24 not taken.
✓ Branch 23 → 36 taken 20 times.
✗ Branch 25 → 26 not taken.
✗ Branch 25 → 37 not taken.
✗ Branch 26 → 27 not taken.
✗ Branch 26 → 28 not taken.
✗ Branch 29 → 30 not taken.
✗ Branch 29 → 31 not taken.
✗ Branch 32 → 33 not taken.
✗ Branch 32 → 34 not taken.
✗ Branch 37 → 38 not taken.
✓ Branch 37 → 50 taken 20 times.
✗ Branch 39 → 40 not taken.
✗ Branch 39 → 51 not taken.
✗ Branch 40 → 41 not taken.
✗ Branch 40 → 42 not taken.
✗ Branch 43 → 44 not taken.
✗ Branch 43 → 45 not taken.
✗ Branch 46 → 47 not taken.
✗ Branch 46 → 48 not taken.
__m_mform_constraint_boundary_MOD___final_m_mform_constraint_boundary_Mformdirichlet:
✗ Branch 3 → 4 not taken.
✓ Branch 3 → 7 taken 132 times.
✗ Branch 4 → 5 not taken.
✗ Branch 4 → 6 not taken.
✓ Branch 8 → 9 taken 132 times.
✓ Branch 8 → 56 taken 132 times.
✗ Branch 10 → 11 not taken.
✓ Branch 10 → 12 taken 132 times.
✓ Branch 12 → 13 taken 132 times.
✗ Branch 12 → 15 not taken.
✓ Branch 13 → 14 taken 132 times.
✗ Branch 13 → 16 not taken.
✗ Branch 15 → 16 not taken.
✗ Branch 15 → 55 not taken.
✗ Branch 16 → 17 not taken.
✓ Branch 16 → 28 taken 132 times.
✗ Branch 18 → 19 not taken.
✗ Branch 18 → 28 not taken.
✗ Branch 19 → 20 not taken.
✗ Branch 19 → 21 not taken.
✗ Branch 21 → 22 not taken.
✗ Branch 21 → 23 not taken.
✗ Branch 23 → 24 not taken.
✗ Branch 23 → 25 not taken.
✗ Branch 25 → 26 not taken.
✗ Branch 25 → 27 not taken.
✗ Branch 28 → 29 not taken.
✓ Branch 28 → 30 taken 132 times.
✗ Branch 30 → 31 not taken.
✓ Branch 30 → 40 taken 132 times.
✗ Branch 32 → 33 not taken.
✗ Branch 32 → 40 not taken.
✗ Branch 33 → 34 not taken.
✗ Branch 33 → 35 not taken.
✗ Branch 35 → 36 not taken.
✗ Branch 35 → 37 not taken.
✗ Branch 37 → 38 not taken.
✗ Branch 37 → 39 not taken.
✗ Branch 40 → 41 not taken.
✓ Branch 40 → 42 taken 132 times.
✗ Branch 42 → 43 not taken.
✓ Branch 42 → 52 taken 132 times.
✗ Branch 44 → 45 not taken.
✗ Branch 44 → 52 not taken.
✗ Branch 45 → 46 not taken.
✗ Branch 45 → 47 not taken.
✗ Branch 47 → 48 not taken.
✗ Branch 47 → 49 not taken.
✗ Branch 49 → 50 not taken.
✗ Branch 49 → 51 not taken.
✗ Branch 52 → 53 not taken.
✓ Branch 52 → 54 taken 132 times.
__m_mform_constraint_boundary_MOD___final_m_mform_constraint_boundary_Mformneumann:
✗ Branch 3 → 4 not taken.
✓ Branch 3 → 7 taken 105 times.
✗ Branch 4 → 5 not taken.
✗ Branch 4 → 6 not taken.
✓ Branch 8 → 9 taken 105 times.
✓ Branch 8 → 56 taken 105 times.
✗ Branch 10 → 11 not taken.
✓ Branch 10 → 12 taken 105 times.
✓ Branch 12 → 13 taken 105 times.
✗ Branch 12 → 15 not taken.
✓ Branch 13 → 14 taken 105 times.
✗ Branch 13 → 16 not taken.
✗ Branch 15 → 16 not taken.
✗ Branch 15 → 55 not taken.
✗ Branch 16 → 17 not taken.
✓ Branch 16 → 28 taken 105 times.
✗ Branch 18 → 19 not taken.
✗ Branch 18 → 28 not taken.
✗ Branch 19 → 20 not taken.
✗ Branch 19 → 21 not taken.
✗ Branch 21 → 22 not taken.
✗ Branch 21 → 23 not taken.
✗ Branch 23 → 24 not taken.
✗ Branch 23 → 25 not taken.
✗ Branch 25 → 26 not taken.
✗ Branch 25 → 27 not taken.
✗ Branch 28 → 29 not taken.
✓ Branch 28 → 30 taken 105 times.
✗ Branch 30 → 31 not taken.
✓ Branch 30 → 40 taken 105 times.
✗ Branch 32 → 33 not taken.
✗ Branch 32 → 40 not taken.
✗ Branch 33 → 34 not taken.
✗ Branch 33 → 35 not taken.
✗ Branch 35 → 36 not taken.
✗ Branch 35 → 37 not taken.
✗ Branch 37 → 38 not taken.
✗ Branch 37 → 39 not taken.
✗ Branch 40 → 41 not taken.
✓ Branch 40 → 42 taken 105 times.
✗ Branch 42 → 43 not taken.
✓ Branch 42 → 52 taken 105 times.
✗ Branch 44 → 45 not taken.
✗ Branch 44 → 52 not taken.
✗ Branch 45 → 46 not taken.
✗ Branch 45 → 47 not taken.
✗ Branch 47 → 48 not taken.
✗ Branch 47 → 49 not taken.
✗ Branch 49 → 50 not taken.
✗ Branch 49 → 51 not taken.
✗ Branch 52 → 53 not taken.
✓ Branch 52 → 54 taken 105 times.
__m_mform_constraint_boundary_MOD___final_m_mform_constraint_boundary_Mformuniqueness:
✗ Branch 3 → 4 not taken.
✓ Branch 3 → 7 taken 20 times.
✗ Branch 4 → 5 not taken.
✗ Branch 4 → 6 not taken.
✓ Branch 8 → 9 taken 20 times.
✓ Branch 8 → 56 taken 20 times.
✗ Branch 10 → 11 not taken.
✓ Branch 10 → 12 taken 20 times.
✓ Branch 12 → 13 taken 20 times.
✗ Branch 12 → 15 not taken.
✓ Branch 13 → 14 taken 20 times.
✗ Branch 13 → 16 not taken.
✗ Branch 15 → 16 not taken.
✗ Branch 15 → 55 not taken.
✗ Branch 16 → 17 not taken.
✓ Branch 16 → 28 taken 20 times.
✗ Branch 18 → 19 not taken.
✗ Branch 18 → 28 not taken.
✗ Branch 19 → 20 not taken.
✗ Branch 19 → 21 not taken.
✗ Branch 21 → 22 not taken.
✗ Branch 21 → 23 not taken.
✗ Branch 23 → 24 not taken.
✗ Branch 23 → 25 not taken.
✗ Branch 25 → 26 not taken.
✗ Branch 25 → 27 not taken.
✗ Branch 28 → 29 not taken.
✓ Branch 28 → 30 taken 20 times.
✗ Branch 30 → 31 not taken.
✓ Branch 30 → 40 taken 20 times.
✗ Branch 32 → 33 not taken.
✗ Branch 32 → 40 not taken.
✗ Branch 33 → 34 not taken.
✗ Branch 33 → 35 not taken.
✗ Branch 35 → 36 not taken.
✗ Branch 35 → 37 not taken.
✗ Branch 37 → 38 not taken.
✗ Branch 37 → 39 not taken.
✗ Branch 40 → 41 not taken.
✓ Branch 40 → 42 taken 20 times.
✗ Branch 42 → 43 not taken.
✓ Branch 42 → 52 taken 20 times.
✗ Branch 44 → 45 not taken.
✗ Branch 44 → 52 not taken.
✗ Branch 45 → 46 not taken.
✗ Branch 45 → 47 not taken.
✗ Branch 47 → 48 not taken.
✗ Branch 47 → 49 not taken.
✗ Branch 49 → 50 not taken.
✗ Branch 49 → 51 not taken.
✗ Branch 52 → 53 not taken.
✓ Branch 52 → 54 taken 20 times.
|
771 | end module m_mform_constraint_boundary |
| 527 |