Sasha Kaurov lue
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
| import glob | |
| import os | |
| curdir = os.getcwd() | |
| dirs = [ name for name in os.listdir('./') if os.path.isdir(os.path.join('./', name)) ] | |
| for j in range(len(dirs)): | |
| os.chdir(curdir) | |
| os.chdir(dirs[j]) |
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
| 2 | |
| +---------------------------------------------+------------+------------+ | |
| | Total wallclock time elapsed since start | 71.7s | | | |
| | | | | | |
| | Section | no. calls | wall time | % of total | | |
| +---------------------------------+-----------+------------+------------+ | |
| | Assembly | 33 | 64.3s | 90% | | |
| | Output | 8 | 1.57s | 2.2% | | |
| | Solve | 33 | 2.3s | 3.2% | | |
| | setup | 1 | 1.88s | 2.6% | |
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
| # | |
| # Pick the MPI C compiler | |
| CC = mpicc | |
| # | |
| # | |
| # cc flag to enable OPENMP | |
| # | |
| OPENMP = -fopenmp | |
| # |
lue
/ plot_mod.py
Created
February 14, 2016 21:17
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
| import numpy as np | |
| import matplotlib.pyplot as plt | |
| import matplotlib.patches as mpatches | |
| dat=np.genfromtxt('dos_ocd_ioeu_eur_sr-qc2.hg19.ch.fl.all.out.dosage.gz.noqc.Indels-removed.5132ind.0.6info_0.8cert.pruned2.evec', dtype=None, skiprows=1) | |
| fam=np.genfromtxt('dos_ocd_ioeu_eur_sr-qc2.hg19.ch.fl.all.out.dosage.gz.noqc.Indels-removed.5132ind.0.6info_0.8cert.fam', dtype=None) | |
| sex=fam['f4'] | |
| stat=fam['f5'] |
lue
/ paper_plot.py
Last active
October 2, 2015 16:58
Difference in 21cm power spectrum with and without inner structure
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
| from prettyfigure.style import * | |
| define_figure_style() | |
| data = np.load('temp.npz') | |
| pkdata = data['pkdata'] | |
| pkdata_corrected = data['pkdata_corrected'] | |
| pkdata_other = data['pkdata_other'] | |
| pkdata_b = data['pkdata_b'] | |
| slice1 = data['slice1'] | |
| slice2 = data['slice2'] |