Skip to content

Instantly share code, notes, and snippets.

View ChrisRackauckas's full-sized avatar
🎯
Focusing

Christopher Rackauckas ChrisRackauckas

🎯
Focusing
View GitHub Profile
; Function muladd2!
; Location: C:\Users\Chris Rackauckas\OneDrive\Computer\Desktop\test.jl:21
; Function Attrs: uwtable
define nonnull %jl_value_t addrspace(10)* @"japi1_muladd2!_35788"(%jl_value_t addrspace(10)*, %jl_va
lue_t addrspace(10)**, i32) #0 {
top:
%gcframe = alloca %jl_value_t addrspace(10)*, i32 3 %3 = bitcast %jl_value_t addrspace(10)** %gcframe to i8*
call void @llvm.memset.p0i8.i32(i8* %3, i8 0, i32 24, i32 0, i1 false)
%4 = alloca %jl_value_t addrspace(10)**, align 8
store volatile %jl_value_t addrspace(10)** %1, %jl_value_t addrspace(10)*** %4, align 8
; Function muladd1!
; Location: C:\Users\Chris Rackauckas\OneDrive\Computer\Desktop\test.jl:16
; Function Attrs: uwtable
define nonnull %jl_value_t addrspace(10)* @"japi1_muladd1!_35787"(%jl_value_t addrspace(10)*, %jl_va
lue_t addrspace(10)**, i32) #0 {
top:
%gcframe = alloca %jl_value_t addrspace(10)*, i32 3 %3 = bitcast %jl_value_t addrspace(10)** %gcframe to i8*
call void @llvm.memset.p0i8.i32(i8* %3, i8 0, i32 24, i32 0, i1 false)
%4 = alloca %jl_value_t addrspace(10)**, align 8
store volatile %jl_value_t addrspace(10)** %1, %jl_value_t addrspace(10)*** %4, align 8
# N-body simulation
# Based off of the Computer Language Benchmarks Game
using Printf
using LinearAlgebra
# Constants
const SOLAR_MASS = 4*pi*pi
const DAYS_PER_YEAR = 365.24
Count File Line Function
1 .\abstractarray.jl 83 axes(::Array{Any,1})
3 .\abstractarray.jl 83 axes
6 .\abstractarray.jl 831 copymutable
1 .\abstractarray.jl 654 copyto!(::Array{Int64,1}, ::Core.Compiler.BitSet)
1 .\abstractarray.j
Note:
┌ Warning: The profile data buffer is full; profiling probably terminated
│ before your program finished. To profile for longer runs, call
│ `Profile.init()` with a larger buffer and/or larger delay.
└ @ Profile C:\cygwin\home\Administrator\buildbot\worker\package_win64\build\usr\share\julia\stdlib\v0.7\Profile\src\Profile.jl:315
was thrown when printing. Profile results:
Count File Line Function
@ChrisRackauckas
ChrisRackauckas / rks.jl
Last active March 7, 2018 10:50
Just showing some RKs, no biggie.
const dz = 0.1
const z = 10
const space_steps = Int(z / dz)
const time_steps = 365 * 24 * 60
const ts = linspace(0, 12, time_steps)
const dt = (365 * 24 * 60 * 60) / time_steps
const alpha = 2e-6
const starting_temp = 30.0
using ApproxFun, OrdinaryDiffEq, Sundials, BenchmarkTools
using Plots; gr()
S = Fourier()
n = 100 # n = 1000 takes forever
T = ApproxFun.plan_transform(S, n)
Ti = ApproxFun.plan_itransform(S, n)
x = points(S, n)
function fast_gm(t,u,p,du)
a,α,Abar,β,D1,D2 = p
@inbounds for i in 2:N-1, j in 2:N-1
du[i,j,1] = D1*(u[i-1,j,1] + u[i+1,j,1] + u[i,j+1,1] + u[i,j-1,1] - 4u[i,j,1]) +
a*u[i,j,1]^2/u[i,j,2] + Abar - α*u[i,j,1]
end
@inbounds for i in 2:N-1, j in 2:N-1
du[i,j,2] = D2*(u[i-1,j,2] + u[i+1,j,2] + u[i,j+1,2] + u[i,j-1,2] - 4u[i,j,2]) +
Console
vals = [0.0, 0.0, 0.0]
nx.value = 1.261707777416795e-13
ηi = [0.0, 0.0, 0.0]
nx = 1.2042517040202507e-13
vals = [-188.765, -148.286, 7.57139]
nx.value = 6.778965875676743e68
vals = [-3.01714, -2.37014, 0.121018]
rm: refusing to remove '.' or '..' directory: skipping '/home/jrun/pkgdefault/.'
rm: refusing to remove '.' or '..' directory: skipping '/home/jrun/pkgdefault/..'
INFO: is_add = true and is_clone = true
INFO: Packages to add : SubString{String}["DifferentialEquations"]
INFO: Initializing package repository /home/jrun/pkgdefault/.julia/v0.6
INFO: Cloning METADATA from https://github.com/JuliaLang/METADATA.jl
INFO: Cloning cache of AlgebraicDiffEq from https://github.com/JuliaDiffEq/AlgebraicDiffEq.jl.git
INFO: Cloning cache of BandedMatrices from https://github.com/JuliaMatrices/BandedMatrices.jl.git
INFO: Cloning cache of BinDeps from https://github.com/JuliaLang/BinDeps.jl.git
INFO: Cloning cache of BoundaryValueDiffEq from https://github.com/JuliaDiffEq/BoundaryValueDiffEq.jl.git