Skip to content
New issue

Have a question about this project? # for a free GitHub account to open an issue and contact its maintainers and the community.

#

By clicking “#”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? # to your account

DEU tutorial fix #298

Merged
merged 4 commits into from
Aug 7, 2020
Merged

DEU tutorial fix #298

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
665f175
update DEU deps
Aug 7, 2020
4bb5999
clean-up DEU
Aug 7, 2020
cd6e581
Update Project.toml
ChrisRackauckas Aug 7, 2020
108d530
Update tutorials/DiffEqUncertainty/02-AD_and_optimization.jmd
ChrisRackauckas Aug 7, 2020
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
30 changes: 1 addition & 29 deletions tutorials/DiffEqUncertainty/01-expectation_introduction.jmd
View file
Original file line number Diff line number Diff line change
@@ -358,32 +358,4 @@ The performance gains realized by leveraging batch GPU processing is problem dep
```{julia; echo=false; skip="notebook"}
using SciMLTutorials
SciMLTutorials.tutorial_footer(WEAVE_ARGS[:folder],WEAVE_ARGS[:file])
```


<!--
## Batch + nout
```julia
using DiffEqGPU

function f(du, u,p,t)
@inbounds begin
du[1] = p[1]*u[1];
end
nothing
end

u0 = Float32[10.0]
p = Float32[-0.3]
tspan = (0.0f0,10.0f0)
prob = ODEProblem(f,u0,tspan,p)

g(sol) = [sol(4.0)[1], sol(6.0)[1]]

u0_dist = [truncated(Normal(3.0f0,2.0f0),-5f0,11f0)]
p_dist = [truncated(Normal(-.7f0, .1f0), -1f0,0f0)]

expectation(g, prob, u0_dist, p_dist, Koopman(), Tsit5(), EnsembleGPUArray();
quadalg = CubaSUAVE(), batch=1000, nout=2)[1]
```
-->
```
9 changes: 2 additions & 7 deletions tutorials/DiffEqUncertainty/02-AD_and_optimization.jmd
View file
Original file line number Diff line number Diff line change
@@ -277,7 +277,7 @@ The probability of impacting the wall is now
```
We can check if this is within tolerance by
```julia
isapprox(λ, 0.01, atol=1e-5)
abs(λ - 0.01) <= 1e-5
```

Again, we plot some Monte Carlo simulations from this result as follows
@@ -302,12 +302,7 @@ begin
end
```



## Conclusion


```{julia; echo=false; skip="notebook"}
using SciMLTutorials
SciMLTutorials.tutorial_footer(WEAVE_ARGS[:folder],WEAVE_ARGS[:file])
```
```
4 changes: 4 additions & 0 deletions tutorials/DiffEqUncertainty/Project.toml
View file
Original file line number Diff line number Diff line change
@@ -5,9 +5,11 @@ DiffEqSensitivity = "41bf760c-e81c-5289-8e54-58b1f1f8abe2"
DiffEqUncertainty = "ef61062a-5684-51dc-bb67-a0fcdec5c97d"
DifferentialEquations = "0c46a032-eb83-5123-abaf-570d42b7fbaa"
Distributions = "31c24e10-a181-5473-b8eb-7969acd0382f"
ForwardDiff = "f6369f11-7733-5829-9624-2563aa707210"
NLopt = "76087f3c-5699-56af-9a33-bf431cd00edd"
OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
Quadrature = "67601950-bd08-11e9-3c89-fd23fb4432d2"

[compat]
BenchmarkTools = "0.5"
@@ -16,6 +18,8 @@ DiffEqSensitivity = "6.28"
DiffEqUncertainty = "1.5"
DifferentialEquations = "6.15"
Distributions = "0.23"
ForwardDiff = "0.10"
NLopt = "0.6"
OrdinaryDiffEq = "5.42"
Plots = "1.5"
Quadrature = "1.3"