# Notebook from ventolab/HGDA
Path: scATACseq_females_3.ipynb

## ATAC + MULTIOME FEMALES (april 2021)_____no_output_____### Convert to Seurat objects for CICERO_____no_output_____

import numpy as np
import scanpy as sc
import pandas as pd
import anndata
/opt/conda/envs/atac_env/lib/python3.8/site-packages/joblib/_multiprocessing_helpers.py:45: UserWarning: [Errno 28] No space left on device.  joblib will operate in serial mode
  warnings.warn('%s.  joblib will operate in serial mode' % (e,))
import rpy2.rinterface_lib.callbacks
import logging_____no_output_____# Ignore R warning messages
#Note: this can be commented out to get more verbose R output
rpy2.rinterface_lib.callbacks.logger.setLevel(logging.ERROR)_____no_output_____import anndata2ri
anndata2ri.activate()
sc.logging.print_versions()WARNING: If you miss a compact list, please try `print_header`!
print(anndata2ri.__version__)
print(anndata.__version__)1.0.6
0.7.5
%load_ext rpy2.ipython_____no_output_____%%R

library(Signac)
library(Seurat)
library(SeuratDisk)
library(GenomeInfoDb)
library(EnsDb.Hsapiens.v86)
library(BSgenome.Hsapiens.UCSC.hg38)
library(ensembldb)
library(JASPAR2018)
library(TFBSTools)
library(patchwork)
library(universalmotif)
library(MotifDb)
library(TFBSTools)
library(patchwork)
library(chromVAR)
library(motifmatchr)_____no_output_____

For future reference most packages were installed here: */tmp/RtmpGjkK1o/downloaded_packages/*_____no_output_____**Set directories**_____no_output_____

# Define variables
outdir = "/nfs/team292/vl6/my_MULTIOME_dir/females_apr2021/"
experiment_prefix = 'females_'_____no_output_____

#### Convert the h5ad created with cisTopic into a Seurat object with peaks assay_____no_output_____

ATAC_scanpy = sc.read_h5ad(outdir + experiment_prefix + "_harmony_ATACwCisTopic.h5ad")
ATAC_scanpy_____no_output_____ATAC_scanpy.X_____no_output_____binary_raw = ATAC_scanpy.layers['binary_raw']
binary_raw_____no_output_____binary_raw.shape_____no_output_____

Step 1: from sparse matrix (python format) to dgcMatrix (R format)

I found this function here: <https://github.com/theislab/anndata2ri/issues/8>_____no_output_____

import numpy as np
from scipy import sparse

import rpy2.robjects as ro
from rpy2.robjects import pandas2ri, numpy2ri
from rpy2.robjects.conversion import localconverter

ro.r("library(Matrix)")

def csr_to_dgc(csr):
    """Convert (and transpose) a csr matrix from python to a R dgCMatrix (not sure if type is consistent)
    """
    numeric = ro.r("as.numeric")
    with localconverter(ro.default_converter + ro.numpy2ri.converter):
        X = ro.r("sparseMatrix")(
            i=numeric(csr.indices),
            p=numeric(csr.indptr),
            x=numeric(csr.data),
            dims=list(csr.shape[::-1]),
            index1=False
        )
    return X_____no_output_____ATAC_R = csr_to_dgc(binary_raw)
ATAC_R_____no_output_____%%R -i ATAC_R -o ATAC_dgT

ATAC_dgT <- as(ATAC_R, "dgTMatrix")
print(dim(ATAC_dgT)) # Check the shape of the matrix [1] 201343  67202

Step 2: save metadata per cell and per peak and add colnames and rownames to the peak matrix_____no_output_____

metadata_cell = ATAC_scanpy.obs
metadata_peak = ATAC_scanpy.var
layers = ATAC_scanpy.layers
reduction_cistopic = ATAC_scanpy.obsm['X_cistopic_40']
reduction_umap = ATAC_scanpy.obsm['X_umap']
reduction_harmony = ATAC_scanpy.obsm['X_harmony']_____no_output_____del ATAC_scanpy.obs['cellatac_clusters']_____no_output_____%%R -i ATAC_dgT -i metadata_cell -i metadata_peak

rownames(ATAC_dgT) <- rownames(metadata_peak) # Add rownames to peak matrix
colnames(ATAC_dgT) <- rownames(metadata_cell) # Add colnames to peak matrix _____no_output_____%%R
print(head(rownames(ATAC_dgT)))
print(head(colnames(ATAC_dgT)))
print(names(ATAC_dgT))[1] "chr1:816836-817441" "chr1:817738-818526" "chr1:825667-826063"
[4] "chr1:826596-827809" "chr1:858468-859047" "chr1:869401-870252"
[1] "01-AAACGAAAGGCTCCTG-1" "01-AAACGAAAGTCCAGAG-1" "01-AAACGAACACCCATCT-1"
[4] "01-AAACGAACAGGTCTGC-1" "01-AAACGAACATAGCCAT-1" "01-AAACGAAGTAGCATGG-1"
NULL

Step 3: build ChromatinAssay with peaks x cells matrix _____no_output_____

%%R -o chrom_assay

 chrom_assay <- CreateChromatinAssay(
    counts = ATAC_dgT,
    sep = c(":", "-"),
    genome = "hg38",
    ranges = NULL,
 )
chrom_assayChromatinAssay data with 201343 features for 67202 cells
Variable features: 0
Genome: hg38
Annotation present: FALSE
Motifs present: FALSE
Fragment files: 0

Step 4: build Seurat object from chromatin assay _____no_output_____

%%R -o ATAC_Seurat -i metadata_cell

ATAC_Seurat <- CreateSeuratObject(
    counts = chrom_assay,
    assay = "peaks",
    meta.data = metadata_cell,
  )
ATAC_SeuratAn object of class Seurat
201343 features across 67202 samples within 1 assay
Active assay: peaks (201343 features, 0 variable features)
%%R
print(table(ATAC_Seurat@meta.data$individual))
print(table(ATAC_Seurat@meta.data$leiden))
  F81 Hrv11 Hrv12 Hrv13 Hrv39 Hrv49 Hrv50 Hrv58 Hrv59
 8553  1163  3391  4674 10451  8586 20240  4912  5232

   0    1    2    3    4    5    6    7    8    9   10   11   12   13   14   15
7042 6881 6815 5322 5262 4591 3988 3595 3496 3383 3165 2856 2529 1598  895  840
  16   17   18   19   20   21   22   23
 818  806  778  732  615  517  420  258

Step 5: Add dimensionality reductions to Seurat object _____no_output_____

%%R -i reduction_cistopic -i metadata_cell -o reduction_cistopic_2

reduction_cistopic_2 <- data.matrix(reduction_cistopic)
key <- "cisTopic_"
colnames(reduction_cistopic_2) <- paste0(key, 1:ncol(x = reduction_cistopic_2))
rownames(reduction_cistopic_2) <- rownames(metadata_cell)
#print(head(reduction_cistopic_2))_____no_output_____%%R -i reduction_umap -i metadata_cell -o reduction_umap_2

reduction_umap_2 <- data.matrix(reduction_umap)
key2 <- "UMAP_"
colnames(reduction_umap_2) <- paste0(key2, 1:ncol(x = reduction_umap_2))
rownames(reduction_umap_2) <- rownames(metadata_cell)
#print(head(reduction_umap_2))_____no_output_____%%R -i reduction_harmony -i metadata_cell -o reduction_harmony_2

reduction_harmony_2 <- data.matrix(reduction_harmony)
key3 <- "Harmony_"
colnames(reduction_harmony_2) <- paste0(key3, 1:ncol(x = reduction_harmony_2))
rownames(reduction_harmony_2) <- rownames(metadata_cell)
#print(head(reduction_harmony_2))_____no_output_____%%R

ATAC_Seurat@reductions$cisTopics <- CreateDimReducObject(embeddings = reduction_cistopic_2, assay = "peaks")
ATAC_Seurat@reductions$umap <- CreateDimReducObject(embeddings = reduction_umap_2, assay = "peaks")
ATAC_Seurat@reductions$harmony <- CreateDimReducObject(embeddings = reduction_harmony_2, assay = "peaks")
ATAC_SeuratAn object of class Seurat
201343 features across 67202 samples within 1 assay
Active assay: peaks (201343 features, 0 variable features)
 3 dimensional reductions calculated: cisTopics, umap, harmony

### Save RDS object of females _____no_output_____

%%R -i outdir -i experiment_prefix

saveRDS(ATAC_Seurat, file = paste0(outdir, experiment_prefix, "_full.rds"))_____no_output_____%%R -i outdir -i experiment_prefix

ATAC_Seurat <- readRDS(file = paste0(outdir, experiment_prefix, "_full.rds"))_____no_output_____%%R

gonads_colors = c('#c466b0', '#ebd3c3',  '#edb7b7','#d96884','#e64e74','#e08b8b', '#d9abb7', '#dba781',
     '#d4db81','#779eed', '#7b9e99',
 '#70ccbe',   '#71a2c7','#91bd80',  '#bdb380','#60bddb',
    '#aad3f2')
DimPlot(object = ATAC_Seurat, group.by = "cell_type", label = TRUE, repel = TRUE,
        label.size = 5, cols = gonads_colors)_____no_output_____

## TF motifs with ChromVar_____no_output_____#### Add motifs_____no_output_____

%%R -o pfm

# Get a list of motif position frequency matrices from the JASPAR, HOCOMOCO, HOMER and SwissRegulon databases
pfm <- getMatrixSet(
  x = JASPAR2018,
  opts = list(species = 9606, all_versions = FALSE)
)

tf_motif_hoco=subset(MotifDb, organism=='Hsapiens' & dataSource=='HOCOMOCOv10')
pfmHoco = convert_motifs(tf_motif_hoco, 'TFBSTools-PWMatrix')
pfmHoco_names =  sapply(pfmHoco, function(x) x@ID)
names(tf_motif_hoco) = pfmHoco_names
for (i in 1:length(pfmHoco) ){
    # Convert into a PCM
    x = t(pfmHoco[[i]]@profileMatrix) - apply(pfmHoco[[i]]@profileMatrix, 2, min) # normalize per column
    pfmHoco[[i]]@profileMatrix = round(t(x / apply(x, 1, max))*100) # mulitply x100 and round
    class(pfmHoco[[i]]@profileMatrix) <- 'integer'
    pfm[[pfmHoco_names[i]]] <- pfmHoco[[i]]
}
# add data from SwissRegulon
tf_motif_SwissRegulon=subset(MotifDb, organism=='Hsapiens' & dataSource=='SwissRegulon')
pfmSwissRegulon = convert_motifs(tf_motif_SwissRegulon, 'TFBSTools-PWMatrix')
pfmSwissRegulon_names =  sapply(pfmSwissRegulon, function(x) x@ID)
names(tf_motif_SwissRegulon) = pfmSwissRegulon_names
for (i in 1:length(pfmSwissRegulon) ){
    # Convert into a PCM
    x = t(pfmSwissRegulon[[i]]@profileMatrix) - apply(pfmSwissRegulon[[i]]@profileMatrix, 2, min) # normalize per column
    pfmSwissRegulon[[i]]@profileMatrix = round(t(x / apply(x, 1, max))*100) # mulitply x100 and round
    class(pfmSwissRegulon[[i]]@profileMatrix) <- 'integer'
    pfm[[pfmSwissRegulon_names[i]]] <- pfmSwissRegulon[[i]]
}

# add data from HOMER
tf_motif_HOMER=subset(MotifDb, dataSource=='HOMER')
pfmHOMER = convert_motifs(tf_motif_HOMER, 'TFBSTools-PWMatrix')
pfmHOMER_names =  sapply(pfmHOMER, function(x) x@ID)
names(pfmHOMER) = paste0( sapply( strsplit(pfmHOMER_names, '\\('), head, 1), '.HOMER')
#print(names(pfmHOMER))
names(pfmHOMER)[93] = "FOXL2.HOMER"

pfmHOMER = pfmHOMER[ grep('[a-z]',  names(pfmHOMER), invert = T) ] # remove non-human (lowercase)
#print(names(pfmHOMER))

for (i in 1:length(pfmHOMER) ){
    # Convert into a PCM
    x = t(pfmHOMER[[i]]@profileMatrix) - apply(pfmHOMER[[i]]@profileMatrix, 2, min) # normalize per column
    pfmHOMER[[i]]@profileMatrix = round(t(x / apply(x, 1, max))*100) # mulitply x100 and round
    class(pfmHOMER[[i]]@profileMatrix) <- 'integer'
    pfmHOMER[[i]]@name = gsub('.HOMER', '', names(pfmHOMER)[i])
    pfm[[names(pfmHOMER)[i]]] <- pfmHOMER[[i]]
}

# Add name in id
idx =  grep('^MA', names(pfm) )
for (i in idx ){
    pfm[[i]]@ID = paste0(pfm[[i]]@name, '.', pfm[[i]]@ID)
    names(pfm)[i] = pfm[[i]]@ID
}

# remove pcm with NAs
pfm = pfm[ ! sapply(pfm, function(x) any(is.na(x@profileMatrix))) ]

length(pfm)
table(unlist(sapply(pfm, function(x) x@tags$dataSource)))[1] 638 147 683

#### Add motifs_____no_output_____

%%R -o ATAC_Seurat_withMotifs

ATAC_Seurat_withMotifs <- AddMotifs(object = ATAC_Seurat,
                         genome= BSgenome.Hsapiens.UCSC.hg38,
                         pfm = pfm,
                         assay = "peaks",
                         verbose = TRUE)

ATAC_Seurat_withMotifsAn object of class Seurat
201343 features across 67202 samples within 1 assay
Active assay: peaks (201343 features, 0 variable features)
 3 dimensional reductions calculated: cisTopics, umap, harmony

#### Compute per-cell motif activity scores with chromVar_____no_output_____

%%R -o ATAC_Seurat_withChromVar

library(BiocParallel)
register(SerialParam())
ATAC_Seurat_withChromVar <- RunChromVAR(
  object = ATAC_Seurat_withMotifs,
  genome = BSgenome.Hsapiens.UCSC.hg38,
  assay="peaks"
)
ATAC_Seurat_withChromVarAn object of class Seurat
203263 features across 67202 samples within 2 assays
Active assay: peaks (201343 features, 0 variable features)
 1 other assay present: chromvar
 3 dimensional reductions calculated: cisTopics, umap, harmony
%%R -i outdir -i experiment_prefix

saveRDS(ATAC_Seurat_withChromVar, file = paste0(outdir, experiment_prefix, "_chromVar_binary.rds"))_____no_output_____%%R -i pfm -o motif2names

library(dplyr)
motif2names <- lapply(pfm, function(x) x@name)
names(motif2names) <- sapply(pfm, function(x) x@ID) %>% gsub('_', '-', .)_____no_output_____%%R

DefaultAssay(ATAC_Seurat_withChromVar) <- "chromvar" # Change default assay _____no_output_____%%R

Idents(ATAC_Seurat_withChromVar) <- ATAC_Seurat_withChromVar@meta.data$cell_type
print(table(Idents(ATAC_Seurat_withChromVar)))
     Germ cells         coelEpi     ovarianSurf           sPAX8    preGC_I_OSR1
           1618             895            3496             615            6881
 preGC_II_KITLG  preGC_III_GJA1              Oi              Gi              PV
           7300            4591            6815            5262            2529
M_MullDuct_LGR5           M_MGP     M_prog_ISL1      Epithelial          Immune
           5322           14131             806            1598             818
    Endothelial          Neural
           3588             937
%%R -i outdir -o all_da_peaks

all_da_peaks = FindAllMarkers(ATAC_Seurat_withChromVar, test.use = 'LR', verbose = T, only.pos = F, random.seed = 1, logfc.threshold = 0, min.pct = 0, return.thresh = 1)
write.csv(all_da_peaks, file = paste0(outdir, 'TFs/DARs/all_DA_peaks_LR.csv'))  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~27m 28s       |++                                                | 2 % ~28m 24s       |++                                                | 3 % ~27m 04s       |+++                                               | 4 % ~25m 30s       |+++                                               | 5 % ~24m 59s       |++++                                              | 6 % ~23m 52s       |++++                                              | 7 % ~23m 36s       |+++++                                             | 8 % ~22m 59s       |+++++                                             | 9 % ~22m 34s       |++++++                                            | 10% ~22m 21s       |++++++                                            | 11% ~22m 01s       |+++++++                                           | 12% ~21m 43s       |+++++++                                           | 14% ~21m 18s       |++++++++                                          | 15% ~20m 54s       |++++++++                                          | 16% ~20m 47s       |+++++++++                                         | 17% ~20m 20s       |+++++++++                                         | 18% ~20m 09s       |++++++++++                                        | 19% ~19m 46s       |++++++++++                                        | 20% ~19m 26s       |+++++++++++                                       | 21% ~19m 09s       |+++++++++++                                       | 22% ~18m 57s       |++++++++++++                                      | 23% ~18m 38s       |++++++++++++                                      | 24% ~18m 24s       |+++++++++++++                                     | 25% ~18m 13s       |++++++++++++++                                    | 26% ~17m 59s       |++++++++++++++                                    | 27% ~17m 46s       |+++++++++++++++                                   | 28% ~17m 33s       |+++++++++++++++                                   | 29% ~17m 16s       |++++++++++++++++                                  | 30% ~16m 58s       |++++++++++++++++                                  | 31% ~16m 49s       |+++++++++++++++++                                 | 32% ~16m 35s       |+++++++++++++++++                                 | 33% ~16m 21s       |++++++++++++++++++                                | 34% ~15m 59s       |++++++++++++++++++                                | 35% ~15m 43s       |+++++++++++++++++++                               | 36% ~15m 27s       |+++++++++++++++++++                               | 38% ~15m 13s       |++++++++++++++++++++                              | 39% ~14m 58s       |++++++++++++++++++++                              | 40% ~14m 41s       |+++++++++++++++++++++                             | 41% ~14m 25s       |+++++++++++++++++++++                             | 42% ~14m 09s       |++++++++++++++++++++++                            | 43% ~13m 52s       |++++++++++++++++++++++                            | 44% ~13m 35s       |+++++++++++++++++++++++                           | 45% ~13m 17s       |+++++++++++++++++++++++                           | 46% ~13m 00s       |++++++++++++++++++++++++                          | 47% ~12m 44s       |++++++++++++++++++++++++                          | 48% ~12m 30s       |+++++++++++++++++++++++++                         | 49% ~12m 14s       |+++++++++++++++++++++++++                         | 50% ~11m 57s       |++++++++++++++++++++++++++                        | 51% ~11m 42s       |+++++++++++++++++++++++++++                       | 52% ~11m 26s       |+++++++++++++++++++++++++++                       | 53% ~11m 12s       |++++++++++++++++++++++++++++                      | 54% ~10m 58s       |++++++++++++++++++++++++++++                      | 55% ~10m 43s       |+++++++++++++++++++++++++++++                     | 56% ~10m 28s       |+++++++++++++++++++++++++++++                     | 57% ~10m 12s       |++++++++++++++++++++++++++++++                    | 58% ~09m 57s       |++++++++++++++++++++++++++++++                    | 59% ~09m 42s       |+++++++++++++++++++++++++++++++                   | 60% ~09m 27s       |+++++++++++++++++++++++++++++++                   | 61% ~09m 13s       |++++++++++++++++++++++++++++++++                  | 62% ~08m 57s       |++++++++++++++++++++++++++++++++                  | 64% ~08m 42s       |+++++++++++++++++++++++++++++++++                 | 65% ~08m 27s       |+++++++++++++++++++++++++++++++++                 | 66% ~08m 12s       |++++++++++++++++++++++++++++++++++                | 67% ~07m 58s       |++++++++++++++++++++++++++++++++++                | 68% ~07m 43s       |+++++++++++++++++++++++++++++++++++               | 69% ~07m 29s       |+++++++++++++++++++++++++++++++++++               | 70% ~07m 14s       |++++++++++++++++++++++++++++++++++++              | 71% ~06m 59s       |++++++++++++++++++++++++++++++++++++              | 72% ~06m 45s       |+++++++++++++++++++++++++++++++++++++             | 73% ~06m 29s       |+++++++++++++++++++++++++++++++++++++             | 74% ~06m 14s       |++++++++++++++++++++++++++++++++++++++            | 75% ~05m 59s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~05m 43s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~05m 29s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~05m 14s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~04m 59s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~04m 44s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~04m 30s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~04m 14s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 59s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~03m 44s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~03m 30s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~03m 15s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 59s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 44s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 29s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~02m 14s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 59s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 44s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 29s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 14s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~60s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~45s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~30s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~15s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=23m 48s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~25m 11s       |++                                                | 2 % ~24m 47s       |++                                                | 3 % ~24m 57s       |+++                                               | 4 % ~25m 02s       |+++                                               | 5 % ~25m 40s       |++++                                              | 6 % ~25m 34s       |++++                                              | 7 % ~25m 33s       |+++++                                             | 8 % ~25m 08s       |+++++                                             | 9 % ~24m 44s       |++++++                                            | 10% ~24m 15s       |++++++                                            | 11% ~23m 46s       |+++++++                                           | 12% ~23m 29s       |+++++++                                           | 14% ~23m 09s       |++++++++                                          | 15% ~22m 51s       |++++++++                                          | 16% ~22m 51s       |+++++++++                                         | 17% ~22m 37s       |+++++++++                                         | 18% ~22m 10s       |++++++++++                                        | 19% ~21m 49s       |++++++++++                                        | 20% ~21m 29s       |+++++++++++                                       | 21% ~21m 18s       |+++++++++++                                       | 22% ~21m 02s       |++++++++++++                                      | 23% ~20m 44s       |++++++++++++                                      | 24% ~20m 24s       |+++++++++++++                                     | 25% ~20m 06s       |++++++++++++++                                    | 26% ~19m 45s       |++++++++++++++                                    | 27% ~19m 25s       |+++++++++++++++                                   | 28% ~19m 10s       |+++++++++++++++                                   | 29% ~18m 57s       |++++++++++++++++                                  | 30% ~18m 38s       |++++++++++++++++                                  | 31% ~18m 18s       |+++++++++++++++++                                 | 32% ~18m 01s       |+++++++++++++++++                                 | 33% ~17m 46s       |++++++++++++++++++                                | 34% ~17m 36s       |++++++++++++++++++                                | 35% ~17m 16s       |+++++++++++++++++++                               | 36% ~16m 56s       |+++++++++++++++++++                               | 38% ~16m 45s       |++++++++++++++++++++                              | 39% ~16m 27s       |++++++++++++++++++++                              | 40% ~16m 10s       |+++++++++++++++++++++                             | 41% ~15m 55s       |+++++++++++++++++++++                             | 42% ~15m 39s       |++++++++++++++++++++++                            | 43% ~15m 21s       |++++++++++++++++++++++                            | 44% ~15m 04s       |+++++++++++++++++++++++                           | 45% ~14m 50s       |+++++++++++++++++++++++                           | 46% ~14m 33s       |++++++++++++++++++++++++                          | 47% ~14m 15s       |++++++++++++++++++++++++                          | 48% ~13m 57s       |+++++++++++++++++++++++++                         | 49% ~13m 39s       |+++++++++++++++++++++++++                         | 50% ~13m 22s       |++++++++++++++++++++++++++                        | 51% ~13m 08s       |+++++++++++++++++++++++++++                       | 52% ~12m 51s       |+++++++++++++++++++++++++++                       | 53% ~12m 33s       |++++++++++++++++++++++++++++                      | 54% ~12m 15s       |++++++++++++++++++++++++++++                      | 55% ~11m 59s       |+++++++++++++++++++++++++++++                     | 56% ~11m 45s       |+++++++++++++++++++++++++++++                     | 57% ~11m 30s       |++++++++++++++++++++++++++++++                    | 58% ~11m 15s       |++++++++++++++++++++++++++++++                    | 59% ~10m 59s       |+++++++++++++++++++++++++++++++                   | 60% ~10m 45s       |+++++++++++++++++++++++++++++++                   | 61% ~10m 29s       |++++++++++++++++++++++++++++++++                  | 62% ~10m 11s       |++++++++++++++++++++++++++++++++                  | 64% ~09m 54s       |+++++++++++++++++++++++++++++++++                 | 65% ~09m 36s       |+++++++++++++++++++++++++++++++++                 | 66% ~09m 19s       |++++++++++++++++++++++++++++++++++                | 67% ~09m 02s       |++++++++++++++++++++++++++++++++++                | 68% ~08m 45s       |+++++++++++++++++++++++++++++++++++               | 69% ~08m 28s       |+++++++++++++++++++++++++++++++++++               | 70% ~08m 12s       |++++++++++++++++++++++++++++++++++++              | 71% ~07m 55s       |++++++++++++++++++++++++++++++++++++              | 72% ~07m 37s       |+++++++++++++++++++++++++++++++++++++             | 73% ~07m 20s       |+++++++++++++++++++++++++++++++++++++             | 74% ~07m 04s       |++++++++++++++++++++++++++++++++++++++            | 75% ~06m 46s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~06m 29s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~06m 12s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~05m 54s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~05m 37s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~05m 20s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~05m 04s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~04m 47s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~04m 30s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~04m 13s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~03m 56s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~03m 39s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~03m 23s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~03m 06s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 49s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~02m 32s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~02m 15s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 58s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 41s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 24s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~01m 07s       |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~50s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~34s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~17s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=26m 49s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~20m 33s       |++                                                | 2 % ~20m 46s       |++                                                | 3 % ~20m 07s       |+++                                               | 4 % ~19m 49s       |+++                                               | 5 % ~19m 45s       |++++                                              | 6 % ~19m 44s       |++++                                              | 7 % ~19m 56s       |+++++                                             | 8 % ~20m 07s       |+++++                                             | 9 % ~20m 02s       |++++++                                            | 10% ~19m 55s       |++++++                                            | 11% ~19m 36s       |+++++++                                           | 12% ~19m 23s       |+++++++                                           | 14% ~19m 01s       |++++++++                                          | 15% ~18m 46s       |++++++++                                          | 16% ~18m 29s       |+++++++++                                         | 17% ~18m 00s       |+++++++++                                         | 18% ~17m 47s       |++++++++++                                        | 19% ~17m 32s       |++++++++++                                        | 20% ~17m 34s       |+++++++++++                                       | 21% ~17m 21s       |+++++++++++                                       | 22% ~17m 05s       |++++++++++++                                      | 23% ~16m 52s       |++++++++++++                                      | 24% ~16m 40s       |+++++++++++++                                     | 25% ~16m 27s       |++++++++++++++                                    | 26% ~16m 12s       |++++++++++++++                                    | 27% ~15m 59s       |+++++++++++++++                                   | 28% ~15m 41s       |+++++++++++++++                                   | 29% ~15m 22s       |++++++++++++++++                                  | 30% ~15m 09s       |++++++++++++++++                                  | 31% ~14m 51s       |+++++++++++++++++                                 | 32% ~14m 37s       |+++++++++++++++++                                 | 33% ~14m 29s       |++++++++++++++++++                                | 34% ~14m 15s       |++++++++++++++++++                                | 35% ~13m 57s       |+++++++++++++++++++                               | 36% ~13m 41s       |+++++++++++++++++++                               | 38% ~13m 26s       |++++++++++++++++++++                              | 39% ~13m 12s       |++++++++++++++++++++                              | 40% ~13m 00s       |+++++++++++++++++++++                             | 41% ~12m 48s       |+++++++++++++++++++++                             | 42% ~12m 34s       |++++++++++++++++++++++                            | 43% ~12m 22s       |++++++++++++++++++++++                            | 44% ~12m 09s       |+++++++++++++++++++++++                           | 45% ~11m 57s       |+++++++++++++++++++++++                           | 46% ~11m 45s       |++++++++++++++++++++++++                          | 47% ~11m 30s       |++++++++++++++++++++++++                          | 48% ~11m 17s       |+++++++++++++++++++++++++                         | 49% ~11m 06s       |+++++++++++++++++++++++++                         | 50% ~10m 54s       |++++++++++++++++++++++++++                        | 51% ~10m 41s       |+++++++++++++++++++++++++++                       | 52% ~10m 28s       |+++++++++++++++++++++++++++                       | 53% ~10m 14s       |++++++++++++++++++++++++++++                      | 54% ~09m 60s       |++++++++++++++++++++++++++++                      | 55% ~09m 47s       |+++++++++++++++++++++++++++++                     | 56% ~09m 33s       |+++++++++++++++++++++++++++++                     | 57% ~09m 18s       |++++++++++++++++++++++++++++++                    | 58% ~09m 05s       |++++++++++++++++++++++++++++++                    | 59% ~08m 50s       |+++++++++++++++++++++++++++++++                   | 60% ~08m 37s       |+++++++++++++++++++++++++++++++                   | 61% ~08m 24s       |++++++++++++++++++++++++++++++++                  | 62% ~08m 10s       |++++++++++++++++++++++++++++++++                  | 64% ~07m 57s       |+++++++++++++++++++++++++++++++++                 | 65% ~07m 43s       |+++++++++++++++++++++++++++++++++                 | 66% ~07m 30s       |++++++++++++++++++++++++++++++++++                | 67% ~07m 16s       |++++++++++++++++++++++++++++++++++                | 68% ~07m 02s       |+++++++++++++++++++++++++++++++++++               | 69% ~06m 49s       |+++++++++++++++++++++++++++++++++++               | 70% ~06m 36s       |++++++++++++++++++++++++++++++++++++              | 71% ~06m 22s       |++++++++++++++++++++++++++++++++++++              | 72% ~06m 09s       |+++++++++++++++++++++++++++++++++++++             | 73% ~05m 55s       |+++++++++++++++++++++++++++++++++++++             | 74% ~05m 41s       |++++++++++++++++++++++++++++++++++++++            | 75% ~05m 28s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~05m 14s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~05m 00s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~04m 46s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~04m 33s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~04m 19s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~04m 05s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 51s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 37s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~03m 23s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~03m 09s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 55s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 42s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 28s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 15s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~02m 01s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 47s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 34s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 21s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 07s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~54s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~40s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~27s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~13s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=21m 27s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~25m 21s       |++                                                | 2 % ~26m 27s       |++                                                | 3 % ~26m 25s       |+++                                               | 4 % ~25m 57s       |+++                                               | 5 % ~25m 02s       |++++                                              | 6 % ~24m 43s       |++++                                              | 7 % ~24m 27s       |+++++                                             | 8 % ~24m 06s       |+++++                                             | 9 % ~23m 51s       |++++++                                            | 10% ~23m 33s       |++++++                                            | 11% ~23m 16s       |+++++++                                           | 12% ~23m 01s       |+++++++                                           | 14% ~22m 47s       |++++++++                                          | 15% ~22m 40s       |++++++++                                          | 16% ~22m 32s       |+++++++++                                         | 17% ~22m 14s       |+++++++++                                         | 18% ~22m 00s       |++++++++++                                        | 19% ~21m 45s       |++++++++++                                        | 20% ~21m 31s       |+++++++++++                                       | 21% ~21m 18s       |+++++++++++                                       | 22% ~21m 04s       |++++++++++++                                      | 23% ~20m 45s       |++++++++++++                                      | 24% ~20m 30s       |+++++++++++++                                     | 25% ~20m 09s       |++++++++++++++                                    | 26% ~19m 52s       |++++++++++++++                                    | 27% ~19m 35s       |+++++++++++++++                                   | 28% ~19m 16s       |+++++++++++++++                                   | 29% ~19m 01s       |++++++++++++++++                                  | 30% ~18m 42s       |++++++++++++++++                                  | 31% ~18m 23s       |+++++++++++++++++                                 | 32% ~18m 07s       |+++++++++++++++++                                 | 33% ~17m 53s       |++++++++++++++++++                                | 34% ~17m 38s       |++++++++++++++++++                                | 35% ~17m 21s       |+++++++++++++++++++                               | 36% ~17m 03s       |+++++++++++++++++++                               | 38% ~16m 48s       |++++++++++++++++++++                              | 39% ~16m 31s       |++++++++++++++++++++                              | 40% ~16m 14s       |+++++++++++++++++++++                             | 41% ~15m 58s       |+++++++++++++++++++++                             | 42% ~15m 43s       |++++++++++++++++++++++                            | 43% ~15m 26s       |++++++++++++++++++++++                            | 44% ~15m 08s       |+++++++++++++++++++++++                           | 45% ~14m 50s       |+++++++++++++++++++++++                           | 46% ~14m 32s       |++++++++++++++++++++++++                          | 47% ~14m 14s       |++++++++++++++++++++++++                          | 48% ~13m 58s       |+++++++++++++++++++++++++                         | 49% ~13m 43s       |+++++++++++++++++++++++++                         | 50% ~13m 29s       |++++++++++++++++++++++++++                        | 51% ~13m 13s       |+++++++++++++++++++++++++++                       | 52% ~12m 56s       |+++++++++++++++++++++++++++                       | 53% ~12m 39s       |++++++++++++++++++++++++++++                      | 54% ~12m 22s       |++++++++++++++++++++++++++++                      | 55% ~12m 05s       |+++++++++++++++++++++++++++++                     | 56% ~11m 47s       |+++++++++++++++++++++++++++++                     | 57% ~11m 32s       |++++++++++++++++++++++++++++++                    | 58% ~11m 16s       |++++++++++++++++++++++++++++++                    | 59% ~10m 58s       |+++++++++++++++++++++++++++++++                   | 60% ~10m 40s       |+++++++++++++++++++++++++++++++                   | 61% ~10m 24s       |++++++++++++++++++++++++++++++++                  | 62% ~10m 07s       |++++++++++++++++++++++++++++++++                  | 64% ~09m 50s       |+++++++++++++++++++++++++++++++++                 | 65% ~09m 33s       |+++++++++++++++++++++++++++++++++                 | 66% ~09m 16s       |++++++++++++++++++++++++++++++++++                | 67% ~08m 59s       |++++++++++++++++++++++++++++++++++                | 68% ~08m 42s       |+++++++++++++++++++++++++++++++++++               | 69% ~08m 26s       |+++++++++++++++++++++++++++++++++++               | 70% ~08m 09s       |++++++++++++++++++++++++++++++++++++              | 71% ~07m 52s       |++++++++++++++++++++++++++++++++++++              | 72% ~07m 35s       |+++++++++++++++++++++++++++++++++++++             | 73% ~07m 18s       |+++++++++++++++++++++++++++++++++++++             | 74% ~07m 02s       |++++++++++++++++++++++++++++++++++++++            | 75% ~06m 46s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~06m 28s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~06m 12s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~05m 55s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~05m 37s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~05m 21s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~05m 04s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~04m 47s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~04m 30s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~04m 13s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~03m 56s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~03m 39s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~03m 23s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~03m 06s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 49s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~02m 33s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~02m 16s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 59s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 42s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 25s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~01m 08s       |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~51s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~34s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~17s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=27m 07s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~17m 59s       |++                                                | 2 % ~17m 41s       |++                                                | 3 % ~17m 22s       |+++                                               | 4 % ~17m 11s       |+++                                               | 5 % ~17m 57s       |++++                                              | 6 % ~17m 51s       |++++                                              | 7 % ~17m 35s       |+++++                                             | 8 % ~17m 15s       |+++++                                             | 9 % ~17m 33s       |++++++                                            | 10% ~17m 13s       |++++++                                            | 11% ~16m 57s       |+++++++                                           | 12% ~16m 56s       |+++++++                                           | 14% ~16m 38s       |++++++++                                          | 15% ~16m 22s       |++++++++                                          | 16% ~16m 03s       |+++++++++                                         | 17% ~15m 51s       |+++++++++                                         | 18% ~15m 37s       |++++++++++                                        | 19% ~15m 23s       |++++++++++                                        | 20% ~15m 10s       |+++++++++++                                       | 21% ~14m 56s       |+++++++++++                                       | 22% ~14m 46s       |++++++++++++                                      | 23% ~14m 38s       |++++++++++++                                      | 24% ~14m 23s       |+++++++++++++                                     | 25% ~14m 09s       |++++++++++++++                                    | 26% ~13m 56s       |++++++++++++++                                    | 27% ~13m 45s       |+++++++++++++++                                   | 28% ~13m 37s       |+++++++++++++++                                   | 29% ~13m 25s       |++++++++++++++++                                  | 30% ~13m 11s       |++++++++++++++++                                  | 31% ~12m 60s       |+++++++++++++++++                                 | 32% ~12m 48s       |+++++++++++++++++                                 | 33% ~12m 38s       |++++++++++++++++++                                | 34% ~12m 23s       |++++++++++++++++++                                | 35% ~12m 10s       |+++++++++++++++++++                               | 36% ~11m 58s       |+++++++++++++++++++                               | 38% ~11m 47s       |++++++++++++++++++++                              | 39% ~11m 35s       |++++++++++++++++++++                              | 40% ~11m 23s       |+++++++++++++++++++++                             | 41% ~11m 12s       |+++++++++++++++++++++                             | 42% ~10m 58s       |++++++++++++++++++++++                            | 43% ~10m 46s       |++++++++++++++++++++++                            | 44% ~10m 34s       |+++++++++++++++++++++++                           | 45% ~10m 19s       |+++++++++++++++++++++++                           | 46% ~10m 09s       |++++++++++++++++++++++++                          | 47% ~09m 56s       |++++++++++++++++++++++++                          | 48% ~09m 44s       |+++++++++++++++++++++++++                         | 49% ~09m 32s       |+++++++++++++++++++++++++                         | 50% ~09m 20s       |++++++++++++++++++++++++++                        | 51% ~09m 07s       |+++++++++++++++++++++++++++                       | 52% ~08m 54s       |+++++++++++++++++++++++++++                       | 53% ~08m 42s       |++++++++++++++++++++++++++++                      | 54% ~08m 30s       |++++++++++++++++++++++++++++                      | 55% ~08m 19s       |+++++++++++++++++++++++++++++                     | 56% ~08m 08s       |+++++++++++++++++++++++++++++                     | 57% ~07m 56s       |++++++++++++++++++++++++++++++                    | 58% ~07m 44s       |++++++++++++++++++++++++++++++                    | 59% ~07m 32s       |+++++++++++++++++++++++++++++++                   | 60% ~07m 20s       |+++++++++++++++++++++++++++++++                   | 61% ~07m 08s       |++++++++++++++++++++++++++++++++                  | 62% ~06m 56s       |++++++++++++++++++++++++++++++++                  | 64% ~06m 45s       |+++++++++++++++++++++++++++++++++                 | 65% ~06m 33s       |+++++++++++++++++++++++++++++++++                 | 66% ~06m 21s       |++++++++++++++++++++++++++++++++++                | 67% ~06m 08s       |++++++++++++++++++++++++++++++++++                | 68% ~05m 58s       |+++++++++++++++++++++++++++++++++++               | 69% ~05m 47s       |+++++++++++++++++++++++++++++++++++               | 70% ~05m 34s       |++++++++++++++++++++++++++++++++++++              | 71% ~05m 23s       |++++++++++++++++++++++++++++++++++++              | 72% ~05m 10s       |+++++++++++++++++++++++++++++++++++++             | 73% ~04m 59s       |+++++++++++++++++++++++++++++++++++++             | 74% ~04m 47s       |++++++++++++++++++++++++++++++++++++++            | 75% ~04m 36s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~04m 24s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~04m 13s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~04m 01s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~03m 49s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~03m 38s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~03m 26s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 14s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 02s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~02m 51s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~02m 40s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 28s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 17s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 06s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~01m 55s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~01m 43s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 32s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 20s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 09s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~57s           |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~46s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~34s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~23s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~11s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=18m 19s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~18m 24s       |++                                                | 2 % ~18m 09s       |++                                                | 3 % ~18m 36s       |+++                                               | 4 % ~18m 37s       |+++                                               | 5 % ~18m 29s       |++++                                              | 6 % ~18m 12s       |++++                                              | 7 % ~18m 03s       |+++++                                             | 8 % ~17m 48s       |+++++                                             | 9 % ~17m 39s       |++++++                                            | 10% ~17m 20s       |++++++                                            | 11% ~16m 58s       |+++++++                                           | 12% ~16m 50s       |+++++++                                           | 14% ~16m 37s       |++++++++                                          | 15% ~16m 19s       |++++++++                                          | 16% ~16m 14s       |+++++++++                                         | 17% ~15m 59s       |+++++++++                                         | 18% ~15m 49s       |++++++++++                                        | 19% ~15m 36s       |++++++++++                                        | 20% ~15m 22s       |+++++++++++                                       | 21% ~15m 07s       |+++++++++++                                       | 22% ~14m 53s       |++++++++++++                                      | 23% ~14m 39s       |++++++++++++                                      | 24% ~14m 25s       |+++++++++++++                                     | 25% ~14m 12s       |++++++++++++++                                    | 26% ~14m 00s       |++++++++++++++                                    | 27% ~13m 47s       |+++++++++++++++                                   | 28% ~13m 36s       |+++++++++++++++                                   | 29% ~13m 23s       |++++++++++++++++                                  | 30% ~13m 11s       |++++++++++++++++                                  | 31% ~12m 59s       |+++++++++++++++++                                 | 32% ~12m 43s       |+++++++++++++++++                                 | 33% ~12m 30s       |++++++++++++++++++                                | 34% ~12m 17s       |++++++++++++++++++                                | 35% ~12m 07s       |+++++++++++++++++++                               | 36% ~11m 54s       |+++++++++++++++++++                               | 38% ~11m 42s       |++++++++++++++++++++                              | 39% ~11m 31s       |++++++++++++++++++++                              | 40% ~11m 19s       |+++++++++++++++++++++                             | 41% ~11m 10s       |+++++++++++++++++++++                             | 42% ~10m 58s       |++++++++++++++++++++++                            | 43% ~10m 47s       |++++++++++++++++++++++                            | 44% ~10m 35s       |+++++++++++++++++++++++                           | 45% ~10m 24s       |+++++++++++++++++++++++                           | 46% ~10m 13s       |++++++++++++++++++++++++                          | 47% ~10m 01s       |++++++++++++++++++++++++                          | 48% ~09m 49s       |+++++++++++++++++++++++++                         | 49% ~09m 37s       |+++++++++++++++++++++++++                         | 50% ~09m 26s       |++++++++++++++++++++++++++                        | 51% ~09m 13s       |+++++++++++++++++++++++++++                       | 52% ~09m 01s       |+++++++++++++++++++++++++++                       | 53% ~08m 49s       |++++++++++++++++++++++++++++                      | 54% ~08m 39s       |++++++++++++++++++++++++++++                      | 55% ~08m 26s       |+++++++++++++++++++++++++++++                     | 56% ~08m 14s       |+++++++++++++++++++++++++++++                     | 57% ~08m 01s       |++++++++++++++++++++++++++++++                    | 58% ~07m 49s       |++++++++++++++++++++++++++++++                    | 59% ~07m 39s       |+++++++++++++++++++++++++++++++                   | 60% ~07m 27s       |+++++++++++++++++++++++++++++++                   | 61% ~07m 15s       |++++++++++++++++++++++++++++++++                  | 62% ~07m 03s       |++++++++++++++++++++++++++++++++                  | 64% ~06m 52s       |+++++++++++++++++++++++++++++++++                 | 65% ~06m 41s       |+++++++++++++++++++++++++++++++++                 | 66% ~06m 29s       |++++++++++++++++++++++++++++++++++                | 67% ~06m 18s       |++++++++++++++++++++++++++++++++++                | 68% ~06m 05s       |+++++++++++++++++++++++++++++++++++               | 69% ~05m 53s       |+++++++++++++++++++++++++++++++++++               | 70% ~05m 42s       |++++++++++++++++++++++++++++++++++++              | 71% ~05m 30s       |++++++++++++++++++++++++++++++++++++              | 72% ~05m 17s       |+++++++++++++++++++++++++++++++++++++             | 73% ~05m 05s       |+++++++++++++++++++++++++++++++++++++             | 74% ~04m 53s       |++++++++++++++++++++++++++++++++++++++            | 75% ~04m 41s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~04m 29s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~04m 18s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~04m 07s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~03m 55s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~03m 43s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~03m 31s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 19s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 07s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~02m 55s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~02m 43s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 32s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 20s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 09s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~01m 57s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~01m 45s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 33s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 21s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 10s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~58s           |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~46s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~35s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~23s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~12s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=18m 34s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~21m 05s       |++                                                | 2 % ~20m 23s       |++                                                | 3 % ~20m 20s       |+++                                               | 4 % ~20m 15s       |+++                                               | 5 % ~20m 09s       |++++                                              | 6 % ~19m 36s       |++++                                              | 7 % ~19m 11s       |+++++                                             | 8 % ~19m 06s       |+++++                                             | 9 % ~18m 59s       |++++++                                            | 10% ~18m 39s       |++++++                                            | 11% ~18m 34s       |+++++++                                           | 12% ~18m 19s       |+++++++                                           | 14% ~18m 05s       |++++++++                                          | 15% ~17m 57s       |++++++++                                          | 16% ~17m 42s       |+++++++++                                         | 17% ~17m 28s       |+++++++++                                         | 18% ~17m 14s       |++++++++++                                        | 19% ~16m 60s       |++++++++++                                        | 20% ~16m 48s       |+++++++++++                                       | 21% ~16m 31s       |+++++++++++                                       | 22% ~16m 17s       |++++++++++++                                      | 23% ~15m 59s       |++++++++++++                                      | 24% ~15m 42s       |+++++++++++++                                     | 25% ~15m 24s       |++++++++++++++                                    | 26% ~15m 08s       |++++++++++++++                                    | 27% ~14m 52s       |+++++++++++++++                                   | 28% ~14m 39s       |+++++++++++++++                                   | 29% ~14m 28s       |++++++++++++++++                                  | 30% ~14m 15s       |++++++++++++++++                                  | 31% ~14m 06s       |+++++++++++++++++                                 | 32% ~13m 54s       |+++++++++++++++++                                 | 33% ~13m 41s       |++++++++++++++++++                                | 34% ~13m 26s       |++++++++++++++++++                                | 35% ~13m 13s       |+++++++++++++++++++                               | 36% ~12m 58s       |+++++++++++++++++++                               | 38% ~12m 42s       |++++++++++++++++++++                              | 39% ~12m 29s       |++++++++++++++++++++                              | 40% ~12m 15s       |+++++++++++++++++++++                             | 41% ~12m 00s       |+++++++++++++++++++++                             | 42% ~11m 47s       |++++++++++++++++++++++                            | 43% ~11m 33s       |++++++++++++++++++++++                            | 44% ~11m 22s       |+++++++++++++++++++++++                           | 45% ~11m 09s       |+++++++++++++++++++++++                           | 46% ~10m 57s       |++++++++++++++++++++++++                          | 47% ~10m 46s       |++++++++++++++++++++++++                          | 48% ~10m 34s       |+++++++++++++++++++++++++                         | 49% ~10m 22s       |+++++++++++++++++++++++++                         | 50% ~10m 08s       |++++++++++++++++++++++++++                        | 51% ~09m 55s       |+++++++++++++++++++++++++++                       | 52% ~09m 43s       |+++++++++++++++++++++++++++                       | 53% ~09m 29s       |++++++++++++++++++++++++++++                      | 54% ~09m 15s       |++++++++++++++++++++++++++++                      | 55% ~09m 02s       |+++++++++++++++++++++++++++++                     | 56% ~08m 50s       |+++++++++++++++++++++++++++++                     | 57% ~08m 36s       |++++++++++++++++++++++++++++++                    | 58% ~08m 22s       |++++++++++++++++++++++++++++++                    | 59% ~08m 09s       |+++++++++++++++++++++++++++++++                   | 60% ~07m 55s       |+++++++++++++++++++++++++++++++                   | 61% ~07m 43s       |++++++++++++++++++++++++++++++++                  | 62% ~07m 29s       |++++++++++++++++++++++++++++++++                  | 64% ~07m 16s       |+++++++++++++++++++++++++++++++++                 | 65% ~07m 04s       |+++++++++++++++++++++++++++++++++                 | 66% ~06m 52s       |++++++++++++++++++++++++++++++++++                | 67% ~06m 39s       |++++++++++++++++++++++++++++++++++                | 68% ~06m 27s       |+++++++++++++++++++++++++++++++++++               | 69% ~06m 14s       |+++++++++++++++++++++++++++++++++++               | 70% ~06m 02s       |++++++++++++++++++++++++++++++++++++              | 71% ~05m 49s       |++++++++++++++++++++++++++++++++++++              | 72% ~05m 37s       |+++++++++++++++++++++++++++++++++++++             | 73% ~05m 24s       |+++++++++++++++++++++++++++++++++++++             | 74% ~05m 11s       |++++++++++++++++++++++++++++++++++++++            | 75% ~04m 58s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~04m 45s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~04m 32s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~04m 20s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~04m 07s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~03m 55s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~03m 42s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 30s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 17s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~03m 05s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~02m 52s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 40s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 27s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 15s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 03s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~01m 50s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 38s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 26s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 14s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 01s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~49s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~37s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~25s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~12s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=19m 40s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~19m 16s       |++                                                | 2 % ~18m 23s       |++                                                | 3 % ~18m 54s       |+++                                               | 4 % ~18m 15s       |+++                                               | 5 % ~17m 59s       |++++                                              | 6 % ~17m 26s       |++++                                              | 7 % ~17m 12s       |+++++                                             | 8 % ~16m 46s       |+++++                                             | 9 % ~16m 36s       |++++++                                            | 10% ~16m 52s       |++++++                                            | 11% ~16m 30s       |+++++++                                           | 12% ~16m 08s       |+++++++                                           | 14% ~15m 53s       |++++++++                                          | 15% ~15m 45s       |++++++++                                          | 16% ~15m 31s       |+++++++++                                         | 17% ~15m 14s       |+++++++++                                         | 18% ~15m 09s       |++++++++++                                        | 19% ~14m 52s       |++++++++++                                        | 20% ~14m 36s       |+++++++++++                                       | 21% ~14m 20s       |+++++++++++                                       | 22% ~14m 10s       |++++++++++++                                      | 23% ~13m 58s       |++++++++++++                                      | 24% ~13m 45s       |+++++++++++++                                     | 25% ~13m 31s       |++++++++++++++                                    | 26% ~13m 19s       |++++++++++++++                                    | 27% ~13m 05s       |+++++++++++++++                                   | 28% ~12m 49s       |+++++++++++++++                                   | 29% ~12m 34s       |++++++++++++++++                                  | 30% ~12m 26s       |++++++++++++++++                                  | 31% ~12m 19s       |+++++++++++++++++                                 | 32% ~12m 06s       |+++++++++++++++++                                 | 33% ~11m 55s       |++++++++++++++++++                                | 34% ~11m 41s       |++++++++++++++++++                                | 35% ~11m 27s       |+++++++++++++++++++                               | 36% ~11m 14s       |+++++++++++++++++++                               | 38% ~11m 04s       |++++++++++++++++++++                              | 39% ~10m 50s       |++++++++++++++++++++                              | 40% ~10m 39s       |+++++++++++++++++++++                             | 41% ~10m 30s       |+++++++++++++++++++++                             | 42% ~10m 19s       |++++++++++++++++++++++                            | 43% ~10m 11s       |++++++++++++++++++++++                            | 44% ~10m 00s       |+++++++++++++++++++++++                           | 45% ~09m 48s       |+++++++++++++++++++++++                           | 46% ~09m 35s       |++++++++++++++++++++++++                          | 47% ~09m 21s       |++++++++++++++++++++++++                          | 48% ~09m 08s       |+++++++++++++++++++++++++                         | 49% ~08m 58s       |+++++++++++++++++++++++++                         | 50% ~08m 47s       |++++++++++++++++++++++++++                        | 51% ~08m 36s       |+++++++++++++++++++++++++++                       | 52% ~08m 24s       |+++++++++++++++++++++++++++                       | 53% ~08m 14s       |++++++++++++++++++++++++++++                      | 54% ~08m 04s       |++++++++++++++++++++++++++++                      | 55% ~07m 53s       |+++++++++++++++++++++++++++++                     | 56% ~07m 42s       |+++++++++++++++++++++++++++++                     | 57% ~07m 30s       |++++++++++++++++++++++++++++++                    | 58% ~07m 18s       |++++++++++++++++++++++++++++++                    | 59% ~07m 06s       |+++++++++++++++++++++++++++++++                   | 60% ~06m 55s       |+++++++++++++++++++++++++++++++                   | 61% ~06m 44s       |++++++++++++++++++++++++++++++++                  | 62% ~06m 33s       |++++++++++++++++++++++++++++++++                  | 64% ~06m 21s       |+++++++++++++++++++++++++++++++++                 | 65% ~06m 09s       |+++++++++++++++++++++++++++++++++                 | 66% ~05m 60s       |++++++++++++++++++++++++++++++++++                | 67% ~05m 48s       |++++++++++++++++++++++++++++++++++                | 68% ~05m 37s       |+++++++++++++++++++++++++++++++++++               | 69% ~05m 25s       |+++++++++++++++++++++++++++++++++++               | 70% ~05m 14s       |++++++++++++++++++++++++++++++++++++              | 71% ~05m 03s       |++++++++++++++++++++++++++++++++++++              | 72% ~04m 53s       |+++++++++++++++++++++++++++++++++++++             | 73% ~04m 42s       |+++++++++++++++++++++++++++++++++++++             | 74% ~04m 31s       |++++++++++++++++++++++++++++++++++++++            | 75% ~04m 20s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~04m 09s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~03m 59s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~03m 49s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~03m 37s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~03m 26s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~03m 15s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 05s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~02m 54s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~02m 43s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~02m 32s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 20s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 10s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~01m 59s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~01m 48s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~01m 37s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 26s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 16s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 05s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~54s           |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~43s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~32s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~22s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~11s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=17m 12s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~18m 15s       |++                                                | 2 % ~16m 59s       |++                                                | 3 % ~16m 45s       |+++                                               | 4 % ~16m 23s       |+++                                               | 5 % ~17m 40s       |++++                                              | 6 % ~17m 35s       |++++                                              | 7 % ~17m 45s       |+++++                                             | 8 % ~17m 20s       |+++++                                             | 9 % ~17m 06s       |++++++                                            | 10% ~17m 07s       |++++++                                            | 11% ~17m 08s       |+++++++                                           | 12% ~17m 00s       |+++++++                                           | 14% ~16m 45s       |++++++++                                          | 15% ~16m 29s       |++++++++                                          | 16% ~16m 14s       |+++++++++                                         | 17% ~16m 02s       |+++++++++                                         | 18% ~15m 51s       |++++++++++                                        | 19% ~15m 45s       |++++++++++                                        | 20% ~15m 35s       |+++++++++++                                       | 21% ~15m 30s       |+++++++++++                                       | 22% ~15m 15s       |++++++++++++                                      | 23% ~14m 57s       |++++++++++++                                      | 24% ~14m 52s       |+++++++++++++                                     | 25% ~14m 37s       |++++++++++++++                                    | 26% ~14m 20s       |++++++++++++++                                    | 27% ~14m 05s       |+++++++++++++++                                   | 28% ~13m 51s       |+++++++++++++++                                   | 29% ~13m 36s       |++++++++++++++++                                  | 30% ~13m 22s       |++++++++++++++++                                  | 31% ~13m 06s       |+++++++++++++++++                                 | 32% ~12m 51s       |+++++++++++++++++                                 | 33% ~12m 37s       |++++++++++++++++++                                | 34% ~12m 25s       |++++++++++++++++++                                | 35% ~12m 13s       |+++++++++++++++++++                               | 36% ~12m 05s       |+++++++++++++++++++                               | 38% ~11m 51s       |++++++++++++++++++++                              | 39% ~11m 36s       |++++++++++++++++++++                              | 40% ~11m 26s       |+++++++++++++++++++++                             | 41% ~11m 14s       |+++++++++++++++++++++                             | 42% ~11m 01s       |++++++++++++++++++++++                            | 43% ~10m 47s       |++++++++++++++++++++++                            | 44% ~10m 34s       |+++++++++++++++++++++++                           | 45% ~10m 21s       |+++++++++++++++++++++++                           | 46% ~10m 08s       |++++++++++++++++++++++++                          | 47% ~09m 57s       |++++++++++++++++++++++++                          | 48% ~09m 45s       |+++++++++++++++++++++++++                         | 49% ~09m 32s       |+++++++++++++++++++++++++                         | 50% ~09m 23s       |++++++++++++++++++++++++++                        | 51% ~09m 12s       |+++++++++++++++++++++++++++                       | 52% ~08m 59s       |+++++++++++++++++++++++++++                       | 53% ~08m 47s       |++++++++++++++++++++++++++++                      | 54% ~08m 35s       |++++++++++++++++++++++++++++                      | 55% ~08m 22s       |+++++++++++++++++++++++++++++                     | 56% ~08m 09s       |+++++++++++++++++++++++++++++                     | 57% ~07m 56s       |++++++++++++++++++++++++++++++                    | 58% ~07m 45s       |++++++++++++++++++++++++++++++                    | 59% ~07m 32s       |+++++++++++++++++++++++++++++++                   | 60% ~07m 19s       |+++++++++++++++++++++++++++++++                   | 61% ~07m 09s       |++++++++++++++++++++++++++++++++                  | 62% ~06m 57s       |++++++++++++++++++++++++++++++++                  | 64% ~06m 47s       |+++++++++++++++++++++++++++++++++                 | 65% ~06m 36s       |+++++++++++++++++++++++++++++++++                 | 66% ~06m 27s       |++++++++++++++++++++++++++++++++++                | 67% ~06m 16s       |++++++++++++++++++++++++++++++++++                | 68% ~06m 05s       |+++++++++++++++++++++++++++++++++++               | 69% ~05m 53s       |+++++++++++++++++++++++++++++++++++               | 70% ~05m 42s       |++++++++++++++++++++++++++++++++++++              | 71% ~05m 31s       |++++++++++++++++++++++++++++++++++++              | 72% ~05m 20s       |+++++++++++++++++++++++++++++++++++++             | 73% ~05m 08s       |+++++++++++++++++++++++++++++++++++++             | 74% ~04m 57s       |++++++++++++++++++++++++++++++++++++++            | 75% ~04m 45s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~04m 34s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~04m 22s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~04m 10s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~03m 58s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~03m 46s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~03m 34s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 22s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 10s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~02m 58s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~02m 46s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 34s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 22s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 10s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~01m 58s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~01m 46s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 34s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 23s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 11s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~59s           |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~47s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~35s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~24s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~12s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=18m 51s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~18m 54s       |++                                                | 2 % ~19m 40s       |++                                                | 3 % ~19m 19s       |+++                                               | 4 % ~18m 49s       |+++                                               | 5 % ~18m 37s       |++++                                              | 6 % ~17m 43s       |++++                                              | 7 % ~16m 36s       |+++++                                             | 8 % ~15m 41s       |+++++                                             | 9 % ~15m 42s       |++++++                                            | 10% ~15m 46s       |++++++                                            | 11% ~15m 38s       |+++++++                                           | 12% ~15m 47s       |+++++++                                           | 14% ~15m 38s       |++++++++                                          | 15% ~15m 32s       |++++++++                                          | 16% ~15m 35s       |+++++++++                                         | 17% ~15m 45s       |+++++++++                                         | 18% ~15m 57s       |++++++++++                                        | 19% ~16m 07s       |++++++++++                                        | 20% ~16m 08s       |+++++++++++                                       | 21% ~16m 07s       |+++++++++++                                       | 22% ~16m 10s       |++++++++++++                                      | 23% ~16m 02s       |++++++++++++                                      | 24% ~15m 55s       |+++++++++++++                                     | 25% ~15m 46s       |++++++++++++++                                    | 26% ~15m 39s       |++++++++++++++                                    | 27% ~15m 29s       |+++++++++++++++                                   | 28% ~15m 15s       |+++++++++++++++                                   | 29% ~14m 59s       |++++++++++++++++                                  | 30% ~14m 50s       |++++++++++++++++                                  | 31% ~14m 36s       |+++++++++++++++++                                 | 32% ~14m 26s       |+++++++++++++++++                                 | 33% ~14m 11s       |++++++++++++++++++                                | 34% ~13m 56s       |++++++++++++++++++                                | 35% ~13m 44s       |+++++++++++++++++++                               | 36% ~13m 32s       |+++++++++++++++++++                               | 38% ~13m 18s       |++++++++++++++++++++                              | 39% ~13m 03s       |++++++++++++++++++++                              | 40% ~12m 50s       |+++++++++++++++++++++                             | 41% ~12m 38s       |+++++++++++++++++++++                             | 42% ~12m 25s       |++++++++++++++++++++++                            | 43% ~12m 12s       |++++++++++++++++++++++                            | 44% ~11m 58s       |+++++++++++++++++++++++                           | 45% ~11m 48s       |+++++++++++++++++++++++                           | 46% ~11m 37s       |++++++++++++++++++++++++                          | 47% ~11m 24s       |++++++++++++++++++++++++                          | 48% ~11m 10s       |+++++++++++++++++++++++++                         | 49% ~10m 57s       |+++++++++++++++++++++++++                         | 50% ~10m 43s       |++++++++++++++++++++++++++                        | 51% ~10m 31s       |+++++++++++++++++++++++++++                       | 52% ~10m 17s       |+++++++++++++++++++++++++++                       | 53% ~10m 04s       |++++++++++++++++++++++++++++                      | 54% ~09m 51s       |++++++++++++++++++++++++++++                      | 55% ~09m 37s       |+++++++++++++++++++++++++++++                     | 56% ~09m 24s       |+++++++++++++++++++++++++++++                     | 57% ~09m 10s       |++++++++++++++++++++++++++++++                    | 58% ~08m 56s       |++++++++++++++++++++++++++++++                    | 59% ~08m 44s       |+++++++++++++++++++++++++++++++                   | 60% ~08m 33s       |+++++++++++++++++++++++++++++++                   | 61% ~08m 20s       |++++++++++++++++++++++++++++++++                  | 62% ~08m 09s       |++++++++++++++++++++++++++++++++                  | 64% ~07m 57s       |+++++++++++++++++++++++++++++++++                 | 65% ~07m 45s       |+++++++++++++++++++++++++++++++++                 | 66% ~07m 31s       |++++++++++++++++++++++++++++++++++                | 67% ~07m 18s       |++++++++++++++++++++++++++++++++++                | 68% ~07m 04s       |+++++++++++++++++++++++++++++++++++               | 69% ~06m 51s       |+++++++++++++++++++++++++++++++++++               | 70% ~06m 38s       |++++++++++++++++++++++++++++++++++++              | 71% ~06m 25s       |++++++++++++++++++++++++++++++++++++              | 72% ~06m 12s       |+++++++++++++++++++++++++++++++++++++             | 73% ~05m 59s       |+++++++++++++++++++++++++++++++++++++             | 74% ~05m 45s       |++++++++++++++++++++++++++++++++++++++            | 75% ~05m 31s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~05m 16s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~05m 03s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~04m 49s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~04m 35s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~04m 21s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~04m 08s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 54s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 40s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~03m 26s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~03m 12s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 58s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 44s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 30s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 16s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~02m 03s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 49s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 35s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 22s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 08s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~55s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~41s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~27s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~14s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=21m 45s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~17m 27s       |++                                                | 2 % ~18m 05s       |++                                                | 3 % ~18m 06s       |+++                                               | 4 % ~17m 33s       |+++                                               | 5 % ~18m 32s       |++++                                              | 6 % ~18m 45s       |++++                                              | 7 % ~19m 06s       |+++++                                             | 8 % ~19m 08s       |+++++                                             | 9 % ~18m 55s       |++++++                                            | 10% ~18m 56s       |++++++                                            | 11% ~18m 42s       |+++++++                                           | 12% ~18m 29s       |+++++++                                           | 14% ~18m 16s       |++++++++                                          | 15% ~17m 56s       |++++++++                                          | 16% ~17m 36s       |+++++++++                                         | 17% ~17m 16s       |+++++++++                                         | 18% ~17m 17s       |++++++++++                                        | 19% ~17m 01s       |++++++++++                                        | 20% ~16m 44s       |+++++++++++                                       | 21% ~16m 25s       |+++++++++++                                       | 22% ~16m 06s       |++++++++++++                                      | 23% ~15m 50s       |++++++++++++                                      | 24% ~15m 39s       |+++++++++++++                                     | 25% ~15m 22s       |++++++++++++++                                    | 26% ~15m 08s       |++++++++++++++                                    | 27% ~14m 50s       |+++++++++++++++                                   | 28% ~14m 36s       |+++++++++++++++                                   | 29% ~14m 21s       |++++++++++++++++                                  | 30% ~14m 04s       |++++++++++++++++                                  | 31% ~13m 46s       |+++++++++++++++++                                 | 32% ~13m 34s       |+++++++++++++++++                                 | 33% ~13m 19s       |++++++++++++++++++                                | 34% ~13m 04s       |++++++++++++++++++                                | 35% ~12m 50s       |+++++++++++++++++++                               | 36% ~12m 35s       |+++++++++++++++++++                               | 38% ~12m 23s       |++++++++++++++++++++                              | 39% ~12m 12s       |++++++++++++++++++++                              | 40% ~11m 58s       |+++++++++++++++++++++                             | 41% ~11m 44s       |+++++++++++++++++++++                             | 42% ~11m 30s       |++++++++++++++++++++++                            | 43% ~11m 17s       |++++++++++++++++++++++                            | 44% ~11m 03s       |+++++++++++++++++++++++                           | 45% ~10m 50s       |+++++++++++++++++++++++                           | 46% ~10m 36s       |++++++++++++++++++++++++                          | 47% ~10m 27s       |++++++++++++++++++++++++                          | 48% ~10m 16s       |+++++++++++++++++++++++++                         | 49% ~10m 03s       |+++++++++++++++++++++++++                         | 50% ~09m 50s       |++++++++++++++++++++++++++                        | 51% ~09m 38s       |+++++++++++++++++++++++++++                       | 52% ~09m 24s       |+++++++++++++++++++++++++++                       | 53% ~09m 11s       |++++++++++++++++++++++++++++                      | 54% ~08m 58s       |++++++++++++++++++++++++++++                      | 55% ~08m 46s       |+++++++++++++++++++++++++++++                     | 56% ~08m 32s       |+++++++++++++++++++++++++++++                     | 57% ~08m 20s       |++++++++++++++++++++++++++++++                    | 58% ~08m 10s       |++++++++++++++++++++++++++++++                    | 59% ~07m 58s       |+++++++++++++++++++++++++++++++                   | 60% ~07m 48s       |+++++++++++++++++++++++++++++++                   | 61% ~07m 35s       |++++++++++++++++++++++++++++++++                  | 62% ~07m 22s       |++++++++++++++++++++++++++++++++                  | 64% ~07m 09s       |+++++++++++++++++++++++++++++++++                 | 65% ~06m 56s       |+++++++++++++++++++++++++++++++++                 | 66% ~06m 45s       |++++++++++++++++++++++++++++++++++                | 67% ~06m 32s       |++++++++++++++++++++++++++++++++++                | 68% ~06m 19s       |+++++++++++++++++++++++++++++++++++               | 69% ~06m 08s       |+++++++++++++++++++++++++++++++++++               | 70% ~05m 56s       |++++++++++++++++++++++++++++++++++++              | 71% ~05m 45s       |++++++++++++++++++++++++++++++++++++              | 72% ~05m 32s       |+++++++++++++++++++++++++++++++++++++             | 73% ~05m 21s       |+++++++++++++++++++++++++++++++++++++             | 74% ~05m 08s       |++++++++++++++++++++++++++++++++++++++            | 75% ~04m 56s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~04m 43s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~04m 31s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~04m 19s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~04m 06s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~03m 53s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~03m 41s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 28s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 16s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~03m 04s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~02m 51s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 39s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 27s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 15s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 03s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~01m 50s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 38s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 26s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 14s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 01s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~49s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~37s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~25s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~12s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=19m 40s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~19m 09s       |++                                                | 2 % ~17m 35s       |++                                                | 3 % ~16m 47s       |+++                                               | 4 % ~16m 33s       |+++                                               | 5 % ~16m 10s       |++++                                              | 6 % ~15m 44s       |++++                                              | 7 % ~15m 26s       |+++++                                             | 8 % ~15m 12s       |+++++                                             | 9 % ~15m 18s       |++++++                                            | 10% ~15m 18s       |++++++                                            | 11% ~15m 19s       |+++++++                                           | 12% ~15m 15s       |+++++++                                           | 14% ~15m 06s       |++++++++                                          | 15% ~14m 56s       |++++++++                                          | 16% ~14m 45s       |+++++++++                                         | 17% ~14m 32s       |+++++++++                                         | 18% ~14m 19s       |++++++++++                                        | 19% ~14m 07s       |++++++++++                                        | 20% ~13m 54s       |+++++++++++                                       | 21% ~13m 43s       |+++++++++++                                       | 22% ~13m 31s       |++++++++++++                                      | 23% ~13m 23s       |++++++++++++                                      | 24% ~13m 16s       |+++++++++++++                                     | 25% ~13m 06s       |++++++++++++++                                    | 26% ~12m 56s       |++++++++++++++                                    | 27% ~12m 49s       |+++++++++++++++                                   | 28% ~12m 40s       |+++++++++++++++                                   | 29% ~12m 30s       |++++++++++++++++                                  | 30% ~12m 19s       |++++++++++++++++                                  | 31% ~12m 08s       |+++++++++++++++++                                 | 32% ~11m 53s       |+++++++++++++++++                                 | 33% ~11m 38s       |++++++++++++++++++                                | 34% ~11m 24s       |++++++++++++++++++                                | 35% ~11m 11s       |+++++++++++++++++++                               | 36% ~10m 57s       |+++++++++++++++++++                               | 38% ~10m 47s       |++++++++++++++++++++                              | 39% ~10m 35s       |++++++++++++++++++++                              | 40% ~10m 21s       |+++++++++++++++++++++                             | 41% ~10m 09s       |+++++++++++++++++++++                             | 42% ~09m 57s       |++++++++++++++++++++++                            | 43% ~09m 47s       |++++++++++++++++++++++                            | 44% ~09m 35s       |+++++++++++++++++++++++                           | 45% ~09m 25s       |+++++++++++++++++++++++                           | 46% ~09m 13s       |++++++++++++++++++++++++                          | 47% ~09m 02s       |++++++++++++++++++++++++                          | 48% ~08m 51s       |+++++++++++++++++++++++++                         | 49% ~08m 39s       |+++++++++++++++++++++++++                         | 50% ~08m 28s       |++++++++++++++++++++++++++                        | 51% ~08m 16s       |+++++++++++++++++++++++++++                       | 52% ~08m 04s       |+++++++++++++++++++++++++++                       | 53% ~07m 53s       |++++++++++++++++++++++++++++                      | 54% ~07m 41s       |++++++++++++++++++++++++++++                      | 55% ~07m 30s       |+++++++++++++++++++++++++++++                     | 56% ~07m 20s       |+++++++++++++++++++++++++++++                     | 57% ~07m 09s       |++++++++++++++++++++++++++++++                    | 58% ~06m 59s       |++++++++++++++++++++++++++++++                    | 59% ~06m 50s       |+++++++++++++++++++++++++++++++                   | 60% ~06m 40s       |+++++++++++++++++++++++++++++++                   | 61% ~06m 30s       |++++++++++++++++++++++++++++++++                  | 62% ~06m 20s       |++++++++++++++++++++++++++++++++                  | 64% ~06m 11s       |+++++++++++++++++++++++++++++++++                 | 65% ~06m 01s       |+++++++++++++++++++++++++++++++++                 | 66% ~05m 51s       |++++++++++++++++++++++++++++++++++                | 67% ~05m 42s       |++++++++++++++++++++++++++++++++++                | 68% ~05m 32s       |+++++++++++++++++++++++++++++++++++               | 69% ~05m 21s       |+++++++++++++++++++++++++++++++++++               | 70% ~05m 11s       |++++++++++++++++++++++++++++++++++++              | 71% ~04m 60s       |++++++++++++++++++++++++++++++++++++              | 72% ~04m 49s       |+++++++++++++++++++++++++++++++++++++             | 73% ~04m 38s       |+++++++++++++++++++++++++++++++++++++             | 74% ~04m 28s       |++++++++++++++++++++++++++++++++++++++            | 75% ~04m 17s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~04m 07s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~03m 56s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~03m 45s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~03m 34s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~03m 24s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~03m 13s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 02s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~02m 51s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~02m 40s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~02m 30s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 19s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 08s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~01m 57s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~01m 47s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~01m 36s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 25s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 14s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 04s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~53s           |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~42s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~32s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~21s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~11s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=16m 52s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~21m 38s       |++                                                | 2 % ~20m 35s       |++                                                | 3 % ~20m 42s       |+++                                               | 4 % ~21m 05s       |+++                                               | 5 % ~20m 34s       |++++                                              | 6 % ~20m 44s       |++++                                              | 7 % ~20m 49s       |+++++                                             | 8 % ~20m 20s       |+++++                                             | 9 % ~19m 58s       |++++++                                            | 10% ~19m 38s       |++++++                                            | 11% ~19m 20s       |+++++++                                           | 12% ~19m 03s       |+++++++                                           | 14% ~18m 46s       |++++++++                                          | 15% ~18m 33s       |++++++++                                          | 16% ~18m 18s       |+++++++++                                         | 17% ~18m 06s       |+++++++++                                         | 18% ~17m 51s       |++++++++++                                        | 19% ~17m 38s       |++++++++++                                        | 20% ~17m 26s       |+++++++++++                                       | 21% ~17m 26s       |+++++++++++                                       | 22% ~17m 33s       |++++++++++++                                      | 23% ~17m 32s       |++++++++++++                                      | 24% ~17m 21s       |+++++++++++++                                     | 25% ~17m 16s       |++++++++++++++                                    | 26% ~17m 12s       |++++++++++++++                                    | 27% ~17m 08s       |+++++++++++++++                                   | 28% ~16m 57s       |+++++++++++++++                                   | 29% ~16m 48s       |++++++++++++++++                                  | 30% ~16m 34s       |++++++++++++++++                                  | 31% ~16m 28s       |+++++++++++++++++                                 | 32% ~16m 13s       |+++++++++++++++++                                 | 33% ~16m 03s       |++++++++++++++++++                                | 34% ~15m 46s       |++++++++++++++++++                                | 35% ~15m 37s       |+++++++++++++++++++                               | 36% ~15m 25s       |+++++++++++++++++++                               | 38% ~15m 12s       |++++++++++++++++++++                              | 39% ~14m 56s       |++++++++++++++++++++                              | 40% ~14m 37s       |+++++++++++++++++++++                             | 41% ~14m 19s       |+++++++++++++++++++++                             | 42% ~14m 02s       |++++++++++++++++++++++                            | 43% ~13m 44s       |++++++++++++++++++++++                            | 44% ~13m 27s       |+++++++++++++++++++++++                           | 45% ~13m 10s       |+++++++++++++++++++++++                           | 46% ~12m 54s       |++++++++++++++++++++++++                          | 47% ~12m 38s       |++++++++++++++++++++++++                          | 48% ~12m 22s       |+++++++++++++++++++++++++                         | 49% ~12m 05s       |+++++++++++++++++++++++++                         | 50% ~11m 45s       |++++++++++++++++++++++++++                        | 51% ~11m 28s       |+++++++++++++++++++++++++++                       | 52% ~11m 12s       |+++++++++++++++++++++++++++                       | 53% ~10m 57s       |++++++++++++++++++++++++++++                      | 54% ~10m 41s       |++++++++++++++++++++++++++++                      | 55% ~10m 27s       |+++++++++++++++++++++++++++++                     | 56% ~10m 12s       |+++++++++++++++++++++++++++++                     | 57% ~09m 53s       |++++++++++++++++++++++++++++++                    | 58% ~09m 34s       |++++++++++++++++++++++++++++++                    | 59% ~09m 18s       |+++++++++++++++++++++++++++++++                   | 60% ~09m 02s       |+++++++++++++++++++++++++++++++                   | 61% ~08m 48s       |++++++++++++++++++++++++++++++++                  | 62% ~08m 33s       |++++++++++++++++++++++++++++++++                  | 64% ~08m 21s       |+++++++++++++++++++++++++++++++++                 | 65% ~08m 08s       |+++++++++++++++++++++++++++++++++                 | 66% ~07m 56s       |++++++++++++++++++++++++++++++++++                | 67% ~07m 43s       |++++++++++++++++++++++++++++++++++                | 68% ~07m 30s       |+++++++++++++++++++++++++++++++++++               | 69% ~07m 17s       |+++++++++++++++++++++++++++++++++++               | 70% ~07m 05s       |++++++++++++++++++++++++++++++++++++              | 71% ~06m 52s       |++++++++++++++++++++++++++++++++++++              | 72% ~06m 39s       |+++++++++++++++++++++++++++++++++++++             | 73% ~06m 24s       |+++++++++++++++++++++++++++++++++++++             | 74% ~06m 10s       |++++++++++++++++++++++++++++++++++++++            | 75% ~05m 55s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~05m 41s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~05m 27s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~05m 12s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~04m 57s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~04m 42s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~04m 27s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~04m 13s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 58s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~03m 43s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~03m 28s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~03m 13s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 58s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 43s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 28s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~02m 13s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 58s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 43s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 28s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 14s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~59s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~44s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~29s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~15s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=23m 30s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~20m 06s       |++                                                | 2 % ~18m 14s       |++                                                | 3 % ~18m 17s       |+++                                               | 4 % ~19m 30s       |+++                                               | 5 % ~20m 26s       |++++                                              | 6 % ~20m 23s       |++++                                              | 7 % ~20m 46s       |+++++                                             | 8 % ~20m 31s       |+++++                                             | 9 % ~20m 16s       |++++++                                            | 10% ~20m 03s       |++++++                                            | 11% ~19m 42s       |+++++++                                           | 12% ~19m 15s       |+++++++                                           | 14% ~18m 54s       |++++++++                                          | 15% ~18m 44s       |++++++++                                          | 16% ~18m 28s       |+++++++++                                         | 17% ~18m 07s       |+++++++++                                         | 18% ~17m 47s       |++++++++++                                        | 19% ~17m 33s       |++++++++++                                        | 20% ~17m 16s       |+++++++++++                                       | 21% ~16m 59s       |+++++++++++                                       | 22% ~16m 49s       |++++++++++++                                      | 23% ~16m 32s       |++++++++++++                                      | 24% ~16m 20s       |+++++++++++++                                     | 25% ~16m 05s       |++++++++++++++                                    | 26% ~15m 48s       |++++++++++++++                                    | 27% ~15m 34s       |+++++++++++++++                                   | 28% ~15m 18s       |+++++++++++++++                                   | 29% ~15m 10s       |++++++++++++++++                                  | 30% ~14m 53s       |++++++++++++++++                                  | 31% ~14m 37s       |+++++++++++++++++                                 | 32% ~14m 22s       |+++++++++++++++++                                 | 33% ~14m 07s       |++++++++++++++++++                                | 34% ~13m 55s       |++++++++++++++++++                                | 35% ~13m 40s       |+++++++++++++++++++                               | 36% ~13m 24s       |+++++++++++++++++++                               | 38% ~13m 09s       |++++++++++++++++++++                              | 39% ~12m 59s       |++++++++++++++++++++                              | 40% ~12m 49s       |+++++++++++++++++++++                             | 41% ~12m 39s       |+++++++++++++++++++++                             | 42% ~12m 29s       |++++++++++++++++++++++                            | 43% ~12m 21s       |++++++++++++++++++++++                            | 44% ~12m 12s       |+++++++++++++++++++++++                           | 45% ~12m 04s       |+++++++++++++++++++++++                           | 46% ~11m 53s       |++++++++++++++++++++++++                          | 47% ~11m 43s       |++++++++++++++++++++++++                          | 48% ~11m 30s       |+++++++++++++++++++++++++                         | 49% ~11m 18s       |+++++++++++++++++++++++++                         | 50% ~11m 07s       |++++++++++++++++++++++++++                        | 51% ~10m 57s       |+++++++++++++++++++++++++++                       | 52% ~10m 45s       |+++++++++++++++++++++++++++                       | 53% ~10m 34s       |++++++++++++++++++++++++++++                      | 54% ~10m 22s       |++++++++++++++++++++++++++++                      | 55% ~10m 08s       |+++++++++++++++++++++++++++++                     | 56% ~09m 53s       |+++++++++++++++++++++++++++++                     | 57% ~09m 39s       |++++++++++++++++++++++++++++++                    | 58% ~09m 25s       |++++++++++++++++++++++++++++++                    | 59% ~09m 11s       |+++++++++++++++++++++++++++++++                   | 60% ~08m 57s       |+++++++++++++++++++++++++++++++                   | 61% ~08m 43s       |++++++++++++++++++++++++++++++++                  | 62% ~08m 29s       |++++++++++++++++++++++++++++++++                  | 64% ~08m 16s       |+++++++++++++++++++++++++++++++++                 | 65% ~08m 03s       |+++++++++++++++++++++++++++++++++                 | 66% ~07m 49s       |++++++++++++++++++++++++++++++++++                | 67% ~07m 34s       |++++++++++++++++++++++++++++++++++                | 68% ~07m 20s       |+++++++++++++++++++++++++++++++++++               | 69% ~07m 07s       |+++++++++++++++++++++++++++++++++++               | 70% ~06m 52s       |++++++++++++++++++++++++++++++++++++              | 71% ~06m 39s       |++++++++++++++++++++++++++++++++++++              | 72% ~06m 24s       |+++++++++++++++++++++++++++++++++++++             | 73% ~06m 09s       |+++++++++++++++++++++++++++++++++++++             | 74% ~05m 55s       |++++++++++++++++++++++++++++++++++++++            | 75% ~05m 41s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~05m 27s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~05m 13s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~04m 58s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~04m 44s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~04m 29s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~04m 15s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~04m 00s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 47s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~03m 33s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~03m 19s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~03m 05s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 51s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 37s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 23s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~02m 08s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 54s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 40s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 26s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 12s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~57s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~43s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~29s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~14s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=22m 57s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~25m 29s       |++                                                | 2 % ~24m 28s       |++                                                | 3 % ~23m 54s       |+++                                               | 4 % ~23m 34s       |+++                                               | 5 % ~24m 02s       |++++                                              | 6 % ~23m 58s       |++++                                              | 7 % ~23m 49s       |+++++                                             | 8 % ~23m 36s       |+++++                                             | 9 % ~23m 42s       |++++++                                            | 10% ~23m 15s       |++++++                                            | 11% ~22m 56s       |+++++++                                           | 12% ~22m 30s       |+++++++                                           | 14% ~22m 06s       |++++++++                                          | 15% ~21m 44s       |++++++++                                          | 16% ~21m 22s       |+++++++++                                         | 17% ~21m 00s       |+++++++++                                         | 18% ~20m 40s       |++++++++++                                        | 19% ~20m 24s       |++++++++++                                        | 20% ~20m 08s       |+++++++++++                                       | 21% ~19m 53s       |+++++++++++                                       | 22% ~19m 32s       |++++++++++++                                      | 23% ~19m 10s       |++++++++++++                                      | 24% ~18m 55s       |+++++++++++++                                     | 25% ~18m 38s       |++++++++++++++                                    | 26% ~18m 24s       |++++++++++++++                                    | 27% ~18m 13s       |+++++++++++++++                                   | 28% ~17m 53s       |+++++++++++++++                                   | 29% ~17m 36s       |++++++++++++++++                                  | 30% ~17m 22s       |++++++++++++++++                                  | 31% ~17m 09s       |+++++++++++++++++                                 | 32% ~16m 56s       |+++++++++++++++++                                 | 33% ~16m 46s       |++++++++++++++++++                                | 34% ~16m 30s       |++++++++++++++++++                                | 35% ~16m 12s       |+++++++++++++++++++                               | 36% ~15m 55s       |+++++++++++++++++++                               | 38% ~15m 37s       |++++++++++++++++++++                              | 39% ~15m 21s       |++++++++++++++++++++                              | 40% ~15m 04s       |+++++++++++++++++++++                             | 41% ~14m 47s       |+++++++++++++++++++++                             | 42% ~14m 32s       |++++++++++++++++++++++                            | 43% ~14m 16s       |++++++++++++++++++++++                            | 44% ~14m 00s       |+++++++++++++++++++++++                           | 45% ~13m 42s       |+++++++++++++++++++++++                           | 46% ~13m 25s       |++++++++++++++++++++++++                          | 47% ~13m 08s       |++++++++++++++++++++++++                          | 48% ~12m 51s       |+++++++++++++++++++++++++                         | 49% ~12m 36s       |+++++++++++++++++++++++++                         | 50% ~12m 21s       |++++++++++++++++++++++++++                        | 51% ~12m 06s       |+++++++++++++++++++++++++++                       | 52% ~11m 49s       |+++++++++++++++++++++++++++                       | 53% ~11m 32s       |++++++++++++++++++++++++++++                      | 54% ~11m 15s       |++++++++++++++++++++++++++++                      | 55% ~10m 59s       |+++++++++++++++++++++++++++++                     | 56% ~10m 42s       |+++++++++++++++++++++++++++++                     | 57% ~10m 26s       |++++++++++++++++++++++++++++++                    | 58% ~10m 15s       |++++++++++++++++++++++++++++++                    | 59% ~09m 59s       |+++++++++++++++++++++++++++++++                   | 60% ~09m 43s       |+++++++++++++++++++++++++++++++                   | 61% ~09m 28s       |++++++++++++++++++++++++++++++++                  | 62% ~09m 12s       |++++++++++++++++++++++++++++++++                  | 64% ~08m 56s       |+++++++++++++++++++++++++++++++++                 | 65% ~08m 40s       |+++++++++++++++++++++++++++++++++                 | 66% ~08m 25s       |++++++++++++++++++++++++++++++++++                | 67% ~08m 10s       |++++++++++++++++++++++++++++++++++                | 68% ~07m 55s       |+++++++++++++++++++++++++++++++++++               | 69% ~07m 40s       |+++++++++++++++++++++++++++++++++++               | 70% ~07m 25s       |++++++++++++++++++++++++++++++++++++              | 71% ~07m 09s       |++++++++++++++++++++++++++++++++++++              | 72% ~06m 54s       |+++++++++++++++++++++++++++++++++++++             | 73% ~06m 40s       |+++++++++++++++++++++++++++++++++++++             | 74% ~06m 25s       |++++++++++++++++++++++++++++++++++++++            | 75% ~06m 10s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~05m 55s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~05m 40s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~05m 24s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~05m 08s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~04m 53s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~04m 38s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~04m 22s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~04m 07s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~03m 52s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~03m 36s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~03m 21s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~03m 06s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 50s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 35s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~02m 19s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~02m 04s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 48s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 33s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 17s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~01m 02s       |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~46s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~31s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~15s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=24m 44s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~19m 47s       |++                                                | 2 % ~19m 21s       |++                                                | 3 % ~18m 20s       |+++                                               | 4 % ~17m 47s       |+++                                               | 5 % ~17m 47s       |++++                                              | 6 % ~17m 50s       |++++                                              | 7 % ~18m 04s       |+++++                                             | 8 % ~17m 36s       |+++++                                             | 9 % ~17m 17s       |++++++                                            | 10% ~17m 13s       |++++++                                            | 11% ~16m 59s       |+++++++                                           | 12% ~16m 45s       |+++++++                                           | 14% ~16m 27s       |++++++++                                          | 15% ~16m 17s       |++++++++                                          | 16% ~16m 19s       |+++++++++                                         | 17% ~16m 10s       |+++++++++                                         | 18% ~16m 01s       |++++++++++                                        | 19% ~15m 56s       |++++++++++                                        | 20% ~15m 47s       |+++++++++++                                       | 21% ~15m 39s       |+++++++++++                                       | 22% ~15m 45s       |++++++++++++                                      | 23% ~15m 39s       |++++++++++++                                      | 24% ~15m 24s       |+++++++++++++                                     | 25% ~15m 33s       |++++++++++++++                                    | 26% ~15m 21s       |++++++++++++++                                    | 27% ~15m 09s       |+++++++++++++++                                   | 28% ~14m 58s       |+++++++++++++++                                   | 29% ~14m 44s       |++++++++++++++++                                  | 30% ~14m 30s       |++++++++++++++++                                  | 31% ~14m 21s       |+++++++++++++++++                                 | 32% ~14m 06s       |+++++++++++++++++                                 | 33% ~13m 49s       |++++++++++++++++++                                | 34% ~13m 33s       |++++++++++++++++++                                | 35% ~13m 16s       |+++++++++++++++++++                               | 36% ~13m 00s       |+++++++++++++++++++                               | 38% ~12m 52s       |++++++++++++++++++++                              | 39% ~12m 40s       |++++++++++++++++++++                              | 40% ~12m 25s       |+++++++++++++++++++++                             | 41% ~12m 10s       |+++++++++++++++++++++                             | 42% ~11m 57s       |++++++++++++++++++++++                            | 43% ~11m 42s       |++++++++++++++++++++++                            | 44% ~11m 27s       |+++++++++++++++++++++++                           | 45% ~11m 12s       |+++++++++++++++++++++++                           | 46% ~10m 57s       |++++++++++++++++++++++++                          | 47% ~10m 43s       |++++++++++++++++++++++++                          | 48% ~10m 29s       |+++++++++++++++++++++++++                         | 49% ~10m 15s       |+++++++++++++++++++++++++                         | 50% ~10m 00s       |++++++++++++++++++++++++++                        | 51% ~09m 47s       |+++++++++++++++++++++++++++                       | 52% ~09m 35s       |+++++++++++++++++++++++++++                       | 53% ~09m 25s       |++++++++++++++++++++++++++++                      | 54% ~09m 11s       |++++++++++++++++++++++++++++                      | 55% ~08m 58s       |+++++++++++++++++++++++++++++                     | 56% ~08m 44s       |+++++++++++++++++++++++++++++                     | 57% ~08m 32s       |++++++++++++++++++++++++++++++                    | 58% ~08m 19s       |++++++++++++++++++++++++++++++                    | 59% ~08m 06s       |+++++++++++++++++++++++++++++++                   | 60% ~07m 54s       |+++++++++++++++++++++++++++++++                   | 61% ~07m 41s       |++++++++++++++++++++++++++++++++                  | 62% ~07m 29s       |++++++++++++++++++++++++++++++++                  | 64% ~07m 17s       |+++++++++++++++++++++++++++++++++                 | 65% ~07m 05s       |+++++++++++++++++++++++++++++++++                 | 66% ~06m 53s       |++++++++++++++++++++++++++++++++++                | 67% ~06m 41s       |++++++++++++++++++++++++++++++++++                | 68% ~06m 29s       |+++++++++++++++++++++++++++++++++++               | 69% ~06m 17s       |+++++++++++++++++++++++++++++++++++               | 70% ~06m 05s       |++++++++++++++++++++++++++++++++++++              | 71% ~05m 53s       |++++++++++++++++++++++++++++++++++++              | 72% ~05m 40s       |+++++++++++++++++++++++++++++++++++++             | 73% ~05m 28s       |+++++++++++++++++++++++++++++++++++++             | 74% ~05m 15s       |++++++++++++++++++++++++++++++++++++++            | 75% ~05m 03s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~04m 50s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~04m 37s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~04m 25s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~04m 12s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~03m 59s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~03m 46s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~03m 34s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~03m 21s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~03m 08s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~02m 56s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~02m 43s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~02m 31s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 18s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 05s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~01m 53s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~01m 40s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 28s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 15s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 03s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~50s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~38s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~25s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~13s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=20m 02s
  |                                                  | 0 % ~calculating   |+                                                 | 1 % ~23m 10s       |++                                                | 2 % ~27m 27s       |++                                                | 3 % ~26m 19s       |+++                                               | 4 % ~25m 47s       |+++                                               | 5 % ~24m 35s       |++++                                              | 6 % ~23m 55s       |++++                                              | 7 % ~23m 43s       |+++++                                             | 8 % ~23m 16s       |+++++                                             | 9 % ~22m 33s       |++++++                                            | 10% ~22m 35s       |++++++                                            | 11% ~22m 28s       |+++++++                                           | 12% ~22m 30s       |+++++++                                           | 14% ~22m 21s       |++++++++                                          | 15% ~22m 27s       |++++++++                                          | 16% ~22m 07s       |+++++++++                                         | 17% ~21m 43s       |+++++++++                                         | 18% ~21m 21s       |++++++++++                                        | 19% ~21m 03s       |++++++++++                                        | 20% ~20m 52s       |+++++++++++                                       | 21% ~20m 37s       |+++++++++++                                       | 22% ~20m 19s       |++++++++++++                                      | 23% ~20m 05s       |++++++++++++                                      | 24% ~19m 51s       |+++++++++++++                                     | 25% ~19m 34s       |++++++++++++++                                    | 26% ~19m 20s       |++++++++++++++                                    | 27% ~18m 60s       |+++++++++++++++                                   | 28% ~18m 39s       |+++++++++++++++                                   | 29% ~18m 20s       |++++++++++++++++                                  | 30% ~18m 07s       |++++++++++++++++                                  | 31% ~17m 53s       |+++++++++++++++++                                 | 32% ~17m 34s       |+++++++++++++++++                                 | 33% ~17m 16s       |++++++++++++++++++                                | 34% ~17m 00s       |++++++++++++++++++                                | 35% ~16m 44s       |+++++++++++++++++++                               | 36% ~16m 26s       |+++++++++++++++++++                               | 38% ~16m 09s       |++++++++++++++++++++                              | 39% ~15m 53s       |++++++++++++++++++++                              | 40% ~15m 34s       |+++++++++++++++++++++                             | 41% ~15m 14s       |+++++++++++++++++++++                             | 42% ~14m 58s       |++++++++++++++++++++++                            | 43% ~14m 44s       |++++++++++++++++++++++                            | 44% ~14m 27s       |+++++++++++++++++++++++                           | 45% ~14m 14s       |+++++++++++++++++++++++                           | 46% ~14m 01s       |++++++++++++++++++++++++                          | 47% ~13m 45s       |++++++++++++++++++++++++                          | 48% ~13m 27s       |+++++++++++++++++++++++++                         | 49% ~13m 11s       |+++++++++++++++++++++++++                         | 50% ~12m 56s       |++++++++++++++++++++++++++                        | 51% ~12m 43s       |+++++++++++++++++++++++++++                       | 52% ~12m 28s       |+++++++++++++++++++++++++++                       | 53% ~12m 12s       |++++++++++++++++++++++++++++                      | 54% ~11m 56s       |++++++++++++++++++++++++++++                      | 55% ~11m 39s       |+++++++++++++++++++++++++++++                     | 56% ~11m 22s       |+++++++++++++++++++++++++++++                     | 57% ~11m 05s       |++++++++++++++++++++++++++++++                    | 58% ~10m 49s       |++++++++++++++++++++++++++++++                    | 59% ~10m 32s       |+++++++++++++++++++++++++++++++                   | 60% ~10m 15s       |+++++++++++++++++++++++++++++++                   | 61% ~09m 57s       |++++++++++++++++++++++++++++++++                  | 62% ~09m 40s       |++++++++++++++++++++++++++++++++                  | 64% ~09m 24s       |+++++++++++++++++++++++++++++++++                 | 65% ~09m 08s       |+++++++++++++++++++++++++++++++++                 | 66% ~08m 52s       |++++++++++++++++++++++++++++++++++                | 67% ~08m 36s       |++++++++++++++++++++++++++++++++++                | 68% ~08m 20s       |+++++++++++++++++++++++++++++++++++               | 69% ~08m 05s       |+++++++++++++++++++++++++++++++++++               | 70% ~07m 48s       |++++++++++++++++++++++++++++++++++++              | 71% ~07m 31s       |++++++++++++++++++++++++++++++++++++              | 72% ~07m 15s       |+++++++++++++++++++++++++++++++++++++             | 73% ~06m 58s       |+++++++++++++++++++++++++++++++++++++             | 74% ~06m 41s       |++++++++++++++++++++++++++++++++++++++            | 75% ~06m 24s       |+++++++++++++++++++++++++++++++++++++++           | 76% ~06m 08s       |+++++++++++++++++++++++++++++++++++++++           | 77% ~05m 52s       |++++++++++++++++++++++++++++++++++++++++          | 78% ~05m 36s       |++++++++++++++++++++++++++++++++++++++++          | 79% ~05m 21s       |+++++++++++++++++++++++++++++++++++++++++         | 80% ~05m 05s       |+++++++++++++++++++++++++++++++++++++++++         | 81% ~04m 49s       |++++++++++++++++++++++++++++++++++++++++++        | 82% ~04m 33s       |++++++++++++++++++++++++++++++++++++++++++        | 83% ~04m 17s       |+++++++++++++++++++++++++++++++++++++++++++       | 84% ~04m 00s       |+++++++++++++++++++++++++++++++++++++++++++       | 85% ~03m 44s       |++++++++++++++++++++++++++++++++++++++++++++      | 86% ~03m 28s       |++++++++++++++++++++++++++++++++++++++++++++      | 88% ~03m 11s       |+++++++++++++++++++++++++++++++++++++++++++++     | 89% ~02m 55s       |+++++++++++++++++++++++++++++++++++++++++++++     | 90% ~02m 39s       |++++++++++++++++++++++++++++++++++++++++++++++    | 91% ~02m 23s       |++++++++++++++++++++++++++++++++++++++++++++++    | 92% ~02m 08s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 93% ~01m 52s       |+++++++++++++++++++++++++++++++++++++++++++++++   | 94% ~01m 36s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 95% ~01m 20s       |++++++++++++++++++++++++++++++++++++++++++++++++  | 96% ~01m 04s       |+++++++++++++++++++++++++++++++++++++++++++++++++ | 97% ~48s           |+++++++++++++++++++++++++++++++++++++++++++++++++ | 98% ~32s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 99% ~16s           |++++++++++++++++++++++++++++++++++++++++++++++++++| 100% elapsed=25m 43s

#### z-scores per cell type annotation_____no_output_____

%%R -i outdir -i experiment_prefix

ATAC_Seurat_withChromVar <- readRDS(file = paste0(outdir, experiment_prefix, "_chromVar_binary.rds"))
DefaultAssay(ATAC_Seurat_withChromVar) <- "chromvar" # Change default assay
ATAC_Seurat_withChromVarAn object of class Seurat
203263 features across 67202 samples within 2 assays
Active assay: chromvar (1920 features, 0 variable features)
 1 other assay present: peaks
 3 dimensional reductions calculated: cisTopics, umap, harmony
%%R -o chromvar_assay

chromvar_assay <- as.data.frame(t(x = as.matrix(x = GetAssayData(object = ATAC_Seurat_withChromVar))))_____no_output_____chromvar_assay.shape_____no_output_____# Mapping of cells to cell types
ATAC_scanpy.obs['barcode'] = ATAC_scanpy.obs_names
cells_celltypes = ATAC_scanpy.obs[['barcode', 'cell_type']]
cells_celltypes.head()_____no_output_____cells_mapping_clusters = cells_celltypes.groupby('cell_type').groups_____no_output_____# Mean by rows
for key in cells_mapping_clusters:
    print(key)
    chromvar_assay.loc[key] = chromvar_assay.loc[cells_mapping_clusters[key]].mean(axis = 0)Germ cells
coelEpi
sPAX8
preGC_I_OSR1
ovarianSurf
preGC_II_KITLG
preGC_III_GJA1
Oi
Gi
M_prog_ISL1
M_MGP
M_MullDuct_LGR5
PV
Epithelial
Endothelial
Immune
Neural
rows_celltypes = list(np.unique(cells_celltypes['cell_type']))
print(rows_celltypes)
# Select only aggregated rows
chromvar_assay = chromvar_assay.loc[rows_celltypes, ]['Endothelial', 'Epithelial', 'Germ cells', 'Gi', 'Immune', 'M_MGP', 'M_MullDuct_LGR5', 'M_prog_ISL1', 'Neural', 'Oi', 'PV', 'coelEpi', 'ovarianSurf', 'preGC_III_GJA1', 'preGC_II_KITLG', 'preGC_I_OSR1', 'sPAX8']
# Normalize by column (TF)
from scipy.stats import zscore
chromvar_assay_zscore = chromvar_assay.apply(zscore)_____no_output_____chromvar_assay_zscore.to_csv(outdir + experiment_prefix + "_chromvar_celltype_zscore.csv")_____no_output_____chromvar_assay_zscore = chromvar_assay_zscore.reindex(['Germ cells',
                                                 'coelEpi',  'sPAX8', 'preGC_I_OSR1',
                                                 'ovarianSurf','preGC_II_KITLG', 'preGC_III_GJA1',
                                                  'Oi',  'Gi',
                                                 'M_prog_ISL1', 'M_MGP', 'M_MullDuct_LGR5',
                                                 'PV','Epithelial', 'Endothelial', 'Immune', 'Neural'
                                          ])
_____no_output_____import seaborn as sns
my_cmap = sns.diverging_palette(240, 10, s = 100, l = 40,n=10, as_cmap=True)
_____no_output_____chromvar_assay_zscore.shape_____no_output_____heatmap_TFs = sns.clustermap(chromvar_assay_zscore, row_cluster=False, annot=False, xticklabels = False, yticklabels = "auto",  figsize=(8, 8),
                         cmap = my_cmap)
heatmap_TFs.savefig("/home/jovyan/MULTIOME_april2021/figures_females/heatmap_females_TFs.pdf")
/opt/conda/envs/atac_env/lib/python3.8/site-packages/seaborn/matrix.py:649: UserWarning: Clustering large matrix with scipy. Installing `fastcluster` may give better performance.
  warnings.warn(msg)

### Select top 3 TFs per cell type _____no_output_____

top_tfs = []
for c in chromvar_assay_zscore.index:
    df_c = chromvar_assay_zscore.loc[c]
    df_c = df_c.sort_values(axis = 0, ascending = False)
    top_4 = df_c[0:5].index
    top_tfs.extend(top_4)_____no_output_____chromvar_assay_zscore_top4 = chromvar_assay_zscore[top_tfs]_____no_output_____chromvar_assay_zscore_top4.columns_____no_output_____to_plot = ['DMRT6.HOMER', 'EGR2.SwissRegulon', 'PAX1.MA0779.1', 'FOXF2.MA0030.1','MNX1.MA0707.1',
            'MESP1.SwissRegulon',  'TCF3.MA0522.2',  'MSC.SwissRegulon','OLIG1.SwissRegulon','HOXC13.MA0907.1',
             'HSF2.MA0770.1', 'TWIST1.MA1123.1','EBF3.SwissRegulon',
           'HNF1B.MA0153.2', 'FLI1.SwissRegulon', 'SPI1.MA0080.4', 'RFX3.SwissRegulon',
           ]
len(to_plot)_____no_output_____chromvar_assay_zscore_top4 = chromvar_assay_zscore_top4[to_plot]
chromvar_assay_zscore_top4.columns_____no_output_____chromvar_assay_zscore_top4.columns = [col.split(".")[0] for col in chromvar_assay_zscore_top4.columns]_____no_output_____import matplotlib
divnorm = matplotlib.colors.DivergingNorm(vmin=chromvar_assay_zscore.values.min(), vcenter=0, vmax=chromvar_assay_zscore.values.max())
<ipython-input-91-63a9665b3b7d>:2: MatplotlibDeprecationWarning:
The DivergingNorm class was deprecated in Matplotlib 3.2 and will be removed two minor releases later. Use TwoSlopeNorm instead.
  divnorm = matplotlib.colors.DivergingNorm(vmin=chromvar_assay_zscore.values.min(), vcenter=0, vmax=chromvar_assay_zscore.values.max())
sns.set(font_scale=1.8)
heatmap_TFs = sns.clustermap(chromvar_assay_zscore_top4, row_cluster=False, col_cluster = False, annot=False, xticklabels = True, yticklabels = "auto",  figsize=(10, 10),
                         cmap = my_cmap, norm=divnorm)
heatmap_TFs.savefig("/home/jovyan/MULTIOME_april2021/figures_females/heatmap_females_TFs_top1_curated.pdf")
_____no_output_____

### Plot TFs associated with DSD_____no_output_____

dsd_tfs = ['ARX', 'ATRX', 'BMP15', 'CBX2', 'DAX1', 'DHH', 'DHX37', 'DMRT1', 'EMX2', 'ESR2', 'FGF9', 'FGFR2', 'FOXL2',
          'GATA4', 'HHAT', 'MAP3K1', 'NR2F2', 'NR5A1', 'NUP107', 'RSPO1', 'SOX3', 'SOX8', 'SOX9', 'SOX10', 'SRY',
          'STARD8', 'TSPYL1', 'WNT4', 'WT1', 'WWOX', 'ZFPM2', 'ZNRF3']_____no_output_____df = pd.DataFrame(chromvar_assay_zscore.columns)_____no_output_____df[0][df[0].str.contains('|'.join(dsd_tfs))]
_____no_output_____dsd_tfs_present = ['ARX.SwissRegulon', 'DMRT1.HOMER', 'EMX2.SwissRegulon', 'ESR2.SwissRegulon',
                  'GATA4.SwissRegulon', 'FOXL2.HOMER',
                  'NR2F2.MA1111.1', 'SOX3.SwissRegulon',
                   'SOX8-HUMAN.H10MO.D', 'SOX9.SwissRegulon',
                  'SOX10.SwissRegulon', 'SRY.SwissRegulon',
                  'WT1-HUMAN.H10MO.D']
_____no_output_____chromvar_assay_zscore_DSD = chromvar_assay_zscore[dsd_tfs_present]
chromvar_assay_zscore_DSD.columns_____no_output_____chromvar_assay_zscore_DSD.columns = [col.split(".")[0] for col in chromvar_assay_zscore_DSD.columns]
chromvar_assay_zscore_DSD.columns = [col.split("-")[0] for col in chromvar_assay_zscore_DSD.columns]
_____no_output_____sns.set(font_scale=1.8)
heatmap_TFs = sns.clustermap(chromvar_assay_zscore_DSD, row_cluster=False, col_cluster = False, annot=False, xticklabels = True, yticklabels = "auto",  figsize=(14, 14),
                         cmap = my_cmap, norm=divnorm)
heatmap_TFs.savefig("/home/jovyan/MULTIOME_april2021/figures_females/DSD_TFs_heatmap.pdf")
_____no_output_____

#### End of notebook_____no_output_____