This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | #include <iostream> | |
| #include <fstream> | |
| #include <sstream> | |
| #include <vector> | |
| #include <cmath> | |
| #include "Eigen/Core" | |
| #include "Eigen/Eigenvalues" | |
| #include <cstdlib> | |
| //global val | 
  
    
      This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | #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); | 
  
    
      This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | #include <iostream> | |
| #include <fstream> | |
| #include <sstream> | |
| #include <string> | |
| #include <vector> | |
| #include <cassert> | |
| #include <cmath> | |
| #include <cstdlib> | |
| #define INLINE __attribute__((always_inline)) | 
  
    
      This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | #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) | 
  
    
      This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | 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" | 
  
    
      This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | #include <cmath> | |
| #include <iostream> | |
| #include <fstream> | |
| #include <string> | |
| #include <cassert> | |
| #include <cstdlib> | |
| #include <vector> | |
| #include <sstream> | |
| #define INLINE __attribute__((always_inline)) | 
  
    
      This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | 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="") | 
  
    
      This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | #include <iostream> | |
| #include <fstream> | |
| #include <string> | |
| #include <vector> | |
| #include <algorithm> | |
| #include <cmath> | |
| #include <cassert> | |
| struct set_XY{ | |
| double q[2]; | 
  
    
      This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | #include <cmath> | |
| #include <iostream> | |
| #include <fstream> | |
| #include <string> | |
| #include <cassert> | |
| #include <cstdlib> | |
| #include <vector> | |
| #include <sstream> | |
| #define INLINE __attribute__((always_inline)) | 
  
    
      This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
      Learn more about bidirectional Unicode characters
    
  
  
    
  | #!/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]) |