thegedge

# -----------------------------------------------------------------------------
# Copyright (C) 2014 Jason Gedge
#
# Permission is hereby granted, free of charge, to any person obtaining a copy
# of this software and associated documentation files (the "Software"), to deal
# in the Software without restriction, including without limitation the rights
# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
# copies of the Software, and to permit persons to whom the Software is
# furnished to do so, subject to the following conditions:
#

thegedge

merge :: Ord a => [a] -> [a] -> [a]
merge [] l = l
merge (a:l1) (b:l2)
  | a <= b    = a : merge l1 (b:l2)
  | otherwise = b : merge l2 (a:l1)

mergeSort :: Ord a => [a] -> [a]
mergeSort [] = []
mergeSort [a] = [a]
mergeSort l = 

thegedge

#!/usr/bin/env python
import timeit

def dot(a, b):
	"""Dot product"""
	return a[0]*b[0] + a[1]*b[1] + a[2]*b[2]

def rs_intersect(R0, d, C, r):
	"""Ray/sphere intersection"""
	R0 = (R0[0] - C[0], R0[1] - C[1], R0[2] - C[2])