Solve the equation system with just the penalization terms to find an implicit update for the IMEX timestepping procedure.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| type(atl_material_type), | intent(in) | :: | material |
Definition of the material, which directly describes the penalization. We expect the conductivity sigma to be defined in materialdat(:,:,3) and the permittivity epsilon in materialdat(:,:,2). |
||
| real(kind=rk), | intent(in) | :: | weighted_dt |
Timestep which is already weighted by the time integration scheme. |
||
| integer, | intent(in) | :: | ndims |
Number of dimensions, the equation system is computed in (2 or 3). |
||
| type(ply_poly_project_type), | intent(inout) | :: | poly_proj |
Description of the projection for the material. |
||
| real(kind=rk), | intent(inout) | :: | state(:,:,:) |
The state variables of the equation system, they will be updated to the solution of the implicit computation for penalization. |
||
| real(kind=rk), | intent(out) | :: | timestep_rhs(:,:,:) |
Right hand side contribution by the implicit calculation. |