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