Skip to content

Instantly share code, notes, and snippets.

@albop
albop / gist:7525675
Created Nov 18, 2013
zero_based_indexing
View gist:7525675
rewrite(x::Number) = x
rewrite(x::Symbol) = x
rewrite(x::String) = x
function rewrite(expr)
if expr.head == :ref
return Expr(expr.head, expr.args[1], [rewrite(:( 1 + $i)) for i in expr.args[2:end]]...)
else
return Expr(expr.head, [rewrite(i) for i in expr.args]...)
end
@albop
albop / gist:10281951
Created Apr 9, 2014
Browser with tabs
View gist:10281951
#!/usr/bin/python
# -*- coding: utf-8 -*-
from __future__ import print_function
import sys
from PySide.QtCore import *
from PySide.QtGui import *
from PySide import QtGui
from PySide.QtWebKit import *
from PySide import QtWebKit
View gist:4b82eb9cdc09f82bbbd4
######################################################################
using PyPlot
using Grid
using Dierckx
using NumericalMath
using splines
#using ApproXD
meshgrid(x::Vector, y::Vector) = (repmat(x, 1, length(y))',
repmat(y, 1, length(x)))
@albop
albop / gist:ab49de476b6bcbdd689c
Last active Aug 29, 2015
Functions needed for vfi with dolo
View gist:ab49de476b6bcbdd689c
from dolo import yaml_import
model = yaml_import('examples/models/rbc.yaml')
s = model.calibration['states']
x = model.calibration['controls']
y = model.calibration['auxiliaries']
e = model.calibration['shocks']
v = model.calibration['values']
p = model.calibration['parameters']
View watcher.py
"""
Filename: watcher.py
Authors: Pablo Winant
Time long computation and send a message when finished.
"""
import contextlib
@contextlib.contextmanager
def watcher(task_name=None, email=None):
@albop
albop / index.html
Last active Aug 29, 2015
mc_sample_path: numba annotate
View index.html
<html>
<head>
<style>
.annotation_table {
color: #000000;
font-family: monospace;
margin: 5px;
@albop
albop / interpolation_speed.jl
Created Jan 7, 2016
splines.jl vs interpolations.jl
View interpolation_speed.jl
d = 3 # number of dimensions
K = 50 # number of points along each dimension
N = 100000 # number of points at which to interpolate
A = rand([K for i = 1:d]...) # filtered coefficients
B = rand(N,d) # points at which to evaluate the splines
max(B, minimum(A)+0.01)
min(B, maximum(A)-0.01)
@albop
albop / eval_cubic_cuda.py
Last active Jan 19, 2016
interpolation with numba.cuda
View eval_cubic_cuda.py
from __future__ import division
from numba import double, int64
from numba import jit, njit
import numpy as np
from numpy import zeros, array
@albop
albop / generated_interp.py
Last active Jan 21, 2016
Half working generated interpolation example
View generated_interp.py
from math import floor
from numba import njit
####
#### Working
####
@njit
def native_index_1d(mat, vec):
return mat[vec[0]]
View erfc_libm vs erfc_python
from cffi import FFI
ffi = FFI()
ffi.cdef('double erfc(double x);')
libm = ffi.dlopen("m")
erfc = libm.erfc
from math import erfc as p_erfc
from numba import njit
from numpy import linspace