Skip to content

Commit

Permalink
use OrdinaryDiffEqRosenbrock
Browse files Browse the repository at this point in the history
  • Loading branch information
JoshuaLampert committed Dec 4, 2024
1 parent bf86e29 commit edf0572
Show file tree
Hide file tree
Showing 11 changed files with 14 additions and 11 deletions.
3 changes: 2 additions & 1 deletion README.md
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,8 @@ julia> Pkg.add("Plots")

To create special node sets, you might also want to install [QuasiMonteCarlo.jl](https://github.com/SciML/QuasiMonteCarlo.jl) or
[Meshes.jl](https://github.com/JuliaGeometry/Meshes.jl) and for solving time-dependent partial differential equations
[OrdinaryDiffEq.jl](https://github.com/SciML/OrdinaryDiffEq.jl) in a similar way as above for Plots.jl. See the documentation for more
[OrdinaryDiffEq.jl](https://github.com/SciML/OrdinaryDiffEq.jl) in a similar way as above for Plots.jl. Consider using the subpackage
OrdinaryDiffEqRosenbrock.jl as it contains the most relevant time integration schemes for this package. See the documentation for more
examples on how to use these packages in combination with KernelInterpolation.jl.

## Usage
Expand Down
2 changes: 1 addition & 1 deletion docs/src/development.md
Original file line number Diff line number Diff line change
Expand Up @@ -19,7 +19,7 @@ cd KernelInterpolation.jl
mkdir run
cd run
julia --project=. -e 'using Pkg; Pkg.develop(PackageSpec(path=".."))' # Install local KernelInterpolation.jl clone
julia --project=. -e 'using Pkg; Pkg.add(["Plots", "QuasiMonteCarlo", "Meshes", "OrdinaryDiffEq"])' # Install additional packages
julia --project=. -e 'using Pkg; Pkg.add(["Plots", "QuasiMonteCarlo", "Meshes", "OrdinaryDiffEqRosenbrock"])' # Install additional packages
```

If you use other packages for executing KernelInterpolation.jl, you can add them to the project in the `run`
Expand Down
3 changes: 2 additions & 1 deletion docs/src/index.md
Original file line number Diff line number Diff line change
Expand Up @@ -38,7 +38,8 @@ julia> Pkg.add("Plots")

To create special node sets, you might also want to install [QuasiMonteCarlo.jl](https://github.com/SciML/QuasiMonteCarlo.jl) or
[Meshes.jl](https://github.com/JuliaGeometry/Meshes.jl) and for solving time-dependent partial differential equations
[OrdinaryDiffEq.jl](https://github.com/SciML/OrdinaryDiffEq.jl) in a similar way as above for Plots.jl. See the documentation for more
[OrdinaryDiffEq.jl](https://github.com/SciML/OrdinaryDiffEq.jl) in a similar way as above for Plots.jl. Consider using the subpackage
OrdinaryDiffEqRosenbrock.jl as it contains the most relevant time integration schemes for this package. See the documentation for more
examples on how to use these packages in combination with KernelInterpolation.jl.

## Usage
Expand Down
2 changes: 2 additions & 0 deletions docs/src/pdes.md
Original file line number Diff line number Diff line change
Expand Up @@ -226,6 +226,8 @@ which already provides a wide range of time integration methods. Note that this
system, which is more difficult to solve than a simple ODE system. Thus, we are restricted to specialized time integration
methods, which can handle DAEs. We recommend using the `Rodas5P` method, which is a Rosenbrock method for stiff DAEs. See
also the [documentation of OrdinaryDiffEq.jl](https://docs.sciml.ai/DiffEqDocs/latest/tutorials/dae_example/) for more information.
`Rodas5P` along with other Rosenbrock methods are implemented in the subpackage `OrdinaryDiffEqRosenbrock.jl`. Therefore, we
recommend using this package in combination with KernelInterpolation.jl to reduce precompile time.

In KernelInterpolation.jl, you can use the constructor of a [`Semidiscretization`](@ref) in a very similar way as
[`SpatialDiscretization`](@ref), but with the additional initial condition. This can be turned into an `ODEProblem` object
Expand Down
2 changes: 1 addition & 1 deletion examples/PDEs/advection_1d_basic.jl
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
using KernelInterpolation
using OrdinaryDiffEq
using OrdinaryDiffEqRosenbrock
using Plots

# source term of advection equation
Expand Down
2 changes: 1 addition & 1 deletion examples/PDEs/advection_3d_basic.jl
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
using KernelInterpolation
using OrdinaryDiffEq
using OrdinaryDiffEqRosenbrock
using LinearAlgebra: norm
using WriteVTK: WriteVTK, paraview_collection

Expand Down
2 changes: 1 addition & 1 deletion examples/PDEs/advection_diffusion_2d_basic.jl
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
using KernelInterpolation
using OrdinaryDiffEq
using OrdinaryDiffEqRosenbrock
using LinearAlgebra: norm
using Plots

Expand Down
2 changes: 1 addition & 1 deletion examples/PDEs/heat_2d_basic.jl
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
using KernelInterpolation
using OrdinaryDiffEq
using OrdinaryDiffEqRosenbrock
using Plots

# right-hand-side of Heat equation
Expand Down
2 changes: 1 addition & 1 deletion examples/PDEs/heat_2d_manufactured.jl
Original file line number Diff line number Diff line change
@@ -1,5 +1,5 @@
using KernelInterpolation
using OrdinaryDiffEq
using OrdinaryDiffEqRosenbrock
using Plots

# right-hand-side of heat equation
Expand Down
3 changes: 1 addition & 2 deletions test/Project.toml
Original file line number Diff line number Diff line change
Expand Up @@ -3,7 +3,7 @@ Aqua = "4c88cf16-eb10-579e-8560-4a9242c79595"
ExplicitImports = "7d51a73a-1435-4ff3-83d9-f097790105c7"
LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
Meshes = "eacbb407-ea5a-433e-ab97-5258b1ca43fa"
OrdinaryDiffEq = "1dea7af3-3e70-54e6-95c3-0bf5283fa5ed"
OrdinaryDiffEqRosenbrock = "43230ef6-c299-4910-a778-202eb28ce4ce"
Plots = "91a5bcdd-55d7-5caf-9e0b-520d859cae80"
QuasiMonteCarlo = "8a4e6c94-4038-4cdc-81c3-7e6ffdb2a71b"
Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
Expand All @@ -18,7 +18,6 @@ Aqua = "0.8.3"
ExplicitImports = "1.0.1"
LinearAlgebra = "1"
Meshes = "0.52.1"
OrdinaryDiffEq = "6.68"
Plots = "1.25.11"
QuasiMonteCarlo = "0.3.1"
Random = "1"
Expand Down
2 changes: 1 addition & 1 deletion test/runtests.jl
Original file line number Diff line number Diff line change
Expand Up @@ -9,7 +9,7 @@ end

@testsnippet AdditionalImports begin
using LinearAlgebra: norm, Symmetric, I
using OrdinaryDiffEq: solve, Rodas5P
using OrdinaryDiffEqRosenbrock: solve, Rodas5P
using StaticArrays: MVector
using Meshes: Meshes, Sphere, Point, PointSet, RegularSampling
end

0 comments on commit edf0572

Please sign in to comment.