Skip to content

Instantly share code, notes, and snippets.

View sritchie's full-sized avatar
🎯
Focusing

Sam Ritchie sritchie

🎯
Focusing
View GitHub Profile
@sritchie
sritchie / problem8.clj
Created September 6, 2012 05:00 — forked from johnynek/problem8.clj
Problem 8
(def bignumber "7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379583319528532088055111254069874715852386305071569329096329522744304355766896648950445244523161731856403098711121722383113622298934233803081353362766142828064444866452387493035890729629049156044077239071381051585930796086670172427121883998797908792274921901699720888093776657273330010533678812202354218097512545405947522435258490771167055601360483958644670632441572215539753697817977846174064955149290862569321978468622482839722413756570560574902614079729686524145351004748216637048440319989000889524345065854122758866688116427171479924442928230863465674813919123162824586178664583591245665294765456828489128831426076900422421902267105562632111110937054421750694165896040807198403850962455444362981230987879927244284909188845801561660979191338754992005240636899125607176060588611646710940507754100225698315520005593572972571636269561882670428252483600823257530420752963450")
(def
@sritchie
sritchie / problem8.clj
Created September 6, 2012 04:59 — forked from johnynek/problem8.clj
Problem 8
(def bignumber "7316717653133062491922511967442657474235534919493496983520312774506326239578318016984801869478851843858615607891129494954595017379583319528532088055111254069874715852386305071569329096329522744304355766896648950445244523161731856403098711121722383113622298934233803081353362766142828064444866452387493035890729629049156044077239071381051585930796086670172427121883998797908792274921901699720888093776657273330010533678812202354218097512545405947522435258490771167055601360483958644670632441572215539753697817977846174064955149290862569321978468622482839722413756570560574902614079729686524145351004748216637048440319989000889524345065854122758866688116427171479924442928230863465674813919123162824586178664583591245665294765456828489128831426076900422421902267105562632111110937054421750694165896040807198403850962455444362981230987879927244284909188845801561660979191338754992005240636899125607176060588611646710940507754100225698315520005593572972571636269561882670428252483600823257530420752963450")
(def
(ns ybot.kwphrase-tests
(:use [ybot.analytics.ga kwphrase]
[ybot datastores]
[midje sweet cascalog]))
(let [ga-data [["20121025"
"http://dopeness.org"
"United States"
"Oregon"
"Portland"
@sritchie
sritchie / sam.clj
Created January 25, 2012 18:02 — forked from robinkraft/sam.clj
first macro!
(defmacro sam [expr pos zero neg]
(if (< 0 (eval expr))
pos
(if (> 0 (eval expr))
neg
zero)))
@sritchie
sritchie / gist:1667190
Created January 24, 2012 01:39 — forked from clizzin/gist:1667186
moment function for cascalog.math
(def moment
(<- [!val !k :> !moment]
(c/avg !val :> !mean)
(- !val !mean :> !dev)
(expt !dev !pow :> !pow-dev)
(c/sum !pow-dev :> !pow-dev-sum)
(div !pow-dev-sum !count :> !moment)))
@sritchie
sritchie / gist:1627929
Created January 17, 2012 18:16 — forked from danhammer/gist:1627819
timing matrix multiplication
import numpy as np
def test_my_mult(n):
A = np.random.rand(n*23).reshape(n,23)
At = A.T
t0 = time.time()
res = np.dot(A.T, A)
print time.time() - t0
print np.shape(res)
@sritchie
sritchie / matrix.py
Created January 17, 2012 18:11 — forked from danhammer/gist:1627807
timing matrix multiplication
(defn feature-vec [n]
(map (partial cons 1)
(for [x (range n)]
(take 22 (repeatedly rand)))))
(defn dot-product [x y]
(reduce + (map * x y)))
(defn transpose
"returns the transposition of a `coll` of vectors"
(defn drop-every-nth
[n coll]
(keep-indexed (fn [idx val]
(when (pos? (mod (inc idx) n))
val))
coll))
(defn top-n-seq [sample-name n]
(let [aa-seq-fields ["?aa" "?v" "?j" "?sn" "?nc" "?c" "?nf" "?rf" "?cdr3"]
sample-name-seqs (<- aa-seq-fields
(aa-seqs :>> aa-seq-fields)
(= ?sn sample-name))]
(<- aa-seq-fields ((c/first-n sample-name-seqs n :sort ["?nc"] :reverse true)
:>> aa-seq-fields))))
(defn top-n-seq [sample-name n]
(let [aa-seq-fields ["?aa" "?v" "?j" "?sn" "?nc" "?c" "?nf" "?rf" "?cdr3"]
sample-name-seqs (<- aa-seq-fields
(aa-seqs :>> aa-seq-fields)
(= ?sn sample-name))]
(<- aa-seq-fields ((c/first-n sample-name-seqs n :sort ["?nc"] :reverse true)
:>> aa-seq-fields))))
(let [aa-seq-fields ["?aa" "?v" "?j" "?sn" "?nc" "?c" "?nf" "?rf" "?cdr3"]]
(?<- (stdout) aa-seq-fields