Skip to content

Instantly share code, notes, and snippets.

@jgoldfar
Created December 4, 2014 23:29
Show Gist options
  • Save jgoldfar/cb0b0acaf49f69dd03fe to your computer and use it in GitHub Desktop.
Save jgoldfar/cb0b0acaf49f69dd03fe to your computer and use it in GitHub Desktop.
segfault->error
cd(dirname(@__FILE__))
include(joinpath("src", "ISP.jl"))
function doISP()
include("ex3.jl")
ed=ISP.ExData(a,b,c,s,u,ux,ut,sp,γf)
sd=ISP.SolvData(ed,10)
const tgrid = sd.grid.tgrid
v_actual = convert(Vector{Float64}, [sd.grid.xbdy, [sd.g(t) for t in tgrid]])
v0 = copy(v_actual)
len = length(v0)
v0[1:(len/2)] = mean(v0[1:(len/2)])
v0[(len/2+1):len] = mean(v0[(len/2+1):len])
mu_actual = Float64[u(s(t), t) for t in tgrid]
nu_actual = Float64[u(zero(t), t) for t in tgrid]
delta = 0.5
l = 5.0
R = 20.0
func=ISP.SimpleTest(sd, mu_actual, nu_actual, v0, delta, l, R)
end
doISP()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment