Skip to content

Instantly share code, notes, and snippets.

@anyuzx
anyuzx / save_sparse_matrix.py
Created Oct 3, 2016
save spicy sparse matrix/array object to file
View save_sparse_matrix.py
from scipy import sparse, io
m = sparse.csr_matrix([[0,0,0],[1,0,0],[0,1,0]])
m # <3x3 sparse matrix of type '<type 'numpy.int64'>' with 2 stored elements in Compressed Sparse Row format>
io.mmwrite("test.mtx", m)
del m
newm = io.mmread("test.mtx")
newm # <3x3 sparse matrix of type '<type 'numpy.int32'>' with 2 stored elements in COOrdinate format>
@anyuzx
anyuzx / acf.py
Created Sep 28, 2016
autocorrelation in python
View acf.py
# Since there is no real autocorrelation compute function in python and numpy
# I write my own one
import numpy as np
# Suppose data array represent a time series data
# each element represent data at given time
# Assume time are equally spaced
def acf(data):
mean = np.mean(data)
@anyuzx
anyuzx / lammps_compute_gyration.md
Last active Feb 28, 2017
walk through LAMMPS compute_gyration.h and compute_gyration.cpp
View lammps_compute_gyration.md

Notes about LAMMPS source codes, compute_gyration.h and compute_gyration.cpp.

compute_gyration.h

#ifdef COMPUTE_CLASS

ComputeStyle(gyration,ComputeGyration)

#else
@anyuzx
anyuzx / 0_reuse_code.js
Created Feb 24, 2016
Here are some things you can do with Gists in GistBox.
View 0_reuse_code.js
// Use Gists to store code you would like to remember later on
console.log(window); // log the "window" object to the console