aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
...
* mg2d: make the refinement depth/exact solve size configurableAnton Khirnov2019-03-11
| | | | API bump
* mg2d: allocate separate data arrays for the multigrid layerAnton Khirnov2019-03-11
| | | | | | | Do not use those from the topmost level directly. This introduces an additional copy, but allows us to decouple level allocation from the solver allocation.
* mg2d: remove unused variableAnton Khirnov2019-03-11
|
* mg2d: refactor levels allocationAnton Khirnov2019-03-04
| | | | This should simplify future improvements.
* mg2d: be less strict for converging to machine epsilonAnton Khirnov2019-03-04
|
* mg2d: fail if the residual grows too much from the original oneAnton Khirnov2019-03-04
|
* mg2d: take into account the factor in front of ∂_xx when computing the CFL ↵Anton Khirnov2019-03-04
| | | | factor
* ell_grid_solve: initialize diff_coeffs to 0Anton Khirnov2019-03-04
| | | | Prevents the uninitialized items from affecting interpolation.
* Implement C/r falloff boundary condition.Anton Khirnov2019-03-04
| | | | API bump.
* ell_grid_solve: do not read from outside of allocated boundary valuesAnton Khirnov2019-02-24
|
* mg2d: print a more explicit error message when we reach maxiterAnton Khirnov2019-02-24
|
* ell_grid_solve: do not read from possibly-uninitialized boundary layerAnton Khirnov2019-02-08
|
* Rename fixdiff boundary condition to reflect.Anton Khirnov2019-02-08
| | | | | | | This is what it actually does, we do not support setting arbitrary derivative values. API bump.
* ell_grid_solve: use the new boundary API to simplify codeAnton Khirnov2019-02-06
|
* mg2d_test: use the new boundary API to simplify codeAnton Khirnov2019-02-06
|
* relax_test: use the new boundary API to simplify codeAnton Khirnov2019-02-06
|
* boundary: add public functions for querying boundary informationAnton Khirnov2019-02-06
| | | | | | Should allow treating boundaries in a generic manner. API bump
* Change the meaning of boundary location to make more sense.Anton Khirnov2019-02-05
| | | | API bump.
* Add forgotten header.Anton Khirnov2019-02-05
|
* meson.build: link-depend on the symbol visibility linker scriptAnton Khirnov2019-02-05
|
* residual_calc.asm: use the correct coefficients for y derivativesAnton Khirnov2019-02-02
|
* mg2d_test: reduce the symmetry of the test problemAnton Khirnov2019-02-02
| | | | Should make for a more robust test.
* mg2d: relax the condition on residual convergenceAnton Khirnov2019-02-02
| | | | | Apparently in certain cases the residual can briefly grow, even though the process converges after.
* Add forgotten boundary.c file.Anton Khirnov2019-02-02
|
* mg2d: print stats properly for exact solvesAnton Khirnov2019-02-02
|
* ell_grid_solve: split residual computation into its own fileAnton Khirnov2019-01-30
|
* mg2d_test: add a test for fixed-derivative boundary conditionAnton Khirnov2019-01-29
| | | | Also test more terms in the equation.
* mg2d: account for effect of the ~u term on the relaxation factorAnton Khirnov2019-01-29
| | | | | The maximum allowed time step in the presence of a -K * u term goes like 2 (dx ** 2) / (4 + K (dx ** 2))
* relax_test: allow running just one iterationAnton Khirnov2019-01-29
|
* Fix callback type mismatchAnton Khirnov2019-01-29
|
* ell_grid_solve: remove unused variablesAnton Khirnov2019-01-29
|
* ell_grid_solve: add missing headerAnton Khirnov2019-01-29
|
* mg2d: fail when a full multigrid iteration does not reduce the residualAnton Khirnov2019-01-29
|
* Solve the discretized system exactly on the coarsest level.Anton Khirnov2019-01-26
|
* ell_relax -> ell_grid_solveAnton Khirnov2019-01-26
| | | | | Generalize the API to allow for multiple solver types. This is done in preparation for the exact linear system inversion solver.
* mg2d: factor out the boundary condition-related APIAnton Khirnov2019-01-25
|
* Stop duplicating some constants between ell_relax and mg2d.Anton Khirnov2019-01-25
|
* mg2d: split the constant definitions into their own headerAnton Khirnov2019-01-25
| | | | This will allow sharing them with ell_relax.
* ell_relax: implement correct handling of the boundary cornersAnton Khirnov2019-01-25
|
* mg2d: log the convergence factors after relax/correctAnton Khirnov2019-01-22
|
* mg2d: fix typoAnton Khirnov2019-01-21
| | | | | Do nb_relax_post steps in the post-correct relaxation phase, not nb_relax_pre.
* mg2d: failer earlier when relaxation does not convergeAnton Khirnov2019-01-21
|
* mg2d: check and propagate the error from mg_solve_subgrid()Anton Khirnov2019-01-21
|
* mg2d: factor out the relaxation step callAnton Khirnov2019-01-16
|
* mg2d: add a context variable to control the log levelAnton Khirnov2019-01-16
| | | | | | Simpler to use than overriding the log callback. API bump.
* log: allow specifying a NULL callback to disable loggingAnton Khirnov2019-01-16
|
* mg2d: add a function for getting maximum supported fd_stencilAnton Khirnov2019-01-16
| | | | | | | This is needed for the caller to know the number of boundary zones that can be filled. API bump.
* ell_relax: compute the residual norm in residual_calc()Anton Khirnov2019-01-13
| | | | It is cheap and avoids an extra step in mg2d.
* ell_relax: do not calculate the residual at the fixed-val boundariesAnton Khirnov2019-01-13
| | | | | Keep it at zero there to begin with, rather than enforce this later. This will be useful in the following commits.
* residual_calc.asm: fix partial storesAnton Khirnov2019-01-13
| | | | .store1 and .store3 were switched