sglyon

sglyon

sglyon

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [

sglyon

sglyon

using Roots
using ODE

@doc doc"""
ode_events{TY}(f::Function, y0::TY, tspan, odeint::Function; kwargs...)

### Arguments

* `f`: function such that dy/dt = f(t, y)
* `y0`: initial condition for y

sglyon

sglyon

Notice that when running it on a Linux VM on OSX I get about a 4x speedup relative to running it natively on OSX.

Also notice from the profile that the Linux version spends about 90% of its time in the myexp call and only about 3% in the line @inbounds d_i2_j += ((points[k, i] - data[k, j])^2). (I know these aren't actually times, but rather times when the profiler registered that line. But it gives an approximation to time spent).

The OSX version spends about 57% of its time (so about 5.3 seconds on average) in myexp and about 33% of its time (about 3 seconds on average) in the line @inbounds d_i2_j += ((points[k, i] - data[k, j])^2).

I can't explain why, on the same hardware, the OSX version is so much slower here.

sglyon

function [density, pts_PCn, di_min] = eval_density(data, points)
[n,d] = size(data);
n_points = size(points,1);

Datan = (data-ones(n,1)*mean(data))./(ones(n,1)*std(data));
[U,S,V] = svd(Datan,0);
PC = Datan*V;
PCn = PC./(ones(n,1)*std(PC));

Pointsn = (points-ones(n_points,1)*mean(data))./(ones(n_points,1)*std(data));

sglyon

sglyon

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [