Skip to content

Instantly share code, notes, and snippets.

View kohnakagawa's full-sized avatar
🌴
On vacation

tsunekoh kohnakagawa

🌴
On vacation
View GitHub Profile
#include <iostream>
#include <fstream>
#include <sstream>
#include <vector>
#include <cmath>
#include "Eigen/Core"
#include "Eigen/Eigenvalues"
#include <cstdlib>
//global val
#include <iostream>
#include <cstdlib>
#include <sys/time.h>
#define BLOCK_SIZE 256
#define GRID_SIZE 500
#define SYS_SIZE (BLOCK_SIZE * GRID_SIZE)
double Uniform(){
return ((double)rand()+1.0)/((double)RAND_MAX+2.0);
#include <iostream>
#include <fstream>
#include <sstream>
#include <string>
#include <vector>
#include <cassert>
#include <cmath>
#include <cstdlib>
#define INLINE __attribute__((always_inline))
#small sample for GMM
#install.packages('mclust') #installing Mclust package
args <- commandArgs(trailingOnly = T)
infile <- args[1]
outfile <- args[2]
time <- args[3]
library(mclust)
dat <- read.table(infile)
dat <- data.matrix(dat)
result.dat <- Mclust(dat)
set dgrid3d 30,30
set hidden3d
set cntrparam levels incremental 1.2,0.1,1.8
set contour
set view 0,0
unset surface
set format z ""
set key outside
splot "chemeq.txt" lw 4 title "surface density"
set xlabel "{/Symbol s}/R" offset 4.,0.8 font "Helvetica, 20"
#include <cmath>
#include <iostream>
#include <fstream>
#include <string>
#include <cassert>
#include <cstdlib>
#include <vector>
#include <sstream>
#define INLINE __attribute__((always_inline))
options(scipen=2)
btime <- 1240000
etime <- 3960000
stptime <- 40000
time <- seq(btime,etime,by=stptime)
sepval <- 32.0
tardir <- "./cutof0.65/chem5e-5inv0.65h0.20o0.052"
for(i in time){
#infile <- paste(tardir,"/z_cord_hyphil_time",i,".txt",sep="")
infile <- paste(tardir,"/z_cord_hyphob_time",i,".txt",sep="")
#include <iostream>
#include <fstream>
#include <string>
#include <vector>
#include <algorithm>
#include <cmath>
#include <cassert>
struct set_XY{
double q[2];
#include <cmath>
#include <iostream>
#include <fstream>
#include <string>
#include <cassert>
#include <cstdlib>
#include <vector>
#include <sstream>
#define INLINE __attribute__((always_inline))
#!/usr/bin/python
import os
import sys
argvs = sys.argv
argc = len(argvs)
if (argc != 4):
print "Usage $ python %s beg_time end_time pud_digit" % argvs[0]
quit()
beg_time = int(argvs[1])
end_time = int(argvs[2])