1 TODO

Remove MSstats logging. – done
Make explicit spearman correlations between methods. – done
- Do both for all data and for the top-50/bottom-50 – done, but weird.
Make 100% certain that the samples are annotated correctly. – done.
Limma/MSstats/EdgeR venn disagram for all up in CF
- Repeat in all down CF
- Repeat with a cutoff, top-50
Individual plots for P/PE proteins.

2 Analyzing data from openMS and friends.

In preprocessing_comet_highres.Rmd, I used the openMS tutorials and supplemental materials from a couple papers to hopefully correctly perform the various preprocessing tasks required to extract intensity data from DIA/SWATH transitions.

The final steps of that process combined the transition intensities from every sample into a metadata file (results/tric/HCD_meta.tsv), an intensity matrix (results/tric/HCD_outmatrix.tsv), and a feature aligned output matrix (results/tric/aligned_comet_HCD.tsv).

My reading of the SWATH2stats and MSstats source code suggests to me that the log2(intensities) of the feature aligned data are our final proxy for protein abundance. At first glance, this suggests to me that these data might follow a distribution similar to RNASeq data (negative binomial, but perhaps with a bigger tail?). In addition, by the time we use tric on the data, we have a count matrix and sample annotation data frames which look remarkably similar to those used in a RNASeq expressionset. Indeed, by the end of the MSstats processing, it creates a MSnSet class of its own which uses fData/exprs/pData.

For the curious, my reasoning for saying that the log intensities are our proxy for abundance comes from MSstats/R/DataProcess.R in a clause which looks like:

if (logTrans == 2) {
  work[["ABUNDANCE"]] <- log2(work[["ABUNDANCE"]])
} else if (logTrans == 10) {
  work[["ABUNDANCE"]] <- log10(work[["ABUNDANCE"]])
} else {
  ## Above there was a check for only log 2 and 10, but we can do e if we want.
  ## I might go back up there and remove that check. Long live e! 2.718282 rules!
  work[["ABUNDANCE"]] <- log(work[["ABUNDANCE"]]) / log(logTrans)
}

(Note: I added the natural log to the set of conditions, but otherwise the logic is unchanged.)

With that in mind, I want to use some tools with which I am familiar in order to try to understand these data. Therefore I will first attempt to coerce my tric aligned data and annotations into a ‘normal’ expressionset. Then I want to do some diagnostic plots which, if I am wrong and these distributions are not as expected, will be conceptually incorrect (I don’t yet think I am wrong).

2.1 Sample annotation via SWATH2stats

I am using the SWATH2stats vignette as my primary source of information. Thus I see that it uses the OpenSWATH_SM3_GoldStandardAutomatedResults_human_peakgroups.txt which has a format nearly identical to the tric output matrix. Thus for the moment I will assume that the proper input for SWATH2stats is ‘results/tric/comet_HCD.tsv’ and not the metadata nor output matrix. In addition, SWATH2stats provides a similar conversion function which takes the tric output and coerces it into a similar matrix after performing its various filtering tasks. Therefore I will use that.

I keep a sample sheet of all the DIA samples used in this analysis in ‘sample_sheets/dia_samples.xlsx’ It should contain all the other required data with one important caveat, I removed 1 sample by ‘commenting’ it (e.g. prefixing it with ‘##’ – which is an admittedly dumb thing to do in an excel file.

One last caveat: I hacked the SWATH2stats sample_annotation() function to add a couple columns in an attempt to make it a little more robust when faced with sample sheets with differently named columns.

In addition, SWATH2stats provides some nice filtering and combination functions which should be considered when generating various expressionset data structures later.

2.1.1 Creating a swath2stats experiment using our comet-derived library data

our_tric_data <- read.csv("results/tric/20180611/whole_8mz/comet_HCD.tsv", sep="\t")

sample_annot <- openxlsx::read.xlsx("sample_sheets/Mtb_dia_samples.xlsx")
rownames(sample_annot) <- make.names(sample_annot[["sampleid"]], unique=TRUE)
## Drop samples starting with comments
keep_idx <- ! grepl(pattern="##", x=sample_annot[["sampleid"]])
sample_annot <- sample_annot[keep_idx, ]
expt_idx <- sample_annot[["expt_id"]] == "may2018" | sample_annot[["expt_id"]] == "mar2018"
expt_idx[is.na(sample_annot[["expt_id"]])] <- FALSE
sample_annot <- sample_annot[expt_idx, ]
mz_idx <- sample_annot[["windowsize"]] == "8"
sample_annot <- sample_annot[mz_idx, ]
## Set the mzXML column to match the filename column in the data.
loaded <- sm(devtools::load_all("~/scratch/git/SWATH2stats"))
## s2s, my witty way of shortening SWATH2stats...
our_s2s_exp <- sample_annotation(data=our_tric_data,
                                 sample_annotation=sample_annot,
                                 run_column="run",
                                 fullpeptidename_column="fullunimodpeptidename")

## Found the same mzXML files in the annotations and data.

2.1.2 Creating a swath2stats experiment using the tuberculist-derived library data

This just repeats the previous stanza, but uses the tuberculist transition-library swath data as its input rather than our own comet data.

There is one important caveat in the following block: I used a regex to remove the second half of geneID_geneName so that later when I merge in the annotation data I have it will match.

tb_tric_data <- read.csv("results/tric/20180611/whole_8mz_tuberculist/comet_HCD.tsv", sep="\t")

tb_tric_data[["ProteinName"]] <- gsub(pattern="^(.*)_.*$", replacement="\\1",
                                      x=tb_tric_data[["ProteinName"]])
tb_s2s_exp <- sample_annotation(data=tb_tric_data,
                                sample_annotation=sample_annot,
                                fullpeptidename_column="fullunimodpeptidename")

## Found the same mzXML files in the annotations and data.

Now I have a couple data structures which should prove useful for the metrics provided by SWATH2stats, MSstats, and my own hpgltools.

3 SWATH2stats continued

The various metrics and filters provided by SWATH2stats seem quite reasonable to me. The only thing that really bothers me is that they are all case sensitive and I found that the most recent tric changed the capitalization of a column, causing these to all fall down. Therefore I went in and made everything case insensitive in a fashion similar to that done by MSstats (except I hate capital letters, so I used tolower() rather than toupper()).

3.1 Run filters on our comet data

## Get correlations on a sample by sample basis
pp(file="images/20180611_our_swath2stats_sample_cor.png")

## Going to write the image to: images/20180611_our_swath2stats_sample_cor.png when dev.off() is called.

sample_cor <- plot_correlation_between_samples(our_s2s_exp, size=2,
                                               fun.aggregate=sum,
                                               column.values="intensity")
dev.off()

## png 
##   2

sample_cond_rep_cor <- plot_correlation_between_samples(our_s2s_exp, size=2,
                                                        comparison=transition_group_id ~
                                                          condition + bioreplicate + run,
                                                        fun.aggregate=sum,
                                                        column.values="intensity")

decoy_lists <- assess_decoy_rate(our_s2s_exp)

## Number of non-decoy peptides: 3259
## Number of decoy peptides: 182
## Decoy rate: 0.0558

## This seems a bit high to me, yesno?
fdr_overall <- assess_fdr_overall(our_s2s_exp, output="Rconsole", plot=TRUE)

byrun_fdr <- assess_fdr_byrun(our_s2s_exp, FFT=0.7, plot=TRUE, output="Rconsole")

## The average FDR by run on assay level is 0.007

## The average FDR by run on peptide level is 0.007

## The average FDR by run on protein level is 0.019

chosen_mscore <- mscore4assayfdr(our_s2s_exp, FFT=0.7, fdr_target=0.02)

## Target assay FDR: 0.02

## Required overall m-score cutoff: 0.0019953
## achieving assay FDR: 0.0189

prot_score <- mscore4protfdr(our_s2s_exp, FFT=0.7, fdr_target=0.02)

## Target protein FDR: 0.02

## Required overall m-score cutoff: 0.00050119
## achieving protein FDR: 0.0181

our_mscore_filtered <- filter_mscore(our_s2s_exp, chosen_mscore)

## Original dimension: 47362, new dimension: 45761, difference: 1601.

our_freq_mscore <- filter_mscore_freqobs(our_s2s_exp, 0.01, 0.8, rm.decoy=FALSE)

## Peptides need to have been quantified in more conditions than: 18.4 in order to pass this percentage-based threshold.

## Fraction of peptides selected: 1

## Original dimension: 47841, new dimension: 47841, difference: 0.

our_data_filtered_fdr <- filter_mscore_fdr(our_mscore_filtered, FFT=0.7,
                                           overall_protein_fdr_target=prot_score,
                                           upper_overall_peptide_fdr_limit=0.05)

## Target protein FDR: 0.000501187233627273

## Required overall m-score cutoff: 0.01
## achieving protein FDR: 0

## filter_mscore_fdr is filtering the data...

## finding m-score cutoff to achieve desired protein FDR in protein master list..

## finding m-score cutoff to achieve desired global peptide FDR..

## Target peptide FDR: 0.05

## Required overall m-score cutoff: 0.01

## Achieving peptide FDR: 0

## Proteins selected: 
## Total proteins selected: 945
## Final target proteins: 945
## Final decoy proteins: 0

## Peptides mapping to these protein entries selected:
## Total mapping peptides: 3203
## Final target peptides: 3203
## Final decoy peptides: 0

## Total peptides selected from:
## Total peptides: 3203
## Final target peptides: 3203
## Final decoy peptides: 0

## Individual run FDR quality of the peptides was not calculated
## as not every run contains a decoy.

## The decoys have been removed from the returned data.

our_only_proteotypic <- filter_proteotypic_peptides(our_data_filtered_fdr)

## Number of proteins detected: 947
## Protein identifiers: Rv1908c, Rv0242c, Rv3224, Rv0667, Rv1133c, Rv3036c

## Number of proteins detected that are supported by a proteotypic peptide: 930
## Number of proteotypic peptides detected: 3158

our_all_filtered <- filter_all_peptides(our_only_proteotypic)

## Number of proteins detected: 930

## First 6 protein identifiers: Rv1908c, Rv0242c, Rv3224, Rv0667, Rv1133c, Rv3036c

our_only_strong <- filter_on_max_peptides(data=our_all_filtered, n_peptides=10)

## Before filtering: 
##   Number of proteins: 930
##   Number of peptides: 3158
## 
## Percentage of peptides removed: 8.61%
## 
## After filtering: 
##   Number of proteins: 928
##   Number of peptides: 2886

our_only_minimum <- filter_on_min_peptides(data=our_only_strong, n_peptides=3)

## Before filtering:
##   Number of proteins: 928
##   Number of peptides: 2886
## 
## Percentage of peptides removed: 0%
## 
## After filtering:
##   Number of proteins: 895
##   Number of peptides: 2886

## I think these matrixes are probably smarter to use than the raw outmatrix from tric.
## But I am not a fan of rerwriting the sample column names.
matrix_prefix <- file.path("results", "swath2stats", ver)
if (!file.exists(matrix_prefix)) {
  dir.create(matrix_prefix)
}
protein_matrix_all <- write_matrix_proteins(
  our_s2s_exp, write.csv=TRUE,
  filename=file.path(matrix_prefix, "our_protein_all.csv"))

## Protein overview matrix results/swath2stats/20180611/our_protein_all.csv written to working folder.

dim(protein_matrix_all)

## [1] 1110   24

protein_matrix_mscore <- write_matrix_proteins(
  our_mscore_filtered, write.csv=TRUE,
  filename=file.path(matrix_prefix, "our_protein_matrix_mscore.csv"))

## Protein overview matrix results/swath2stats/20180611/our_protein_matrix_mscore.csv written to working folder.

dim(protein_matrix_mscore)

## [1] 945  24

peptide_matrix_mscore <- write_matrix_peptides(
  our_mscore_filtered, write.csv=TRUE,
  filename=file.path(matrix_prefix, "our_peptide_matrix_mscore.csv"))

## Peptide overview matrix results/swath2stats/20180611/our_peptide_matrix_mscore.csv written to working folder.

dim(peptide_matrix_mscore)

## [1] 3203   24

protein_matrix_minimum <- write_matrix_proteins(
  our_only_minimum, write.csv=TRUE,
  filename=file.path(matrix_prefix, "our_protein_matrix_minimum.csv"))

## Protein overview matrix results/swath2stats/20180611/our_protein_matrix_minimum.csv written to working folder.

dim(protein_matrix_minimum)

## [1] 895  24

peptide_matrix_minimum <- write_matrix_peptides(
  our_only_minimum, write.csv=TRUE,
  filename=file.path(matrix_prefix, "our_peptide_matrix_minimum.csv"))

## Peptide overview matrix results/swath2stats/20180611/our_peptide_matrix_minimum.csv written to working folder.

dim(peptide_matrix_minimum)

## [1] 38491    24

rt_cor <- plot_correlation_between_samples(
  our_only_minimum, column.values="intensity", fun.aggregate=sum, size=2)

## I have no effing clue what this plot means.
variation <- plot_variation(our_only_minimum, fun.aggregate=sum)

cols <- colnames(our_only_minimum)
our_disaggregated <- disaggregate(our_only_minimum, all.columns=TRUE)

## The library contains6transitions per precursor.
## The data table was transformed into a table containing one row per transition.

our_msstats_input <- convert_for_msstats(our_disaggregated)

## One or several columns required by MSstats were not in the data. The columns were created and filled with NAs.
## Missing columns: fragmention, productcharge, isotopelabeltype

## isotopelabeltype was filled with light.

3.2 Run filters on the tuberculist

Repeat the above stanza using the tuberculist data.

## Get correlations on a sample by sample basis
pp(file="images/20180611_tb_swath2stats_sample_cor.png")

## Going to write the image to: images/20180611_tb_swath2stats_sample_cor.png when dev.off() is called.

sample_cor <- plot_correlation_between_samples(tb_s2s_exp, size=2,
                                               fun.aggregate=sum,
                                               column.values="intensity")
dev.off()

## png 
##   2

sample_cond_rep_cor <- plot_correlation_between_samples(tb_s2s_exp, size=2,
                                                        comparison=transition_group_id ~
                                                          condition + bioreplicate + run,
                                                        fun.aggregate=sum,
                                                        column.values="intensity")

decoy_lists <- assess_decoy_rate(tb_s2s_exp)

## Number of non-decoy peptides: 23975
## Number of decoy peptides: 1762
## Decoy rate: 0.0735

## This seems a bit high to me, yesno?
fdr_overall <- assess_fdr_overall(tb_s2s_exp, output="Rconsole", plot=TRUE)

byrun_fdr <- assess_fdr_byrun(tb_s2s_exp, FFT=0.7, plot=TRUE, output="Rconsole")

## The average FDR by run on assay level is 0.008

## The average FDR by run on peptide level is 0.009

## The average FDR by run on protein level is 0.044

chosen_mscore <- mscore4assayfdr(tb_s2s_exp, FFT=0.7, fdr_target=0.02)

## Target assay FDR: 0.02

## Required overall m-score cutoff: 0.0056234
## achieving assay FDR: 0.019

prot_score <- mscore4protfdr(tb_s2s_exp, FFT=0.7, fdr_target=0.02)

## Target protein FDR: 0.02

## Required overall m-score cutoff: 0.00070795
## achieving protein FDR: 0.0198

tb_mscore_filtered <- filter_mscore(tb_s2s_exp, chosen_mscore)

## Original dimension: 252491, new dimension: 231348, difference: 21143.

tb_freq_mscore <- filter_mscore_freqobs(tb_s2s_exp, 0.01, 0.8, rm.decoy=FALSE)

## Peptides need to have been quantified in more conditions than: 18.4 in order to pass this percentage-based threshold.

## Fraction of peptides selected: 1

## Original dimension: 256996, new dimension: 256996, difference: 0.

tb_data_filtered_fdr <- filter_mscore_fdr(tb_mscore_filtered, FFT=0.7,
                                           overall_protein_fdr_target=prot_score,
                                           upper_overall_peptide_fdr_limit=0.05)

## Target protein FDR: 0.000707945784384137

## Required overall m-score cutoff: 0.01
## achieving protein FDR: 0

## filter_mscore_fdr is filtering the data...

## finding m-score cutoff to achieve desired protein FDR in protein master list..

## finding m-score cutoff to achieve desired global peptide FDR..

## Target peptide FDR: 0.05

## Required overall m-score cutoff: 0.01

## Achieving peptide FDR: 0

## Proteins selected: 
## Total proteins selected: 3087
## Final target proteins: 3087
## Final decoy proteins: 0

## Peptides mapping to these protein entries selected:
## Total mapping peptides: 22671
## Final target peptides: 22671
## Final decoy peptides: 0

## Total peptides selected from:
## Total peptides: 22671
## Final target peptides: 22671
## Final decoy peptides: 0

## Individual run FDR quality of the peptides was not calculated
## as not every run contains a decoy.

## The decoys have been removed from the returned data.

tb_only_proteotypic <- filter_proteotypic_peptides(tb_data_filtered_fdr)

## Number of proteins detected: 3101
## Protein identifiers: Rv1270c, Rv0161, Rv1613, Rv2540c, Rv0440, Rv1327c

## Number of proteins detected that are supported by a proteotypic peptide: 2964
## Number of proteotypic peptides detected: 22510

tb_all_filtered <- filter_all_peptides(tb_only_proteotypic)

## Number of proteins detected: 2966

## First 6 protein identifiers: Rv1270c, Rv0161, Rv1613, Rv2540c, Rv0440, Rv1327c

tb_only_strong <- filter_on_max_peptides(data=tb_all_filtered, n_peptides=10)

## Before filtering: 
##   Number of proteins: 2964
##   Number of peptides: 22510
## 
## Percentage of peptides removed: 23.15%
## 
## After filtering: 
##   Number of proteins: 2936
##   Number of peptides: 17299

tb_only_minimum <- filter_on_min_peptides(data=tb_only_strong, n_peptides=3)

## Before filtering:
##   Number of proteins: 2936
##   Number of peptides: 17299
## 
## Percentage of peptides removed: 0.03%
## 
## After filtering:
##   Number of proteins: 2679
##   Number of peptides: 17293

## I think these matrixes are probably smarter to use than the raw outmatrix from tric.
## But I am not a fan of rerwriting the sample column names.
matrix_prefix <- file.path("results", "swath2stats", ver)
if (!file.exists(matrix_prefix)) {
  dir.create(matrix_prefix)
}
protein_matrix_all <- write_matrix_proteins(
  tb_s2s_exp, write.csv=TRUE,
  filename=file.path(matrix_prefix, "tb_protein_all.csv"))

## Protein overview matrix results/swath2stats/20180611/tb_protein_all.csv written to working folder.

dim(protein_matrix_all)

## [1] 4570   24

protein_matrix_mscore <- write_matrix_proteins(
  tb_mscore_filtered, write.csv=TRUE,
  filename=file.path(matrix_prefix, "tb_protein_matrix_mscore.csv"))

## Protein overview matrix results/swath2stats/20180611/tb_protein_matrix_mscore.csv written to working folder.

dim(protein_matrix_mscore)

## [1] 3087   24

peptide_matrix_mscore <- write_matrix_peptides(
  tb_mscore_filtered, write.csv=TRUE,
  filename=file.path(matrix_prefix, "tb_peptide_matrix_mscore.csv"))

## Peptide overview matrix results/swath2stats/20180611/tb_peptide_matrix_mscore.csv written to working folder.

dim(peptide_matrix_mscore)

## [1] 22671    24

protein_matrix_minimum <- write_matrix_proteins(
  tb_only_minimum, write.csv=TRUE,
  filename=file.path(matrix_prefix, "tb_protein_matrix_minimum.csv"))

## Protein overview matrix results/swath2stats/20180611/tb_protein_matrix_minimum.csv written to working folder.

dim(protein_matrix_minimum)

## [1] 2679   24

peptide_matrix_minimum <- write_matrix_peptides(
  tb_only_minimum, write.csv=TRUE,
  filename=file.path(matrix_prefix, "tb_peptide_matrix_minimum.csv"))

## Peptide overview matrix results/swath2stats/20180611/tb_peptide_matrix_minimum.csv written to working folder.

dim(peptide_matrix_minimum)

## [1] 164565     24

rt_cor <- plot_correlation_between_samples(size=2,
  tb_only_minimum, column.values="intensity", fun.aggregate=sum)

## I have no effing clue what this plot means.
variation <- plot_variation(tb_only_minimum, fun.aggregate=sum)

cols <- colnames(tb_only_minimum)
tb_disaggregated <- disaggregate(tb_only_minimum, all.columns=TRUE)

## The library contains5transitions per precursor.
## The data table was transformed into a table containing one row per transition.

tb_msstats_input <- convert_for_msstats(tb_disaggregated)

## One or several columns required by MSstats were not in the data. The columns were created and filled with NAs.
## Missing columns: fragmention, productcharge, isotopelabeltype

## isotopelabeltype was filled with light.

3.3 Some new plots

In response to some interesting queries from Yan, I made a few little functions which query and plot data from the scored data provided by openswath/pyprophet. Let us look at their results here.

3.3.1 Using the comet data

our_pyprophet_fun <- sm(extract_pyprophet_data(metadata="sample_sheets/Mtb_dia_samples.xlsx"))

our_mass_plot <- sm(plot_pyprophet_distribution(our_pyprophet_fun, column="mass"))
our_mass_plot[["boxplot"]]

our_mass_plot[["density"]]

our_deltart_plot_all <- sm(plot_pyprophet_distribution(
  our_pyprophet_fun, column="delta_rt"))
our_deltart_plot_all[["boxplot"]]

our_deltart_plot_real <- sm(plot_pyprophet_distribution(
  our_pyprophet_fun,
  column="delta_rt", keep_decoys=FALSE))
our_deltart_plot_real[["boxplot"]]

our_deltart_plot_decoys <- sm(plot_pyprophet_distribution(
  our_pyprophet_fun,
  column="delta_rt", keep_real=FALSE))
our_deltart_plot_decoys[["boxplot"]]

our_deltart_plot_decoys[["density"]]

our_widthvsmass <- sm(plot_pyprophet_data(our_pyprophet_fun, legend=FALSE))
our_widthvsmass$plot

3.3.2 The tuberculist data

tb_pyprophet_fun <- sm(extract_pyprophet_data(metadata="sample_sheets/Mtb_dia_samples.xlsx",
                                              pyprophet_column="tuberculistscored"))

tb_mass_plot <- sm(plot_pyprophet_distribution(tb_pyprophet_fun, column="mass"))
tb_mass_plot[["boxplot"]]

tb_mass_plot[["density"]]

tb_deltart_plot_all <- sm(plot_pyprophet_distribution(
  tb_pyprophet_fun, column="delta_rt"))
tb_deltart_plot_all[["boxplot"]]

tb_deltart_plot_real <- sm(plot_pyprophet_distribution(
  tb_pyprophet_fun,
  column="delta_rt", keep_decoys=FALSE))
tb_deltart_plot_real[["boxplot"]]

tb_deltart_plot_decoys <- sm(plot_pyprophet_distribution(
  tb_pyprophet_fun,
  column="delta_rt", keep_real=FALSE))
tb_deltart_plot_decoys[["boxplot"]]

tb_deltart_plot_decoys[["density"]]

tb_widthvsmass <- sm(plot_pyprophet_data(tb_pyprophet_fun, legend=FALSE))
tb_widthvsmass$plot

3.4 MSstats

msstats.org seems to provide a complete solution for performing reasonable metrics of this data.

I am currently reading: http://msstats.org/wp-content/uploads/2017/01/MSstats_v3.7.3_manual.pdf

I made some moderately intrusive changes to MSstats to make it clearer, as well.

3.4.1 Comet derived msstats

devtools::load_all("~/scratch/git/MSstats")

## Loading MSstats

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

## Warning: character(0)

our_msstats_quant <- sm(dataProcess(our_msstats_input))
our_msstats_plots <- sm(dataProcessPlots(our_msstats_quant, type="QCPLOT"))

my_levels <- levels(as.factor(our_msstats_input$condition))
my_levels

## [1] "comp_cf"     "comp_whole"  "delta_cf"    "delta_whole" "wt_cf"      
## [6] "wt_whole"

comparisons <- ghetto_contrast_matrix(
  numerators=c("wt_cf", "delta_cf", "comp_cf",
               "delta_cf", "comp_cf", "delta_whole",
               "comp_whole"),
  denominators=c("wt_whole", "delta_whole", "comp_whole",
                 "wt_cf", "wt_cf", "wt_whole",
                 "wt_whole"))
our_results <- list()
for (c in 1:length(rownames(comparisons))) {
  name <- rownames(comparisons)[c]
  message("Starting ", name)
  comp <- comparisons[c, ]
  comparison <- t(as.matrix(comp))
  rownames(comparison) <- name
  our_results[name] <- sm(MSstats::groupComparison(contrast.matrix=comparison,
                                                   data=our_msstats_quant))
}

## Starting wt_cf_vs_wt_whole

## Warning in our_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting delta_cf_vs_delta_whole

## Warning in our_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting comp_cf_vs_comp_whole

## Warning in our_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting delta_cf_vs_wt_cf

## Warning in our_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting comp_cf_vs_wt_cf

## Warning in our_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting delta_whole_vs_wt_whole

## Warning in our_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting comp_whole_vs_wt_whole

## Warning in our_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

3.4.2 Tuberculist derived msstats

tb_msstats_quant <- sm(dataProcess(tb_msstats_input))
##tb_msstats_plots <- sm(dataProcessPlots(tb_msstats_quant, type="QCPLOT"))

my_levels <- levels(as.factor(tb_msstats_input$condition))
my_levels

## [1] "comp_cf"     "comp_whole"  "delta_cf"    "delta_whole" "wt_cf"      
## [6] "wt_whole"

comparisons <- ghetto_contrast_matrix(
  numerators=c("wt_cf", "delta_cf", "comp_cf",
               "delta_cf", "comp_cf", "delta_whole",
               "comp_whole"),
  denominators=c("wt_whole", "delta_whole", "comp_whole",
                 "wt_cf", "wt_cf", "wt_whole",
                 "wt_whole"))
tb_results <- list()
for (c in 1:length(rownames(comparisons))) {
  name <- rownames(comparisons)[c]
  message("Starting ", name)
  comp <- comparisons[c, ]
  comparison <- t(as.matrix(comp))
  rownames(comparison) <- name
  tb_results[name] <- sm(MSstats::groupComparison(contrast.matrix=comparison,
                                                   data=tb_msstats_quant))
}

## Starting wt_cf_vs_wt_whole

## Warning in tb_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting delta_cf_vs_delta_whole

## Warning in tb_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting comp_cf_vs_comp_whole

## Warning in tb_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting delta_cf_vs_wt_cf

## Warning in tb_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting comp_cf_vs_wt_cf

## Warning in tb_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting delta_whole_vs_wt_whole

## Warning in tb_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

## Starting comp_whole_vs_wt_whole

## Warning in tb_results[name] <- sm(MSstats::groupComparison(contrast.matrix
## = comparison, : number of items to replace is not a multiple of replacement
## length

3.4.3 P/PE protein QC plots for Yan

Yan asked for the p/pe protein qc plots. ok. I changed the dataProcessPlots to return something useful, so that should be possible now.

pe_genes <- read.table("reference/annotated_pe_genes.txt")[[1]]

## Unfortunately, the names did not get set in my changed version of dataProcessPlots...
plotlst <- our_msstats_plots$QCPLOT
available_plots <- gsub(pattern="^1/", replacement="",
                        x=levels(our_msstats_quant$ProcessedData$PROTEIN))
names(plotlst) <- available_plots

pe_in_avail_idx <- pe_genes %in% available_plots
pe_in_avail <- pe_genes[pe_in_avail_idx]
pe_plots <- plotlst[pe_in_avail]
pdf(file="pe_qc_plots.pdf")
for (p in 1:length(pe_plots)) {
  plot(pe_plots[[p]])
}

## Warning: Removed 132 rows containing non-finite values (stat_boxplot).

## Warning: Removed 220 rows containing non-finite values (stat_boxplot).

## Warning: Removed 418 rows containing non-finite values (stat_boxplot).

## Warning: Removed 330 rows containing non-finite values (stat_boxplot).

## Warning: Removed 132 rows containing non-finite values (stat_boxplot).

## Warning: Removed 1452 rows containing non-finite values (stat_boxplot).

## Warning: Removed 1052 rows containing non-finite values (stat_boxplot).

## Warning: Removed 858 rows containing non-finite values (stat_boxplot).

## Warning: Removed 264 rows containing non-finite values (stat_boxplot).

## Warning: Removed 154 rows containing non-finite values (stat_boxplot).

## Warning: Removed 132 rows containing non-finite values (stat_boxplot).

## Warning: Removed 352 rows containing non-finite values (stat_boxplot).

## Warning: Removed 616 rows containing non-finite values (stat_boxplot).

## Warning: Removed 638 rows containing non-finite values (stat_boxplot).

## Warning: Removed 550 rows containing non-finite values (stat_boxplot).

## Warning: Removed 440 rows containing non-finite values (stat_boxplot).

## Warning: Removed 814 rows containing non-finite values (stat_boxplot).

## Warning: Removed 462 rows containing non-finite values (stat_boxplot).

## Warning: Removed 242 rows containing non-finite values (stat_boxplot).

dev.off()

## png 
##   2

length(pe_plots)

## [1] 19

4 Create hpgltools expressionset

Since I am not certain I understand these data, I will take the intensities from SWATH2stats, metadata, and annotation data; attempt to create a ‘normal’ expressionset; poke at it to see what I can learn.

4.1 Massaging the metadata

I want to use the same metadata as were used for MSstats. It has a few important differences from the requirements of hpgltools: pretty much only that I do not allow rownames/sampleIDs to start with a number.

4.2 Massaging the intensity matrix

I do not want the \1 before the protein names, I already merged them into one entry per gene vis SWATH2stats.

4.2.1 Make matrices with our libraries

our_prot_mtrx <- read.csv(file.path("results", "swath2stats", ver, "our_protein_matrix_minimum.csv"))
rownames(our_prot_mtrx) <- gsub(pattern="^1\\/", replacement="", x=our_prot_mtrx[["proteinname"]])
our_prot_mtrx <- our_prot_mtrx[, -1]
## Important question: Did SWATH2stats reorder my data?
colnames(our_prot_mtrx) <- gsub(pattern="^X", replacement="s", x=colnames(our_prot_mtrx))

4.2.2 Make matrices with the tb libraries

tb_prot_mtrx <- read.csv(file.path("results", "swath2stats", ver, "tb_protein_matrix_minimum.csv"))
rownames(tb_prot_mtrx) <- gsub(pattern="^1\\/", replacement="", x=tb_prot_mtrx[["proteinname"]])
tb_prot_mtrx <- tb_prot_mtrx[, -1]
## Important question: Did SWATH2stats reorder my data?
colnames(tb_prot_mtrx) <- gsub(pattern="^X", replacement="s", x=colnames(tb_prot_mtrx))

4.3 Merge the pieces

Now we should have sufficient pieces to make an expressionset.

While here, I will also split the data into a cf and whole-cell pair of data structures.

4.3.1 Our data

## Drop the metadata not in the protein matrix:
## And ensure that they are the same order.
reordered <- colnames(our_prot_mtrx)
metadata <- sample_annot[reordered, ]

our_protein_expt <- sm(create_expt(metadata=metadata,
                                   count_dataframe=our_prot_mtrx,
                                   gene_info=mtb_annotations))

our_whole_expt <- subset_expt(our_protein_expt, subset="collectiontype=='whole'")

## There were 23, now there are 11 samples.

our_cf_expt <- subset_expt(our_protein_expt, subset="collectiontype=='cf'")

## There were 23, now there are 12 samples.

4.3.2 Tb data

## Drop the metadata not in the protein matrix:
## And ensure that they are the same order.
reordered <- colnames(tb_prot_mtrx)
metadata <- sample_annot[reordered, ]

tb_protein_expt <- sm(create_expt(metadata,
                                  count_dataframe=tb_prot_mtrx,
                                  gene_info=mtb_annotations))
tb_whole_expt <- subset_expt(tb_protein_expt, subset="collectiontype=='whole'")

## There were 23, now there are 11 samples.

tb_cf_expt <- subset_expt(tb_protein_expt, subset="collectiontype=='cf'")

## There were 23, now there are 12 samples.

4.4 Metrics of the full data set

4.4.1 Our libraries

our_protein_metrics <- sm(graph_metrics(our_protein_expt))

our_protein_norm <- sm(normalize_expt(our_protein_expt, transform="log2", convert="cpm",
                                      norm="quant", filter=TRUE))
our_protein_norm_metrics <- sm(graph_metrics(our_protein_norm))

our_protein_fsva <- sm(normalize_expt(our_protein_expt, transform="log2",
                                   batch="fsva", filter=TRUE))
our_protein_fsva_metrics <- sm(graph_metrics(our_protein_fsva))

4.4.2 Tb libraries

tb_protein_metrics <- sm(graph_metrics(tb_protein_expt))

tb_protein_norm <- sm(normalize_expt(tb_protein_expt, transform="log2", convert="cpm",
                                      norm="quant", filter=TRUE))
tb_protein_norm_metrics <- sm(graph_metrics(tb_protein_norm))

tb_protein_fsva <- sm(normalize_expt(tb_protein_expt, transform="log2", convert="cpm",
                                      batch="fsva", filter=TRUE))
tb_protein_fsva_metrics <- sm(graph_metrics(tb_protein_fsva))

4.5 Metrics of the whole-cell data set

4.5.1 Our libraries

our_whole_metrics <- sm(graph_metrics(our_whole_expt))

our_whole_norm <- sm(normalize_expt(our_whole_expt, transform="log2", convert="cpm",
                                    norm="quant", filter=TRUE))
our_whole_norm_metrics <- sm(graph_metrics(our_whole_norm))

our_whole_fsva <- sm(normalize_expt(our_whole_expt, transform="log2", convert="cpm",
                                    batch="fsva", filter=TRUE))
our_whole_fsva_metrics <- sm(graph_metrics(our_whole_fsva))

4.5.2 Tb libraries

tb_whole_metrics <- sm(graph_metrics(tb_whole_expt))

tb_whole_norm <- sm(normalize_expt(tb_whole_expt, transform="log2", convert="cpm",
                                    norm="quant", filter=TRUE))
tb_whole_norm_metrics <- sm(graph_metrics(tb_whole_norm))

tb_whole_fsva <- sm(normalize_expt(tb_whole_expt, transform="log2", convert="cpm",
                                    batch="fsva", filter=TRUE))
tb_whole_fsva_metrics <- sm(graph_metrics(tb_whole_fsva))

4.6 Metrics of the filtrate data set

4.6.1 Our libraries

our_cf_metrics <- sm(graph_metrics(our_cf_expt))

our_cf_norm <- sm(normalize_expt(our_cf_expt, transform="log2", convert="cpm",
                                 norm="quant", filter=TRUE))
our_cf_norm_metrics <- sm(graph_metrics(our_cf_norm))

our_cf_fsva <- sm(normalize_expt(our_cf_expt, transform="log2", convert="cpm",
                                 batch="fsva", filter=TRUE))
our_cf_fsva_metrics <- sm(graph_metrics(our_cf_fsva))

4.6.2 Tb libraries

tb_cf_metrics <- sm(graph_metrics(tb_cf_expt))

tb_cf_norm <- sm(normalize_expt(tb_cf_expt, transform="log2", convert="cpm",
                                 norm="quant", filter=TRUE))
tb_cf_norm_metrics <- sm(graph_metrics(tb_cf_norm))

tb_cf_fsva <- sm(normalize_expt(tb_cf_expt, transform="log2", convert="cpm",
                                 batch="fsva", filter=TRUE))
tb_cf_fsva_metrics <- sm(graph_metrics(tb_cf_fsva))

4.7 plot some metrics

4.7.1 Our libraries

pp(image=our_protein_metrics$libsize,
   file=file.path("images", paste0(ver, "_our_libsize.png")))

## Writing the image to: images/20180611_our_libsize.png and calling dev.off().

## It seems to me that the scale of the data is all within an order of magnitude or two.
## I cannot get used to these absurdly large numbers though.
pp(image=our_protein_norm_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_our_norm_pca.png")))

## Writing the image to: images/20180611_our_norm_pca.png and calling dev.off().

## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

## There appears to be a nice split in the data, however the un-assayable batch
## effect is a problem.
pp(image=our_protein_fsva_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_fsva_pca.png")))

## Writing the image to: images/20180611_fsva_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

## fsva seems to get some handle on the data, but I don't think we should rely
## upon it.
pp(image=our_protein_norm_metrics$corheat,
   file=file.path("images", paste0(ver, "_norm_corheat.png")))

## Writing the image to: images/20180611_norm_corheat.png and calling dev.off().

## Once again, the whole-cell/culture-filtrate split is very large.
pp(image=our_protein_metrics$density,
   file=file.path("images", paste0(ver, "_raw_density.png")))

## Writing the image to: images/20180611_raw_density.png and calling dev.off().

## There are two obvious distributions in the data, once again split between types.
pp(image=our_protein_metrics$boxplot,
   file=file.path("images", paste0(ver, "_boxplot.png")))

## Writing the image to: images/20180611_boxplot.png and calling dev.off().

## This recapitulates the previous plot.

pp(image=our_whole_metrics$libsize,
   file=file.path("images", paste0(ver, "_whole_libsize.png")))

## Writing the image to: images/20180611_whole_libsize.png and calling dev.off().

pp(image=our_whole_norm_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_whole_norm_pca.png")))

## Writing the image to: images/20180611_whole_norm_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

pp(image=our_whole_fsva_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_whole_fsva_pca.png")))

## Writing the image to: images/20180611_whole_fsva_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

pp(image=our_whole_norm_metrics$corheat,
   file=file.path("images", paste0(ver, "_whole_norm_corheat.png")))

## Writing the image to: images/20180611_whole_norm_corheat.png and calling dev.off().

pp(image=our_whole_metrics$density,
   file=file.path("images", paste0(ver, "_whole_raw_density.png")))

## Writing the image to: images/20180611_whole_raw_density.png and calling dev.off().

pp(image=our_whole_metrics$boxplot,
   file=file.path("images", paste0(ver, "_whole_boxplot.png")))

## Writing the image to: images/20180611_whole_boxplot.png and calling dev.off().

pp(image=our_cf_metrics$libsize,
   file=file.path("images", paste0(ver, "_libsize.png")))

## Writing the image to: images/20180611_libsize.png and calling dev.off().

pp(image=our_cf_norm_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_norm_pca.png")))

## Writing the image to: images/20180611_norm_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

pp(image=our_cf_fsva_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_fsva_pca.png")))

## Writing the image to: images/20180611_fsva_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

pp(image=our_cf_norm_metrics$corheat,
   file=file.path("images", paste0(ver, "_norm_corheat.png")))

## Writing the image to: images/20180611_norm_corheat.png and calling dev.off().

pp(image=our_cf_metrics$density,
   file=file.path("images", paste0(ver, "_raw_density.png")))

## Writing the image to: images/20180611_raw_density.png and calling dev.off().

pp(image=our_cf_metrics$boxplot,
   file=file.path("images", paste0(ver, "_boxplot.png")))

## Writing the image to: images/20180611_boxplot.png and calling dev.off().

4.7.2 Tb libraries

pp(image=tb_protein_metrics$libsize,
   file=file.path("images", paste0(ver, "_tb_libsize.png")))

## Writing the image to: images/20180611_tb_libsize.png and calling dev.off().

## It seems to me that the scale of the data is all within an order of magnitude or two.
## I cannot get used to these absurdly large numbers though.
pp(image=tb_protein_norm_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_tb_norm_pca.png")))

## Writing the image to: images/20180611_tb_norm_pca.png and calling dev.off().

## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

## There appears to be a nice split in the data, however the un-assayable batch
## effect is a problem.
pp(image=tb_protein_fsva_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_fsva_pca.png")))

## Writing the image to: images/20180611_fsva_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

## fsva seems to get some handle on the data, but I don't think we should rely
## upon it.
pp(image=tb_protein_norm_metrics$corheat,
   file=file.path("images", paste0(ver, "_norm_corheat.png")))

## Writing the image to: images/20180611_norm_corheat.png and calling dev.off().

## Once again, the whole-cell/culture-filtrate split is very large.
pp(image=tb_protein_metrics$density,
   file=file.path("images", paste0(ver, "_raw_density.png")))

## Writing the image to: images/20180611_raw_density.png and calling dev.off().

## There are two obvious distributions in the data, once again split between types.
pp(image=tb_protein_metrics$boxplot,
   file=file.path("images", paste0(ver, "_boxplot.png")))

## Writing the image to: images/20180611_boxplot.png and calling dev.off().

## This recapitulates the previous plot.

pp(image=tb_whole_metrics$libsize,
   file=file.path("images", paste0(ver, "_whole_libsize.png")))

## Writing the image to: images/20180611_whole_libsize.png and calling dev.off().

pp(image=tb_whole_norm_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_whole_norm_pca.png")))

## Writing the image to: images/20180611_whole_norm_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

pp(image=tb_whole_fsva_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_whole_fsva_pca.png")))

## Writing the image to: images/20180611_whole_fsva_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

pp(image=tb_whole_norm_metrics$corheat,
   file=file.path("images", paste0(ver, "_whole_norm_corheat.png")))

## Writing the image to: images/20180611_whole_norm_corheat.png and calling dev.off().

pp(image=tb_whole_metrics$density,
   file=file.path("images", paste0(ver, "_whole_raw_density.png")))

## Writing the image to: images/20180611_whole_raw_density.png and calling dev.off().

pp(image=tb_whole_metrics$boxplot,
   file=file.path("images", paste0(ver, "_whole_boxplot.png")))

## Writing the image to: images/20180611_whole_boxplot.png and calling dev.off().

pp(image=tb_cf_metrics$libsize,
   file=file.path("images", paste0(ver, "_libsize.png")))

## Writing the image to: images/20180611_libsize.png and calling dev.off().

pp(image=tb_cf_norm_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_norm_pca.png")))

## Writing the image to: images/20180611_norm_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

pp(image=tb_cf_fsva_metrics$pcaplot,
   file=file.path("images", paste0(ver, "_fsva_pca.png")))

## Writing the image to: images/20180611_fsva_pca.png and calling dev.off().
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse
## Too few points to calculate an ellipse

pp(image=tb_cf_norm_metrics$corheat,
   file=file.path("images", paste0(ver, "_norm_corheat.png")))

## Writing the image to: images/20180611_norm_corheat.png and calling dev.off().

pp(image=tb_cf_metrics$density,
   file=file.path("images", paste0(ver, "_raw_density.png")))

## Writing the image to: images/20180611_raw_density.png and calling dev.off().

pp(image=tb_cf_metrics$boxplot,
   file=file.path("images", paste0(ver, "_boxplot.png")))

## Writing the image to: images/20180611_boxplot.png and calling dev.off().

5 Attempt some quantification comparisons?

5.1 Our libraries

our_pairwise_filt <- sm(normalize_expt(our_protein_expt, filter=TRUE))
our_pairwise_comp <- sm(all_pairwise(our_pairwise_filt, model_batch="fsva", force=TRUE))

our_pairwise_nobatch <- sm(all_pairwise(our_pairwise_filt, model_batch=FALSE, force=TRUE))

5.2 Tb libraries

tb_pairwise_filt <- sm(normalize_expt(tb_protein_expt, filter=TRUE))
tb_pairwise_comp <- sm(all_pairwise(tb_pairwise_filt, model_batch="fsva", force=TRUE))

tb_pairwise_nobatch <- sm(all_pairwise(tb_pairwise_filt, model_batch=FALSE, force=TRUE))

5.2.1 Show a few metrics from the hpgltools pairwise comparisons

our_pairwise_comp$comparison$heat

tb_pairwise_comp$comparison$heat

6 For each msstats run, do a DE table

6.1 wt_cf vs wt_whole

6.1.1 Our libraries

keepers <- list(
  "wtcf_vs_wtwhole" = c("wt_cf", "wt_whole"))
droppers <- c("undefined")
names(droppers) <- "log2fc"

## Make sure to set the rownames so it will merge into the excel file.
set_name <- "wt_cf_vs_wt_whole"
rownames(our_results[[set_name]]) <- our_results[[set_name]][["Protein"]]
wtcf_wtwhole_tables <- sm(combine_de_tables(
  our_pairwise_comp, keepers=keepers, extra_annot=our_results[[set_name]],
  excludes=droppers,
  excel=paste0("excel/our_wtcf_vs_wtwhole_tables-v", ver, ".xlsx")))

wtcf_nobatch_wtwhole_tables <- sm(combine_de_tables(
  our_pairwise_nobatch, keepers=keepers, extra_annot=our_results[[set_name]],
  excludes=droppers,
  excel=paste0("excel/our_wtcf_vs_wtwhole_nobatch_tables-v", ver, ".xlsx")))

comp_table <- wtcf_wtwhole_tables$data[[set_name]]
cor.test(comp_table$log2fc, comp_table$deseq_logfc, method="spearman")

## Warning in cor.test.default(comp_table$log2fc, comp_table$deseq_logfc,
## method = "spearman"): Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  comp_table$log2fc and comp_table$deseq_logfc
## S = 6.7e+07, p-value <2e-16
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##    rho 
## 0.3996

comp_table <- wtcf_nobatch_wtwhole_tables$data[[set_name]]
cor.test(comp_table$log2fc, comp_table$deseq_logfc, method="spearman")

## Warning in cor.test.default(comp_table$log2fc, comp_table$deseq_logfc,
## method = "spearman"): Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  comp_table$log2fc and comp_table$deseq_logfc
## S = 5.3e+07, p-value <2e-16
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##    rho 
## 0.5197

6.1.2 Tb libraries

keepers <- list(
  "wtcf_vs_wtwhole" = c("wt_cf", "wt_whole"))
droppers <- c("undefined")
names(droppers) <- "log2fc"

## Make sure to set the rownames so it will merge into the excel file.
set_name <- "wt_cf_vs_wt_whole"
rownames(tb_results[[set_name]]) <- tb_results[[set_name]][["Protein"]]
wtcf_wtwhole_tables <- sm(combine_de_tables(
  tb_pairwise_comp, keepers=keepers, extra_annot=tb_results[[set_name]],
  excludes=droppers,
  excel=paste0("excel/tb_wtcf_vs_wtwhole_tables-v", ver, ".xlsx")))

wtcf_nobatch_wtwhole_tables <- sm(combine_de_tables(
  tb_pairwise_nobatch, keepers=keepers, extra_annot=tb_results[[set_name]],
  excludes=droppers,
  excel=paste0("excel/tb_wtcf_vs_wtwhole_nobatch_tables-v", ver, ".xlsx")))

comp_table <- wtcf_wtwhole_tables$data[[set_name]]
cor.test(comp_table$log2fc, comp_table$deseq_logfc, method="spearman")

## Warning in cor.test.default(comp_table$log2fc, comp_table$deseq_logfc,
## method = "spearman"): Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  comp_table$log2fc and comp_table$deseq_logfc
## S = 1.1e+09, p-value <2e-16
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##    rho 
## 0.6209

comp_table <- wtcf_nobatch_wtwhole_tables$data[[set_name]]
cor.test(comp_table$log2fc, comp_table$deseq_logfc, method="spearman")

## Warning in cor.test.default(comp_table$log2fc, comp_table$deseq_logfc,
## method = "spearman"): Cannot compute exact p-value with ties

## 
##  Spearman's rank correlation rho
## 
## data:  comp_table$log2fc and comp_table$deseq_logfc
## S = 9.5e+08, p-value <2e-16
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##    rho 
## 0.6581

6.2 delta_cf vs delta_whole

6.2.1 Our libraries

keepers <- list(
  "deltacf_vs_deltawhole" = c("delta_cf", "delta_whole"))
set_name <- "delta_cf_vs_delta_whole"
rownames(our_results[[set_name]]) <- our_results[[set_name]][["Protein"]]
deltacf_deltawhole_tables <- sm(combine_de_tables(
  our_pairwise_comp, keepers=keepers, extra_annot=our_results[[set_name]],
  excel=paste0("excel/our_deltacf_vs_deltawhole_tables-v", ver, ".xlsx")))

6.2.2 Tb libraries

keepers <- list(
  "deltacf_vs_deltawhole" = c("delta_cf", "delta_whole"))
set_name <- "delta_cf_vs_delta_whole"
rownames(tb_results[[set_name]]) <- tb_results[[set_name]][["Protein"]]
deltacf_deltawhole_tables <- sm(combine_de_tables(
  tb_pairwise_comp, keepers=keepers, extra_annot=tb_results[[set_name]],
  excel=paste0("excel/tb_deltacf_vs_deltawhole_tables-v", ver, ".xlsx")))

6.3 comp_cf vs comp_whole

6.3.1 Our libraries

keepers <- list(
  "compcf_vs_compwhole" = c("comp_cf", "comp_whole"))
set_name <- "comp_cf_vs_comp_whole"
rownames(our_results[[set_name]]) <- our_results[[set_name]][["Protein"]]
compcf_compwhole_tables <- sm(combine_de_tables(
  our_pairwise_comp, keepers=keepers, extra_annot=our_results[[set_name]],
  excel=paste0("excel/our_compcf_vs_compwhole_tables-v", ver, ".xlsx")))

6.3.2 Tb libraries

keepers <- list(
  "compcf_vs_compwhole" = c("comp_cf", "comp_whole"))
set_name <- "comp_cf_vs_comp_whole"
rownames(tb_results[[set_name]]) <- tb_results[[set_name]][["Protein"]]
compcf_compwhole_tables <- sm(combine_de_tables(
  tb_pairwise_comp, keepers=keepers, extra_annot=tb_results[[set_name]],
  excel=paste0("excel/tb_compcf_vs_compwhole_tables-v", ver, ".xlsx")))

6.4 delta_cf vs wt_cf

6.4.1 Our libraries

keepers <- list(
  "deltacf_vs_wtcf" = c("delta_cf", "wt_cf"))
set_name <- "delta_cf_vs_wt_cf"
rownames(our_results[[set_name]]) <- our_results[[set_name]][["Protein"]]
deltacf_wtcf_tables <- sm(combine_de_tables(
  our_pairwise_comp, keepers=keepers, extra_annot=our_results[[set_name]],
  excel=paste0("excel/our_deltacf_vs_wtcf_tables-v", ver, ".xlsx")))

6.4.2 Tb libraries

keepers <- list(
  "deltacf_vs_wtcf" = c("delta_cf", "wt_cf"))
set_name <- "delta_cf_vs_wt_cf"
rownames(tb_results[[set_name]]) <- tb_results[[set_name]][["Protein"]]
deltacf_wtcf_tables <- sm(combine_de_tables(
  tb_pairwise_comp, keepers=keepers, extra_annot=tb_results[[set_name]],
  excel=paste0("excel/tb_deltacf_vs_wtcf_tables-v", ver, ".xlsx")))

6.5 comp_cf vs wt_cf

6.5.1 Our libraries

keepers <- list(
  "compcf_vs_wtcf" = c("comp_cf", "wt_cf"))
set_name <- "comp_cf_vs_wt_cf"
rownames(our_results[[set_name]]) <- our_results[[set_name]][["Protein"]]
compcf_wtcf_tables <- sm(combine_de_tables(
  our_pairwise_comp, keepers=keepers, extra_annot=our_results[[set_name]],
  excel=paste0("excel/our_compcf_vs_wtcf_tables-v", ver, ".xlsx")))

6.5.2 Tb libraries

keepers <- list(
  "compcf_vs_wtcf" = c("comp_cf", "wt_cf"))
set_name <- "comp_cf_vs_wt_cf"
rownames(tb_results[[set_name]]) <- tb_results[[set_name]][["Protein"]]
compcf_wtcf_tables <- sm(combine_de_tables(
  tb_pairwise_comp, keepers=keepers, extra_annot=tb_results[[set_name]],
  excel=paste0("excel/tb_compcf_vs_wtcf_tables-v", ver, ".xlsx")))

6.6 delta_whole vs wt_whole

6.6.1 Our libraries

keepers <- list(
  "wtcf_vs_wtwhole" = c("wt_cf", "wt_whole"))
set_name <- "delta_whole_vs_wt_whole"
rownames(our_results[[set_name]]) <- our_results[[set_name]][["Protein"]]
wtcf_wtwhole_tables <- sm(combine_de_tables(
  our_pairwise_comp, keepers=keepers, extra_annot=our_results[[set_name]],
  excel=paste0("excel/our_deltawhole_vs_wtwhole_tables-v", ver, ".xlsx")))

6.6.2 Tb libraries

keepers <- list(
  "wtcf_vs_wtwhole" = c("wt_cf", "wt_whole"))
set_name <- "delta_whole_vs_wt_whole"
rownames(tb_results[[set_name]]) <- tb_results[[set_name]][["Protein"]]
wtcf_wtwhole_tables <- sm(combine_de_tables(
  tb_pairwise_comp, keepers=keepers, extra_annot=tb_results[[set_name]],
  excel=paste0("excel/tb_deltawhole_vs_wtwhole_tables-v", ver, ".xlsx")))

6.7 comp_whole vs wt_whole

6.7.1 Our libraries

keepers <- list(
  "compwhole_vs_wtwhole" = c("comp_whole", "wt_whole"))
set_name <- "comp_whole_vs_wt_whole"
rownames(our_results[[set_name]]) <- our_results[[set_name]][["Protein"]]
compwhole_wtwhole_tables <- sm(combine_de_tables(
  our_pairwise_comp, keepers=keepers, extra_annot=our_results[[set_name]],
  excel=paste0("excel/our_compwhole_vs_wtwhole_tables-v", ver, ".xlsx")))

6.7.2 Tb libraries

keepers <- list(
  "compwhole_vs_wtwhole" = c("comp_whole", "wt_whole"))
set_name <- "comp_whole_vs_wt_whole"
rownames(tb_results[[set_name]]) <- tb_results[[set_name]][["Protein"]]
compwhole_wtwhole_tables <- sm(combine_de_tables(
  tb_pairwise_comp, keepers=keepers, extra_annot=tb_results[[set_name]],
  excel=paste0("excel/tb_compwhole_vs_wtwhole_tables-v", ver, ".xlsx")))

our_table <- our_pairwise_comp$deseq$all_tables$wt_whole_vs_wt_cf
tb_table <- tb_pairwise_comp$deseq$all_tables$wt_whole_vs_wt_cf
our_order <- order(our_table[["limma_logfc"]])

## Error in order(our_table[["limma_logfc"]]): argument 1 is not a vector

rownames(tb_table) <- gsub(pattern="^(.*)_.*", replacement="\\1", x=rownames(tb_table))

test_table <- merge(our_table, tb_table, by="row.names")
cor.test(test_table$logFC.x, test_table$logFC.y)

## 
##  Pearson's product-moment correlation
## 
## data:  test_table$logFC.x and test_table$logFC.y
## t = 15, df = 890, p-value <2e-16
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.4021 0.5063
## sample estimates:
##    cor 
## 0.4558

7 TODO

2018-04-10: Make sure my invocations of SWATH2stats/MSstats are correct.

if (!isTRUE(get0("skip_load"))) {
  message(paste0("This is hpgltools commit: ", get_git_commit()))
  this_save <- paste0(gsub(pattern="\\.Rmd", replace="", x=rmd_file), "-v", ver, ".rda.xz")
  message(paste0("Saving to ", this_save))
  tmp <- sm(saveme(filename=this_save))
  pander::pander(sessionInfo())
}

LS0tCnRpdGxlOiAiTS50dWJlcmN1bG9zaXMgMjAxODA2MTE6IEFuYWx5emluZyBkYXRhIGZyb20gT3BlblN3YXRoV29ya0Zsb3cvVFJJQy4iCmF1dGhvcjogImF0YiBhYmVsZXdAZ21haWwuY29tIgpkYXRlOiAiYHIgU3lzLkRhdGUoKWAiCm91dHB1dDoKIGh0bWxfZG9jdW1lbnQ6CiAgY29kZV9kb3dubG9hZDogdHJ1ZQogIGNvZGVfZm9sZGluZzogc2hvdwogIGZpZ19jYXB0aW9uOiB0cnVlCiAgZmlnX2hlaWdodDogNwogIGZpZ193aWR0aDogNwogIGhpZ2hsaWdodDogZGVmYXVsdAogIGtlZXBfbWQ6IGZhbHNlCiAgbW9kZTogc2VsZmNvbnRhaW5lZAogIG51bWJlcl9zZWN0aW9uczogdHJ1ZQogIHNlbGZfY29udGFpbmVkOiB0cnVlCiAgdGhlbWU6IHJlYWRhYmxlCiAgdG9jOiB0cnVlCiAgdG9jX2Zsb2F0OgogICAgY29sbGFwc2VkOiBmYWxzZQogICAgc21vb3RoX3Njcm9sbDogZmFsc2UKLS0tCgo8c3R5bGU+CiAgYm9keSAubWFpbi1jb250YWluZXIgewogICAgbWF4LXdpZHRoOiAxNjAwcHg7CiAgfQo8L3N0eWxlPgoKYGBge3Igb3B0aW9ucywgaW5jbHVkZT1GQUxTRX0KaWYgKCFpc1RSVUUoZ2V0MCgic2tpcF9sb2FkIikpKSB7CiAgbGlicmFyeShocGdsdG9vbHMpCiAgdHQgPC0gZGV2dG9vbHM6OmxvYWRfYWxsKCJ+L2hwZ2x0b29scyIpCiAga25pdHI6Om9wdHNfa25pdCRzZXQocHJvZ3Jlc3M9VFJVRSwKICAgICAgICAgICAgICAgICAgICAgICB2ZXJib3NlPVRSVUUsCiAgICAgICAgICAgICAgICAgICAgICAgd2lkdGg9OTAsCiAgICAgICAgICAgICAgICAgICAgICAgZWNobz1UUlVFKQogIGtuaXRyOjpvcHRzX2NodW5rJHNldChlcnJvcj1UUlVFLAogICAgICAgICAgICAgICAgICAgICAgICBmaWcud2lkdGg9OCwKICAgICAgICAgICAgICAgICAgICAgICAgZmlnLmhlaWdodD04LAogICAgICAgICAgICAgICAgICAgICAgICBkcGk9OTYpCiAgb2xkX29wdGlvbnMgPC0gb3B0aW9ucyhkaWdpdHM9NCwKICAgICAgICAgICAgICAgICAgICAgICAgIHN0cmluZ3NBc0ZhY3RvcnM9RkFMU0UsCiAgICAgICAgICAgICAgICAgICAgICAgICBrbml0ci5kdXBsaWNhdGUubGFiZWw9ImFsbG93IikKICBnZ3Bsb3QyOjp0aGVtZV9zZXQoZ2dwbG90Mjo6dGhlbWVfYncoYmFzZV9zaXplPTEwKSkKICB2ZXIgPC0gIjIwMTgwNjExIgogIHByZXZpb3VzX2ZpbGUgPC0gIjAxX3ByZXByb2Nlc3NpbmdfY29tZXRfaGlnaHJlcy5SbWQiCgogIHRtcCA8LSB0cnkoc20obG9hZG1lKGZpbGVuYW1lPXBhc3RlMChnc3ViKHBhdHRlcm49IlxcLlJtZCIsIHJlcGxhY2U9IiIsIHg9cHJldmlvdXNfZmlsZSksICItdiIsIHZlciwgIi5yZGEueHoiKSkpKQogIHByZXZpb3VzX2ZpbGUgPC0gIjAxX2Fubm90YXRpb25fMjAxODA2MTEuUm1kIgogIHRtcCA8LSB0cnkoc20obG9hZG1lKGZpbGVuYW1lPXBhc3RlMChnc3ViKHBhdHRlcm49IlxcLlJtZCIsIHJlcGxhY2U9IiIsIHg9cHJldmlvdXNfZmlsZSksICItdiIsIHZlciwgIi5yZGEueHoiKSkpKQoKICBybWRfZmlsZSA8LSBwYXN0ZTAoIjAyX3N3YXRoMnN0YXRzXyIsIHZlciwgIi5SbWQiKQp9CmBgYAoKIyBUT0RPCgoqIFJlbW92ZSBNU3N0YXRzIGxvZ2dpbmcuIC0tIGRvbmUKKiBNYWtlIGV4cGxpY2l0IHNwZWFybWFuIGNvcnJlbGF0aW9ucyBiZXR3ZWVuIG1ldGhvZHMuIC0tIGRvbmUKICAgICogRG8gYm90aCBmb3IgYWxsIGRhdGEgYW5kIGZvciB0aGUgdG9wLTUwL2JvdHRvbS01MCAtLSBkb25lLCBidXQgd2VpcmQuCiogTWFrZSAxMDAlIGNlcnRhaW4gdGhhdCB0aGUgc2FtcGxlcyBhcmUgYW5ub3RhdGVkIGNvcnJlY3RseS4gLS0gZG9uZS4KKiBMaW1tYS9NU3N0YXRzL0VkZ2VSIHZlbm4gZGlzYWdyYW0gZm9yIGFsbCB1cCBpbiBDRgogICAqIFJlcGVhdCBpbiBhbGwgZG93biBDRgogICAqIFJlcGVhdCB3aXRoIGEgY3V0b2ZmLCB0b3AtNTAKKiBJbmRpdmlkdWFsIHBsb3RzIGZvciBQL1BFIHByb3RlaW5zLgoKQW5hbHl6aW5nIGRhdGEgZnJvbSBvcGVuTVMgYW5kIGZyaWVuZHMuCj09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PQoKSW4gcHJlcHJvY2Vzc2luZ19jb21ldF9oaWdocmVzLlJtZCwgSSB1c2VkIHRoZSBvcGVuTVMgdHV0b3JpYWxzIGFuZCBzdXBwbGVtZW50YWwKbWF0ZXJpYWxzIGZyb20gYSBjb3VwbGUgcGFwZXJzIHRvIGhvcGVmdWxseSBjb3JyZWN0bHkgcGVyZm9ybSB0aGUgdmFyaW91cwpwcmVwcm9jZXNzaW5nIHRhc2tzIHJlcXVpcmVkIHRvIGV4dHJhY3QgaW50ZW5zaXR5IGRhdGEgZnJvbSBESUEvU1dBVEgKdHJhbnNpdGlvbnMuCgpUaGUgZmluYWwgc3RlcHMgb2YgdGhhdCBwcm9jZXNzIGNvbWJpbmVkIHRoZSB0cmFuc2l0aW9uIGludGVuc2l0aWVzIGZyb20gZXZlcnkKc2FtcGxlIGludG8gYSBtZXRhZGF0YSBmaWxlIChyZXN1bHRzL3RyaWMvSENEX21ldGEudHN2KSwgYW4gaW50ZW5zaXR5IG1hdHJpeAoocmVzdWx0cy90cmljL0hDRF9vdXRtYXRyaXgudHN2KSwgYW5kIGEgZmVhdHVyZSBhbGlnbmVkIG91dHB1dCBtYXRyaXgKKHJlc3VsdHMvdHJpYy9hbGlnbmVkX2NvbWV0X0hDRC50c3YpLgoKTXkgcmVhZGluZyBvZiB0aGUgU1dBVEgyc3RhdHMgYW5kIE1Tc3RhdHMgc291cmNlIGNvZGUgc3VnZ2VzdHMgdG8gbWUgdGhhdCB0aGUKbG9nMihpbnRlbnNpdGllcykgb2YgdGhlIGZlYXR1cmUgYWxpZ25lZCBkYXRhIGFyZSBvdXIgZmluYWwgcHJveHkgZm9yIHByb3RlaW4KYWJ1bmRhbmNlLiAgQXQgZmlyc3QgZ2xhbmNlLCB0aGlzIHN1Z2dlc3RzIHRvIG1lIHRoYXQgdGhlc2UgZGF0YSBtaWdodCBmb2xsb3cgYQpkaXN0cmlidXRpb24gc2ltaWxhciB0byBSTkFTZXEgZGF0YSAobmVnYXRpdmUgYmlub21pYWwsIGJ1dCBwZXJoYXBzIHdpdGggYQpiaWdnZXIgdGFpbD8pLiAgSW4gYWRkaXRpb24sIGJ5IHRoZSB0aW1lIHdlIHVzZSB0cmljIG9uIHRoZSBkYXRhLCB3ZSBoYXZlIGEKY291bnQgbWF0cml4IGFuZCBzYW1wbGUgYW5ub3RhdGlvbiBkYXRhIGZyYW1lcyB3aGljaCBsb29rIHJlbWFya2FibHkgc2ltaWxhciB0bwp0aG9zZSB1c2VkIGluIGEgUk5BU2VxIGV4cHJlc3Npb25zZXQuICBJbmRlZWQsIGJ5IHRoZSBlbmQgb2YgdGhlIE1Tc3RhdHMKcHJvY2Vzc2luZywgaXQgY3JlYXRlcyBhIE1TblNldCBjbGFzcyBvZiBpdHMgb3duIHdoaWNoIHVzZXMgZkRhdGEvZXhwcnMvcERhdGEuCgpGb3IgdGhlIGN1cmlvdXMsIG15IHJlYXNvbmluZyBmb3Igc2F5aW5nIHRoYXQgdGhlIGxvZyBpbnRlbnNpdGllcyBhcmUgb3VyIHByb3h5CmZvciBhYnVuZGFuY2UgY29tZXMgZnJvbSBNU3N0YXRzL1IvRGF0YVByb2Nlc3MuUiBpbiBhIGNsYXVzZSB3aGljaCBsb29rcyBsaWtlOgoKYGBge3Igc25pcHBldCwgZXZhbD1GQUxTRX0KaWYgKGxvZ1RyYW5zID09IDIpIHsKICB3b3JrW1siQUJVTkRBTkNFIl1dIDwtIGxvZzIod29ya1tbIkFCVU5EQU5DRSJdXSkKfSBlbHNlIGlmIChsb2dUcmFucyA9PSAxMCkgewogIHdvcmtbWyJBQlVOREFOQ0UiXV0gPC0gbG9nMTAod29ya1tbIkFCVU5EQU5DRSJdXSkKfSBlbHNlIHsKICAjIyBBYm92ZSB0aGVyZSB3YXMgYSBjaGVjayBmb3Igb25seSBsb2cgMiBhbmQgMTAsIGJ1dCB3ZSBjYW4gZG8gZSBpZiB3ZSB3YW50LgogICMjIEkgbWlnaHQgZ28gYmFjayB1cCB0aGVyZSBhbmQgcmVtb3ZlIHRoYXQgY2hlY2suIExvbmcgbGl2ZSBlISAyLjcxODI4MiBydWxlcyEKICB3b3JrW1siQUJVTkRBTkNFIl1dIDwtIGxvZyh3b3JrW1siQUJVTkRBTkNFIl1dKSAvIGxvZyhsb2dUcmFucykKfQpgYGAKCihOb3RlOiBJIGFkZGVkIHRoZSBuYXR1cmFsIGxvZyB0byB0aGUgc2V0IG9mIGNvbmRpdGlvbnMsIGJ1dCBvdGhlcndpc2UgdGhlIGxvZ2ljCmlzIHVuY2hhbmdlZC4pCgpXaXRoIHRoYXQgaW4gbWluZCwgSSB3YW50IHRvIHVzZSBzb21lIHRvb2xzIHdpdGggd2hpY2ggSSBhbSBmYW1pbGlhciBpbiBvcmRlciB0bwp0cnkgdG8gdW5kZXJzdGFuZCB0aGVzZSBkYXRhLiAgVGhlcmVmb3JlIEkgd2lsbCBmaXJzdCBhdHRlbXB0IHRvIGNvZXJjZSBteSB0cmljCmFsaWduZWQgZGF0YSBhbmQgYW5ub3RhdGlvbnMgaW50byBhICdub3JtYWwnIGV4cHJlc3Npb25zZXQuICBUaGVuIEkgd2FudCB0byBkbwpzb21lIGRpYWdub3N0aWMgcGxvdHMgd2hpY2gsIGlmIEkgYW0gd3JvbmcgYW5kIHRoZXNlIGRpc3RyaWJ1dGlvbnMgYXJlIG5vdCBhcwpleHBlY3RlZCwgd2lsbCBiZSBjb25jZXB0dWFsbHkgaW5jb3JyZWN0IChJIGRvbid0IHlldCB0aGluayBJIGFtIHdyb25nKS4KCiMjIFNhbXBsZSBhbm5vdGF0aW9uIHZpYSBTV0FUSDJzdGF0cwoKSSBhbSB1c2luZyB0aGUgU1dBVEgyc3RhdHMgdmlnbmV0dGUgYXMgbXkgcHJpbWFyeSBzb3VyY2Ugb2YgaW5mb3JtYXRpb24uICBUaHVzIEkKc2VlIHRoYXQgaXQgdXNlcyB0aGUKT3BlblNXQVRIX1NNM19Hb2xkU3RhbmRhcmRBdXRvbWF0ZWRSZXN1bHRzX2h1bWFuX3BlYWtncm91cHMudHh0IHdoaWNoIGhhcyBhCmZvcm1hdCBuZWFybHkgaWRlbnRpY2FsIHRvIHRoZSB0cmljIG91dHB1dCBtYXRyaXguICBUaHVzIGZvciB0aGUgbW9tZW50IEkgd2lsbAphc3N1bWUgdGhhdCB0aGUgcHJvcGVyIGlucHV0IGZvciBTV0FUSDJzdGF0cyBpcyAncmVzdWx0cy90cmljL2NvbWV0X0hDRC50c3YnIGFuZApub3QgdGhlIG1ldGFkYXRhIG5vciBvdXRwdXQgbWF0cml4LiAgSW4gYWRkaXRpb24sIFNXQVRIMnN0YXRzIHByb3ZpZGVzIGEgc2ltaWxhcgpjb252ZXJzaW9uIGZ1bmN0aW9uIHdoaWNoIHRha2VzIHRoZSB0cmljIG91dHB1dCBhbmQgY29lcmNlcyBpdCBpbnRvIGEgc2ltaWxhcgptYXRyaXggYWZ0ZXIgcGVyZm9ybWluZyBpdHMgdmFyaW91cyBmaWx0ZXJpbmcgdGFza3MuICBUaGVyZWZvcmUgSSB3aWxsIHVzZSB0aGF0LgoKSSBrZWVwIGEgc2FtcGxlIHNoZWV0IG9mIGFsbCB0aGUgRElBIHNhbXBsZXMgdXNlZCBpbiB0aGlzIGFuYWx5c2lzIGluCidzYW1wbGVfc2hlZXRzL2RpYV9zYW1wbGVzLnhsc3gnICBJdCBzaG91bGQgY29udGFpbiBhbGwgdGhlIG90aGVyIHJlcXVpcmVkIGRhdGEKd2l0aCBvbmUgaW1wb3J0YW50IGNhdmVhdCwgSSByZW1vdmVkIDEgc2FtcGxlIGJ5ICdjb21tZW50aW5nJyBpdCAoZS5nLiBwcmVmaXhpbmcKaXQgd2l0aCAnIyMnIC0tIHdoaWNoIGlzIGFuIGFkbWl0dGVkbHkgZHVtYiB0aGluZyB0byBkbyBpbiBhbiBleGNlbCBmaWxlLgoKT25lIGxhc3QgY2F2ZWF0OiBJIGhhY2tlZCB0aGUgU1dBVEgyc3RhdHMgc2FtcGxlX2Fubm90YXRpb24oKSBmdW5jdGlvbiB0byBhZGQgYQpjb3VwbGUgY29sdW1ucyBpbiBhbiBhdHRlbXB0IHRvIG1ha2UgaXQgYSBsaXR0bGUgbW9yZSByb2J1c3Qgd2hlbiBmYWNlZCB3aXRoCnNhbXBsZSBzaGVldHMgd2l0aCBkaWZmZXJlbnRseSBuYW1lZCBjb2x1bW5zLgoKSW4gYWRkaXRpb24sIFNXQVRIMnN0YXRzIHByb3ZpZGVzIHNvbWUgbmljZSBmaWx0ZXJpbmcgYW5kIGNvbWJpbmF0aW9uCmZ1bmN0aW9ucyB3aGljaCBzaG91bGQgYmUgY29uc2lkZXJlZCB3aGVuIGdlbmVyYXRpbmcgdmFyaW91cyBleHByZXNzaW9uc2V0IGRhdGEKc3RydWN0dXJlcyBsYXRlci4KCiMjIyBDcmVhdGluZyBhIHN3YXRoMnN0YXRzIGV4cGVyaW1lbnQgdXNpbmcgb3VyIGNvbWV0LWRlcml2ZWQgbGlicmFyeSBkYXRhCgpgYGB7ciBvdXJfc3dhdGgyc3RhdHNfaW5pdGlhbH0Kb3VyX3RyaWNfZGF0YSA8LSByZWFkLmNzdigicmVzdWx0cy90cmljLzIwMTgwNjExL3dob2xlXzhtei9jb21ldF9IQ0QudHN2Iiwgc2VwPSJcdCIpCgpzYW1wbGVfYW5ub3QgPC0gb3Blbnhsc3g6OnJlYWQueGxzeCgic2FtcGxlX3NoZWV0cy9NdGJfZGlhX3NhbXBsZXMueGxzeCIpCnJvd25hbWVzKHNhbXBsZV9hbm5vdCkgPC0gbWFrZS5uYW1lcyhzYW1wbGVfYW5ub3RbWyJzYW1wbGVpZCJdXSwgdW5pcXVlPVRSVUUpCiMjIERyb3Agc2FtcGxlcyBzdGFydGluZyB3aXRoIGNvbW1lbnRzCmtlZXBfaWR4IDwtICEgZ3JlcGwocGF0dGVybj0iIyMiLCB4PXNhbXBsZV9hbm5vdFtbInNhbXBsZWlkIl1dKQpzYW1wbGVfYW5ub3QgPC0gc2FtcGxlX2Fubm90W2tlZXBfaWR4LCBdCmV4cHRfaWR4IDwtIHNhbXBsZV9hbm5vdFtbImV4cHRfaWQiXV0gPT0gIm1heTIwMTgiIHwgc2FtcGxlX2Fubm90W1siZXhwdF9pZCJdXSA9PSAibWFyMjAxOCIKZXhwdF9pZHhbaXMubmEoc2FtcGxlX2Fubm90W1siZXhwdF9pZCJdXSldIDwtIEZBTFNFCnNhbXBsZV9hbm5vdCA8LSBzYW1wbGVfYW5ub3RbZXhwdF9pZHgsIF0KbXpfaWR4IDwtIHNhbXBsZV9hbm5vdFtbIndpbmRvd3NpemUiXV0gPT0gIjgiCnNhbXBsZV9hbm5vdCA8LSBzYW1wbGVfYW5ub3RbbXpfaWR4LCBdCiMjIFNldCB0aGUgbXpYTUwgY29sdW1uIHRvIG1hdGNoIHRoZSBmaWxlbmFtZSBjb2x1bW4gaW4gdGhlIGRhdGEuCmxvYWRlZCA8LSBzbShkZXZ0b29sczo6bG9hZF9hbGwoIn4vc2NyYXRjaC9naXQvU1dBVEgyc3RhdHMiKSkKIyMgczJzLCBteSB3aXR0eSB3YXkgb2Ygc2hvcnRlbmluZyBTV0FUSDJzdGF0cy4uLgpvdXJfczJzX2V4cCA8LSBzYW1wbGVfYW5ub3RhdGlvbihkYXRhPW91cl90cmljX2RhdGEsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHNhbXBsZV9hbm5vdGF0aW9uPXNhbXBsZV9hbm5vdCwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgcnVuX2NvbHVtbj0icnVuIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZnVsbHBlcHRpZGVuYW1lX2NvbHVtbj0iZnVsbHVuaW1vZHBlcHRpZGVuYW1lIikKYGBgCgojIyMgQ3JlYXRpbmcgYSBzd2F0aDJzdGF0cyBleHBlcmltZW50IHVzaW5nIHRoZSB0dWJlcmN1bGlzdC1kZXJpdmVkIGxpYnJhcnkgZGF0YQoKVGhpcyBqdXN0IHJlcGVhdHMgdGhlIHByZXZpb3VzIHN0YW56YSwgYnV0IHVzZXMgdGhlIHR1YmVyY3VsaXN0CnRyYW5zaXRpb24tbGlicmFyeSBzd2F0aCBkYXRhIGFzIGl0cyBpbnB1dCByYXRoZXIgdGhhbiBvdXIgb3duIGNvbWV0CmRhdGEuCgpUaGVyZSBpcyBvbmUgaW1wb3J0YW50IGNhdmVhdCBpbiB0aGUgZm9sbG93aW5nIGJsb2NrOiBJIHVzZWQgYSByZWdleCB0byByZW1vdmUKdGhlIHNlY29uZCBoYWxmIG9mIGdlbmVJRF9nZW5lTmFtZSBzbyB0aGF0IGxhdGVyIHdoZW4gSSBtZXJnZSBpbiB0aGUgYW5ub3RhdGlvbgpkYXRhIEkgaGF2ZSBpdCB3aWxsIG1hdGNoLgoKYGBge3IgdGJfc3dhdGgyc3RhdHNfaW5pdGlhbH0KdGJfdHJpY19kYXRhIDwtIHJlYWQuY3N2KCJyZXN1bHRzL3RyaWMvMjAxODA2MTEvd2hvbGVfOG16X3R1YmVyY3VsaXN0L2NvbWV0X0hDRC50c3YiLCBzZXA9Ilx0IikKCnRiX3RyaWNfZGF0YVtbIlByb3RlaW5OYW1lIl1dIDwtIGdzdWIocGF0dGVybj0iXiguKilfLiokIiwgcmVwbGFjZW1lbnQ9IlxcMSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgeD10Yl90cmljX2RhdGFbWyJQcm90ZWluTmFtZSJdXSkKdGJfczJzX2V4cCA8LSBzYW1wbGVfYW5ub3RhdGlvbihkYXRhPXRiX3RyaWNfZGF0YSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBzYW1wbGVfYW5ub3RhdGlvbj1zYW1wbGVfYW5ub3QsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZnVsbHBlcHRpZGVuYW1lX2NvbHVtbj0iZnVsbHVuaW1vZHBlcHRpZGVuYW1lIikKYGBgCgpOb3cgSSBoYXZlIGEgY291cGxlIGRhdGEgc3RydWN0dXJlcyB3aGljaCBzaG91bGQgcHJvdmUgdXNlZnVsIGZvciB0aGUgbWV0cmljcwpwcm92aWRlZCBieSBTV0FUSDJzdGF0cywgTVNzdGF0cywgYW5kIG15IG93biBocGdsdG9vbHMuCgojIFNXQVRIMnN0YXRzIGNvbnRpbnVlZAoKVGhlIHZhcmlvdXMgbWV0cmljcyBhbmQgZmlsdGVycyBwcm92aWRlZCBieSBTV0FUSDJzdGF0cyBzZWVtIHF1aXRlIHJlYXNvbmFibGUgdG8KbWUuICBUaGUgb25seSB0aGluZyB0aGF0IHJlYWxseSBib3RoZXJzIG1lIGlzIHRoYXQgdGhleSBhcmUgYWxsIGNhc2Ugc2Vuc2l0aXZlCmFuZCBJIGZvdW5kIHRoYXQgdGhlIG1vc3QgcmVjZW50IHRyaWMgY2hhbmdlZCB0aGUgY2FwaXRhbGl6YXRpb24gb2YgYSBjb2x1bW4sCmNhdXNpbmcgdGhlc2UgdG8gYWxsIGZhbGwgZG93bi4gIFRoZXJlZm9yZSBJIHdlbnQgaW4gYW5kIG1hZGUgZXZlcnl0aGluZyBjYXNlCmluc2Vuc2l0aXZlIGluIGEgZmFzaGlvbiBzaW1pbGFyIHRvIHRoYXQgZG9uZSBieSBNU3N0YXRzIChleGNlcHQgSSBoYXRlIGNhcGl0YWwKbGV0dGVycywgc28gSSB1c2VkIHRvbG93ZXIoKSByYXRoZXIgdGhhbiB0b3VwcGVyKCkpLgoKIyMgUnVuIGZpbHRlcnMgb24gb3VyIGNvbWV0IGRhdGEKCmBgYHtyIG91cl9zd2F0aDJzdGF0c19wcm9jZXNzaW5nfQojIyBHZXQgY29ycmVsYXRpb25zIG9uIGEgc2FtcGxlIGJ5IHNhbXBsZSBiYXNpcwpwcChmaWxlPSJpbWFnZXMvMjAxODA2MTFfb3VyX3N3YXRoMnN0YXRzX3NhbXBsZV9jb3IucG5nIikKc2FtcGxlX2NvciA8LSBwbG90X2NvcnJlbGF0aW9uX2JldHdlZW5fc2FtcGxlcyhvdXJfczJzX2V4cCwgc2l6ZT0yLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGZ1bi5hZ2dyZWdhdGU9c3VtLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbHVtbi52YWx1ZXM9ImludGVuc2l0eSIpCmRldi5vZmYoKQpzYW1wbGVfY29uZF9yZXBfY29yIDwtIHBsb3RfY29ycmVsYXRpb25fYmV0d2Vlbl9zYW1wbGVzKG91cl9zMnNfZXhwLCBzaXplPTIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29tcGFyaXNvbj10cmFuc2l0aW9uX2dyb3VwX2lkIH4KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbmRpdGlvbiArIGJpb3JlcGxpY2F0ZSArIHJ1biwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBmdW4uYWdncmVnYXRlPXN1bSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb2x1bW4udmFsdWVzPSJpbnRlbnNpdHkiKQoKZGVjb3lfbGlzdHMgPC0gYXNzZXNzX2RlY295X3JhdGUob3VyX3Myc19leHApCiMjIFRoaXMgc2VlbXMgYSBiaXQgaGlnaCB0byBtZSwgeWVzbm8/CmZkcl9vdmVyYWxsIDwtIGFzc2Vzc19mZHJfb3ZlcmFsbChvdXJfczJzX2V4cCwgb3V0cHV0PSJSY29uc29sZSIsIHBsb3Q9VFJVRSkKCmJ5cnVuX2ZkciA8LSBhc3Nlc3NfZmRyX2J5cnVuKG91cl9zMnNfZXhwLCBGRlQ9MC43LCBwbG90PVRSVUUsIG91dHB1dD0iUmNvbnNvbGUiKQpjaG9zZW5fbXNjb3JlIDwtIG1zY29yZTRhc3NheWZkcihvdXJfczJzX2V4cCwgRkZUPTAuNywgZmRyX3RhcmdldD0wLjAyKQpwcm90X3Njb3JlIDwtIG1zY29yZTRwcm90ZmRyKG91cl9zMnNfZXhwLCBGRlQ9MC43LCBmZHJfdGFyZ2V0PTAuMDIpCgpvdXJfbXNjb3JlX2ZpbHRlcmVkIDwtIGZpbHRlcl9tc2NvcmUob3VyX3Myc19leHAsIGNob3Nlbl9tc2NvcmUpCm91cl9mcmVxX21zY29yZSA8LSBmaWx0ZXJfbXNjb3JlX2ZyZXFvYnMob3VyX3Myc19leHAsIDAuMDEsIDAuOCwgcm0uZGVjb3k9RkFMU0UpCm91cl9kYXRhX2ZpbHRlcmVkX2ZkciA8LSBmaWx0ZXJfbXNjb3JlX2ZkcihvdXJfbXNjb3JlX2ZpbHRlcmVkLCBGRlQ9MC43LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgb3ZlcmFsbF9wcm90ZWluX2Zkcl90YXJnZXQ9cHJvdF9zY29yZSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHVwcGVyX292ZXJhbGxfcGVwdGlkZV9mZHJfbGltaXQ9MC4wNSkKb3VyX29ubHlfcHJvdGVvdHlwaWMgPC0gZmlsdGVyX3Byb3Rlb3R5cGljX3BlcHRpZGVzKG91cl9kYXRhX2ZpbHRlcmVkX2ZkcikKb3VyX2FsbF9maWx0ZXJlZCA8LSBmaWx0ZXJfYWxsX3BlcHRpZGVzKG91cl9vbmx5X3Byb3Rlb3R5cGljKQpvdXJfb25seV9zdHJvbmcgPC0gZmlsdGVyX29uX21heF9wZXB0aWRlcyhkYXRhPW91cl9hbGxfZmlsdGVyZWQsIG5fcGVwdGlkZXM9MTApCm91cl9vbmx5X21pbmltdW0gPC0gZmlsdGVyX29uX21pbl9wZXB0aWRlcyhkYXRhPW91cl9vbmx5X3N0cm9uZywgbl9wZXB0aWRlcz0zKQoKIyMgSSB0aGluayB0aGVzZSBtYXRyaXhlcyBhcmUgcHJvYmFibHkgc21hcnRlciB0byB1c2UgdGhhbiB0aGUgcmF3IG91dG1hdHJpeCBmcm9tIHRyaWMuCiMjIEJ1dCBJIGFtIG5vdCBhIGZhbiBvZiByZXJ3cml0aW5nIHRoZSBzYW1wbGUgY29sdW1uIG5hbWVzLgptYXRyaXhfcHJlZml4IDwtIGZpbGUucGF0aCgicmVzdWx0cyIsICJzd2F0aDJzdGF0cyIsIHZlcikKaWYgKCFmaWxlLmV4aXN0cyhtYXRyaXhfcHJlZml4KSkgewogIGRpci5jcmVhdGUobWF0cml4X3ByZWZpeCkKfQpwcm90ZWluX21hdHJpeF9hbGwgPC0gd3JpdGVfbWF0cml4X3Byb3RlaW5zKAogIG91cl9zMnNfZXhwLCB3cml0ZS5jc3Y9VFJVRSwKICBmaWxlbmFtZT1maWxlLnBhdGgobWF0cml4X3ByZWZpeCwgIm91cl9wcm90ZWluX2FsbC5jc3YiKSkKZGltKHByb3RlaW5fbWF0cml4X2FsbCkKcHJvdGVpbl9tYXRyaXhfbXNjb3JlIDwtIHdyaXRlX21hdHJpeF9wcm90ZWlucygKICBvdXJfbXNjb3JlX2ZpbHRlcmVkLCB3cml0ZS5jc3Y9VFJVRSwKICBmaWxlbmFtZT1maWxlLnBhdGgobWF0cml4X3ByZWZpeCwgIm91cl9wcm90ZWluX21hdHJpeF9tc2NvcmUuY3N2IikpCmRpbShwcm90ZWluX21hdHJpeF9tc2NvcmUpCnBlcHRpZGVfbWF0cml4X21zY29yZSA8LSB3cml0ZV9tYXRyaXhfcGVwdGlkZXMoCiAgb3VyX21zY29yZV9maWx0ZXJlZCwgd3JpdGUuY3N2PVRSVUUsCiAgZmlsZW5hbWU9ZmlsZS5wYXRoKG1hdHJpeF9wcmVmaXgsICJvdXJfcGVwdGlkZV9tYXRyaXhfbXNjb3JlLmNzdiIpKQpkaW0ocGVwdGlkZV9tYXRyaXhfbXNjb3JlKQpwcm90ZWluX21hdHJpeF9taW5pbXVtIDwtIHdyaXRlX21hdHJpeF9wcm90ZWlucygKICBvdXJfb25seV9taW5pbXVtLCB3cml0ZS5jc3Y9VFJVRSwKICBmaWxlbmFtZT1maWxlLnBhdGgobWF0cml4X3ByZWZpeCwgIm91cl9wcm90ZWluX21hdHJpeF9taW5pbXVtLmNzdiIpKQpkaW0ocHJvdGVpbl9tYXRyaXhfbWluaW11bSkKcGVwdGlkZV9tYXRyaXhfbWluaW11bSA8LSB3cml0ZV9tYXRyaXhfcGVwdGlkZXMoCiAgb3VyX29ubHlfbWluaW11bSwgd3JpdGUuY3N2PVRSVUUsCiAgZmlsZW5hbWU9ZmlsZS5wYXRoKG1hdHJpeF9wcmVmaXgsICJvdXJfcGVwdGlkZV9tYXRyaXhfbWluaW11bS5jc3YiKSkKZGltKHBlcHRpZGVfbWF0cml4X21pbmltdW0pCgpydF9jb3IgPC0gcGxvdF9jb3JyZWxhdGlvbl9iZXR3ZWVuX3NhbXBsZXMoCiAgb3VyX29ubHlfbWluaW11bSwgY29sdW1uLnZhbHVlcz0iaW50ZW5zaXR5IiwgZnVuLmFnZ3JlZ2F0ZT1zdW0sIHNpemU9MikKIyMgSSBoYXZlIG5vIGVmZmluZyBjbHVlIHdoYXQgdGhpcyBwbG90IG1lYW5zLgp2YXJpYXRpb24gPC0gcGxvdF92YXJpYXRpb24ob3VyX29ubHlfbWluaW11bSwgZnVuLmFnZ3JlZ2F0ZT1zdW0pCgpjb2xzIDwtIGNvbG5hbWVzKG91cl9vbmx5X21pbmltdW0pCm91cl9kaXNhZ2dyZWdhdGVkIDwtIGRpc2FnZ3JlZ2F0ZShvdXJfb25seV9taW5pbXVtLCBhbGwuY29sdW1ucz1UUlVFKQpvdXJfbXNzdGF0c19pbnB1dCA8LSBjb252ZXJ0X2Zvcl9tc3N0YXRzKG91cl9kaXNhZ2dyZWdhdGVkKQpgYGAKCiMjIFJ1biBmaWx0ZXJzIG9uIHRoZSB0dWJlcmN1bGlzdAoKUmVwZWF0IHRoZSBhYm92ZSBzdGFuemEgdXNpbmcgdGhlIHR1YmVyY3VsaXN0IGRhdGEuCgpgYGB7ciB0Yl9zd2F0aDJzdGF0c19wcm9jZXNzaW5nfQojIyBHZXQgY29ycmVsYXRpb25zIG9uIGEgc2FtcGxlIGJ5IHNhbXBsZSBiYXNpcwpwcChmaWxlPSJpbWFnZXMvMjAxODA2MTFfdGJfc3dhdGgyc3RhdHNfc2FtcGxlX2Nvci5wbmciKQpzYW1wbGVfY29yIDwtIHBsb3RfY29ycmVsYXRpb25fYmV0d2Vlbl9zYW1wbGVzKHRiX3Myc19leHAsIHNpemU9MiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBmdW4uYWdncmVnYXRlPXN1bSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb2x1bW4udmFsdWVzPSJpbnRlbnNpdHkiKQpkZXYub2ZmKCkKc2FtcGxlX2NvbmRfcmVwX2NvciA8LSBwbG90X2NvcnJlbGF0aW9uX2JldHdlZW5fc2FtcGxlcyh0Yl9zMnNfZXhwLCBzaXplPTIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29tcGFyaXNvbj10cmFuc2l0aW9uX2dyb3VwX2lkIH4KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvbmRpdGlvbiArIGJpb3JlcGxpY2F0ZSArIHJ1biwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBmdW4uYWdncmVnYXRlPXN1bSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBjb2x1bW4udmFsdWVzPSJpbnRlbnNpdHkiKQoKZGVjb3lfbGlzdHMgPC0gYXNzZXNzX2RlY295X3JhdGUodGJfczJzX2V4cCkKIyMgVGhpcyBzZWVtcyBhIGJpdCBoaWdoIHRvIG1lLCB5ZXNubz8KZmRyX292ZXJhbGwgPC0gYXNzZXNzX2Zkcl9vdmVyYWxsKHRiX3Myc19leHAsIG91dHB1dD0iUmNvbnNvbGUiLCBwbG90PVRSVUUpCgpieXJ1bl9mZHIgPC0gYXNzZXNzX2Zkcl9ieXJ1bih0Yl9zMnNfZXhwLCBGRlQ9MC43LCBwbG90PVRSVUUsIG91dHB1dD0iUmNvbnNvbGUiKQpjaG9zZW5fbXNjb3JlIDwtIG1zY29yZTRhc3NheWZkcih0Yl9zMnNfZXhwLCBGRlQ9MC43LCBmZHJfdGFyZ2V0PTAuMDIpCnByb3Rfc2NvcmUgPC0gbXNjb3JlNHByb3RmZHIodGJfczJzX2V4cCwgRkZUPTAuNywgZmRyX3RhcmdldD0wLjAyKQoKdGJfbXNjb3JlX2ZpbHRlcmVkIDwtIGZpbHRlcl9tc2NvcmUodGJfczJzX2V4cCwgY2hvc2VuX21zY29yZSkKdGJfZnJlcV9tc2NvcmUgPC0gZmlsdGVyX21zY29yZV9mcmVxb2JzKHRiX3Myc19leHAsIDAuMDEsIDAuOCwgcm0uZGVjb3k9RkFMU0UpCnRiX2RhdGFfZmlsdGVyZWRfZmRyIDwtIGZpbHRlcl9tc2NvcmVfZmRyKHRiX21zY29yZV9maWx0ZXJlZCwgRkZUPTAuNywKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG92ZXJhbGxfcHJvdGVpbl9mZHJfdGFyZ2V0PXByb3Rfc2NvcmUsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB1cHBlcl9vdmVyYWxsX3BlcHRpZGVfZmRyX2xpbWl0PTAuMDUpCnRiX29ubHlfcHJvdGVvdHlwaWMgPC0gZmlsdGVyX3Byb3Rlb3R5cGljX3BlcHRpZGVzKHRiX2RhdGFfZmlsdGVyZWRfZmRyKQp0Yl9hbGxfZmlsdGVyZWQgPC0gZmlsdGVyX2FsbF9wZXB0aWRlcyh0Yl9vbmx5X3Byb3Rlb3R5cGljKQp0Yl9vbmx5X3N0cm9uZyA8LSBmaWx0ZXJfb25fbWF4X3BlcHRpZGVzKGRhdGE9dGJfYWxsX2ZpbHRlcmVkLCBuX3BlcHRpZGVzPTEwKQp0Yl9vbmx5X21pbmltdW0gPC0gZmlsdGVyX29uX21pbl9wZXB0aWRlcyhkYXRhPXRiX29ubHlfc3Ryb25nLCBuX3BlcHRpZGVzPTMpCgojIyBJIHRoaW5rIHRoZXNlIG1hdHJpeGVzIGFyZSBwcm9iYWJseSBzbWFydGVyIHRvIHVzZSB0aGFuIHRoZSByYXcgb3V0bWF0cml4IGZyb20gdHJpYy4KIyMgQnV0IEkgYW0gbm90IGEgZmFuIG9mIHJlcndyaXRpbmcgdGhlIHNhbXBsZSBjb2x1bW4gbmFtZXMuCm1hdHJpeF9wcmVmaXggPC0gZmlsZS5wYXRoKCJyZXN1bHRzIiwgInN3YXRoMnN0YXRzIiwgdmVyKQppZiAoIWZpbGUuZXhpc3RzKG1hdHJpeF9wcmVmaXgpKSB7CiAgZGlyLmNyZWF0ZShtYXRyaXhfcHJlZml4KQp9CnByb3RlaW5fbWF0cml4X2FsbCA8LSB3cml0ZV9tYXRyaXhfcHJvdGVpbnMoCiAgdGJfczJzX2V4cCwgd3JpdGUuY3N2PVRSVUUsCiAgZmlsZW5hbWU9ZmlsZS5wYXRoKG1hdHJpeF9wcmVmaXgsICJ0Yl9wcm90ZWluX2FsbC5jc3YiKSkKZGltKHByb3RlaW5fbWF0cml4X2FsbCkKcHJvdGVpbl9tYXRyaXhfbXNjb3JlIDwtIHdyaXRlX21hdHJpeF9wcm90ZWlucygKICB0Yl9tc2NvcmVfZmlsdGVyZWQsIHdyaXRlLmNzdj1UUlVFLAogIGZpbGVuYW1lPWZpbGUucGF0aChtYXRyaXhfcHJlZml4LCAidGJfcHJvdGVpbl9tYXRyaXhfbXNjb3JlLmNzdiIpKQpkaW0ocHJvdGVpbl9tYXRyaXhfbXNjb3JlKQpwZXB0aWRlX21hdHJpeF9tc2NvcmUgPC0gd3JpdGVfbWF0cml4X3BlcHRpZGVzKAogIHRiX21zY29yZV9maWx0ZXJlZCwgd3JpdGUuY3N2PVRSVUUsCiAgZmlsZW5hbWU9ZmlsZS5wYXRoKG1hdHJpeF9wcmVmaXgsICJ0Yl9wZXB0aWRlX21hdHJpeF9tc2NvcmUuY3N2IikpCmRpbShwZXB0aWRlX21hdHJpeF9tc2NvcmUpCnByb3RlaW5fbWF0cml4X21pbmltdW0gPC0gd3JpdGVfbWF0cml4X3Byb3RlaW5zKAogIHRiX29ubHlfbWluaW11bSwgd3JpdGUuY3N2PVRSVUUsCiAgZmlsZW5hbWU9ZmlsZS5wYXRoKG1hdHJpeF9wcmVmaXgsICJ0Yl9wcm90ZWluX21hdHJpeF9taW5pbXVtLmNzdiIpKQpkaW0ocHJvdGVpbl9tYXRyaXhfbWluaW11bSkKcGVwdGlkZV9tYXRyaXhfbWluaW11bSA8LSB3cml0ZV9tYXRyaXhfcGVwdGlkZXMoCiAgdGJfb25seV9taW5pbXVtLCB3cml0ZS5jc3Y9VFJVRSwKICBmaWxlbmFtZT1maWxlLnBhdGgobWF0cml4X3ByZWZpeCwgInRiX3BlcHRpZGVfbWF0cml4X21pbmltdW0uY3N2IikpCmRpbShwZXB0aWRlX21hdHJpeF9taW5pbXVtKQoKcnRfY29yIDwtIHBsb3RfY29ycmVsYXRpb25fYmV0d2Vlbl9zYW1wbGVzKHNpemU9MiwKICB0Yl9vbmx5X21pbmltdW0sIGNvbHVtbi52YWx1ZXM9ImludGVuc2l0eSIsIGZ1bi5hZ2dyZWdhdGU9c3VtKQojIyBJIGhhdmUgbm8gZWZmaW5nIGNsdWUgd2hhdCB0aGlzIHBsb3QgbWVhbnMuCnZhcmlhdGlvbiA8LSBwbG90X3ZhcmlhdGlvbih0Yl9vbmx5X21pbmltdW0sIGZ1bi5hZ2dyZWdhdGU9c3VtKQoKY29scyA8LSBjb2xuYW1lcyh0Yl9vbmx5X21pbmltdW0pCnRiX2Rpc2FnZ3JlZ2F0ZWQgPC0gZGlzYWdncmVnYXRlKHRiX29ubHlfbWluaW11bSwgYWxsLmNvbHVtbnM9VFJVRSkKdGJfbXNzdGF0c19pbnB1dCA8LSBjb252ZXJ0X2Zvcl9tc3N0YXRzKHRiX2Rpc2FnZ3JlZ2F0ZWQpCmBgYAoKIyMgU29tZSBuZXcgcGxvdHMKCkluIHJlc3BvbnNlIHRvIHNvbWUgaW50ZXJlc3RpbmcgcXVlcmllcyBmcm9tIFlhbiwgSSBtYWRlIGEgZmV3IGxpdHRsZSBmdW5jdGlvbnMKd2hpY2ggcXVlcnkgYW5kIHBsb3QgZGF0YSBmcm9tIHRoZSBzY29yZWQgZGF0YSBwcm92aWRlZCBieSBvcGVuc3dhdGgvcHlwcm9waGV0LgpMZXQgdXMgbG9vayBhdCB0aGVpciByZXN1bHRzIGhlcmUuCgojIyMgVXNpbmcgdGhlIGNvbWV0IGRhdGEKCmBgYHtyIG91cl9weXByb3BoZXRfcGxvdHN9Cm91cl9weXByb3BoZXRfZnVuIDwtIHNtKGV4dHJhY3RfcHlwcm9waGV0X2RhdGEobWV0YWRhdGE9InNhbXBsZV9zaGVldHMvTXRiX2RpYV9zYW1wbGVzLnhsc3giKSkKCm91cl9tYXNzX3Bsb3QgPC0gc20ocGxvdF9weXByb3BoZXRfZGlzdHJpYnV0aW9uKG91cl9weXByb3BoZXRfZnVuLCBjb2x1bW49Im1hc3MiKSkKb3VyX21hc3NfcGxvdFtbImJveHBsb3QiXV0Kb3VyX21hc3NfcGxvdFtbImRlbnNpdHkiXV0KCm91cl9kZWx0YXJ0X3Bsb3RfYWxsIDwtIHNtKHBsb3RfcHlwcm9waGV0X2Rpc3RyaWJ1dGlvbigKICBvdXJfcHlwcm9waGV0X2Z1biwgY29sdW1uPSJkZWx0YV9ydCIpKQpvdXJfZGVsdGFydF9wbG90X2FsbFtbImJveHBsb3QiXV0KCm91cl9kZWx0YXJ0X3Bsb3RfcmVhbCA8LSBzbShwbG90X3B5cHJvcGhldF9kaXN0cmlidXRpb24oCiAgb3VyX3B5cHJvcGhldF9mdW4sCiAgY29sdW1uPSJkZWx0YV9ydCIsIGtlZXBfZGVjb3lzPUZBTFNFKSkKb3VyX2RlbHRhcnRfcGxvdF9yZWFsW1siYm94cGxvdCJdXQoKb3VyX2RlbHRhcnRfcGxvdF9kZWNveXMgPC0gc20ocGxvdF9weXByb3BoZXRfZGlzdHJpYnV0aW9uKAogIG91cl9weXByb3BoZXRfZnVuLAogIGNvbHVtbj0iZGVsdGFfcnQiLCBrZWVwX3JlYWw9RkFMU0UpKQpvdXJfZGVsdGFydF9wbG90X2RlY295c1tbImJveHBsb3QiXV0Kb3VyX2RlbHRhcnRfcGxvdF9kZWNveXNbWyJkZW5zaXR5Il1dCgpvdXJfd2lkdGh2c21hc3MgPC0gc20ocGxvdF9weXByb3BoZXRfZGF0YShvdXJfcHlwcm9waGV0X2Z1biwgbGVnZW5kPUZBTFNFKSkKb3VyX3dpZHRodnNtYXNzJHBsb3QKYGBgCgojIyMgVGhlIHR1YmVyY3VsaXN0IGRhdGEKCmBgYHtyIHRiX3B5cHJvcGhldF9wbG90c30KdGJfcHlwcm9waGV0X2Z1biA8LSBzbShleHRyYWN0X3B5cHJvcGhldF9kYXRhKG1ldGFkYXRhPSJzYW1wbGVfc2hlZXRzL010Yl9kaWFfc2FtcGxlcy54bHN4IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHB5cHJvcGhldF9jb2x1bW49InR1YmVyY3VsaXN0c2NvcmVkIikpCgp0Yl9tYXNzX3Bsb3QgPC0gc20ocGxvdF9weXByb3BoZXRfZGlzdHJpYnV0aW9uKHRiX3B5cHJvcGhldF9mdW4sIGNvbHVtbj0ibWFzcyIpKQp0Yl9tYXNzX3Bsb3RbWyJib3hwbG90Il1dCnRiX21hc3NfcGxvdFtbImRlbnNpdHkiXV0KCnRiX2RlbHRhcnRfcGxvdF9hbGwgPC0gc20ocGxvdF9weXByb3BoZXRfZGlzdHJpYnV0aW9uKAogIHRiX3B5cHJvcGhldF9mdW4sIGNvbHVtbj0iZGVsdGFfcnQiKSkKdGJfZGVsdGFydF9wbG90X2FsbFtbImJveHBsb3QiXV0KCnRiX2RlbHRhcnRfcGxvdF9yZWFsIDwtIHNtKHBsb3RfcHlwcm9waGV0X2Rpc3RyaWJ1dGlvbigKICB0Yl9weXByb3BoZXRfZnVuLAogIGNvbHVtbj0iZGVsdGFfcnQiLCBrZWVwX2RlY295cz1GQUxTRSkpCnRiX2RlbHRhcnRfcGxvdF9yZWFsW1siYm94cGxvdCJdXQoKdGJfZGVsdGFydF9wbG90X2RlY295cyA8LSBzbShwbG90X3B5cHJvcGhldF9kaXN0cmlidXRpb24oCiAgdGJfcHlwcm9waGV0X2Z1biwKICBjb2x1bW49ImRlbHRhX3J0Iiwga2VlcF9yZWFsPUZBTFNFKSkKdGJfZGVsdGFydF9wbG90X2RlY295c1tbImJveHBsb3QiXV0KdGJfZGVsdGFydF9wbG90X2RlY295c1tbImRlbnNpdHkiXV0KCnRiX3dpZHRodnNtYXNzIDwtIHNtKHBsb3RfcHlwcm9waGV0X2RhdGEodGJfcHlwcm9waGV0X2Z1biwgbGVnZW5kPUZBTFNFKSkKdGJfd2lkdGh2c21hc3MkcGxvdApgYGAKCiMjIE1Tc3RhdHMKCm1zc3RhdHMub3JnIHNlZW1zIHRvIHByb3ZpZGUgYSBjb21wbGV0ZSBzb2x1dGlvbiBmb3IgcGVyZm9ybWluZyByZWFzb25hYmxlIG1ldHJpY3Mgb2YgdGhpcyBkYXRhLgoKSSBhbSBjdXJyZW50bHkgcmVhZGluZzogaHR0cDovL21zc3RhdHMub3JnL3dwLWNvbnRlbnQvdXBsb2Fkcy8yMDE3LzAxL01Tc3RhdHNfdjMuNy4zX21hbnVhbC5wZGYKCkkgbWFkZSBzb21lIG1vZGVyYXRlbHkgaW50cnVzaXZlIGNoYW5nZXMgdG8gTVNzdGF0cyB0byBtYWtlIGl0IGNsZWFyZXIsIGFzIHdlbGwuCgojIyMgQ29tZXQgZGVyaXZlZCBtc3N0YXRzCgpgYGB7ciBvdXJfbXNzdGF0c19xdWFudH0KZGV2dG9vbHM6OmxvYWRfYWxsKCJ+L3NjcmF0Y2gvZ2l0L01Tc3RhdHMiKQpvdXJfbXNzdGF0c19xdWFudCA8LSBzbShkYXRhUHJvY2VzcyhvdXJfbXNzdGF0c19pbnB1dCkpCm91cl9tc3N0YXRzX3Bsb3RzIDwtIHNtKGRhdGFQcm9jZXNzUGxvdHMob3VyX21zc3RhdHNfcXVhbnQsIHR5cGU9IlFDUExPVCIpKQoKbXlfbGV2ZWxzIDwtIGxldmVscyhhcy5mYWN0b3Iob3VyX21zc3RhdHNfaW5wdXQkY29uZGl0aW9uKSkKbXlfbGV2ZWxzCmNvbXBhcmlzb25zIDwtIGdoZXR0b19jb250cmFzdF9tYXRyaXgoCiAgbnVtZXJhdG9ycz1jKCJ3dF9jZiIsICJkZWx0YV9jZiIsICJjb21wX2NmIiwKICAgICAgICAgICAgICAgImRlbHRhX2NmIiwgImNvbXBfY2YiLCAiZGVsdGFfd2hvbGUiLAogICAgICAgICAgICAgICAiY29tcF93aG9sZSIpLAogIGRlbm9taW5hdG9ycz1jKCJ3dF93aG9sZSIsICJkZWx0YV93aG9sZSIsICJjb21wX3dob2xlIiwKICAgICAgICAgICAgICAgICAid3RfY2YiLCAid3RfY2YiLCAid3Rfd2hvbGUiLAogICAgICAgICAgICAgICAgICJ3dF93aG9sZSIpKQpvdXJfcmVzdWx0cyA8LSBsaXN0KCkKZm9yIChjIGluIDE6bGVuZ3RoKHJvd25hbWVzKGNvbXBhcmlzb25zKSkpIHsKICBuYW1lIDwtIHJvd25hbWVzKGNvbXBhcmlzb25zKVtjXQogIG1lc3NhZ2UoIlN0YXJ0aW5nICIsIG5hbWUpCiAgY29tcCA8LSBjb21wYXJpc29uc1tjLCBdCiAgY29tcGFyaXNvbiA8LSB0KGFzLm1hdHJpeChjb21wKSkKICByb3duYW1lcyhjb21wYXJpc29uKSA8LSBuYW1lCiAgb3VyX3Jlc3VsdHNbbmFtZV0gPC0gc20oTVNzdGF0czo6Z3JvdXBDb21wYXJpc29uKGNvbnRyYXN0Lm1hdHJpeD1jb21wYXJpc29uLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkYXRhPW91cl9tc3N0YXRzX3F1YW50KSkKfQpgYGAKCiMjIyBUdWJlcmN1bGlzdCBkZXJpdmVkIG1zc3RhdHMKCmBgYHtyIHRiX21zc3RhdHNfcXVhbnR9CnRiX21zc3RhdHNfcXVhbnQgPC0gc20oZGF0YVByb2Nlc3ModGJfbXNzdGF0c19pbnB1dCkpCiMjdGJfbXNzdGF0c19wbG90cyA8LSBzbShkYXRhUHJvY2Vzc1Bsb3RzKHRiX21zc3RhdHNfcXVhbnQsIHR5cGU9IlFDUExPVCIpKQoKbXlfbGV2ZWxzIDwtIGxldmVscyhhcy5mYWN0b3IodGJfbXNzdGF0c19pbnB1dCRjb25kaXRpb24pKQpteV9sZXZlbHMKY29tcGFyaXNvbnMgPC0gZ2hldHRvX2NvbnRyYXN0X21hdHJpeCgKICBudW1lcmF0b3JzPWMoInd0X2NmIiwgImRlbHRhX2NmIiwgImNvbXBfY2YiLAogICAgICAgICAgICAgICAiZGVsdGFfY2YiLCAiY29tcF9jZiIsICJkZWx0YV93aG9sZSIsCiAgICAgICAgICAgICAgICJjb21wX3dob2xlIiksCiAgZGVub21pbmF0b3JzPWMoInd0X3dob2xlIiwgImRlbHRhX3dob2xlIiwgImNvbXBfd2hvbGUiLAogICAgICAgICAgICAgICAgICJ3dF9jZiIsICJ3dF9jZiIsICJ3dF93aG9sZSIsCiAgICAgICAgICAgICAgICAgInd0X3dob2xlIikpCnRiX3Jlc3VsdHMgPC0gbGlzdCgpCmZvciAoYyBpbiAxOmxlbmd0aChyb3duYW1lcyhjb21wYXJpc29ucykpKSB7CiAgbmFtZSA8LSByb3duYW1lcyhjb21wYXJpc29ucylbY10KICBtZXNzYWdlKCJTdGFydGluZyAiLCBuYW1lKQogIGNvbXAgPC0gY29tcGFyaXNvbnNbYywgXQogIGNvbXBhcmlzb24gPC0gdChhcy5tYXRyaXgoY29tcCkpCiAgcm93bmFtZXMoY29tcGFyaXNvbikgPC0gbmFtZQogIHRiX3Jlc3VsdHNbbmFtZV0gPC0gc20oTVNzdGF0czo6Z3JvdXBDb21wYXJpc29uKGNvbnRyYXN0Lm1hdHJpeD1jb21wYXJpc29uLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBkYXRhPXRiX21zc3RhdHNfcXVhbnQpKQp9CmBgYAoKIyMjIFAvUEUgcHJvdGVpbiBRQyBwbG90cyBmb3IgWWFuCgpZYW4gYXNrZWQgZm9yIHRoZSBwL3BlIHByb3RlaW4gcWMgcGxvdHMuIG9rLiAgSSBjaGFuZ2VkIHRoZSBkYXRhUHJvY2Vzc1Bsb3RzIHRvCnJldHVybiBzb21ldGhpbmcgdXNlZnVsLCBzbyB0aGF0IHNob3VsZCBiZSBwb3NzaWJsZSBub3cuCgpgYGB7ciBwZX0KcGVfZ2VuZXMgPC0gcmVhZC50YWJsZSgicmVmZXJlbmNlL2Fubm90YXRlZF9wZV9nZW5lcy50eHQiKVtbMV1dCgojIyBVbmZvcnR1bmF0ZWx5LCB0aGUgbmFtZXMgZGlkIG5vdCBnZXQgc2V0IGluIG15IGNoYW5nZWQgdmVyc2lvbiBvZiBkYXRhUHJvY2Vzc1Bsb3RzLi4uCnBsb3Rsc3QgPC0gb3VyX21zc3RhdHNfcGxvdHMkUUNQTE9UCmF2YWlsYWJsZV9wbG90cyA8LSBnc3ViKHBhdHRlcm49Il4xLyIsIHJlcGxhY2VtZW50PSIiLAogICAgICAgICAgICAgICAgICAgICAgICB4PWxldmVscyhvdXJfbXNzdGF0c19xdWFudCRQcm9jZXNzZWREYXRhJFBST1RFSU4pKQpuYW1lcyhwbG90bHN0KSA8LSBhdmFpbGFibGVfcGxvdHMKCnBlX2luX2F2YWlsX2lkeCA8LSBwZV9nZW5lcyAlaW4lIGF2YWlsYWJsZV9wbG90cwpwZV9pbl9hdmFpbCA8LSBwZV9nZW5lc1twZV9pbl9hdmFpbF9pZHhdCnBlX3Bsb3RzIDwtIHBsb3Rsc3RbcGVfaW5fYXZhaWxdCnBkZihmaWxlPSJwZV9xY19wbG90cy5wZGYiKQpmb3IgKHAgaW4gMTpsZW5ndGgocGVfcGxvdHMpKSB7CiAgcGxvdChwZV9wbG90c1tbcF1dKQp9CmRldi5vZmYoKQpsZW5ndGgocGVfcGxvdHMpCmBgYAoKIyBDcmVhdGUgaHBnbHRvb2xzIGV4cHJlc3Npb25zZXQKClNpbmNlIEkgYW0gbm90IGNlcnRhaW4gSSB1bmRlcnN0YW5kIHRoZXNlIGRhdGEsIEkgd2lsbCB0YWtlIHRoZSBpbnRlbnNpdGllcyBmcm9tClNXQVRIMnN0YXRzLCBtZXRhZGF0YSwgYW5kIGFubm90YXRpb24gZGF0YTsgIGF0dGVtcHQgdG8gY3JlYXRlIGEgJ25vcm1hbCcKZXhwcmVzc2lvbnNldDsgcG9rZSBhdCBpdCB0byBzZWUgd2hhdCBJIGNhbiBsZWFybi4KCiMjIE1hc3NhZ2luZyB0aGUgbWV0YWRhdGEKCkkgd2FudCB0byB1c2UgdGhlIHNhbWUgbWV0YWRhdGEgYXMgd2VyZSB1c2VkIGZvciBNU3N0YXRzLiAgSXQgaGFzIGEgZmV3CmltcG9ydGFudCBkaWZmZXJlbmNlcyBmcm9tIHRoZSByZXF1aXJlbWVudHMgb2YgaHBnbHRvb2xzOiBwcmV0dHkgbXVjaCBvbmx5IHRoYXQKSSBkbyBub3QgYWxsb3cgcm93bmFtZXMvc2FtcGxlSURzIHRvIHN0YXJ0IHdpdGggYSBudW1iZXIuCgojIyBNYXNzYWdpbmcgdGhlIGludGVuc2l0eSBtYXRyaXgKCkkgZG8gbm90IHdhbnQgdGhlIFwxIGJlZm9yZSB0aGUgcHJvdGVpbiBuYW1lcywgSSBhbHJlYWR5IG1lcmdlZCB0aGVtIGludG8gb25lCmVudHJ5IHBlciBnZW5lIHZpcyBTV0FUSDJzdGF0cy4KCiMjIyBNYWtlIG1hdHJpY2VzIHdpdGggb3VyIGxpYnJhcmllcwoKYGBge3Igb3VyX3Byb3RlaW5fbWF0cml4fQpvdXJfcHJvdF9tdHJ4IDwtIHJlYWQuY3N2KGZpbGUucGF0aCgicmVzdWx0cyIsICJzd2F0aDJzdGF0cyIsIHZlciwgIm91cl9wcm90ZWluX21hdHJpeF9taW5pbXVtLmNzdiIpKQpyb3duYW1lcyhvdXJfcHJvdF9tdHJ4KSA8LSBnc3ViKHBhdHRlcm49Il4xXFwvIiwgcmVwbGFjZW1lbnQ9IiIsIHg9b3VyX3Byb3RfbXRyeFtbInByb3RlaW5uYW1lIl1dKQpvdXJfcHJvdF9tdHJ4IDwtIG91cl9wcm90X210cnhbLCAtMV0KIyMgSW1wb3J0YW50IHF1ZXN0aW9uOiBEaWQgU1dBVEgyc3RhdHMgcmVvcmRlciBteSBkYXRhPwpjb2xuYW1lcyhvdXJfcHJvdF9tdHJ4KSA8LSBnc3ViKHBhdHRlcm49Il5YIiwgcmVwbGFjZW1lbnQ9InMiLCB4PWNvbG5hbWVzKG91cl9wcm90X210cngpKQpgYGAKCiMjIyBNYWtlIG1hdHJpY2VzIHdpdGggdGhlIHRiIGxpYnJhcmllcwoKYGBge3IgdGJfcHJvdGVpbl9tYXRyaXh9CnRiX3Byb3RfbXRyeCA8LSByZWFkLmNzdihmaWxlLnBhdGgoInJlc3VsdHMiLCAic3dhdGgyc3RhdHMiLCB2ZXIsICJ0Yl9wcm90ZWluX21hdHJpeF9taW5pbXVtLmNzdiIpKQpyb3duYW1lcyh0Yl9wcm90X210cngpIDwtIGdzdWIocGF0dGVybj0iXjFcXC8iLCByZXBsYWNlbWVudD0iIiwgeD10Yl9wcm90X210cnhbWyJwcm90ZWlubmFtZSJdXSkKdGJfcHJvdF9tdHJ4IDwtIHRiX3Byb3RfbXRyeFssIC0xXQojIyBJbXBvcnRhbnQgcXVlc3Rpb246IERpZCBTV0FUSDJzdGF0cyByZW9yZGVyIG15IGRhdGE/CmNvbG5hbWVzKHRiX3Byb3RfbXRyeCkgPC0gZ3N1YihwYXR0ZXJuPSJeWCIsIHJlcGxhY2VtZW50PSJzIiwgeD1jb2xuYW1lcyh0Yl9wcm90X210cngpKQpgYGAKCiMjIE1lcmdlIHRoZSBwaWVjZXMKCk5vdyB3ZSBzaG91bGQgaGF2ZSBzdWZmaWNpZW50IHBpZWNlcyB0byBtYWtlIGFuIGV4cHJlc3Npb25zZXQuCgpXaGlsZSBoZXJlLCBJIHdpbGwgYWxzbyBzcGxpdCB0aGUgZGF0YSBpbnRvIGEgY2YgYW5kIHdob2xlLWNlbGwgcGFpciBvZiBkYXRhIHN0cnVjdHVyZXMuCgojIyMgT3VyIGRhdGEKCmBgYHtyIG91cl9leHB0fQojIyBEcm9wIHRoZSBtZXRhZGF0YSBub3QgaW4gdGhlIHByb3RlaW4gbWF0cml4OgojIyBBbmQgZW5zdXJlIHRoYXQgdGhleSBhcmUgdGhlIHNhbWUgb3JkZXIuCnJlb3JkZXJlZCA8LSBjb2xuYW1lcyhvdXJfcHJvdF9tdHJ4KQptZXRhZGF0YSA8LSBzYW1wbGVfYW5ub3RbcmVvcmRlcmVkLCBdCgpvdXJfcHJvdGVpbl9leHB0IDwtIHNtKGNyZWF0ZV9leHB0KG1ldGFkYXRhPW1ldGFkYXRhLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvdW50X2RhdGFmcmFtZT1vdXJfcHJvdF9tdHJ4LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGdlbmVfaW5mbz1tdGJfYW5ub3RhdGlvbnMpKQoKb3VyX3dob2xlX2V4cHQgPC0gc3Vic2V0X2V4cHQob3VyX3Byb3RlaW5fZXhwdCwgc3Vic2V0PSJjb2xsZWN0aW9udHlwZT09J3dob2xlJyIpCm91cl9jZl9leHB0IDwtIHN1YnNldF9leHB0KG91cl9wcm90ZWluX2V4cHQsIHN1YnNldD0iY29sbGVjdGlvbnR5cGU9PSdjZiciKQpgYGAKCiMjIyBUYiBkYXRhCgpgYGB7ciB0Yl9leHB0fQojIyBEcm9wIHRoZSBtZXRhZGF0YSBub3QgaW4gdGhlIHByb3RlaW4gbWF0cml4OgojIyBBbmQgZW5zdXJlIHRoYXQgdGhleSBhcmUgdGhlIHNhbWUgb3JkZXIuCnJlb3JkZXJlZCA8LSBjb2xuYW1lcyh0Yl9wcm90X210cngpCm1ldGFkYXRhIDwtIHNhbXBsZV9hbm5vdFtyZW9yZGVyZWQsIF0KCnRiX3Byb3RlaW5fZXhwdCA8LSBzbShjcmVhdGVfZXhwdChtZXRhZGF0YSwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGNvdW50X2RhdGFmcmFtZT10Yl9wcm90X210cngsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBnZW5lX2luZm89bXRiX2Fubm90YXRpb25zKSkKdGJfd2hvbGVfZXhwdCA8LSBzdWJzZXRfZXhwdCh0Yl9wcm90ZWluX2V4cHQsIHN1YnNldD0iY29sbGVjdGlvbnR5cGU9PSd3aG9sZSciKQp0Yl9jZl9leHB0IDwtIHN1YnNldF9leHB0KHRiX3Byb3RlaW5fZXhwdCwgc3Vic2V0PSJjb2xsZWN0aW9udHlwZT09J2NmJyIpCmBgYAoKIyMgTWV0cmljcyBvZiB0aGUgZnVsbCBkYXRhIHNldAoKIyMjIE91ciBsaWJyYXJpZXMKCmBgYHtyIG91cl9wcm90ZWluX21ldHJpY3MsIGZpZy5zaG93PSdoaWRlJ30Kb3VyX3Byb3RlaW5fbWV0cmljcyA8LSBzbShncmFwaF9tZXRyaWNzKG91cl9wcm90ZWluX2V4cHQpKQpvdXJfcHJvdGVpbl9ub3JtIDwtIHNtKG5vcm1hbGl6ZV9leHB0KG91cl9wcm90ZWluX2V4cHQsIHRyYW5zZm9ybT0ibG9nMiIsIGNvbnZlcnQ9ImNwbSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbm9ybT0icXVhbnQiLCBmaWx0ZXI9VFJVRSkpCm91cl9wcm90ZWluX25vcm1fbWV0cmljcyA8LSBzbShncmFwaF9tZXRyaWNzKG91cl9wcm90ZWluX25vcm0pKQpvdXJfcHJvdGVpbl9mc3ZhIDwtIHNtKG5vcm1hbGl6ZV9leHB0KG91cl9wcm90ZWluX2V4cHQsIHRyYW5zZm9ybT0ibG9nMiIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYmF0Y2g9ImZzdmEiLCBmaWx0ZXI9VFJVRSkpCm91cl9wcm90ZWluX2ZzdmFfbWV0cmljcyA8LSBzbShncmFwaF9tZXRyaWNzKG91cl9wcm90ZWluX2ZzdmEpKQpgYGAKCiMjIyBUYiBsaWJyYXJpZXMKCmBgYHtyIHRiX3Byb3RlaW5fbWV0cmljcywgZmlnLnNob3c9J2hpZGUnfQp0Yl9wcm90ZWluX21ldHJpY3MgPC0gc20oZ3JhcGhfbWV0cmljcyh0Yl9wcm90ZWluX2V4cHQpKQp0Yl9wcm90ZWluX25vcm0gPC0gc20obm9ybWFsaXplX2V4cHQodGJfcHJvdGVpbl9leHB0LCB0cmFuc2Zvcm09ImxvZzIiLCBjb252ZXJ0PSJjcG0iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG5vcm09InF1YW50IiwgZmlsdGVyPVRSVUUpKQp0Yl9wcm90ZWluX25vcm1fbWV0cmljcyA8LSBzbShncmFwaF9tZXRyaWNzKHRiX3Byb3RlaW5fbm9ybSkpCnRiX3Byb3RlaW5fZnN2YSA8LSBzbShub3JtYWxpemVfZXhwdCh0Yl9wcm90ZWluX2V4cHQsIHRyYW5zZm9ybT0ibG9nMiIsIGNvbnZlcnQ9ImNwbSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYmF0Y2g9ImZzdmEiLCBmaWx0ZXI9VFJVRSkpCnRiX3Byb3RlaW5fZnN2YV9tZXRyaWNzIDwtIHNtKGdyYXBoX21ldHJpY3ModGJfcHJvdGVpbl9mc3ZhKSkKYGBgCgojIyBNZXRyaWNzIG9mIHRoZSB3aG9sZS1jZWxsIGRhdGEgc2V0CgojIyMgT3VyIGxpYnJhcmllcwoKYGBge3Igb3VyX3dob2xlX21ldHJpY3MsIGZpZy5zaG93PSdoaWRlJ30Kb3VyX3dob2xlX21ldHJpY3MgPC0gc20oZ3JhcGhfbWV0cmljcyhvdXJfd2hvbGVfZXhwdCkpCm91cl93aG9sZV9ub3JtIDwtIHNtKG5vcm1hbGl6ZV9leHB0KG91cl93aG9sZV9leHB0LCB0cmFuc2Zvcm09ImxvZzIiLCBjb252ZXJ0PSJjcG0iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBub3JtPSJxdWFudCIsIGZpbHRlcj1UUlVFKSkKb3VyX3dob2xlX25vcm1fbWV0cmljcyA8LSBzbShncmFwaF9tZXRyaWNzKG91cl93aG9sZV9ub3JtKSkKb3VyX3dob2xlX2ZzdmEgPC0gc20obm9ybWFsaXplX2V4cHQob3VyX3dob2xlX2V4cHQsIHRyYW5zZm9ybT0ibG9nMiIsIGNvbnZlcnQ9ImNwbSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGJhdGNoPSJmc3ZhIiwgZmlsdGVyPVRSVUUpKQpvdXJfd2hvbGVfZnN2YV9tZXRyaWNzIDwtIHNtKGdyYXBoX21ldHJpY3Mob3VyX3dob2xlX2ZzdmEpKQpgYGAKCiMjIyBUYiBsaWJyYXJpZXMKCmBgYHtyIHRiX3dob2xlX21ldHJpY3MsIGZpZy5zaG93PSdoaWRlJ30KdGJfd2hvbGVfbWV0cmljcyA8LSBzbShncmFwaF9tZXRyaWNzKHRiX3dob2xlX2V4cHQpKQp0Yl93aG9sZV9ub3JtIDwtIHNtKG5vcm1hbGl6ZV9leHB0KHRiX3dob2xlX2V4cHQsIHRyYW5zZm9ybT0ibG9nMiIsIGNvbnZlcnQ9ImNwbSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG5vcm09InF1YW50IiwgZmlsdGVyPVRSVUUpKQp0Yl93aG9sZV9ub3JtX21ldHJpY3MgPC0gc20oZ3JhcGhfbWV0cmljcyh0Yl93aG9sZV9ub3JtKSkKdGJfd2hvbGVfZnN2YSA8LSBzbShub3JtYWxpemVfZXhwdCh0Yl93aG9sZV9leHB0LCB0cmFuc2Zvcm09ImxvZzIiLCBjb252ZXJ0PSJjcG0iLAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBiYXRjaD0iZnN2YSIsIGZpbHRlcj1UUlVFKSkKdGJfd2hvbGVfZnN2YV9tZXRyaWNzIDwtIHNtKGdyYXBoX21ldHJpY3ModGJfd2hvbGVfZnN2YSkpCmBgYAoKIyMgTWV0cmljcyBvZiB0aGUgZmlsdHJhdGUgZGF0YSBzZXQKCiMjIyBPdXIgbGlicmFyaWVzCgpgYGB7ciBvdXJfY2ZfbWV0cmljcywgZmlnLnNob3c9J2hpZGUnfQpvdXJfY2ZfbWV0cmljcyA8LSBzbShncmFwaF9tZXRyaWNzKG91cl9jZl9leHB0KSkKb3VyX2NmX25vcm0gPC0gc20obm9ybWFsaXplX2V4cHQob3VyX2NmX2V4cHQsIHRyYW5zZm9ybT0ibG9nMiIsIGNvbnZlcnQ9ImNwbSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG5vcm09InF1YW50IiwgZmlsdGVyPVRSVUUpKQpvdXJfY2Zfbm9ybV9tZXRyaWNzIDwtIHNtKGdyYXBoX21ldHJpY3Mob3VyX2NmX25vcm0pKQpvdXJfY2ZfZnN2YSA8LSBzbShub3JtYWxpemVfZXhwdChvdXJfY2ZfZXhwdCwgdHJhbnNmb3JtPSJsb2cyIiwgY29udmVydD0iY3BtIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgYmF0Y2g9ImZzdmEiLCBmaWx0ZXI9VFJVRSkpCm91cl9jZl9mc3ZhX21ldHJpY3MgPC0gc20oZ3JhcGhfbWV0cmljcyhvdXJfY2ZfZnN2YSkpCmBgYAoKIyMjIFRiIGxpYnJhcmllcwoKYGBge3IgdGJfY2ZfbWV0cmljcywgZmlnLnNob3c9J2hpZGUnfQp0Yl9jZl9tZXRyaWNzIDwtIHNtKGdyYXBoX21ldHJpY3ModGJfY2ZfZXhwdCkpCnRiX2NmX25vcm0gPC0gc20obm9ybWFsaXplX2V4cHQodGJfY2ZfZXhwdCwgdHJhbnNmb3JtPSJsb2cyIiwgY29udmVydD0iY3BtIiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbm9ybT0icXVhbnQiLCBmaWx0ZXI9VFJVRSkpCnRiX2NmX25vcm1fbWV0cmljcyA8LSBzbShncmFwaF9tZXRyaWNzKHRiX2NmX25vcm0pKQp0Yl9jZl9mc3ZhIDwtIHNtKG5vcm1hbGl6ZV9leHB0KHRiX2NmX2V4cHQsIHRyYW5zZm9ybT0ibG9nMiIsIGNvbnZlcnQ9ImNwbSIsCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIGJhdGNoPSJmc3ZhIiwgZmlsdGVyPVRSVUUpKQp0Yl9jZl9mc3ZhX21ldHJpY3MgPC0gc20oZ3JhcGhfbWV0cmljcyh0Yl9jZl9mc3ZhKSkKYGBgCgojIyBwbG90IHNvbWUgbWV0cmljcwoKIyMjIE91ciBsaWJyYXJpZXMKCmBgYHtyIG91cl9wcmludF9tZXRyaWNzfQpwcChpbWFnZT1vdXJfcHJvdGVpbl9tZXRyaWNzJGxpYnNpemUsCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX291cl9saWJzaXplLnBuZyIpKSkKIyMgSXQgc2VlbXMgdG8gbWUgdGhhdCB0aGUgc2NhbGUgb2YgdGhlIGRhdGEgaXMgYWxsIHdpdGhpbiBhbiBvcmRlciBvZiBtYWduaXR1ZGUgb3IgdHdvLgojIyBJIGNhbm5vdCBnZXQgdXNlZCB0byB0aGVzZSBhYnN1cmRseSBsYXJnZSBudW1iZXJzIHRob3VnaC4KcHAoaW1hZ2U9b3VyX3Byb3RlaW5fbm9ybV9tZXRyaWNzJHBjYXBsb3QsCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX291cl9ub3JtX3BjYS5wbmciKSkpCiMjIFRoZXJlIGFwcGVhcnMgdG8gYmUgYSBuaWNlIHNwbGl0IGluIHRoZSBkYXRhLCBob3dldmVyIHRoZSB1bi1hc3NheWFibGUgYmF0Y2gKIyMgZWZmZWN0IGlzIGEgcHJvYmxlbS4KcHAoaW1hZ2U9b3VyX3Byb3RlaW5fZnN2YV9tZXRyaWNzJHBjYXBsb3QsCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX2ZzdmFfcGNhLnBuZyIpKSkKIyMgZnN2YSBzZWVtcyB0byBnZXQgc29tZSBoYW5kbGUgb24gdGhlIGRhdGEsIGJ1dCBJIGRvbid0IHRoaW5rIHdlIHNob3VsZCByZWx5CiMjIHVwb24gaXQuCnBwKGltYWdlPW91cl9wcm90ZWluX25vcm1fbWV0cmljcyRjb3JoZWF0LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl9ub3JtX2NvcmhlYXQucG5nIikpKQojIyBPbmNlIGFnYWluLCB0aGUgd2hvbGUtY2VsbC9jdWx0dXJlLWZpbHRyYXRlIHNwbGl0IGlzIHZlcnkgbGFyZ2UuCnBwKGltYWdlPW91cl9wcm90ZWluX21ldHJpY3MkZGVuc2l0eSwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfcmF3X2RlbnNpdHkucG5nIikpKQojIyBUaGVyZSBhcmUgdHdvIG9idmlvdXMgZGlzdHJpYnV0aW9ucyBpbiB0aGUgZGF0YSwgb25jZSBhZ2FpbiBzcGxpdCBiZXR3ZWVuIHR5cGVzLgpwcChpbWFnZT1vdXJfcHJvdGVpbl9tZXRyaWNzJGJveHBsb3QsCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX2JveHBsb3QucG5nIikpKQojIyBUaGlzIHJlY2FwaXR1bGF0ZXMgdGhlIHByZXZpb3VzIHBsb3QuCgpwcChpbWFnZT1vdXJfd2hvbGVfbWV0cmljcyRsaWJzaXplLAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl93aG9sZV9saWJzaXplLnBuZyIpKSkKcHAoaW1hZ2U9b3VyX3dob2xlX25vcm1fbWV0cmljcyRwY2FwbG90LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl93aG9sZV9ub3JtX3BjYS5wbmciKSkpCnBwKGltYWdlPW91cl93aG9sZV9mc3ZhX21ldHJpY3MkcGNhcGxvdCwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfd2hvbGVfZnN2YV9wY2EucG5nIikpKQpwcChpbWFnZT1vdXJfd2hvbGVfbm9ybV9tZXRyaWNzJGNvcmhlYXQsCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX3dob2xlX25vcm1fY29yaGVhdC5wbmciKSkpCnBwKGltYWdlPW91cl93aG9sZV9tZXRyaWNzJGRlbnNpdHksCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX3dob2xlX3Jhd19kZW5zaXR5LnBuZyIpKSkKcHAoaW1hZ2U9b3VyX3dob2xlX21ldHJpY3MkYm94cGxvdCwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfd2hvbGVfYm94cGxvdC5wbmciKSkpCgpwcChpbWFnZT1vdXJfY2ZfbWV0cmljcyRsaWJzaXplLAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl9saWJzaXplLnBuZyIpKSkKcHAoaW1hZ2U9b3VyX2NmX25vcm1fbWV0cmljcyRwY2FwbG90LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl9ub3JtX3BjYS5wbmciKSkpCnBwKGltYWdlPW91cl9jZl9mc3ZhX21ldHJpY3MkcGNhcGxvdCwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfZnN2YV9wY2EucG5nIikpKQpwcChpbWFnZT1vdXJfY2Zfbm9ybV9tZXRyaWNzJGNvcmhlYXQsCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX25vcm1fY29yaGVhdC5wbmciKSkpCnBwKGltYWdlPW91cl9jZl9tZXRyaWNzJGRlbnNpdHksCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX3Jhd19kZW5zaXR5LnBuZyIpKSkKcHAoaW1hZ2U9b3VyX2NmX21ldHJpY3MkYm94cGxvdCwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfYm94cGxvdC5wbmciKSkpCmBgYAoKIyMjIFRiIGxpYnJhcmllcwoKYGBge3IgdGJfcHJpbnRfbWV0cmljc30KcHAoaW1hZ2U9dGJfcHJvdGVpbl9tZXRyaWNzJGxpYnNpemUsCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX3RiX2xpYnNpemUucG5nIikpKQojIyBJdCBzZWVtcyB0byBtZSB0aGF0IHRoZSBzY2FsZSBvZiB0aGUgZGF0YSBpcyBhbGwgd2l0aGluIGFuIG9yZGVyIG9mIG1hZ25pdHVkZSBvciB0d28uCiMjIEkgY2Fubm90IGdldCB1c2VkIHRvIHRoZXNlIGFic3VyZGx5IGxhcmdlIG51bWJlcnMgdGhvdWdoLgpwcChpbWFnZT10Yl9wcm90ZWluX25vcm1fbWV0cmljcyRwY2FwbG90LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl90Yl9ub3JtX3BjYS5wbmciKSkpCiMjIFRoZXJlIGFwcGVhcnMgdG8gYmUgYSBuaWNlIHNwbGl0IGluIHRoZSBkYXRhLCBob3dldmVyIHRoZSB1bi1hc3NheWFibGUgYmF0Y2gKIyMgZWZmZWN0IGlzIGEgcHJvYmxlbS4KcHAoaW1hZ2U9dGJfcHJvdGVpbl9mc3ZhX21ldHJpY3MkcGNhcGxvdCwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfZnN2YV9wY2EucG5nIikpKQojIyBmc3ZhIHNlZW1zIHRvIGdldCBzb21lIGhhbmRsZSBvbiB0aGUgZGF0YSwgYnV0IEkgZG9uJ3QgdGhpbmsgd2Ugc2hvdWxkIHJlbHkKIyMgdXBvbiBpdC4KcHAoaW1hZ2U9dGJfcHJvdGVpbl9ub3JtX21ldHJpY3MkY29yaGVhdCwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfbm9ybV9jb3JoZWF0LnBuZyIpKSkKIyMgT25jZSBhZ2FpbiwgdGhlIHdob2xlLWNlbGwvY3VsdHVyZS1maWx0cmF0ZSBzcGxpdCBpcyB2ZXJ5IGxhcmdlLgpwcChpbWFnZT10Yl9wcm90ZWluX21ldHJpY3MkZGVuc2l0eSwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfcmF3X2RlbnNpdHkucG5nIikpKQojIyBUaGVyZSBhcmUgdHdvIG9idmlvdXMgZGlzdHJpYnV0aW9ucyBpbiB0aGUgZGF0YSwgb25jZSBhZ2FpbiBzcGxpdCBiZXR3ZWVuIHR5cGVzLgpwcChpbWFnZT10Yl9wcm90ZWluX21ldHJpY3MkYm94cGxvdCwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfYm94cGxvdC5wbmciKSkpCiMjIFRoaXMgcmVjYXBpdHVsYXRlcyB0aGUgcHJldmlvdXMgcGxvdC4KCnBwKGltYWdlPXRiX3dob2xlX21ldHJpY3MkbGlic2l6ZSwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfd2hvbGVfbGlic2l6ZS5wbmciKSkpCnBwKGltYWdlPXRiX3dob2xlX25vcm1fbWV0cmljcyRwY2FwbG90LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl93aG9sZV9ub3JtX3BjYS5wbmciKSkpCnBwKGltYWdlPXRiX3dob2xlX2ZzdmFfbWV0cmljcyRwY2FwbG90LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl93aG9sZV9mc3ZhX3BjYS5wbmciKSkpCnBwKGltYWdlPXRiX3dob2xlX25vcm1fbWV0cmljcyRjb3JoZWF0LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl93aG9sZV9ub3JtX2NvcmhlYXQucG5nIikpKQpwcChpbWFnZT10Yl93aG9sZV9tZXRyaWNzJGRlbnNpdHksCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX3dob2xlX3Jhd19kZW5zaXR5LnBuZyIpKSkKcHAoaW1hZ2U9dGJfd2hvbGVfbWV0cmljcyRib3hwbG90LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl93aG9sZV9ib3hwbG90LnBuZyIpKSkKCnBwKGltYWdlPXRiX2NmX21ldHJpY3MkbGlic2l6ZSwKICAgZmlsZT1maWxlLnBhdGgoImltYWdlcyIsIHBhc3RlMCh2ZXIsICJfbGlic2l6ZS5wbmciKSkpCnBwKGltYWdlPXRiX2NmX25vcm1fbWV0cmljcyRwY2FwbG90LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl9ub3JtX3BjYS5wbmciKSkpCnBwKGltYWdlPXRiX2NmX2ZzdmFfbWV0cmljcyRwY2FwbG90LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl9mc3ZhX3BjYS5wbmciKSkpCnBwKGltYWdlPXRiX2NmX25vcm1fbWV0cmljcyRjb3JoZWF0LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl9ub3JtX2NvcmhlYXQucG5nIikpKQpwcChpbWFnZT10Yl9jZl9tZXRyaWNzJGRlbnNpdHksCiAgIGZpbGU9ZmlsZS5wYXRoKCJpbWFnZXMiLCBwYXN0ZTAodmVyLCAiX3Jhd19kZW5zaXR5LnBuZyIpKSkKcHAoaW1hZ2U9dGJfY2ZfbWV0cmljcyRib3hwbG90LAogICBmaWxlPWZpbGUucGF0aCgiaW1hZ2VzIiwgcGFzdGUwKHZlciwgIl9ib3hwbG90LnBuZyIpKSkKYGBgCgojIEF0dGVtcHQgc29tZSBxdWFudGlmaWNhdGlvbiBjb21wYXJpc29ucz8KCiMjIE91ciBsaWJyYXJpZXMKCmBgYHtyIG91cl9kZV90ZXN0aW5nLCBmaWcuc2hvdz0naGlkZSd9Cm91cl9wYWlyd2lzZV9maWx0IDwtIHNtKG5vcm1hbGl6ZV9leHB0KG91cl9wcm90ZWluX2V4cHQsIGZpbHRlcj1UUlVFKSkKb3VyX3BhaXJ3aXNlX2NvbXAgPC0gc20oYWxsX3BhaXJ3aXNlKG91cl9wYWlyd2lzZV9maWx0LCBtb2RlbF9iYXRjaD0iZnN2YSIsIGZvcmNlPVRSVUUpKQpvdXJfcGFpcndpc2Vfbm9iYXRjaCA8LSBzbShhbGxfcGFpcndpc2Uob3VyX3BhaXJ3aXNlX2ZpbHQsIG1vZGVsX2JhdGNoPUZBTFNFLCBmb3JjZT1UUlVFKSkKYGBgCgojIyBUYiBsaWJyYXJpZXMKCmBgYHtyIHRiX2RlX3Rlc3RpbmcsIGZpZy5zaG93PSdoaWRlJ30KdGJfcGFpcndpc2VfZmlsdCA8LSBzbShub3JtYWxpemVfZXhwdCh0Yl9wcm90ZWluX2V4cHQsIGZpbHRlcj1UUlVFKSkKdGJfcGFpcndpc2VfY29tcCA8LSBzbShhbGxfcGFpcndpc2UodGJfcGFpcndpc2VfZmlsdCwgbW9kZWxfYmF0Y2g9ImZzdmEiLCBmb3JjZT1UUlVFKSkKdGJfcGFpcndpc2Vfbm9iYXRjaCA8LSBzbShhbGxfcGFpcndpc2UodGJfcGFpcndpc2VfZmlsdCwgbW9kZWxfYmF0Y2g9RkFMU0UsIGZvcmNlPVRSVUUpKQpgYGAKCiMjIyBTaG93IGEgZmV3IG1ldHJpY3MgZnJvbSB0aGUgaHBnbHRvb2xzIHBhaXJ3aXNlIGNvbXBhcmlzb25zCgpgYGB7ciBzaG93X2RlX3Rlc3Rpbmd9Cm91cl9wYWlyd2lzZV9jb21wJGNvbXBhcmlzb24kaGVhdAp0Yl9wYWlyd2lzZV9jb21wJGNvbXBhcmlzb24kaGVhdApgYGAKCiMgRm9yIGVhY2ggbXNzdGF0cyBydW4sIGRvIGEgREUgdGFibGUKCiMjIHd0X2NmIHZzIHd0X3dob2xlCgojIyMgT3VyIGxpYnJhcmllcwoKYGBge3Igb3VyX2NvbWJpbmVfd3RjZl93dHdob2xlLCBmaWcuc2hvdz0naGlkZSd9CmtlZXBlcnMgPC0gbGlzdCgKICAid3RjZl92c193dHdob2xlIiA9IGMoInd0X2NmIiwgInd0X3dob2xlIikpCmRyb3BwZXJzIDwtIGMoInVuZGVmaW5lZCIpCm5hbWVzKGRyb3BwZXJzKSA8LSAibG9nMmZjIgoKIyMgTWFrZSBzdXJlIHRvIHNldCB0aGUgcm93bmFtZXMgc28gaXQgd2lsbCBtZXJnZSBpbnRvIHRoZSBleGNlbCBmaWxlLgpzZXRfbmFtZSA8LSAid3RfY2ZfdnNfd3Rfd2hvbGUiCnJvd25hbWVzKG91cl9yZXN1bHRzW1tzZXRfbmFtZV1dKSA8LSBvdXJfcmVzdWx0c1tbc2V0X25hbWVdXVtbIlByb3RlaW4iXV0Kd3RjZl93dHdob2xlX3RhYmxlcyA8LSBzbShjb21iaW5lX2RlX3RhYmxlcygKICBvdXJfcGFpcndpc2VfY29tcCwga2VlcGVycz1rZWVwZXJzLCBleHRyYV9hbm5vdD1vdXJfcmVzdWx0c1tbc2V0X25hbWVdXSwKICBleGNsdWRlcz1kcm9wcGVycywKICBleGNlbD1wYXN0ZTAoImV4Y2VsL291cl93dGNmX3ZzX3d0d2hvbGVfdGFibGVzLXYiLCB2ZXIsICIueGxzeCIpKSkKCnd0Y2Zfbm9iYXRjaF93dHdob2xlX3RhYmxlcyA8LSBzbShjb21iaW5lX2RlX3RhYmxlcygKICBvdXJfcGFpcndpc2Vfbm9iYXRjaCwga2VlcGVycz1rZWVwZXJzLCBleHRyYV9hbm5vdD1vdXJfcmVzdWx0c1tbc2V0X25hbWVdXSwKICBleGNsdWRlcz1kcm9wcGVycywKICBleGNlbD1wYXN0ZTAoImV4Y2VsL291cl93dGNmX3ZzX3d0d2hvbGVfbm9iYXRjaF90YWJsZXMtdiIsIHZlciwgIi54bHN4IikpKQoKY29tcF90YWJsZSA8LSB3dGNmX3d0d2hvbGVfdGFibGVzJGRhdGFbW3NldF9uYW1lXV0KY29yLnRlc3QoY29tcF90YWJsZSRsb2cyZmMsIGNvbXBfdGFibGUkZGVzZXFfbG9nZmMsIG1ldGhvZD0ic3BlYXJtYW4iKQoKY29tcF90YWJsZSA8LSB3dGNmX25vYmF0Y2hfd3R3aG9sZV90YWJsZXMkZGF0YVtbc2V0X25hbWVdXQpjb3IudGVzdChjb21wX3RhYmxlJGxvZzJmYywgY29tcF90YWJsZSRkZXNlcV9sb2dmYywgbWV0aG9kPSJzcGVhcm1hbiIpCmBgYAoKIyMjIFRiIGxpYnJhcmllcwoKYGBge3IgdGJfY29tYmluZV93dGNmX3d0d2hvbGUsIGZpZy5zaG93PSdoaWRlJ30Ka2VlcGVycyA8LSBsaXN0KAogICJ3dGNmX3ZzX3d0d2hvbGUiID0gYygid3RfY2YiLCAid3Rfd2hvbGUiKSkKZHJvcHBlcnMgPC0gYygidW5kZWZpbmVkIikKbmFtZXMoZHJvcHBlcnMpIDwtICJsb2cyZmMiCgojIyBNYWtlIHN1cmUgdG8gc2V0IHRoZSByb3duYW1lcyBzbyBpdCB3aWxsIG1lcmdlIGludG8gdGhlIGV4Y2VsIGZpbGUuCnNldF9uYW1lIDwtICJ3dF9jZl92c193dF93aG9sZSIKcm93bmFtZXModGJfcmVzdWx0c1tbc2V0X25hbWVdXSkgPC0gdGJfcmVzdWx0c1tbc2V0X25hbWVdXVtbIlByb3RlaW4iXV0Kd3RjZl93dHdob2xlX3RhYmxlcyA8LSBzbShjb21iaW5lX2RlX3RhYmxlcygKICB0Yl9wYWlyd2lzZV9jb21wLCBrZWVwZXJzPWtlZXBlcnMsIGV4dHJhX2Fubm90PXRiX3Jlc3VsdHNbW3NldF9uYW1lXV0sCiAgZXhjbHVkZXM9ZHJvcHBlcnMsCiAgZXhjZWw9cGFzdGUwKCJleGNlbC90Yl93dGNmX3ZzX3d0d2hvbGVfdGFibGVzLXYiLCB2ZXIsICIueGxzeCIpKSkKCnd0Y2Zfbm9iYXRjaF93dHdob2xlX3RhYmxlcyA8LSBzbShjb21iaW5lX2RlX3RhYmxlcygKICB0Yl9wYWlyd2lzZV9ub2JhdGNoLCBrZWVwZXJzPWtlZXBlcnMsIGV4dHJhX2Fubm90PXRiX3Jlc3VsdHNbW3NldF9uYW1lXV0sCiAgZXhjbHVkZXM9ZHJvcHBlcnMsCiAgZXhjZWw9cGFzdGUwKCJleGNlbC90Yl93dGNmX3ZzX3d0d2hvbGVfbm9iYXRjaF90YWJsZXMtdiIsIHZlciwgIi54bHN4IikpKQoKY29tcF90YWJsZSA8LSB3dGNmX3d0d2hvbGVfdGFibGVzJGRhdGFbW3NldF9uYW1lXV0KY29yLnRlc3QoY29tcF90YWJsZSRsb2cyZmMsIGNvbXBfdGFibGUkZGVzZXFfbG9nZmMsIG1ldGhvZD0ic3BlYXJtYW4iKQoKY29tcF90YWJsZSA8LSB3dGNmX25vYmF0Y2hfd3R3aG9sZV90YWJsZXMkZGF0YVtbc2V0X25hbWVdXQpjb3IudGVzdChjb21wX3RhYmxlJGxvZzJmYywgY29tcF90YWJsZSRkZXNlcV9sb2dmYywgbWV0aG9kPSJzcGVhcm1hbiIpCmBgYAoKIyMgZGVsdGFfY2YgdnMgZGVsdGFfd2hvbGUKCiMjIyBPdXIgbGlicmFyaWVzCgpgYGB7ciBvdXJfY29tYmluZV9kZWx0YWNmX2RlbHRhd2hvbGUsIGZpZy5zaG93PSdoaWRlJ30Ka2VlcGVycyA8LSBsaXN0KAogICJkZWx0YWNmX3ZzX2RlbHRhd2hvbGUiID0gYygiZGVsdGFfY2YiLCAiZGVsdGFfd2hvbGUiKSkKc2V0X25hbWUgPC0gImRlbHRhX2NmX3ZzX2RlbHRhX3dob2xlIgpyb3duYW1lcyhvdXJfcmVzdWx0c1tbc2V0X25hbWVdXSkgPC0gb3VyX3Jlc3VsdHNbW3NldF9uYW1lXV1bWyJQcm90ZWluIl1dCmRlbHRhY2ZfZGVsdGF3aG9sZV90YWJsZXMgPC0gc20oY29tYmluZV9kZV90YWJsZXMoCiAgb3VyX3BhaXJ3aXNlX2NvbXAsIGtlZXBlcnM9a2VlcGVycywgZXh0cmFfYW5ub3Q9b3VyX3Jlc3VsdHNbW3NldF9uYW1lXV0sCiAgZXhjZWw9cGFzdGUwKCJleGNlbC9vdXJfZGVsdGFjZl92c19kZWx0YXdob2xlX3RhYmxlcy12IiwgdmVyLCAiLnhsc3giKSkpCmBgYAoKIyMjIFRiIGxpYnJhcmllcwoKYGBge3IgdGJfY29tYmluZV9kZWx0YWNmX2RlbHRhd2hvbGUsIGZpZy5zaG93PSdoaWRlJ30Ka2VlcGVycyA8LSBsaXN0KAogICJkZWx0YWNmX3ZzX2RlbHRhd2hvbGUiID0gYygiZGVsdGFfY2YiLCAiZGVsdGFfd2hvbGUiKSkKc2V0X25hbWUgPC0gImRlbHRhX2NmX3ZzX2RlbHRhX3dob2xlIgpyb3duYW1lcyh0Yl9yZXN1bHRzW1tzZXRfbmFtZV1dKSA8LSB0Yl9yZXN1bHRzW1tzZXRfbmFtZV1dW1siUHJvdGVpbiJdXQpkZWx0YWNmX2RlbHRhd2hvbGVfdGFibGVzIDwtIHNtKGNvbWJpbmVfZGVfdGFibGVzKAogIHRiX3BhaXJ3aXNlX2NvbXAsIGtlZXBlcnM9a2VlcGVycywgZXh0cmFfYW5ub3Q9dGJfcmVzdWx0c1tbc2V0X25hbWVdXSwKICBleGNlbD1wYXN0ZTAoImV4Y2VsL3RiX2RlbHRhY2ZfdnNfZGVsdGF3aG9sZV90YWJsZXMtdiIsIHZlciwgIi54bHN4IikpKQpgYGAKCiMjIGNvbXBfY2YgdnMgY29tcF93aG9sZQoKIyMjIE91ciBsaWJyYXJpZXMKCmBgYHtyIG91cl9jb21iaW5lX2NvbXBjZl9jb21wd2hvbGUsIGZpZy5zaG93PSdoaWRlJ30Ka2VlcGVycyA8LSBsaXN0KAogICJjb21wY2ZfdnNfY29tcHdob2xlIiA9IGMoImNvbXBfY2YiLCAiY29tcF93aG9sZSIpKQpzZXRfbmFtZSA8LSAiY29tcF9jZl92c19jb21wX3dob2xlIgpyb3duYW1lcyhvdXJfcmVzdWx0c1tbc2V0X25hbWVdXSkgPC0gb3VyX3Jlc3VsdHNbW3NldF9uYW1lXV1bWyJQcm90ZWluIl1dCmNvbXBjZl9jb21wd2hvbGVfdGFibGVzIDwtIHNtKGNvbWJpbmVfZGVfdGFibGVzKAogIG91cl9wYWlyd2lzZV9jb21wLCBrZWVwZXJzPWtlZXBlcnMsIGV4dHJhX2Fubm90PW91cl9yZXN1bHRzW1tzZXRfbmFtZV1dLAogIGV4Y2VsPXBhc3RlMCgiZXhjZWwvb3VyX2NvbXBjZl92c19jb21wd2hvbGVfdGFibGVzLXYiLCB2ZXIsICIueGxzeCIpKSkKYGBgCgojIyMgVGIgbGlicmFyaWVzCgpgYGB7ciB0Yl9jb21iaW5lX2NvbXBjZl9jb21wd2hvbGUsIGZpZy5zaG93PSdoaWRlJ30Ka2VlcGVycyA8LSBsaXN0KAogICJjb21wY2ZfdnNfY29tcHdob2xlIiA9IGMoImNvbXBfY2YiLCAiY29tcF93aG9sZSIpKQpzZXRfbmFtZSA8LSAiY29tcF9jZl92c19jb21wX3dob2xlIgpyb3duYW1lcyh0Yl9yZXN1bHRzW1tzZXRfbmFtZV1dKSA8LSB0Yl9yZXN1bHRzW1tzZXRfbmFtZV1dW1siUHJvdGVpbiJdXQpjb21wY2ZfY29tcHdob2xlX3RhYmxlcyA8LSBzbShjb21iaW5lX2RlX3RhYmxlcygKICB0Yl9wYWlyd2lzZV9jb21wLCBrZWVwZXJzPWtlZXBlcnMsIGV4dHJhX2Fubm90PXRiX3Jlc3VsdHNbW3NldF9uYW1lXV0sCiAgZXhjZWw9cGFzdGUwKCJleGNlbC90Yl9jb21wY2ZfdnNfY29tcHdob2xlX3RhYmxlcy12IiwgdmVyLCAiLnhsc3giKSkpCmBgYAoKIyMgZGVsdGFfY2YgdnMgd3RfY2YKCiMjIyBPdXIgbGlicmFyaWVzCgpgYGB7ciBvdXJfY29tYmluZV9kZWx0YWNmX3d0Y2YsIGZpZy5zaG93PSdoaWRlJ30Ka2VlcGVycyA8LSBsaXN0KAogICJkZWx0YWNmX3ZzX3d0Y2YiID0gYygiZGVsdGFfY2YiLCAid3RfY2YiKSkKc2V0X25hbWUgPC0gImRlbHRhX2NmX3ZzX3d0X2NmIgpyb3duYW1lcyhvdXJfcmVzdWx0c1tbc2V0X25hbWVdXSkgPC0gb3VyX3Jlc3VsdHNbW3NldF9uYW1lXV1bWyJQcm90ZWluIl1dCmRlbHRhY2Zfd3RjZl90YWJsZXMgPC0gc20oY29tYmluZV9kZV90YWJsZXMoCiAgb3VyX3BhaXJ3aXNlX2NvbXAsIGtlZXBlcnM9a2VlcGVycywgZXh0cmFfYW5ub3Q9b3VyX3Jlc3VsdHNbW3NldF9uYW1lXV0sCiAgZXhjZWw9cGFzdGUwKCJleGNlbC9vdXJfZGVsdGFjZl92c193dGNmX3RhYmxlcy12IiwgdmVyLCAiLnhsc3giKSkpCmBgYAoKIyMjIFRiIGxpYnJhcmllcwoKYGBge3IgdGJfY29tYmluZV9kZWx0YWNmX3d0Y2YsIGZpZy5zaG93PSdoaWRlJ30Ka2VlcGVycyA8LSBsaXN0KAogICJkZWx0YWNmX3ZzX3d0Y2YiID0gYygiZGVsdGFfY2YiLCAid3RfY2YiKSkKc2V0X25hbWUgPC0gImRlbHRhX2NmX3ZzX3d0X2NmIgpyb3duYW1lcyh0Yl9yZXN1bHRzW1tzZXRfbmFtZV1dKSA8LSB0Yl9yZXN1bHRzW1tzZXRfbmFtZV1dW1siUHJvdGVpbiJdXQpkZWx0YWNmX3d0Y2ZfdGFibGVzIDwtIHNtKGNvbWJpbmVfZGVfdGFibGVzKAogIHRiX3BhaXJ3aXNlX2NvbXAsIGtlZXBlcnM9a2VlcGVycywgZXh0cmFfYW5ub3Q9dGJfcmVzdWx0c1tbc2V0X25hbWVdXSwKICBleGNlbD1wYXN0ZTAoImV4Y2VsL3RiX2RlbHRhY2ZfdnNfd3RjZl90YWJsZXMtdiIsIHZlciwgIi54bHN4IikpKQpgYGAKCiMjIGNvbXBfY2YgdnMgd3RfY2YKCiMjIyBPdXIgbGlicmFyaWVzCgpgYGB7ciBvdXJfY29tYmluZV9jb21wY2Zfd3RjZiwgZmlnLnNob3c9J2hpZGUnfQprZWVwZXJzIDwtIGxpc3QoCiAgImNvbXBjZl92c193dGNmIiA9IGMoImNvbXBfY2YiLCAid3RfY2YiKSkKc2V0X25hbWUgPC0gImNvbXBfY2ZfdnNfd3RfY2YiCnJvd25hbWVzKG91cl9yZXN1bHRzW1tzZXRfbmFtZV1dKSA8LSBvdXJfcmVzdWx0c1tbc2V0X25hbWVdXVtbIlByb3RlaW4iXV0KY29tcGNmX3d0Y2ZfdGFibGVzIDwtIHNtKGNvbWJpbmVfZGVfdGFibGVzKAogIG91cl9wYWlyd2lzZV9jb21wLCBrZWVwZXJzPWtlZXBlcnMsIGV4dHJhX2Fubm90PW91cl9yZXN1bHRzW1tzZXRfbmFtZV1dLAogIGV4Y2VsPXBhc3RlMCgiZXhjZWwvb3VyX2NvbXBjZl92c193dGNmX3RhYmxlcy12IiwgdmVyLCAiLnhsc3giKSkpCmBgYAoKIyMjIFRiIGxpYnJhcmllcwoKYGBge3IgdGJfY29tYmluZV9jb21wY2Zfd3RjZiwgZmlnLnNob3c9J2hpZGUnfQprZWVwZXJzIDwtIGxpc3QoCiAgImNvbXBjZl92c193dGNmIiA9IGMoImNvbXBfY2YiLCAid3RfY2YiKSkKc2V0X25hbWUgPC0gImNvbXBfY2ZfdnNfd3RfY2YiCnJvd25hbWVzKHRiX3Jlc3VsdHNbW3NldF9uYW1lXV0pIDwtIHRiX3Jlc3VsdHNbW3NldF9uYW1lXV1bWyJQcm90ZWluIl1dCmNvbXBjZl93dGNmX3RhYmxlcyA8LSBzbShjb21iaW5lX2RlX3RhYmxlcygKICB0Yl9wYWlyd2lzZV9jb21wLCBrZWVwZXJzPWtlZXBlcnMsIGV4dHJhX2Fubm90PXRiX3Jlc3VsdHNbW3NldF9uYW1lXV0sCiAgZXhjZWw9cGFzdGUwKCJleGNlbC90Yl9jb21wY2ZfdnNfd3RjZl90YWJsZXMtdiIsIHZlciwgIi54bHN4IikpKQpgYGAKCiMjIGRlbHRhX3dob2xlIHZzIHd0X3dob2xlCgojIyMgT3VyIGxpYnJhcmllcwoKYGBge3Igb3VyX2NvbWJpbmVfZGVsYXdob2xlX3d0d2hvbGUsIGZpZy5zaG93PSdoaWRlJ30Ka2VlcGVycyA8LSBsaXN0KAogICJ3dGNmX3ZzX3d0d2hvbGUiID0gYygid3RfY2YiLCAid3Rfd2hvbGUiKSkKc2V0X25hbWUgPC0gImRlbHRhX3dob2xlX3ZzX3d0X3dob2xlIgpyb3duYW1lcyhvdXJfcmVzdWx0c1tbc2V0X25hbWVdXSkgPC0gb3VyX3Jlc3VsdHNbW3NldF9uYW1lXV1bWyJQcm90ZWluIl1dCnd0Y2Zfd3R3aG9sZV90YWJsZXMgPC0gc20oY29tYmluZV9kZV90YWJsZXMoCiAgb3VyX3BhaXJ3aXNlX2NvbXAsIGtlZXBlcnM9a2VlcGVycywgZXh0cmFfYW5ub3Q9b3VyX3Jlc3VsdHNbW3NldF9uYW1lXV0sCiAgZXhjZWw9cGFzdGUwKCJleGNlbC9vdXJfZGVsdGF3aG9sZV92c193dHdob2xlX3RhYmxlcy12IiwgdmVyLCAiLnhsc3giKSkpCmBgYAoKIyMjIFRiIGxpYnJhcmllcwoKYGBge3IgdGJfY29tYmluZV9kZWxhd2hvbGVfd3R3aG9sZSwgZmlnLnNob3c9J2hpZGUnfQprZWVwZXJzIDwtIGxpc3QoCiAgInd0Y2ZfdnNfd3R3aG9sZSIgPSBjKCJ3dF9jZiIsICJ3dF93aG9sZSIpKQpzZXRfbmFtZSA8LSAiZGVsdGFfd2hvbGVfdnNfd3Rfd2hvbGUiCnJvd25hbWVzKHRiX3Jlc3VsdHNbW3NldF9uYW1lXV0pIDwtIHRiX3Jlc3VsdHNbW3NldF9uYW1lXV1bWyJQcm90ZWluIl1dCnd0Y2Zfd3R3aG9sZV90YWJsZXMgPC0gc20oY29tYmluZV9kZV90YWJsZXMoCiAgdGJfcGFpcndpc2VfY29tcCwga2VlcGVycz1rZWVwZXJzLCBleHRyYV9hbm5vdD10Yl9yZXN1bHRzW1tzZXRfbmFtZV1dLAogIGV4Y2VsPXBhc3RlMCgiZXhjZWwvdGJfZGVsdGF3aG9sZV92c193dHdob2xlX3RhYmxlcy12IiwgdmVyLCAiLnhsc3giKSkpCmBgYAoKIyMgY29tcF93aG9sZSB2cyB3dF93aG9sZQoKIyMjIE91ciBsaWJyYXJpZXMKCmBgYHtyIG91cl9jb21iaW5lX2NvbXB3aG9sZV93dHdob2xlLCBmaWcuc2hvdz0naGlkZSd9CmtlZXBlcnMgPC0gbGlzdCgKICAiY29tcHdob2xlX3ZzX3d0d2hvbGUiID0gYygiY29tcF93aG9sZSIsICJ3dF93aG9sZSIpKQpzZXRfbmFtZSA8LSAiY29tcF93aG9sZV92c193dF93aG9sZSIKcm93bmFtZXMob3VyX3Jlc3VsdHNbW3NldF9uYW1lXV0pIDwtIG91cl9yZXN1bHRzW1tzZXRfbmFtZV1dW1siUHJvdGVpbiJdXQpjb21wd2hvbGVfd3R3aG9sZV90YWJsZXMgPC0gc20oY29tYmluZV9kZV90YWJsZXMoCiAgb3VyX3BhaXJ3aXNlX2NvbXAsIGtlZXBlcnM9a2VlcGVycywgZXh0cmFfYW5ub3Q9b3VyX3Jlc3VsdHNbW3NldF9uYW1lXV0sCiAgZXhjZWw9cGFzdGUwKCJleGNlbC9vdXJfY29tcHdob2xlX3ZzX3d0d2hvbGVfdGFibGVzLXYiLCB2ZXIsICIueGxzeCIpKSkKYGBgCgojIyMgVGIgbGlicmFyaWVzCgpgYGB7ciB0Yl9jb21iaW5lX2NvbXB3aG9sZV93dHdob2xlLCBmaWcuc2hvdz0naGlkZSd9CmtlZXBlcnMgPC0gbGlzdCgKICAiY29tcHdob2xlX3ZzX3d0d2hvbGUiID0gYygiY29tcF93aG9sZSIsICJ3dF93aG9sZSIpKQpzZXRfbmFtZSA8LSAiY29tcF93aG9sZV92c193dF93aG9sZSIKcm93bmFtZXModGJfcmVzdWx0c1tbc2V0X25hbWVdXSkgPC0gdGJfcmVzdWx0c1tbc2V0X25hbWVdXVtbIlByb3RlaW4iXV0KY29tcHdob2xlX3d0d2hvbGVfdGFibGVzIDwtIHNtKGNvbWJpbmVfZGVfdGFibGVzKAogIHRiX3BhaXJ3aXNlX2NvbXAsIGtlZXBlcnM9a2VlcGVycywgZXh0cmFfYW5ub3Q9dGJfcmVzdWx0c1tbc2V0X25hbWVdXSwKICBleGNlbD1wYXN0ZTAoImV4Y2VsL3RiX2NvbXB3aG9sZV92c193dHdob2xlX3RhYmxlcy12IiwgdmVyLCAiLnhsc3giKSkpCmBgYAoKYGBge3IgY29tcGFyZV9vdXJfdGJ9Cm91cl90YWJsZSA8LSBvdXJfcGFpcndpc2VfY29tcCRkZXNlcSRhbGxfdGFibGVzJHd0X3dob2xlX3ZzX3d0X2NmCnRiX3RhYmxlIDwtIHRiX3BhaXJ3aXNlX2NvbXAkZGVzZXEkYWxsX3RhYmxlcyR3dF93aG9sZV92c193dF9jZgpvdXJfb3JkZXIgPC0gb3JkZXIob3VyX3RhYmxlW1sibGltbWFfbG9nZmMiXV0pCnJvd25hbWVzKHRiX3RhYmxlKSA8LSBnc3ViKHBhdHRlcm49Il4oLiopXy4qIiwgcmVwbGFjZW1lbnQ9IlxcMSIsIHg9cm93bmFtZXModGJfdGFibGUpKQoKdGVzdF90YWJsZSA8LSBtZXJnZShvdXJfdGFibGUsIHRiX3RhYmxlLCBieT0icm93Lm5hbWVzIikKY29yLnRlc3QodGVzdF90YWJsZSRsb2dGQy54LCB0ZXN0X3RhYmxlJGxvZ0ZDLnkpCmBgYAoKIyBUT0RPCgoqIDIwMTgtMDQtMTA6ICBNYWtlIHN1cmUgbXkgaW52b2NhdGlvbnMgb2YgU1dBVEgyc3RhdHMvTVNzdGF0cyBhcmUgY29ycmVjdC4KCmBgYHtyIHNhdmVtZSwgZXZhbD1GQUxTRX0KaWYgKCFpc1RSVUUoZ2V0MCgic2tpcF9sb2FkIikpKSB7CiAgbWVzc2FnZShwYXN0ZTAoIlRoaXMgaXMgaHBnbHRvb2xzIGNvbW1pdDogIiwgZ2V0X2dpdF9jb21taXQoKSkpCiAgdGhpc19zYXZlIDwtIHBhc3RlMChnc3ViKHBhdHRlcm49IlxcLlJtZCIsIHJlcGxhY2U9IiIsIHg9cm1kX2ZpbGUpLCAiLXYiLCB2ZXIsICIucmRhLnh6IikKICBtZXNzYWdlKHBhc3RlMCgiU2F2aW5nIHRvICIsIHRoaXNfc2F2ZSkpCiAgdG1wIDwtIHNtKHNhdmVtZShmaWxlbmFtZT10aGlzX3NhdmUpKQogIHBhbmRlcjo6cGFuZGVyKHNlc3Npb25JbmZvKCkpCn0KYGBgCg==

M.tuberculosis 20180611: Analyzing data from OpenSwathWorkFlow/TRIC.

atb abelew@gmail.com

2018-06-14