Alexander Fuchs alex-fu27
alex-fu27
/ ErrorMWE.jl
Last active
March 4, 2025 08:08
Periodic boundary conditions cant be applied because assembly changes sparsity pattern?
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| using Ferrite | |
| function assemble_element!(Me :: Matrix, Ke::Matrix, cellvalues::CellValues, q :: Ferrite.Vec{2}) | |
| n_basefuncs = getnbasefunctions(cellvalues) | |
| # Reset to 0 | |
| fill!(Me, 0) | |
| fill!(Ke, 0) | |
| # Loop over quadrature points | |
| for q_point in 1:getnquadpoints(cellvalues) |