@@ -16,7 +16,7 @@ function get_periodic_coupling_info(
1616 b2,
1717 is_opposite:: Function ;
1818 factor_vectordofs = " auto" ,
19- factor_components = " auto" ,
19+ factor_components = " auto"
2020 )
2121
2222 FEType = eltype (FES)
@@ -198,15 +198,15 @@ function interpolate_on_boundaryfaces(
198198end
199199
200200"""
201- get_periodic_coupling_matrix(
202- FES::FESpace,
203- xgrid::ExtendableGrid,
204- b_from,
205- b_to,
206- give_opposite!::Function;
207- mask = :auto,
208- sparsity_tol = 1.0e-12
209- )
201+ get_periodic_coupling_matrix(
202+ FES::FESpace,
203+ xgrid::ExtendableGrid,
204+ b_from,
205+ b_to,
206+ give_opposite!::Function;
207+ mask = :auto,
208+ sparsity_tol = 1.0e-12
209+ )
210210
211211Compute a coupling information for each dof on one boundary as a linear combination of dofs on another boundary
212212
@@ -226,8 +226,8 @@ Example: If b_from is at x[1] = 0 and the opposite boundary is at y[1] = 1, then
226226The return value is a (𝑛 × 𝑛) sparse matrix 𝐴 (𝑛 is the total number of dofs) containing the periodic coupling information.
227227The relation ship between the degrees of freedome is dofᵢ = ∑ⱼ Aⱼᵢ ⋅ dofⱼ.
228228It is guaranteed that
229- i) Aⱼᵢ=0 if dofᵢ is 𝑛𝑜𝑡 on the boundary b_from.
230- ii) Aⱼᵢ=0 if the opposite of dofᵢ is not in the same grid cell as dofⱼ.
229+ i) Aⱼᵢ=0 if dofᵢ is 𝑛𝑜𝑡 on the boundary b_from.
230+ ii) Aⱼᵢ=0 if the opposite of dofᵢ is not in the same grid cell as dofⱼ.
231231Note that A is transposed for efficient col-wise storage.
232232
233233"""
@@ -370,7 +370,9 @@ function get_periodic_coupling_matrix(
370370
371371 # interpolate on the opposite boundary using x_trafo = give_opposite
372372 interpolate! (
373- fe_vector_target[1 ], ON_FACES, eval_point, items = search_areas[face_numbers_of_bfaces[i_boundary_face]],
373+ fe_vector_target[1 ],
374+ ON_FACES, eval_point,
375+ items = search_areas[face_numbers_of_bfaces[i_boundary_face]],
374376 )
375377
376378 # deactivate entry
@@ -512,7 +514,7 @@ function tmul!(
512514 A:: AbstractMatrix{T} ,
513515 x:: AbstractVector{T} ,
514516 α = 1.0 ,
515- β = 0.0 ,
517+ β = 0.0
516518 ) where {T <: AbstractFloat }
517519 return LinearAlgebra. BLAS. gemv! (' T' , α, A, x, β, y)
518520end
0 commit comments