lijas
lijas
Creating the sparsity pattern with affine constraints was extremely slow. ``` grid = generate_grid(Quadrilateral, (300,300)) dh = DofHandler(grid) push!(dh, :u, 2) close!(dh); ch_periodic = ConstraintHandler(dh); periodic = PeriodicDirichlet( :u, ["left"...
I think it is pretty common to want to extract the dofs on some face in order to calculate for example the reaction forces on a boundary. This PR adds...
Also noticed that edges was wrong order compared to VTK, so changed order of that. https://vtk.org/doc/nightly/html/classvtkQuadraticHexahedron.html Also, the order of the faces is wrong, which will affect vtk output if...
I want to create an continuum-shell element with n number of dofs on edges in the thickness direction, and zero edge-dofs inplane. Currently the DofHandler assumes the same number of...
An example with large deformation plasticity. I basically copied the hyperelasticity problem but added an elasto-plastic material model. Don't know if this example is needed since there already exists an...
I have a linear system in the form: A = [K C'; C, 0], (C is relatively small) which requires a large amount of iterations with minres/gmres in order to...
It would be nice to be able to give solver parameters as arguments when defining the solver type. Currently I think one has to give the solver parameters as keywordargument...
I had some plans of using block arrays with iterative solvers, but it seems like matrix-vector multiplications with block arrays is very slow. ``` using BlockArrays, SparseArrays N = 10000...
@kimauth mentioned that the current way of handling shells in the DofHandler interfears with cohesive zone elements. This PR adds a wrapper interpolation for shells that does not mess with...