File tree 5 files changed +33
-25
lines changed 5 files changed +33
-25
lines changed Original file line number Diff line number Diff line change 3
3
This package implements the stochastic Galerkin finite element method for certain model problems,
4
4
including iterative solvers and a posteriori error control in two dimensions.
5
5
6
- ## Installation
7
-
8
- Instantiate this project after cloning or updating via
9
- ```
10
- $ julia --project=.
11
- julia> using Pkg
12
- julia> Pkg.instantiate()
13
- ```
14
-
15
6
## Scripts/Examples
16
7
17
8
- ` scripts/poisson.jl ` : Adaptive stochastic Galerkin FEM for the Poisson problem with stochastic diffusion coefficient
Original file line number Diff line number Diff line change 1
1
using Documenter
2
2
using ExtendableASGFEM
3
3
4
- push! (LOAD_PATH ," ../src/" )
5
-
6
4
makedocs (
7
5
modules = [ExtendableASGFEM],
8
6
sitename = " ExtendableASGFEM.jl" ,
9
7
authors = " Christian Merdon, Martin Eigel" ,
10
- repo = " https://lab.wias-berlin.de/merdon/ASGFEMJulia" ,
11
- # repolink = "https://lab.wias-berlin.de/merdon/ASGFEMJulia",
12
- format = Documenter. HTML (; mathengine = MathJax3 ()),
13
- clean = false ,
14
- source = " src" ,
15
- build = " build" ,
8
+ format = Documenter. HTML (; repolink = " https://github.com/WIAS-PDELib/ExtendableASGFEM.jl" , mathengine = MathJax3 ()),
9
+ clean = false ,
16
10
checkdocs = :all ,
17
- warnonly = true ,
11
+ warnonly = false ,
18
12
doctest = true ,
19
13
pages = [
20
14
" Home" => " index.md"
@@ -37,4 +31,9 @@ makedocs(
37
31
" plots.md"
38
32
]
39
33
],
40
- )
34
+ )
35
+
36
+
37
+ deploydocs (
38
+ repo = " github.com/WIAS-PDELib/ExtendableASGFEM.jl" ,
39
+ )
Original file line number Diff line number Diff line change @@ -13,7 +13,7 @@ the subresiduals that are associated to each multi-index
13
13
14
14
``` @autodocs
15
15
Modules = [ExtendableASGFEM]
16
- Pages = ["estimators .jl"]
16
+ Pages = ["estimate .jl"]
17
17
Order = [:type, :function]
18
18
```
19
19
@@ -29,4 +29,17 @@ Here are some methods that help with enriching the set of multi-indices.
29
29
Modules = [ExtendableASGFEM]
30
30
Pages = ["mopcontrol.jl"]
31
31
Order = [:type, :function]
32
+ ```
33
+
34
+ ## Monte carlo sampling estimator
35
+
36
+ There is also a hierarchical Monte carlo error estimator available that
37
+ compares the solution with a higher order discrete solution for sampled
38
+ deterministic problems. This is merely intended as a way to compute the
39
+ reference error to assess the efficiency of the residual-based error estimator.
40
+
41
+ ``` @autodocs
42
+ Modules = [ExtendableASGFEM]
43
+ Pages = ["sampling_error.jl"]
44
+ Order = [:type, :function]
32
45
```
Original file line number Diff line number Diff line change 2
2
3
3
([source code](SOURCE_URL))
4
4
5
- # runs AFEM loop for stochastic Poisson problem and saves results to data folder
5
+ runs AFEM loop for stochastic Poisson problem, save/load results to/from data folder, plot results
6
+
7
+ main usage:
8
+ - run experiment: run(; problem = problem, kwargs...)
9
+ - load results : show_results(; kwargs...)
10
+ - produce plots : produce_plots(; kwargs...)
6
11
7
12
possible values for problem are
8
- PoissonProblemPrimal = Poisson problem with linear coefficient a
9
- LogTransformedPoissonProblemPrimal = log-transformed Poisson problem with exponential coefficient exp(a)
10
- LogTransformedPoissonProblemDual = dual formulation
13
+ - PoissonProblemPrimal = Poisson problem with linear coefficient a
14
+ - LogTransformedPoissonProblemPrimal = log-transformed Poisson problem with exponential coefficient exp(a)
15
+ - LogTransformedPoissonProblemDual = dual formulation
11
16
12
17
=#
13
18
@@ -64,7 +69,7 @@ default_args = Dict(
64
69
" θ_spatial" => 0.5 ,
65
70
" factor_tail" => 1 ,
66
71
" tail_extension" => [10 ,2 ], # for [0] mode, and all others
67
- " maxdofs" => 1e5 ,
72
+ " maxdofs" => 1e4 ,
68
73
" initial_modes" => [[0 ]],
69
74
" nsamples" => 150 ,
70
75
" use_equilibration_estimator" => false ,
You can’t perform that action at this time.
0 commit comments