acadien

(custom-set-variables
 ;; custom-set-variables was added by Custom.
 ;; If you edit it by hand, you could mess it up, so be careful.
 ;; Your init file should contain only one such instance.
 ;; If there is more than one, they won't work right.
 '(custom-safe-themes (quote ("8aebf25556399b58091e533e455dd50a6a9cba958cc4ebb0aab175863c25b9a4" "a8245b7cc985a0610d71f9852e9f2767ad1b852c2bdea6f4aadc12cce9c4d6d0" default)))
 '(inhibit-startup-screen t))
(custom-set-faces
 ;; custom-set-faces was added by Custom.
 ;; If you edit it by hand, you could mess it up, so be careful.

acadien

#!/usr/bin/python

import numpy as np
from scipy import integrate
import sys

import matplotlib as mpl
from matplotlib import pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
from matplotlib import animation

acadien

// Use Gists to store code you would like to remember later on
console.log(window); // log the "window" object to the console

acadien

#!/usr/bin/python

import sys,subprocess

#Simple code for starting jobs quick! Made for ATPESC 2014

def usage():
    print "%s <directory/name> <# processors>"%(sys.argv[0])

try:

acadien

import matplotlib as mpl
import pylab as pl

def colorHist(pngFile):
  a = pl.imread(pngFile)
  
  b = mpl.colors.rgb_to_hsv(a[...,:3])[...,1].flatten()
  
  vals,bins,dummy = pl.hist(b,256,visible=False)

acadien

#!/usr/bin/python

from numpy import *

x = random.random([10,3])
x[range(10), fabs(x-0.5).argsort()[:,0]]


#Break it down niao
x1 = fabs(x-0.5)     # numbers closest to 0.5 will be closest to 0 now.

acadien

#!/usr/bin/python

from scipy import weave
from scipy.weave import converters
from numpy import *

#Yeah that's right stick your code in a fucking comment cause its a weave hack.
msdCode = """
double c,d;
double *seta,*setb;

acadien

#!/usr/bin/python

import numpy as np
import mayavi.mlab as mlab

x, y, z = np.ogrid[-10:10:20j, -10:10:20j, -10:10:20j]
s = np.sin(x*y*z)/(x*y*z)


mlab.pipeline.image_plane_widget(mlab.pipeline.scalar_field(s),

acadien

#Python Euclidean distance
dist(a,b):
    return sum([(i-j)**2 for i,j in zip(a,b)])**0.5


#Scipy-Weave Euclidean distance (x10 faster than standard Python implementation)
distcode = """
double c=0.0;
for(int i=0;i<64;i++)
  c += (a[i]-b[i])*(a[i]-b[i]);

acadien

#First the standard python list style code:

#Calculate the distance between 2 atoms:
def dist(p1,p2):
    return sum([(a-b)*(a-b) for a,b in zip(p1,p2)])**0.5

#Calculate the angle between 3 atoms:
def ang(p1,p2,p3):
    P12=dist(p1,p2)
    P13=dist(p1,p3)