lue

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])

lue

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% |

lue

#
# Pick the MPI C compiler
CC = mpicc
#

#
# cc flag to enable OPENMP
#
OPENMP = -fopenmp
#

lue

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

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']