High-order, parallel, matrix-free finite element solver for Stokes and Brinkman flow written in Julia.
Details and benchmarks can be found here.
There are two implementations of the software, a parallel Julia implementation and a serial Matlab implementation. For each implementation, there are three versions that correspond to various finite element assembly methods, using assembled matrices, matrix-free methods using matrix-vector products, and matrix-free methods using matrix-matrix products.
From the top-level directory, cd julia-parallel/efficient-operators
. Start julia
and run include(sss.jl)
to start an interactive session that will prompt you for a type of test problem to solve (lid-driven cavity or brinkman flow through obstacles), mesh size, preconditioning steps, etc.