mforets

using ReachabilityAnalysis
using ReachabilityAnalysis.TaylorModels
using ReachabilityAnalysis.TaylorSeries
using Random
using StatsBase

# -----------------
# Numbers
# -----------------

mforets

using Plots
using Unitful
using DifferentialEquations
using GeneralAstrodynamics
using LazySets

# Find a periodic orbit near Earth (orbit, period)
Oₑ, Pₑ = halo(SunEarth; Az=100_000u"km", L=2)

# Find a periodic orbit near Jupiter (orbit, period)

mforets

### Code

```julia
function makie_3d_plot(scene, data::Vector{IntervalBox{3, N}}, color=nothing) where N

    # plotting a single box in 3D doesn't work as expected
    @assert length(data) > 1

    positions = [Point3f0(mid(X)) for X in data]
    scales = [Vec3f0(diam.(X)) for X in data]

mforets

Cursos ofrecidos por el DMA (CURE, Maldonado) para el Diploma de Espec. en Matem.

Nombre del curso	Semestre dictado	Docente encargado
Álgebra Lineal Numérica	1er semestre 2020	Marcelo Forets

mforets

$ ./runTests.sh 
keyfiles =
{
  [1,1] = onsasExample_staticVonMisesTruss_test.m
  [1,2] = onsasExample_uniaxialExtension_test.m
  [1,3] = onsasExample_uniformCurvatureCantilever_test.m
}

totalRuns =  3

mforets

Status

mforets@localhost arch-comp-linear-2020]$ julia
               _
   _       _ _(_)_     |  Documentation: https://docs.julialang.org
  (_)     | (_) (_)    |
   _ _   _| |_  __ _   |  Type "?" for help, "]?" for Pkg help.
  | | | | | | |/ _` |  |
 | | |_| | | | (_| | | Version 1.4.1 (2020-04-14)

mforets

using Plots, BenchmarkTools, StaticArrays
using Revise, ReachabilityAnalysis
const RA = ReachabilityAnalysis
include("/home/mforets/.julia/dev/ReachabilityAnalysis/test/models/hybrid/embrake.jl")
LazySets.deactivate_assertions()

function LazySets.diameter(R::ReachabilityAnalysis.AbstractLazyReachSet; vars::Int)
    overapproximate(project(R, vars=vars), Interval) |> set |> diameter
end

mforets

no pv, no jitter, delta = 1e-7, GLGM06

first_in = findfirst((X->begin
                #= In[14]:11 =#
                set(X) ⊆ P
            end), vsol.Xk) = 888783
all((X->begin
            #= In[14]:14 =#
            set(X) ⊆ P

mforets

#using Plots, BenchmarkTools, StaticArrays #using Revise, ReachabilityAnalysis #const RA = ReachabilityAnalysis include("/home/mforets/.julia/dev/ReachabilityAnalysis/test/models/hybrid/embrake.jl") LazySets.deactivate_assertions()

function LazySets.diameter(R::ReachabilityAnalysis.AbstractLazyReachSet; vars::Int) overapproximate(project(R, vars=vars), Interval) |> set |> diameter end

mforets

```julia
using Plots, BenchmarkTools
using Revise, ReachabilityAnalysis
include("/home/mforets/.julia/dev/ReachabilityAnalysis/test/models/hybrid/embrake.jl")
LazySets.deactivate_assertions()
```

## No pv no jitter

```