EB3I n1 2025 scRNAseq
-
PRE-PROCESSING (I)
-
Load a count matrix, empty droplets & ambient RNA filtering

1 PREAMBLE

1.1 Purpose of this session

This file describes the different steps to perform the first part of the single cell RNAseq data analysis training course for the EB3I n1 2025, covering these steps :

  • Load a single cell raw counts matrix

  • Remove barcodes from empty droplets

  • Estimate and remove ambient RNA contamination (“soup”)

1.2 An important notice

  • These early (but mandatory) steps of the analysis are not covered by the renowned, higher-level analysis frameworks like Seurat nor scran/scater.

  • Thus, we will write/use a certain amount of lines of code, some of which of a higher complexity for R beginners.

  • Newcomer, do not fear !

    • While you are welcome to try writing some by yourself, depending on your R knowledge and self-confidence in it …

    • … you do not have to know how these code work …

    • … thus you should not feel ashamed copy-pasting the pre-written code hereafter …

    • … but in both cases, please try your best understanding what these code do, and even more, why !

    • In this purpose, you may use whatever help you can by searching.

1.3 What to run ?

  • This training presentation, written in Rmarkdown, contains many chunks ( = code blocks)

  • You do not have to run them all !

  • Chunks follow a color code :

1.3.1 Chunks you ARE expected to run :

# chunk_run

## Example of chunk to run
getwd()
Show output 
[1] "/shared/projects/2538_eb3i_n1_2025/atelier_scrnaseq/RMD/Preproc.1"


1.3.2 Chunks you ARE NOT expected to run

Hopefully, we will demo them. You may run them by yourself, young padawan, but if you get late or worst, break something, you’ll be the one to blame !

# chunk_demo

## Example of chunk to NOT run
utils::sessionInfo()
Show output 
R version 4.4.1 (2024-06-14)
Platform: x86_64-conda-linux-gnu
Running under: Ubuntu 22.04.5 LTS

Matrix products: default
BLAS/LAPACK: /shared/ifbstor1/software/miniconda/envs/r-4.4.1/lib/libopenblasp-r0.3.29.so;  LAPACK version 3.12.0

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

time zone: Europe/Paris
tzcode source: system (glibc)

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

loaded via a namespace (and not attached):
 [1] digest_0.6.37     R6_2.6.1          bookdown_0.39     fastmap_1.2.0    
 [5] xfun_0.52         cachem_1.1.0      knitr_1.50        htmltools_0.5.8.1
 [9] rmarkdown_2.29    lifecycle_1.0.4   cli_3.6.5         sass_0.4.10      
[13] jquerylib_0.1.4   compiler_4.4.1    rstudioapi_0.17.1 tools_4.4.1      
[17] evaluate_1.0.3    bslib_0.9.0       rmdformats_1.0.4  yaml_2.3.10      
[21] jsonlite_2.0.0    rlang_1.1.6


1.3.3 Questions

Some questions we would like you to answer during the training course. Please be fair and do not take a look at the answer before we tell you so :)

# chunk_question

## What's the Answer to Life ? The Universe ? Everything ??


1.3.4 Answers

They usually follow questions…

# chunk_answer

cat("The Answer to the Great Question... is ... 
Forty-two. Six by nine : forty-two. That's it. 
That's all there is. (I always thought something was 
fundamentally wrong with the universe.)")
Show output 
The Answer to the Great Question... is ... 
Forty-two. Six by nine : forty-two. That's it. 
That's all there is. (I always thought something was 
fundamentally wrong with the universe.)


1.3.5 Additional exercises / questions

These are beyond the scope of the training course, but you may play it if you feel bored

# chunk_beyond

## Example of chunk to play with but beyond the scope
getwd()
Show output 
[1] "/shared/projects/2538_eb3i_n1_2025/atelier_scrnaseq/RMD/Preproc.1"


2 Start Rstudio

3 Warm-up

  • We now set common parameters as new variables, once and for all for this session :
# setparam


## Set your project name
# WARNING : Do not just copy-paste this ! It's MY project name ! Put YOURS !!
project_name <- "ebaii_sc_teachers"


## Control if the project_name exists on the cluster
cat('PATH CHECK : ', dir.exists(paste0('/shared/projects/', project_name)))
Show output 
PATH CHECK :  TRUE
## Seed for the RNG
my_seed <- 1337L

## Empty droplets max p-value
max_p <- 1E-03

4 Prepare the data structure

  • During the current and remaining training sessions, we will work on different data objects, some retrieved from external resources, some we will create.

  • In order to reuse these for further steps, or other sessions, we need to store them, thus create a directory structure on the cluster.

  • We will create this structure in your project folder (the one you registered at the IFB Core, in which you stored your own data for the tutoring).

  • To ensure that you will always use the same structure for your work, the code hereafter has been writen so that you just have to define a single variable that corresponds to your project name.

    • In MY case, as a SC training teacher, it is ebaii_sc_teachers.

4.1 Main directory

  • Create the training course directory
# maindir

## Prepare the path in a "TD_dir" variable
TD_dir <- paste0("/shared/projects/", project_name, "/SC_TD")

## Create the root directory
dir.create(path = TD_dir, recursive = TRUE)

## Print the root directory on-screen
print(TD_dir)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD"

4.2 Current session

# sessiondir

## Create the session (Preproc.1) directory
session_dir <- paste0(TD_dir, "/01_Preproc.1")
dir.create(path = session_dir, recursive = TRUE)

## Print the session directory on-screen
print(session_dir)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/01_Preproc.1"

4.3 Input data directory

# indir

## Create the INPUT data directory
input_dir <- paste0(session_dir, "/DATA")
dir.create(path = input_dir, recursive = TRUE)

## Print the input directory on-screen
print(input_dir)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/01_Preproc.1/DATA"

4.4 Output results directory

# outdir

## Create the OUTPUT data directory
output_dir <- paste0(session_dir, "/RESULTS")
dir.create(path = output_dir, recursive = TRUE)

## Print the output directory on-screen
print(output_dir)
[1] "/shared/projects/ebaii_sc_teachers/SC_TD/01_Preproc.1/RESULTS"



Note : Actually, for this session we won’t export any result, so creating this RESULTS folder was just for practice purpose.


Now, the we have everything needed to :

  • load data (input_dir)

  • save the results we will generate (output_dir).

5 Additional resources

However, some upcoming steps require a bit complicated methods that are not fun to write, nor copy-paste multiple times. So, we will now prepare code functions that will ease these steps. After executing these functions code, we may call (invoke) them when needed.

5.1 Functions

5.1.1 SoupX helper function

This function was written to ease the use of SoupX on a count matrix It was tested on SoupX=1.6.2 with igraph=1.5.1, Seurat=4.4.0,5.1.0

Parameters :

  • scmat_filt : (sparse)matrix corresponding to an empty droplets -filtered count matrix

  • scmat_raw : (sparse)matrix corresponding to an NON- empty droplets -filtered count matrix

  • soupQuantile, contaminationRange : see ?SoupX::autoEstCont

  • contaminationRange : range of the expected soup proportion

  • soupRange : estimate soup fraction from features with total reads comprised in this range of counts (only used when scmat_raw != NULL)

  • return_object : if TRUE, return the “unsouped” count matrix ; if FALSE, only perform soup estimation and return the estimated proportion value (rho)

  • doPlot : Perform the SoupX estimation plot (rho distribution)

# soupx_func

SoupX_auto <- function(scmat_filt = NULL, scmat_raw = NULL, soupQuantile = 0.9,
                       contaminationRange = c(.01, .8), soupRange = c(0,100),
                       return_object = FALSE, doPlot = FALSE) {
  
  ## Checks
  if(is.null(scmat_filt)) stop('A filtered count matrix is required !')
  
  if(is.null(scmat_raw)) message('No unfiltered raw counts matrix provided. Estimation will be based on filtered matrix only.')
  
  ## If no raw matrix
  if (is.null(scmat_raw)) {
    spChanRaw <- SoupX::SoupChannel(
      tod = scmat_filt, 
      toc = scmat_filt, 
      calcSoupProfile = FALSE)
    sc_rowsum <- sparseMatrixStats::rowSums2(scmat_filt)
    spProf <- data.frame(
      row.names = rownames(scmat_filt), 
      est = sc_rowsum/sum(scmat_filt), 
      counts = sc_rowsum)
    spChan <- SoupX::setSoupProfile(spChanRaw, spProf)
  } else {
    spChan <- SoupX::SoupChannel(
      tod = scmat_raw, 
      toc = scmat_filt, 
      calcSoupProfile = FALSE)
    if (min(spChan$nDropUMIs) > max(soupRange)) stop(
      'Minimum found counts per barcode is : ', 
      min(spChan$nDropUMIs), 
      ', which is smaller than the upper bound of soupRange ! Please increase soupRange max !')
    spChan <- SoupX::estimateSoup(sc = spChan, soupRange = soupRange)
  }
  ## Display Top 20 contributing genes
  if (!return_object) {
    cat('\nSoup-contributing features (Top 20) :\n')
    print(knitr::kable(head(
      spChan$soupProfile[order(spChan$soupProfile$est, decreasing = TRUE), ], 
      n = 20)))
  }
  ## Quick clustering needed
  spClust <- scran::quickCluster(scmat_filt, method = "igraph")
  ## Adding clusters to the SoupChannel object
  spChan <- SoupX::setClusters(sc = spChan, clusters = spClust)
  ## Estimating soup
  sX <- SoupX::autoEstCont(sc = spChan, doPlot = doPlot, tfidfMin = 1, 
                           soupQuantile = soupQuantile, maxMarkers = 100, 
                           contaminationRange = contaminationRange, 
                           rhoMaxFDR = .2, priorRho = .05, priorRhoStdDev = .1, 
                           forceAccept = FALSE)
  
  ## Removing soup (adjusting counts)
  if(return_object) {
    cat('Counts BEFORE SoupX : ', sum(scmat_filt), '\n')
    scmat_soupx <- SoupX::adjustCounts(
      sX, method = 'subtraction', roundToInt = TRUE, 
      tol = .001, pCut = .01)
    cat('Counts AFTER SoupX : ', sum(scmat_soupx), '\n')
    rm(scmat_filt)
    return(scmat_soupx)
  } else return(sX$fit$rhoEst)
}


6 Load a 10X Cell Ranger raw count matrix

6.1 The dataset

  • The data we will work on is hosted on Zenodo. To ease its retrieval, we can define its ID as a variable

    # zid

## This is actually to help me updating this training without hassle
zen_id <- '14034221'

## This is the path to the current EB3I backup
sessionid <- '2538_eb3i_n1_2025'

  • Here will we train ourselves to load into R a single cell RNAseq data produced by 10X Genomics’ software Cell Ranger.

  • We will work with a public dataset provided by the manufacturer, that consists in ~ 10,000 PBMC (peripheral bone marrow cells) from a human donor (downsampled to 1/10th, for smaller computation time).

  • The experiment was performed with the 3’ capture kit v3.

  • The analysis was performed with Cell Ranger v3, mapping on the GRCh38-2020-A manufacturer reference.

  • The data consists into the typical 10x 3-files structure, hosted in a Zenodo respository (Id : 14034221)

    • barcodes.tsv.gz : the list of putative droplet barcodes
    • features.tsv.gz : the list of features (genes)
    • matrix.mtx.gz : the feature x barcode expression count matrix

6.2 Download data

  • The relatively complex code chunk hereafter has been written to :

    • let you download the data files from :

      • Zenodo by default

      • if backup = TRUE, a local backup instead (in case we loose the internet connection)

    • if the requested data are already available locally :

      • load this from the local source

      • if force = TRUE a new download is forced

# dlzen10x

### Named files (will be used later on !)
mtx_file <- "matrix.mtx.gz"
features_file <- "features.tsv.gz"
barcodes_file <- "barcodes.tsv.gz"
summary_file <- "pbmc_10k_v3_summary.html"

## Filename(s) to retrieve
toget_files <- c(mtx_file,
                 features_file,
                 barcodes_file,
                 summary_file)

## Folder to store retrieved files
local_folder <- input_dir

## Use local backup ?
backup <- FALSE
if(backup) message("Using local backup !")

## Force download ?
force <- FALSE
if(force) message("Forcing (re)download !")

### Define remote folder
remote_folder <- if (backup) paste("/shared/projects/", sessionid, "/atelier_scrnaseq/TD/BACKUP/10X/") else paste0("https://zenodo.org/records/", zen_id, "/files/")

### Reconstruct the input paths
remote_path <- paste0(remote_folder, "/", toget_files)

### Reconstruct the output paths
local_path <- paste0(local_folder, "/", toget_files)

## Retrieve files (if they don't exist), in loop
for (tg in seq_along(toget_files)) {
  ## If the file does not locally exist
  if (!file.exists(local_path[tg]) | force) {
    ## Retrieve data
    if(backup) {
      file.copy(from = remote_path[tg],
                to = local_path[tg])
    } else {
      download.file(url = remote_path[tg], 
                    destfile = local_path[tg])
    }
    ## Check if downloaded files exist locally
    if(file.exists(local_path[tg])) message("\tOK")
  } else message(paste0(toget_files[tg], " already downloaded !"))
}
matrix.mtx.gz already downloaded !
features.tsv.gz already downloaded !
barcodes.tsv.gz already downloaded !
pbmc_10k_v3_summary.html already downloaded !

6.3 Load into R

  • Load the 10X data files into R :

    # q_r10X

How ?
    • A good starter

    Show output

     
    https://lmddgtfy.net/?q=load%2010x%20data%20in%20R
    • Back to your local help ressource
    # a_r10X2

## Reading the function help page
?Seurat::Read10X



    We can now apply it to our data :

    # load10x

## Loading unfiltered 10X matrix
scmat <- Seurat::Read10X(data.dir = input_dir)


  • However, we don’t know what the loaded object actually looks like in R…

  • Questions :

    • Do you know a way to …
    # q_whatis

... know what IS the type of object that was created ?


    # a_whatis

## To know the type of an R object : methods::is()
?methods::is
    # gois

## Get the type
methods::is(scmat)

    Show output

     
    [1] "dgCMatrix"     "CsparseMatrix" "dsparseMatrix" "generalMatrix"
[5] "AnyMatrix"     "V3Matrix"      "dMatrix"       "sparseMatrix" 
[9] "Matrix"       



    # q_struc

... observe the STRucture of the resulting object ?


    # a_struc

## To get a basic structure description : utils::str()
?utils::str
    # go_str

## Get the basic structure
utils::str(scmat)

    Show output

     
    Formal class 'dgCMatrix' [package "Matrix"] with 6 slots
  ..@ i       : int [1:1718404] 7922 18295 23126 24791 31764 32866 13690 17058 2511 33306 ...
  ..@ p       : int [1:147457] 0 0 0 0 0 0 6 7 8 8 ...
  ..@ Dim     : int [1:2] 33694 147456
  ..@ Dimnames:List of 2
  .. ..$ : chr [1:33694] "RP11-34P13.3" "FAM138A" "OR4F5" "RP11-34P13.7" ...
  .. ..$ : chr [1:147456] "TCCCGATAGCCTTGAT-1" "AGATCTGAGACTGTAA-1" "GCTTGAAGTGCGAAAC-1" "ATCATGGGTATTCGTG-1" ...
  ..@ x       : num [1:1718404] 1 1 1 1 1 1 1 1 1 1 ...
  ..@ factors : list()


    Not that meaningful, though… At least we know we have a matrix, so one can get its dimensions !

    # scmat_dim

## Dimensions of a >1D object
dim(scmat)

    Show output

     
    [1]  33694 147456


6.4 Features cleanup

Seurat does not like at all underscores (_) in feature names …

We will have to :

  • control if our dataset contains some ? Here, we will use the base::grep() function that looks in character strings x if they contain a provided set of characters as a pattern, and returns their position (when value = FALSE, default) or value itself (value = TRUE)

    # us_check1

## Check feature names with underscore(s)
base::grep(pattern = '_', 
     x = rownames(scmat), 
     value = TRUE)

    Show output

     
     [1] "RP11-442N24__B.1" "RP11-99J16__A.2"  "RP11-59D5__B.2"   "RP11-445L13__B.3"
 [5] "RP11-544L8__B.4"  "XXyac-YX65C7_A.2" "XXyac-YX65C7_A.3" "RP11-524D16__A.3"
 [9] "XX-DJ76P10__A.2"  "RP11-1157N2__B.2" "RP1-213J1P__B.1"  "RP1-213J1P__B.2" 
[13] "RP4-633O19__A.1"  "RP4-754E20__A.5"  "CTA-280A3__B.2"


  • if so, clean it ! Here we will use the base::gsub() function that looks for a pattern in a character string x and substitutes it with a replacement character string.

    # us_clean

## Clean features : replacing underscores by dashes
rownames(scmat) <- base::gsub(
  pattern = "_", 
  replacement = "-", 
  x = rownames(scmat))

  • then, control the efficiency of our cleaning (with a new use of base::grep)

    # us_check2

## Post-cleanup control
base::grep(pattern = "_", 
           x = rownames(scmat), 
           value = TRUE)

    Show output

     
    character(0)


7 Empty droplets

  • While a large part of our barcodes contain no count at all, barcodes with at least one total count may not mandatorily correspond to true cells, due to ambient RNA.

  • Here we want to identify empty droplets, then remove them from the matrix.

7.1 Detection

  • To detect the margins of “true” cells versus empty droplets, we will rely on the so-called double-kneeplot : a plot of UMI counts per barcode, in function of their increasingly ordered rank.


  • To draw such plot, we need to rank the barcodes thanks to their total counts.

7.1.1 Ranking barcodes

  • We will use the DropletUtils package

    # h_barcodeRanks

## Read the function help page
?DropletUtils::barcodeRanks()
    # bcranks

## Generate the rank statistics
bc_rank <- DropletUtils::barcodeRanks(scmat)

  • Check the ranking of barcodes, according to their counts

    # bcrank_top

## Show a summary of the ranks (ordered for display convenience)
print(bc_rank[order(bc_rank$rank),])
    DataFrame with 147456 rows and 3 columns
                        rank     total    fitted
                   <numeric> <integer> <numeric>
ATCATGGCAGCCAATT-1         1     21195        NA
GGGCATCGTAGCTTGT-1         2     18193        NA
AGTGGGATCTTAACCT-1         3     16966        NA
GTCGGGTGTGCAGTAG-1         4     16383        NA
TTAGGCAAGCCGCCTA-1         5     15406        NA
...                      ...       ...       ...
TTGAACGTCAGAGGTG-1    100998         0        NA
ACGATGTAGTTGTCGT-1    100998         0        NA
GTCAAGTTCAGTTGAC-1    100998         0        NA
TCGCGAGGTTGCGCAC-1    100998         0        NA
TATCTCAGTTACCAGT-1    100998         0        NA


7.1.2 Kneeplot

Now, we can draw our kneeplot :

# kneeplot1

## Plot
graphics::plot(x = bc_rank$rank, 
               y = bc_rank$total + 1, 
               log = "xy", 
               xlab = "Barcode rank", 
               ylab = "Total UMIs", 
               col = "black", 
               pch = ".", 
               cex = 5, 
               main = "Kneeplot")


  • Question :

    # q_knee

Looking at the kneeplot, can you roughly predict how many barcodes will be 
kept as cells (ie, as non-empty droplets) ?


    # a_knee

## . The "cliff" in the kneeplot is at a 
##  rank value a bit below ~ 1,000.

Then, we can use these ranks to identify “true” cells :

# h_emptyDrops

## Read the function help page
?DropletUtils::emptyDrops

The function we will use relies on random number generation, so we will fix the RNG seed, for reproducibility

# emptydrops1

## Set the RNG seed
set.seed(my_seed)

## Identify empty droplets
bc_rank2 <- DropletUtils::emptyDrops(scmat)

## Whats is the structure of the returned output ?
str(bc_rank2)
Show output 
Formal class 'DFrame' [package "S4Vectors"] with 6 slots
  ..@ rownames       : chr [1:147456] "TCCCGATAGCCTTGAT-1" "AGATCTGAGACTGTAA-1" "GCTTGAAGTGCGAAAC-1" "ATCATGGGTATTCGTG-1" ...
  ..@ nrows          : int 147456
  ..@ elementType    : chr "ANY"
  ..@ elementMetadata: NULL
  ..@ metadata       :List of 5
  .. ..$ lower  : num 100
  .. ..$ niters : num 10000
  .. ..$ ambient: num [1:16738, 1] 8.43e-07 3.38e-05 9.91e-06 8.43e-07 8.43e-07 ...
  .. .. ..- attr(*, "dimnames")=List of 2
  .. .. .. ..$ : chr [1:16738] "RP11-34P13.7" "FO538757.2" "AP006222.2" "RP4-669L17.10" ...
  .. .. .. ..$ : NULL
  .. ..$ alpha  : num 4171
  .. ..$ retain : num 2611
  ..@ listData       :List of 5
  .. ..$ Total  : int [1:147456] 0 0 0 0 0 6 1 1 0 0 ...
  .. ..$ LogProb: num [1:147456] NA NA NA NA NA NA NA NA NA NA ...
  .. ..$ PValue : num [1:147456] NA NA NA NA NA NA NA NA NA NA ...
  .. ..$ Limited: logi [1:147456] NA NA NA NA NA NA ...
  .. ..$ FDR    : num [1:147456] NA NA NA NA NA NA NA NA NA NA ...

We can show a summary of the ordered barcodes :

# emptysmry

## Show a summary of the qualified ranks (ordered for display convenience)
print(bc_rank2[order(bc_rank2$Total, decreasing = TRUE),])
Show output 
DataFrame with 147456 rows and 5 columns
                       Total   LogProb    PValue   Limited       FDR
                   <integer> <numeric> <numeric> <logical> <numeric>
ATCATGGCAGCCAATT-1     21195  -13101.3 9.999e-05      TRUE         0
GGGCATCGTAGCTTGT-1     18193  -12055.5 9.999e-05      TRUE         0
AGTGGGATCTTAACCT-1     16966  -11594.8 9.999e-05      TRUE         0
GTCGGGTGTGCAGTAG-1     16383  -11994.6 9.999e-05      TRUE         0
TTAGGCAAGCCGCCTA-1     15406  -10687.9 9.999e-05      TRUE         0
...                      ...       ...       ...       ...       ...
TTGAACGTCAGAGGTG-1         0        NA        NA        NA        NA
ACGATGTAGTTGTCGT-1         0        NA        NA        NA        NA
GTCAAGTTCAGTTGAC-1         0        NA        NA        NA        NA
TCGCGAGGTTGCGCAC-1         0        NA        NA        NA        NA
TATCTCAGTTACCAGT-1         0        NA        NA        NA        NA


7.1.3 Selection

Barcodes considered as empty droplets have no p-value (FDR = NA)

# bcr2_na

## Creating a validation column (empty droplets have NA as p-value)
bc_rank2$VALID <- !is.na(bc_rank2$FDR) 

## Quantify "real cells"
table(bc_rank2$VALID)
Show output 

 FALSE   TRUE 
146397   1059


We can control the stringency by setting a FDR-adjusted p-value threshold

# valid_fdr

## Recall the max p-value we set as a variable early on
max_p
Show output 
[1] 0.001
## We've set the max p-value to 1E-03
bc_rank2$VALID[bc_rank2$FDR >= max_p] <- FALSE

## Control
table(bc_rank2$VALID)
Show output 

 FALSE   TRUE 
146581    875


Now, we can display our true cells on the kneeplot :

# kneeplot2

## Plot with highlight of the selected barcodes (red)
graphics::plot(x = bc_rank$rank, 
               y = bc_rank$total+1, 
               log = "xy", 
               xlab = "Barcode rank", 
               ylab = "Total UMIs", 
               col = ifelse(bc_rank2$VALID, "red", "black"), 
               pch = ".", 
               cex = ifelse(bc_rank2$VALID, 7, 3), 
               main = paste0("Kneeplot (", 
                             length(which(bc_rank2$VALID)), 
                             " barcodes kept)"))
Show plot


  • Question :

    # q_kneered

Is there something unexpected for the retained "TRUE" cells (red) in this kneeplot ?


    # a_kneered

## . You may observe that the selection of "true" cells (red)
##   does not strictly correspond to a "cut" in the kneeplot
##   descending curve.
##
## . This is because emptyDrops does not only "follow" the curve
##   to determine a threshold corresponding to a "minimal count" 
##   value to consider a barcode as a cell, but also performs a 
##   statistical analysis for each barcode, considering how much 
##   its expression profile ressembles the one of other cells, 
##   even with a very different (lower) global level of expression.



7.2 Removal

We just need to restrict our count matrix to “true” cells.

  • How many barcodes do we have at the moment ?

    # scdim3

dim(scmat)

    Show output

     
    [1]  33694 147456


  • Filter empty droplets out :

    # edfilter

## Restrict to valid barcodes
scmat_cells <- scmat[, bc_rank2$VALID]

## Control
dim(scmat_cells)

    Show output

     
    [1] 33694   875


8 Ambient RNA

  • The counts measured in this matrix of (now considered “true”) single cells do not reflect the measurement of these cells expression only

  • Due to other over-lysed / dying cells in the emulsion, we also measured fraction of ambient RNA that adds up to the “real” expression

  • We may estimate it by observing counts existing at a minimal level across cells that greatly differ in their expression profile

  • Just in case, due to the fun tool name we will use, SoupX, I’ll write down “soup” as an alias for “ambient RNA”, but both terms have the same meaning.


8.1 Estimate the “soup”

SoupX can estimate the rate and composition of ambient RNA, in several ways :

  • Manual mode : providing a list of features expected to reflect the soup (genes expressed at high levels in a majority of the expected cell types)

  • Automatic modes, using the empty droplets-filtered raw count matrix …

    • … along with the unfiltered matrix (ie, containing all quantified droplets). This is the most efficient mode.

    • … solely : this is less efficient, but useful when the unfiltered matrix is not available.

# h_SoupX_auto

## Reading the SoupX package main help page
?SoupX::SoupX

Let’s run SoupX, just asking for an estimation of the “soup” proportion :

# soupxrhoE

## TIP : In case of error :
# install.packages("irlba", type="source", force=TRUE)

## SoupX for ambient fraction estimation
soup_frac <- SoupX_auto(
  scmat_filt = scmat_cells, 
  scmat_raw = scmat)
Show output 

Soup-contributing features (Top 20) :


|       |       est| counts|
|:------|---------:|------:|
|MALAT1 | 0.0321703|  17473|
|B2M    | 0.0194002|  10537|
|TMSB4X | 0.0165850|   9008|
|EEF1A1 | 0.0129267|   7021|
|RPL21  | 0.0103251|   5608|
|RPS27  | 0.0100269|   5446|
|RPL13  | 0.0091542|   4972|
|RPL13A | 0.0087252|   4739|
|RPL10  | 0.0082852|   4500|
|RPLP1  | 0.0079225|   4303|
|RPL34  | 0.0078801|   4280|
|RPS12  | 0.0076481|   4154|
|RPS6   | 0.0075211|   4085|
|RPS18  | 0.0074530|   4048|
|RPS27A | 0.0073093|   3970|
|RPS2   | 0.0072965|   3963|
|RPS3A  | 0.0072909|   3960|
|RPL32  | 0.0072744|   3951|
|RPL41  | 0.0067865|   3686|
|RPL11  | 0.0063446|   3446|


NOTE : We can observe MALAT1 as the top first gene (as well as multiple riboprotein-coding genes).


What’s the estimated fraction of soup in our counts ?

# soupxfrac

cat("Soup fraction : ", soup_frac)
Show output 
Soup fraction :  0.054


  • Question :

    # q_souprate

Should we try to remove such an amount of ambient RNA ?


    # a_souprate

The best way to know is to try it out =D


8.2 Remove the “soup”

SoupX can remove the soup (subtract counts for identified soup genes to all barcodes/cells) using different methods :

  • Automatically (safer)

  • Specifying a fraction of counts to remove.

    • Pros :

      • By removing X % of reads :

        • if soup was effectively present (at a rate <= X), this will mostly affect soup genes first.

        • If not, it will decrease counts globally, thus have negligible effect in differences between cells / cell types.

    • Cons :

      • Counts are already rare in droplet-based single cell technologies… Removing them “by default” makes them even more rare…

      • If actual soup rate is higher than the defined rate X to remove, this mode will remain uneffective…

# soupxrm

## Run SoupX in "removal" mode
scmat_unsoup <- SoupX_auto(
  scmat_filt = scmat_cells, 
  scmat_raw = scmat, 
  return_object = TRUE)
Show output 
Counts BEFORE SoupX :  3766216
Show output 
Counts AFTER SoupX :  3563040


8.3 Effect of the “soup” removal

  • Now, we want to characterize the effect of the removal of the ambient RNA, by comparing the pre- and post- soupX matrices.

  • An easy way is to describe, then visualize such an effect.



8.3.1 Describe and visualize

8.3.1.1 Before SoupX (DEMO)

We can describe our count matrix thanks to some useful metrics :

  • Total counts
# descBS_totcount

## All counts
base::sum(scmat_cells)
Show output 
[1] 3766216
  • Sparsity (fraction of 0s)
# descBS_sparsity

## Compute sparsity (number of matrix values at 0 / cross-product of matrix dim)
base::sum(sparseMatrixStats::colCounts(
  x = scmat_cells, value = 0)
  ) / base::prod(dim(scmat_cells))
Show output 
[1] 0.9614275
  • Distribution of counts in cells
# descBS_nCount

## Sums up counts for each cell
base::summary(sparseMatrixStats::colSums2(x = scmat_cells, na.rm = TRUE))
Show output 
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
    116    3069    3957    4304    5014   21195
  • Distribution of expressed genes in cells
# descBS_nFeature

## All genes minus 0-count genes
base::summary(nrow(scmat_cells) - sparseMatrixStats::colCounts(
  x = scmat_cells, 
  value = 0, 
  na.rm = TRUE))
Show output 
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
     89    1052    1242    1300    1484    3752

8.3.1.2 After SoupX

We can describe our count matrix thanks to some useful metrics :

  • Total counts
# descAS_totcount

## All counts
base::sum(scmat_unsoup)
Show output 
[1] 3563040
  • Sparsity (fraction of 0s)
# descAS_sparsity

## Compute sparsity (number of matrix values at 0 / cross-product of matrix dim)
base::sum(sparseMatrixStats::colCounts(
  x = scmat_unsoup, value = 0)
  ) / base::prod(dim(scmat_unsoup))
Show output 
[1] 0.9635223
  • Distribution of counts in cells
# descAS_nCount

## Sums up counts for each cell
base::summary(sparseMatrixStats::colSums2(x = scmat_unsoup, na.rm = TRUE))
Show output 
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
    113    2912    3735    4072    4754   19797
  • Distribution of expressed genes in cells
# descAS_nFeature

## Compute sparsity (number of matrix values at 0 / cross-product of matrix dim)
base::summary(nrow(scmat_cells) - sparseMatrixStats::colCounts(
  x = scmat_unsoup, 
  value = 0, 
  na.rm = TRUE))
Show output 
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
     88    1002    1175    1229    1400    3435



  • Questions :

    # q_descdiff

Can you describe the difference (especially for counts and sparsity) ?

Was it expected ?


    # a_descdiff

## . Total counts, max value, counts per cell and number of
##   expressed features per cell, all have decreased.
##
## . The sparsity level has slightly increased (removing the soup
##   obliterated all counts for some features in some cells).
##
## . All of this is absolutely expected.


What are the features most affected by the soup removal ?

# soupfeatures

## Compute the fraction matrix (PRE / POST)
##  NOTE : we may have 0 counts in the divider, 
##  so we increment both matrix by +1 !
cont_frac <- (scmat_cells + 1) / (scmat_unsoup + 1)

## Fraction of counts decreased by soup removal, per feature (ordered)
feat_frac <- sort(sparseMatrixStats::rowMeans2(cont_frac), decreasing = TRUE)

### Display Top 5 features
utils::head(feat_frac, decreasing = TRUE)
Show output 
     LYZ   S100A9   S100A8  HLA-DRA     CD74     CST3 
1.478310 1.453317 1.384649 1.343360 1.279586 1.204137 
## Removing objects to free some RAM
rm(cont_frac, feat_frac)
# qsoupfeat

Does something about these genes strike you ?

Hint : fraction or difference ?

# soupfeaturesdiff

## Compute the DIFFERENCE matrix (PRE / POST)
##  NOTE : we may have 0 counts in the divider, 
##  so we increment both matrix by +1 !
cont_dif <- scmat_cells - scmat_unsoup

## Fraction of counts decreased by soup removal, per feature (ordered)
feat_dif <- sort(sparseMatrixStats::rowMeans2(cont_dif), decreasing = TRUE)

### Display Top 5 features
utils::head(feat_dif, decreasing = TRUE)
Show output 
  MALAT1      B2M   TMSB4X   EEF1A1    RPL21    RPS27 
7.533714 4.528000 3.885714 3.060571 2.435429 2.372571 
## Removing objects to free some RAM
rm(cont_dif, feat_dif)


8.3.2 Plot the “top 1” feature

We want to visualize the expression of the “top” soup gene, LYZ, in our dataset.

To do so, we will generate a projection of the data into a reduced (2-D) mathematical space. This is done by using a series of R function from the Seurat package, that we won’t describe here and now : understanding these commands below is not the scope of our current session, but it will be during the next few days.

Actually, these commands perform all the steps we’ll see in the next few days, in a row, but with default values (thus, almost certainly inadequate to our current data).

In consequence, please just consider the output as a way to quickly and dirtily visualize the data : the depicted topology should be considered as very rough and mostly imperfect.

8.3.2.1 Before SoupX

# bsoupx

## Roll process, using a temporary "sobj_presoupx" Seurat object
sobj_presoupx <- Seurat::CreateSeuratObject(counts = scmat_cells, project = "PBMC10K_BeforeSoupX")
sobj_presoupx <- Seurat::NormalizeData(object = sobj_presoupx, verbose = FALSE)
sobj_presoupx <- Seurat::ScaleData(object = sobj_presoupx, verbose = FALSE)
sobj_presoupx <- Seurat::FindVariableFeatures(object = sobj_presoupx, verbose = FALSE)
sobj_presoupx <- Seurat::RunPCA(object = sobj_presoupx, npcs = 21, verbose = FALSE)
sobj_presoupx <- Seurat::RunUMAP(object = sobj_presoupx, dims = c(1:20), verbose = FALSE)

## Top1 feature plot
fp_pre <- Seurat::FeaturePlot(object = sobj_presoupx, features = 'LYZ') + Seurat::DarkTheme()
print(fp_pre)
Show plot


8.3.2.2 After SoupX

# asoupx

## Roll process, using a temporary "sobj_postsoupx" Seurat object
sobj_postsoupx <- Seurat::CreateSeuratObject(counts = scmat_unsoup, project = "PBMC10K_AfterSoupX")
sobj_postsoupx <- Seurat::NormalizeData(object = sobj_postsoupx, verbose = FALSE)
sobj_postsoupx <- Seurat::ScaleData(object = sobj_postsoupx, verbose = FALSE)
sobj_postsoupx <- Seurat::FindVariableFeatures(object = sobj_postsoupx, verbose = FALSE)
sobj_postsoupx <- Seurat::RunPCA(object = sobj_postsoupx, npcs = 21, verbose = FALSE)
sobj_postsoupx <- Seurat::RunUMAP(object = sobj_postsoupx, dims = c(1:20), verbose = FALSE)

## Top1 feature plot
fp_post <- Seurat::FeaturePlot(object = sobj_postsoupx, features = 'LYZ') + Seurat::DarkTheme()
print(fp_post)
Show plot


Merging plots for ease of use :

# soupx_umaps

## Using the patchwork package to merge plots (and ggplot2 to add titles)
patchwork::wrap_plots(
  list(
    fp_pre & ggplot2::ggtitle(label = "LYZ before SoupX"), 
    fp_post & ggplot2::ggtitle(label = "LYZ after SoupX")), 
  nrow = 1)
Show plot


  • Question

    # q_umapsoupx

Can you compare the two plots :
  - for the LYZ expression across the 2-D space ?
  - about the space topologies ?
  - for the LYZ expression globally ?


    # a_umapsoupx

## . The level of expression measured in cells with smaller expression
##   before SoupX has gone to almost none after : the ambient expression
##   of this gene has efficiently been removed.
##
## . Hopefully, the cluster(s) with higher expression remain(s) high.
##
## . The expression scale upper bound after SoupX is higher than before ?! 
##   But we REMOVED some counts ?! While being counter-intuitive, this is
##   a positive consequence of the ambient removal on the scaling of
##   barcodes expression (see later). 
##
## . The range of both X and Y axes have increased, depicting a better
##   separation of / distance between some of the clusters.
##
## . The global topology of clusters remains close (but not identical).
## 
## . Clusters compacity has increased.


8.4 Conclusion

  • Removing ambient RNA has a positive effect :

    • on the quality of the expression level measurements

    • on the observed topology

    • … despite an estimation of only ~5% !

  • Actually, ambient RNA level and composition is (one of) the major source(s) of the batch effect bias that may alter the integration of different samples.





Beyond :

  1. One can check if the “unsouped” (ie, post-SoupX) matrix retain some soup ?. How would you do it ?

    # b_soupx2

## SoupX for ambient fraction estimation
soup2_frac <- SoupX_auto(
  scmat_filt = scmat_unsoup, 
  scmat_raw = scmat)

## Run SoupX in "removal" mode
scmat_unsoup2 <- SoupX_auto(
  scmat_filt = scmat_unsoup, 
  scmat_raw = scmat, 
  return_object = TRUE)

  2. Save the sobj_presoupx R object on disk.

    # b_save_answer1

## A single object : use RDS format with saveRDS()
saveRDS(object = sobj_presoupx, file = paste0(output_dir, "/sobj_presoupx.RDS"))

  3. Save both the sobj_presoupx and sobj_postsoupx objects in a single archive on disk.

    # b_save_answer2

## Multiple objects : use Rdata format with save()
save(list = c(sobj_presoupx, sobj_postsoupx), file = paste0(output_dir, "/sobjects.RDA"))

  4. Same as 1. but use the bzip2 compression algorithm

    # b_save_answer3

## Compress with bzip2
saveRDS(object = sobj, file = paste0(output_dir, "/sobj.RDS"), compress = "bzip2")





9 Rsession

For reproducibility and context, it is recommended to include in your RMarkdown the list of loaded packages and their version.

# rsession

utils::sessionInfo()
Show output 
R version 4.4.1 (2024-06-14)
Platform: x86_64-conda-linux-gnu
Running under: Ubuntu 22.04.5 LTS

Matrix products: default
BLAS/LAPACK: /shared/ifbstor1/software/miniconda/envs/r-4.4.1/lib/libopenblasp-r0.3.29.so;  LAPACK version 3.12.0

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

time zone: Europe/Paris
tzcode source: system (glibc)

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

loaded via a namespace (and not attached):
  [1] RcppAnnoy_0.0.22            splines_4.4.1              
  [3] later_1.4.2                 tibble_3.2.1               
  [5] R.oo_1.27.0                 polyclip_1.10-7            
  [7] fastDummies_1.7.5           lifecycle_1.0.4            
  [9] edgeR_4.2.2                 globals_0.18.0             
 [11] lattice_0.22-6              MASS_7.3-65                
 [13] magrittr_2.0.3              limma_3.60.6               
 [15] plotly_4.10.4               sass_0.4.10                
 [17] rmarkdown_2.29              jquerylib_0.1.4            
 [19] yaml_2.3.10                 metapod_1.12.0             
 [21] httpuv_1.6.15               Seurat_5.3.0               
 [23] sctransform_0.4.2           spam_2.11-1                
 [25] sp_2.2-0                    spatstat.sparse_3.1-0      
 [27] reticulate_1.42.0           cowplot_1.1.3              
 [29] pbapply_1.7-2               RColorBrewer_1.1-3         
 [31] abind_1.4-8                 zlibbioc_1.50.0            
 [33] Rtsne_0.17                  GenomicRanges_1.56.2       
 [35] purrr_1.0.4                 R.utils_2.13.0             
 [37] BiocGenerics_0.50.0         GenomeInfoDbData_1.2.12    
 [39] IRanges_2.38.1              S4Vectors_0.42.1           
 [41] ggrepel_0.9.6               irlba_2.3.5.1              
 [43] listenv_0.9.1               spatstat.utils_3.1-4       
 [45] goftest_1.2-3               RSpectra_0.16-2            
 [47] spatstat.random_3.4-1       dqrng_0.4.1                
 [49] fitdistrplus_1.2-2          parallelly_1.45.0          
 [51] DelayedMatrixStats_1.26.0   codetools_0.2-20           
 [53] DropletUtils_1.24.0         DelayedArray_0.30.1        
 [55] scuttle_1.14.0              tidyselect_1.2.1           
 [57] UCSC.utils_1.0.0            farver_2.1.2               
 [59] ScaledMatrix_1.12.0         SoupX_1.6.2                
 [61] matrixStats_1.5.0           stats4_4.4.1               
 [63] rmdformats_1.0.4            spatstat.explore_3.4-3     
 [65] jsonlite_2.0.0              BiocNeighbors_1.22.0       
 [67] progressr_0.15.1            ggridges_0.5.6             
 [69] survival_3.7-0              tools_4.4.1                
 [71] ica_1.0-3                   Rcpp_1.0.14                
 [73] glue_1.8.0                  gridExtra_2.3              
 [75] SparseArray_1.4.8           xfun_0.52                  
 [77] MatrixGenerics_1.16.0       GenomeInfoDb_1.40.1        
 [79] dplyr_1.1.4                 HDF5Array_1.32.0           
 [81] withr_3.0.2                 fastmap_1.2.0              
 [83] bluster_1.14.0              rhdf5filters_1.16.0        
 [85] rsvd_1.0.5                  digest_0.6.37              
 [87] R6_2.6.1                    mime_0.13                  
 [89] scattermore_1.2             tensor_1.5                 
 [91] dichromat_2.0-0.1           spatstat.data_3.1-6        
 [93] R.methodsS3_1.8.2           tidyr_1.3.1                
 [95] generics_0.1.4              data.table_1.17.4          
 [97] httr_1.4.7                  htmlwidgets_1.6.4          
 [99] S4Arrays_1.4.1              uwot_0.2.3                 
[101] pkgconfig_2.0.3             gtable_0.3.6               
[103] lmtest_0.9-40               SingleCellExperiment_1.26.0
[105] XVector_0.44.0              htmltools_0.5.8.1          
[107] dotCall64_1.2               bookdown_0.39              
[109] SeuratObject_5.1.0          scales_1.4.0               
[111] Biobase_2.64.0              png_0.1-8                  
[113] spatstat.univar_3.1-3       scran_1.32.0               
[115] knitr_1.50                  rstudioapi_0.17.1          
[117] reshape2_1.4.4              nlme_3.1-165               
[119] cachem_1.1.0                zoo_1.8-14                 
[121] rhdf5_2.48.0                stringr_1.5.1              
[123] KernSmooth_2.23-24          parallel_4.4.1             
[125] miniUI_0.1.2                pillar_1.10.2              
[127] grid_4.4.1                  vctrs_0.6.5                
[129] RANN_2.6.2                  promises_1.3.2             
[131] BiocSingular_1.20.0         beachmat_2.20.0            
[133] xtable_1.8-4                cluster_2.1.6              
[135] evaluate_1.0.3              cli_3.6.5                  
[137] locfit_1.5-9.9              compiler_4.4.1             
[139] rlang_1.1.6                 crayon_1.5.3               
[141] future.apply_1.11.3         labeling_0.4.3             
[143] plyr_1.8.9                  stringi_1.8.7              
[145] viridisLite_0.4.2           deldir_2.0-4               
[147] BiocParallel_1.38.0         lazyeval_0.2.2             
[149] spatstat.geom_3.4-1         Matrix_1.7-3               
[151] RcppHNSW_0.6.0              patchwork_1.3.0            
[153] sparseMatrixStats_1.16.0    future_1.49.0              
[155] ggplot2_3.5.2               Rhdf5lib_1.26.0            
[157] statmod_1.5.0               shiny_1.10.0               
[159] SummarizedExperiment_1.34.0 ROCR_1.0-11                
[161] igraph_2.1.4                bslib_0.9.0
LS0tCnRpdGxlOiAiPENFTlRFUj5FQjNJIG4xIDIwMjUgc2NSTkFzZXE8QlI+LTxCUj48Qj5QUkUtUFJPQ0VTU0lORyAoSSk8L0I+PEJSPi08QlI+TG9hZCBhIGNvdW50IG1hdHJpeCwgZW1wdHkgZHJvcGxldHMgJiBhbWJpZW50IFJOQSBmaWx0ZXJpbmc8L0NFTlRFUj4iCmRhdGU6ICIyMDI1LTE2LTIxLjIyIgphdXRob3I6CiAgLSBuYW1lOiAiRUIzSSBuMSBzY1JOQXNlcSBUZWFtIgogIC0gbmFtZTogIkJhc3RpZW4gSk9CIgogICAgZW1haWw6ICJiYXN0aWVuLmpvYkBndXN0YXZlcm91c3N5LmZyIgogIC0gbmFtZTogIldpbGxpYW0gSkFSQVNTSUVSIgogICAgZW1haWw6ICJ3LmphcmFzc2llckBnbWFpbC5jb20iIApvdXRwdXQ6CiAgcm1kZm9ybWF0czo6cmVhZHRoZWRvd246CiAgICBmaWdfd2lkdGg6IDgKICAgIGZpZ19oZWlnaHQ6IDYKICAgIGhpZ2hsaWdodDogdGFuZ28gICMjIFRoZW1lIGZvciB0aGUgY29kZSBjaHVua3MKICAgIGVtYmVkX2ZvbnRzOiBUUlVFCiAgICBudW1iZXJfc2VjdGlvbnM6IHRydWUgICMjIEFkZHMgbnVtYmVyIHRvIGhlYWRlcnMgKHNlY3Rpb25zKQogICAgdGhlbWU6IGZsYXRseSAgIyMgQ1NTIHRoZW1lIGZvciB0aGUgSFRNTCBwYWdlCiAgICBjb2xsYXBzZWQ6IHRydWUgICMjIEJ5IGRlZmF1bHQsIHRoZSBUT0MgaXMgZm9sZGVkCiAgICB0b2NfZGVwdGg6IDMKICAgIHNtb290aF9zY3JvbGw6IHRydWUgIyMgU21vb3RoIHNjcm9sbCBvZiB0aGUgSFRNTCBwYWdlCiAgICBzZWxmX2NvbnRhaW5lZDogdHJ1ZSAjIyBJbmNsdWRlcyBhbGwgcGxvdHMvaW1hZ2VzIHdpdGhpbiB0aGUgSFRNTAogICAgY29kZV9kb3dubG9hZDogdHJ1ZSAjIyBBZGRzIGEgYnV0dG9uIHRvIGRvd25sb2FkIHRoZSBSbWQKICAgIGNvZGVfZm9sZGluZzogc2hvdwogICAgdGh1bWJuYWlsczogZmFsc2UKICAgIGxpZ2h0Ym94OiB0cnVlCiAgICBmaWdfY2FwdGlvbjogZmFsc2UKICAgIGdhbGxlcnk6IHRydWUKICAgIHVzZV9ib29rZG93bjogdHJ1ZQphbHdheXNfYWxsb3dfaHRtbDogdHJ1ZSAjIyBBbGxvdyBwbGFpbiBIVE1MIGNvZGUgaW4gdGhlIFJtZAplZGl0b3Jfb3B0aW9uczogCiAgbWFya2Rvd246IAogICAgd3JhcDogNzIKLS0tCgo8IS0tIGtuaXQgc2V0dXAgLS0+CgpgYGB7ciBrbml0X3NldHVwLCBlY2hvID0gRkFMU0V9CmtuaXRyOjpvcHRzX2NodW5rJHNldCgKICBlY2hvID0gVFJVRSwgICAgICAgICMgUHJpbnQgdGhlIGNvZGUKICBldmFsID0gVFJVRSwgICAgICAgICMgUnVuIGNvbW1hbmQgbGluZXMKICBtZXNzYWdlID0gRkFMU0UsICAgICMgUHJpbnQgbWVzc2FnZXMKICBwcm9tcHQgPSBGQUxTRSwgICAgICMgRG8gbm90IGRpc3BsYXkgcHJvbXB0CiAgY29tbWVudCA9IE5BLCAgICAgICAjIE5vIGNvbW1lbnRzIG9uIHRoaXMgc2VjdGlvbgogIHdhcm5pbmcgPSBGQUxTRSwgICAgIyBEaXNwbGF5IHdhcm5pbmdzCiAgdGlkeSA9IEZBTFNFLAogIGZpZy5hbGlnbj0iY2VudGVyIiwgCiAgIyByZXN1bHRzID0gJ2hpZGUnLAogIHdpZHRoID0gMTAwICAgICAgICMgTnVtYmVyIG9mIGNoYXJhY3RlcnMgcGVyIGxpbmUKKQpgYGAKCjwhLS0gQ1NTIHRvIGNvbG9yIGNodW5rcyBhbmQgb3V0cHV0cyAtLT4KCmBgYHtjc3MsIGVjaG89RkFMU0V9Ci5ub3RydW4gewogIGJhY2tncm91bmQtY29sb3I6IGxpZ2h0Z3JleSAhaW1wb3J0YW50OwogIGJvcmRlcjogM3B4IHNvbGlkIGJsYWNrICFpbXBvcnRhbnQ7Cn0KLm5vdHJ1bm8gewogIGJhY2tncm91bmQtY29sb3I6IGxpZ2h0Z3JleSAhaW1wb3J0YW50OwogIGNvbG9yIDogYmxhY2sgIWltcG9ydGFudDsKfQoucXVlc3Rpb24gewogIGJhY2tncm91bmQtY29sb3I6IGFxdWFtYXJpbmUgIWltcG9ydGFudDsKICBjb2xvciA6IGJsYWNrICFpbXBvcnRhbnQ7CiAgYm9yZGVyOiAzcHggc29saWQgbGltZWdyZWVuICFpbXBvcnRhbnQ7Cn0KLnF1ZXN0aW9ubyB7CiAgYmFja2dyb3VuZC1jb2xvcjogYXF1YW1hcmluZSAhaW1wb3J0YW50OwogIGNvbG9yIDogYmxhY2sgIWltcG9ydGFudDsKfQouYW5zd2VyIHsKICBiYWNrZ3JvdW5kLWNvbG9yOiBuYXZham93aGl0ZSAhaW1wb3J0YW50OwogIGJvcmRlcjogM3B4IHNvbGlkIGJyb3duICFpbXBvcnRhbnQ7Cn0KLmFuc3dlcm8gewogIGJhY2tncm91bmQtY29sb3I6IG5hdmFqb3doaXRlICFpbXBvcnRhbnQ7CiAgY29sb3IgOiBibGFjayAhaW1wb3J0YW50Owp9Ci5iZXlvbmQgewogIGJhY2tncm91bmQtY29sb3I6IHZpb2xldCAhaW1wb3J0YW50OwogIGJvcmRlcjogM3B4IHNvbGlkIHB1cnBsZSAhaW1wb3J0YW50Owp9Ci5iZXlvbmRvIHsKICBiYWNrZ3JvdW5kLWNvbG9yOiB2aW9sZXQgIWltcG9ydGFudDsKICBjb2xvciA6IGJsYWNrICFpbXBvcnRhbnQ7Cn0KYGBgCgo8IS0tIEhvb2sgdG8gaGFuZGxlIGNvZGUgYmxvY2tzIG91dHB1dCBmb2xkaW5nIC0tPgoKYGBge3Iga25pdF9ob29rLCBlY2hvID0gRkFMU0V9Cmhvb2tzID0ga25pdHI6OmtuaXRfaG9va3MkZ2V0KCkKaG9va19mb2xkYWJsZSA9IGZ1bmN0aW9uKHR5cGUpIHsKICBmb3JjZSh0eXBlKQogIGZ1bmN0aW9uKHgsIG9wdGlvbnMpIHsKICAgIHJlcyA9IGhvb2tzW1t0eXBlXV0oeCwgb3B0aW9ucykKICAgIAogICAgaWYgKGlzRkFMU0Uob3B0aW9uc1tbcGFzdGUwKCJmb2xkLiIsIHR5cGUpXV0pKSByZXR1cm4ocmVzKQogICAgCiAgICBwYXN0ZTAoCiAgICAgICI8ZGV0YWlscz48c3VtbWFyeT5TaG93ICIsIHR5cGUsICI8L3N1bW1hcnk+XG5cbiIsCiAgICAgIHJlcywKICAgICAgIlxuXG48L2RldGFpbHM+IgogICAgKQogIH0KfQprbml0cjo6a25pdF9ob29rcyRzZXQoCiAgb3V0cHV0ID0gaG9va19mb2xkYWJsZSgib3V0cHV0IiksCiAgcGxvdCA9IGhvb2tfZm9sZGFibGUoInBsb3QiKQopCmBgYAoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCjxjZW50ZXI+IVtdKGViM2lfYmFubmVyLnBuZyk8L2NlbnRlcj4KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgojIFBSRUFNQkxFCgojIyBQdXJwb3NlIG9mIHRoaXMgc2Vzc2lvbgoKVGhpcyBmaWxlIGRlc2NyaWJlcyB0aGUgZGlmZmVyZW50IHN0ZXBzIHRvIHBlcmZvcm0gdGhlICoqZmlyc3QgcGFydCoqIG9mCnRoZSAqKnNpbmdsZSBjZWxsIFJOQXNlcSoqIGRhdGEgYW5hbHlzaXMgdHJhaW5pbmcgY291cnNlIGZvciB0aGUgKipFQjNJCm4xIDIwMjUqKiwgY292ZXJpbmcgdGhlc2Ugc3RlcHMgOgoKLSAgICoqTG9hZCoqIGEgc2luZ2xlIGNlbGwgcmF3IGNvdW50cyAqKm1hdHJpeCoqCgotICAgKipSZW1vdmUqKiBiYXJjb2RlcyBmcm9tICoqZW1wdHkgZHJvcGxldHMqKgoKLSAgICoqRXN0aW1hdGUqKiBhbmQgcmVtb3ZlICoqYW1iaWVudCBSTkEqKiBjb250YW1pbmF0aW9uICooInNvdXAiKSoKCiMjIEFuIGltcG9ydGFudCBub3RpY2UKCjxjZW50ZXI+IVtdKGd5cm8ucG5nKTwvY2VudGVyPgoKLSAgIFRoZXNlIGVhcmx5IChidXQgbWFuZGF0b3J5KSBzdGVwcyBvZiB0aGUgYW5hbHlzaXMgYXJlIG5vdCBjb3ZlcmVkIGJ5CiAgICB0aGUgcmVub3duZWQsIGhpZ2hlci1sZXZlbCBhbmFseXNpcyBmcmFtZXdvcmtzIGxpa2UgKipTZXVyYXQqKiBub3IKICAgICoqc2NyYW4vc2NhdGVyKiouCgotICAgVGh1cywgd2Ugd2lsbCB3cml0ZS91c2UgYSBjZXJ0YWluIGFtb3VudCBvZiAqKmxpbmVzIG9mIGNvZGUqKiwgc29tZQogICAgb2Ygd2hpY2ggb2YgYSBoaWdoZXIgY29tcGxleGl0eSBmb3IgUiBiZWdpbm5lcnMuCgotICAgKipOZXdjb21lciwgZG8gbm90IGZlYXIqKiAhCgogICAgLSAgIFdoaWxlIHlvdSBhcmUgKip3ZWxjb21lIHRvIHRyeSoqIHdyaXRpbmcgc29tZSBieSB5b3Vyc2VsZiwKICAgICAgICBkZXBlbmRpbmcgb24geW91ciBSIGtub3dsZWRnZSBhbmQgc2VsZi1jb25maWRlbmNlIGluIGl0IC4uLgoKICAgIC0gICAuLi4geW91IGRvIG5vdCBoYXZlIHRvIGtub3cgWyoqaG93Kipdey51bmRlcmxpbmV9IHRoZXNlIGNvZGUKICAgICAgICB3b3JrIC4uLgoKICAgIC0gICAuLi4gdGh1cyB5b3UgKipzaG91bGQgbm90IGZlZWwgYXNoYW1lZCoqIGNvcHktcGFzdGluZyB0aGUKICAgICAgICBwcmUtd3JpdHRlbiBjb2RlIGhlcmVhZnRlciAuLi4KCiAgICAtICAgLi4uIGJ1dCBpbiBib3RoIGNhc2VzLCBwbGVhc2UgdHJ5ICoqeW91ciBiZXN0IHVuZGVyc3RhbmRpbmcqKgogICAgICAgIFsqKndoYXQqKl17LnVuZGVybGluZX0gdGhlc2UgY29kZSBkbywgYW5kIGV2ZW4gbW9yZSwKICAgICAgICBbKip3aHkqKl17LnVuZGVybGluZX0gIQoKICAgIC0gICBJbiB0aGlzIHB1cnBvc2UsIHlvdSBtYXkgdXNlIHdoYXRldmVyICoqaGVscCoqIHlvdSBjYW4gKipieSoqCiAgICAgICAgKipzZWFyY2hpbmcqKi4KCiMjIFdoYXQgdG8gcnVuID8KCi0gICBUaGlzIHRyYWluaW5nIHByZXNlbnRhdGlvbiwgd3JpdHRlbiBpbiBSbWFya2Rvd24sIGNvbnRhaW5zICoqbWFueSoqCiAgICBjaHVua3MgKCA9IGNvZGUgYmxvY2tzKQoKLSAgIFlvdSAqKmRvIG5vdCBoYXZlIHRvKiogcnVuIHRoZW0gYWxsICEKCi0gICBDaHVua3MgZm9sbG93IGEgKipjb2xvciBjb2RlKiogOgoKIyMjIENodW5rcyB5b3UgKipBUkUqKiBleHBlY3RlZCB0byBydW4gOgoKYGBge3IgY2h1bmtfcnVufQojIGNodW5rX3J1bgoKIyMgRXhhbXBsZSBvZiBjaHVuayB0byBydW4KZ2V0d2QoKQoKYGBgCgo8YnI+CgojIyMgQ2h1bmtzIHlvdSAqKkFSRSBOT1QqKiBleHBlY3RlZCB0byBydW4KCkhvcGVmdWxseSwgd2Ugd2lsbCBkZW1vIHRoZW0uIFlvdSBtYXkgcnVuIHRoZW0gYnkgeW91cnNlbGYsIHlvdW5nCnBhZGF3YW4sIGJ1dCBpZiB5b3UgZ2V0IGxhdGUgb3Igd29yc3QsIGJyZWFrIHNvbWV0aGluZywgKnlvdSonbGwgYmUgdGhlCm9uZSB0byBibGFtZSAhCgpgYGB7ciBjaHVua19kZW1vLCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiMgY2h1bmtfZGVtbwoKIyMgRXhhbXBsZSBvZiBjaHVuayB0byBOT1QgcnVuCnV0aWxzOjpzZXNzaW9uSW5mbygpCgpgYGAKCjxicj4KCiMjIyBRdWVzdGlvbnMKClNvbWUgcXVlc3Rpb25zIHdlIHdvdWxkIGxpa2UgeW91IHRvIGFuc3dlciBkdXJpbmcgdGhlIHRyYWluaW5nIGNvdXJzZS4KUGxlYXNlIGJlIGZhaXIgYW5kIGRvIG5vdCB0YWtlIGEgbG9vayBhdCB0aGUgYW5zd2VyIGJlZm9yZSB3ZSB0ZWxsIHlvdQpzbyA6KQoKYGBge3IgY2h1bmtfcXVlc3Rpb24sIGNsYXNzLnNvdXJjZT0icXVlc3Rpb24iLCBldmFsID0gRkFMU0V9CiMgY2h1bmtfcXVlc3Rpb24KCiMjIFdoYXQncyB0aGUgQW5zd2VyIHRvIExpZmUgPyBUaGUgVW5pdmVyc2UgPyBFdmVyeXRoaW5nID8/CgpgYGAKCjxicj4KCiMjIyBBbnN3ZXJzCgpUaGV5IHVzdWFsbHkgZm9sbG93IHF1ZXN0aW9ucy4uLgoKYGBge3IgY2h1bmtfYW5zd2VyLCBjbGFzcy5zb3VyY2UgPSBjKCJmb2xkLWhpZGUiLCAiYW5zd2VyIiksIGNsYXNzLm91dHB1dD0iYW5zd2VybyJ9CiMgY2h1bmtfYW5zd2VyCgpjYXQoIlRoZSBBbnN3ZXIgdG8gdGhlIEdyZWF0IFF1ZXN0aW9uLi4uIGlzIC4uLiAKRm9ydHktdHdvLiBTaXggYnkgbmluZSA6IGZvcnR5LXR3by4gVGhhdCdzIGl0LiAKVGhhdCdzIGFsbCB0aGVyZSBpcy4gKEkgYWx3YXlzIHRob3VnaHQgc29tZXRoaW5nIHdhcyAKZnVuZGFtZW50YWxseSB3cm9uZyB3aXRoIHRoZSB1bml2ZXJzZS4pIikKCmBgYAoKPGJyPgoKIyMjIEFkZGl0aW9uYWwgZXhlcmNpc2VzIC8gcXVlc3Rpb25zCgpUaGVzZSBhcmUgYmV5b25kIHRoZSBzY29wZSBvZiB0aGUgdHJhaW5pbmcgY291cnNlLCBidXQgeW91IG1heSBwbGF5IGl0CippZiB5b3UgZmVlbCBib3JlZCoKCmBgYHtyIGNodW5rX2JleW9uZCwgY2xhc3Muc291cmNlPSJiZXlvbmQiLCBjbGFzcy5vdXRwdXQ9ImJleW9uZG8ifQojIGNodW5rX2JleW9uZAoKIyMgRXhhbXBsZSBvZiBjaHVuayB0byBwbGF5IHdpdGggYnV0IGJleW9uZCB0aGUgc2NvcGUKZ2V0d2QoKQoKYGBgCgo8YnI+CgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyBTdGFydCBSc3R1ZGlvCgotICAgVXNpbmcgdGhlIFtPcGVuT25EZW1hbmQvUnN0dWRpbyBjaGVhdAogICAgc2hlZXRdKGh0dHBzOi8vbW9vZGxlLmZyYW5jZS1iaW9pbmZvcm1hdGlxdWUuZnIvcGx1Z2luZmlsZS5waHAvMTQ3NS9tb2RfZm9sZGVyL2NvbnRlbnQvMC9Pb0RfUl9Sc3R1ZGlvLmh0bWwpe3RhcmdldD0iX2JsYW5rIn0sCiAgICBjb25uZWN0IHRvIHRoZSBbT3Blbk9uRGVtYW5kCiAgICBwb3J0YWxdKGh0dHBzOi8vb25kZW1hbmQuY2x1c3Rlci5mcmFuY2UtYmlvaW5mb3JtYXRpcXVlLmZyKXt0YXJnZXQ9Il9ibGFuayJ9CiAgICBhbmQgKipjcmVhdGUgYSBSc3R1ZGlvIHNlc3Npb24qKiB3aXRoIHRoZSByaWdodCByZXNvdXJjZQogICAgcmVxdWlyZW1lbnRzLgoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCiMgV2FybS11cAoKLSAgIFdlIG5vdyBzZXQgKipjb21tb24gcGFyYW1ldGVycyoqIGFzIG5ldyB2YXJpYWJsZXMsIG9uY2UgYW5kIGZvciBhbGwKICAgIGZvciB0aGlzIHNlc3Npb24gOgoKYGBge3Igc2V0cGFyYW19CiMgc2V0cGFyYW0KCgojIyBTZXQgeW91ciBwcm9qZWN0IG5hbWUKIyBXQVJOSU5HIDogRG8gbm90IGp1c3QgY29weS1wYXN0ZSB0aGlzICEgSXQncyBNWSBwcm9qZWN0IG5hbWUgISBQdXQgWU9VUlMgISEKcHJvamVjdF9uYW1lIDwtICJlYmFpaV9zY190ZWFjaGVycyIKCgojIyBDb250cm9sIGlmIHRoZSBwcm9qZWN0X25hbWUgZXhpc3RzIG9uIHRoZSBjbHVzdGVyCmNhdCgnUEFUSCBDSEVDSyA6ICcsIGRpci5leGlzdHMocGFzdGUwKCcvc2hhcmVkL3Byb2plY3RzLycsIHByb2plY3RfbmFtZSkpKQoKIyMgU2VlZCBmb3IgdGhlIFJORwpteV9zZWVkIDwtIDEzMzdMCgojIyBFbXB0eSBkcm9wbGV0cyBtYXggcC12YWx1ZQptYXhfcCA8LSAxRS0wMwoKYGBgCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyBQcmVwYXJlIHRoZSBkYXRhIHN0cnVjdHVyZQoKLSAgIER1cmluZyB0aGUgY3VycmVudCBhbmQgcmVtYWluaW5nIHRyYWluaW5nIHNlc3Npb25zLCB3ZSB3aWxsIHdvcmsgb24KICAgIGRpZmZlcmVudCAqKmRhdGEqKiAqKm9iamVjdHMqKiwgc29tZSByZXRyaWV2ZWQgZnJvbSAqKmV4dGVybmFsKioKICAgIHJlc291cmNlcywgc29tZSB3ZSB3aWxsICoqY3JlYXRlKiouCgotICAgSW4gb3JkZXIgdG8gcmV1c2UgdGhlc2UgZm9yIGZ1cnRoZXIgc3RlcHMsIG9yIG90aGVyIHNlc3Npb25zLCB3ZQogICAgbmVlZCB0byAqKnN0b3JlKiogdGhlbSwgdGh1cyBjcmVhdGUgYSAqKmRpcmVjdG9yeSBzdHJ1Y3R1cmUqKiBvbiB0aGUKICAgIGNsdXN0ZXIuCgotICAgV2Ugd2lsbCBjcmVhdGUgdGhpcyBzdHJ1Y3R1cmUgaW4gKip5b3VyIHByb2plY3QqKiBmb2xkZXIgKHRoZSBvbmUKICAgIHlvdSByZWdpc3RlcmVkIGF0IHRoZSBJRkIgQ29yZSwgaW4gd2hpY2ggeW91IHN0b3JlZCB5b3VyIG93biBkYXRhCiAgICBmb3IgdGhlIHR1dG9yaW5nKS4KCi0gICBUbyBlbnN1cmUgdGhhdCAqeW91KiB3aWxsIGFsd2F5cyB1c2UgdGhlIHNhbWUgc3RydWN0dXJlIGZvciAqeW91cioKICAgIHdvcmssIHRoZSBjb2RlIGhlcmVhZnRlciBoYXMgYmVlbiB3cml0ZW4gc28gdGhhdCB5b3UganVzdCBoYXZlIHRvCiAgICBkZWZpbmUgYSBzaW5nbGUgKip2YXJpYWJsZSoqIHRoYXQgY29ycmVzcG9uZHMgdG8geW91ciBwcm9qZWN0IG5hbWUuCgogICAgLSAgIEluICpNWSogY2FzZSwgYXMgYSBTQyB0cmFpbmluZyB0ZWFjaGVyLCBpdCBpcwogICAgICAgICoqYGViYWlpX3NjX3RlYWNoZXJzYCoqLgoKIyMgTWFpbiBkaXJlY3RvcnkKCi0gICBDcmVhdGUgdGhlIHRyYWluaW5nIGNvdXJzZSBkaXJlY3RvcnkKCmBgYHtyIG1haW5kaXIsIGZvbGQub3V0cHV0ID0gRkFMU0V9CiMgbWFpbmRpcgoKIyMgUHJlcGFyZSB0aGUgcGF0aCBpbiBhICJURF9kaXIiIHZhcmlhYmxlClREX2RpciA8LSBwYXN0ZTAoIi9zaGFyZWQvcHJvamVjdHMvIiwgcHJvamVjdF9uYW1lLCAiL1NDX1REIikKCiMjIENyZWF0ZSB0aGUgcm9vdCBkaXJlY3RvcnkKZGlyLmNyZWF0ZShwYXRoID0gVERfZGlyLCByZWN1cnNpdmUgPSBUUlVFKQoKIyMgUHJpbnQgdGhlIHJvb3QgZGlyZWN0b3J5IG9uLXNjcmVlbgpwcmludChURF9kaXIpCgpgYGAKCiMjIEN1cnJlbnQgc2Vzc2lvbgoKYGBge3Igc2Vzc2lvbmRpciwgZm9sZC5vdXRwdXQgPSBGQUxTRX0KIyBzZXNzaW9uZGlyCgojIyBDcmVhdGUgdGhlIHNlc3Npb24gKFByZXByb2MuMSkgZGlyZWN0b3J5CnNlc3Npb25fZGlyIDwtIHBhc3RlMChURF9kaXIsICIvMDFfUHJlcHJvYy4xIikKZGlyLmNyZWF0ZShwYXRoID0gc2Vzc2lvbl9kaXIsIHJlY3Vyc2l2ZSA9IFRSVUUpCgojIyBQcmludCB0aGUgc2Vzc2lvbiBkaXJlY3Rvcnkgb24tc2NyZWVuCnByaW50KHNlc3Npb25fZGlyKQpgYGAKCiMjIElucHV0IGRhdGEgZGlyZWN0b3J5CgpgYGB7ciBpbmRpciwgZm9sZC5vdXRwdXQ9RkFMU0V9CiMgaW5kaXIKCiMjIENyZWF0ZSB0aGUgSU5QVVQgZGF0YSBkaXJlY3RvcnkKaW5wdXRfZGlyIDwtIHBhc3RlMChzZXNzaW9uX2RpciwgIi9EQVRBIikKZGlyLmNyZWF0ZShwYXRoID0gaW5wdXRfZGlyLCByZWN1cnNpdmUgPSBUUlVFKQoKIyMgUHJpbnQgdGhlIGlucHV0IGRpcmVjdG9yeSBvbi1zY3JlZW4KcHJpbnQoaW5wdXRfZGlyKQpgYGAKCiMjIE91dHB1dCByZXN1bHRzIGRpcmVjdG9yeQoKYGBge3Igb3V0ZGlyLCBmb2xkLm91dHB1dD1GQUxTRX0KIyBvdXRkaXIKCiMjIENyZWF0ZSB0aGUgT1VUUFVUIGRhdGEgZGlyZWN0b3J5Cm91dHB1dF9kaXIgPC0gcGFzdGUwKHNlc3Npb25fZGlyLCAiL1JFU1VMVFMiKQpkaXIuY3JlYXRlKHBhdGggPSBvdXRwdXRfZGlyLCByZWN1cnNpdmUgPSBUUlVFKQoKIyMgUHJpbnQgdGhlIG91dHB1dCBkaXJlY3Rvcnkgb24tc2NyZWVuCnByaW50KG91dHB1dF9kaXIpCmBgYAoKPGJyPjxicj4KClsqKk5vdGUgOioqXXsudW5kZXJsaW5lfSAqQWN0dWFsbHksIGZvciB0aGlzIHNlc3Npb24gd2Ugd29uJ3QgZXhwb3J0IGFueQpyZXN1bHQsIHNvIGNyZWF0aW5nIHRoaXMgKipgUkVTVUxUU2AqKiBmb2xkZXIgd2FzIGp1c3QgZm9yIHByYWN0aWNlCnB1cnBvc2UuKgoKPGJyPgoKTm93LCB0aGUgd2UgaGF2ZSBldmVyeXRoaW5nIG5lZWRlZCB0byA6CgotICAgKipsb2FkKiogZGF0YSAoYGlucHV0X2RpcmApCgotICAgKipzYXZlKiogdGhlIHJlc3VsdHMgd2Ugd2lsbCBnZW5lcmF0ZSAoYG91dHB1dF9kaXJgKS4KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgojIEFkZGl0aW9uYWwgcmVzb3VyY2VzCgoqKkhvd2V2ZXIqKiwgc29tZSB1cGNvbWluZyBzdGVwcyByZXF1aXJlIGEgYml0IGNvbXBsaWNhdGVkIG1ldGhvZHMgdGhhdAphcmUgbm90IGZ1biB0byB3cml0ZSwgbm9yIGNvcHktcGFzdGUgbXVsdGlwbGUgdGltZXMuIFNvLCB3ZSB3aWxsIG5vdwpwcmVwYXJlIGNvZGUgKipmdW5jdGlvbnMqKiB0aGF0IHdpbGwgZWFzZSB0aGVzZSBzdGVwcy4gQWZ0ZXIgZXhlY3V0aW5nCnRoZXNlIGZ1bmN0aW9ucyBjb2RlLCB3ZSBtYXkgY2FsbCAoaW52b2tlKSB0aGVtIHdoZW4gbmVlZGVkLgoKIyMgRnVuY3Rpb25zCgojIyMgU291cFggaGVscGVyIGZ1bmN0aW9uCgpUaGlzIGZ1bmN0aW9uIHdhcyB3cml0dGVuIHRvIGVhc2UgdGhlIHVzZSBvZiBTb3VwWCBvbiBhIGNvdW50IG1hdHJpeApJdCB3YXMgdGVzdGVkIG9uIGBTb3VwWGA9MS42LjIgd2l0aCBgaWdyYXBoYD0xLjUuMSwgYFNldXJhdGA9NC40LjAsNS4xLjAKCioqUGFyYW1ldGVycyoqIDoKCi0gICBgc2NtYXRfZmlsdGAgOiAoc3BhcnNlKW1hdHJpeCBjb3JyZXNwb25kaW5nIHRvIGFuIGVtcHR5IGRyb3BsZXRzIC1maWx0ZXJlZCBjb3VudCBtYXRyaXgKCi0gICBgc2NtYXRfcmF3YCA6IChzcGFyc2UpbWF0cml4IGNvcnJlc3BvbmRpbmcgdG8gYW4gTk9OLSBlbXB0eSBkcm9wbGV0cyAtZmlsdGVyZWQgY291bnQgbWF0cml4CgotICAgYHNvdXBRdWFudGlsZWAsIGBjb250YW1pbmF0aW9uUmFuZ2VgIDogc2VlIGA/U291cFg6OmF1dG9Fc3RDb250YAoKLSAgIGBjb250YW1pbmF0aW9uUmFuZ2VgIDogcmFuZ2Ugb2YgdGhlICpleHBlY3RlZCogc291cCBwcm9wb3J0aW9uCgotICAgYHNvdXBSYW5nZWAgOiBlc3RpbWF0ZSBzb3VwIGZyYWN0aW9uIGZyb20gZmVhdHVyZXMgd2l0aCB0b3RhbCByZWFkcwpjb21wcmlzZWQgaW4gdGhpcyByYW5nZSBvZiBjb3VudHMgKG9ubHkgdXNlZCB3aGVuIGBzY21hdF9yYXcgIT0gTlVMTGApCgotICAgYHJldHVybl9vYmplY3RgIDogaWYgYFRSVUVgLCByZXR1cm4gdGhlICJ1bnNvdXBlZCIgY291bnQgbWF0cml4IDsgaWYKYEZBTFNFYCwgb25seSBwZXJmb3JtIHNvdXAgZXN0aW1hdGlvbiBhbmQgcmV0dXJuIHRoZSBlc3RpbWF0ZWQKcHJvcG9ydGlvbiB2YWx1ZSAocmhvKQoKLSAgIGBkb1Bsb3RgIDogUGVyZm9ybSB0aGUgU291cFggZXN0aW1hdGlvbiBwbG90IChyaG8gZGlzdHJpYnV0aW9uKQoKYGBge3Igc291cHhfZnVuY30KIyBzb3VweF9mdW5jCgpTb3VwWF9hdXRvIDwtIGZ1bmN0aW9uKHNjbWF0X2ZpbHQgPSBOVUxMLCBzY21hdF9yYXcgPSBOVUxMLCBzb3VwUXVhbnRpbGUgPSAwLjksCiAgICAgICAgICAgICAgICAgICAgICAgY29udGFtaW5hdGlvblJhbmdlID0gYyguMDEsIC44KSwgc291cFJhbmdlID0gYygwLDEwMCksCiAgICAgICAgICAgICAgICAgICAgICAgcmV0dXJuX29iamVjdCA9IEZBTFNFLCBkb1Bsb3QgPSBGQUxTRSkgewogIAogICMjIENoZWNrcwogIGlmKGlzLm51bGwoc2NtYXRfZmlsdCkpIHN0b3AoJ0EgZmlsdGVyZWQgY291bnQgbWF0cml4IGlzIHJlcXVpcmVkICEnKQogIAogIGlmKGlzLm51bGwoc2NtYXRfcmF3KSkgbWVzc2FnZSgnTm8gdW5maWx0ZXJlZCByYXcgY291bnRzIG1hdHJpeCBwcm92aWRlZC4gRXN0aW1hdGlvbiB3aWxsIGJlIGJhc2VkIG9uIGZpbHRlcmVkIG1hdHJpeCBvbmx5LicpCiAgCiAgIyMgSWYgbm8gcmF3IG1hdHJpeAogIGlmIChpcy5udWxsKHNjbWF0X3JhdykpIHsKICAgIHNwQ2hhblJhdyA8LSBTb3VwWDo6U291cENoYW5uZWwoCiAgICAgIHRvZCA9IHNjbWF0X2ZpbHQsIAogICAgICB0b2MgPSBzY21hdF9maWx0LCAKICAgICAgY2FsY1NvdXBQcm9maWxlID0gRkFMU0UpCiAgICBzY19yb3dzdW0gPC0gc3BhcnNlTWF0cml4U3RhdHM6OnJvd1N1bXMyKHNjbWF0X2ZpbHQpCiAgICBzcFByb2YgPC0gZGF0YS5mcmFtZSgKICAgICAgcm93Lm5hbWVzID0gcm93bmFtZXMoc2NtYXRfZmlsdCksIAogICAgICBlc3QgPSBzY19yb3dzdW0vc3VtKHNjbWF0X2ZpbHQpLCAKICAgICAgY291bnRzID0gc2Nfcm93c3VtKQogICAgc3BDaGFuIDwtIFNvdXBYOjpzZXRTb3VwUHJvZmlsZShzcENoYW5SYXcsIHNwUHJvZikKICB9IGVsc2UgewogICAgc3BDaGFuIDwtIFNvdXBYOjpTb3VwQ2hhbm5lbCgKICAgICAgdG9kID0gc2NtYXRfcmF3LCAKICAgICAgdG9jID0gc2NtYXRfZmlsdCwgCiAgICAgIGNhbGNTb3VwUHJvZmlsZSA9IEZBTFNFKQogICAgaWYgKG1pbihzcENoYW4kbkRyb3BVTUlzKSA+IG1heChzb3VwUmFuZ2UpKSBzdG9wKAogICAgICAnTWluaW11bSBmb3VuZCBjb3VudHMgcGVyIGJhcmNvZGUgaXMgOiAnLCAKICAgICAgbWluKHNwQ2hhbiRuRHJvcFVNSXMpLCAKICAgICAgJywgd2hpY2ggaXMgc21hbGxlciB0aGFuIHRoZSB1cHBlciBib3VuZCBvZiBzb3VwUmFuZ2UgISBQbGVhc2UgaW5jcmVhc2Ugc291cFJhbmdlIG1heCAhJykKICAgIHNwQ2hhbiA8LSBTb3VwWDo6ZXN0aW1hdGVTb3VwKHNjID0gc3BDaGFuLCBzb3VwUmFuZ2UgPSBzb3VwUmFuZ2UpCiAgfQogICMjIERpc3BsYXkgVG9wIDIwIGNvbnRyaWJ1dGluZyBnZW5lcwogIGlmICghcmV0dXJuX29iamVjdCkgewogICAgY2F0KCdcblNvdXAtY29udHJpYnV0aW5nIGZlYXR1cmVzIChUb3AgMjApIDpcbicpCiAgICBwcmludChrbml0cjo6a2FibGUoaGVhZCgKICAgICAgc3BDaGFuJHNvdXBQcm9maWxlW29yZGVyKHNwQ2hhbiRzb3VwUHJvZmlsZSRlc3QsIGRlY3JlYXNpbmcgPSBUUlVFKSwgXSwgCiAgICAgIG4gPSAyMCkpKQogIH0KICAjIyBRdWljayBjbHVzdGVyaW5nIG5lZWRlZAogIHNwQ2x1c3QgPC0gc2NyYW46OnF1aWNrQ2x1c3RlcihzY21hdF9maWx0LCBtZXRob2QgPSAiaWdyYXBoIikKICAjIyBBZGRpbmcgY2x1c3RlcnMgdG8gdGhlIFNvdXBDaGFubmVsIG9iamVjdAogIHNwQ2hhbiA8LSBTb3VwWDo6c2V0Q2x1c3RlcnMoc2MgPSBzcENoYW4sIGNsdXN0ZXJzID0gc3BDbHVzdCkKICAjIyBFc3RpbWF0aW5nIHNvdXAKICBzWCA8LSBTb3VwWDo6YXV0b0VzdENvbnQoc2MgPSBzcENoYW4sIGRvUGxvdCA9IGRvUGxvdCwgdGZpZGZNaW4gPSAxLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgc291cFF1YW50aWxlID0gc291cFF1YW50aWxlLCBtYXhNYXJrZXJzID0gMTAwLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgY29udGFtaW5hdGlvblJhbmdlID0gY29udGFtaW5hdGlvblJhbmdlLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgcmhvTWF4RkRSID0gLjIsIHByaW9yUmhvID0gLjA1LCBwcmlvclJob1N0ZERldiA9IC4xLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgZm9yY2VBY2NlcHQgPSBGQUxTRSkKICAKICAjIyBSZW1vdmluZyBzb3VwIChhZGp1c3RpbmcgY291bnRzKQogIGlmKHJldHVybl9vYmplY3QpIHsKICAgIGNhdCgnQ291bnRzIEJFRk9SRSBTb3VwWCA6ICcsIHN1bShzY21hdF9maWx0KSwgJ1xuJykKICAgIHNjbWF0X3NvdXB4IDwtIFNvdXBYOjphZGp1c3RDb3VudHMoCiAgICAgIHNYLCBtZXRob2QgPSAnc3VidHJhY3Rpb24nLCByb3VuZFRvSW50ID0gVFJVRSwgCiAgICAgIHRvbCA9IC4wMDEsIHBDdXQgPSAuMDEpCiAgICBjYXQoJ0NvdW50cyBBRlRFUiBTb3VwWCA6ICcsIHN1bShzY21hdF9zb3VweCksICdcbicpCiAgICBybShzY21hdF9maWx0KQogICAgcmV0dXJuKHNjbWF0X3NvdXB4KQogIH0gZWxzZSByZXR1cm4oc1gkZml0JHJob0VzdCkKfQoKYGBgCgo8YnI+CgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyBMb2FkIGEgMTBYIENlbGwgUmFuZ2VyIHJhdyBjb3VudCBtYXRyaXgKCiMjIFRoZSBkYXRhc2V0CgotICAgVGhlIGRhdGEgd2Ugd2lsbCB3b3JrIG9uIGlzIGhvc3RlZCBvbgogICAgW1plbm9kb10oaHR0cHM6Ly96ZW5vZG8ub3JnKXt0YXJnZXQ9Il9ibGFuayJ9LiBUbyBlYXNlIGl0cwogICAgcmV0cmlldmFsLCB3ZSBjYW4gZGVmaW5lIGl0cyBJRCBhcyBhIHZhcmlhYmxlCgogICAgYGBge3IgemlkfQogICAgIyB6aWQKICAgIAogICAgIyMgVGhpcyBpcyBhY3R1YWxseSB0byBoZWxwIG1lIHVwZGF0aW5nIHRoaXMgdHJhaW5pbmcgd2l0aG91dCBoYXNzbGUKICAgIHplbl9pZCA8LSAnMTQwMzQyMjEnCiAgICAKICAgICMjIFRoaXMgaXMgdGhlIHBhdGggdG8gdGhlIGN1cnJlbnQgRUIzSSBiYWNrdXAKICAgIHNlc3Npb25pZCA8LSAnMjUzOF9lYjNpX24xXzIwMjUnCiAgICAKICAgIGBgYAoKLSAgIEhlcmUgd2lsbCB3ZSB0cmFpbiBvdXJzZWx2ZXMgdG8gbG9hZCBpbnRvIFIgYSBzaW5nbGUgY2VsbCBSTkFzZXEKICAgIGRhdGEgcHJvZHVjZWQgYnkgMTBYIEdlbm9taWNzJyBzb2Z0d2FyZSAqKmBDZWxsIFJhbmdlcmAqKi4KCi0gICBXZSB3aWxsIHdvcmsgd2l0aCBhIHB1YmxpYyBkYXRhc2V0IHByb3ZpZGVkIGJ5IHRoZSBtYW51ZmFjdHVyZXIsCiAgICB0aGF0IGNvbnNpc3RzIGluIFx+ICoqMTAsMDAwIFBCTUMqKiAocGVyaXBoZXJhbCBib25lIG1hcnJvdyBjZWxscykKICAgIGZyb20gYSBodW1hbiBkb25vciAoZG93bnNhbXBsZWQgdG8gMS8xMHRoLCBmb3Igc21hbGxlciBjb21wdXRhdGlvbiB0aW1lKS4KCi0gICBUaGUgZXhwZXJpbWVudCB3YXMgcGVyZm9ybWVkIHdpdGggdGhlICoqMycgY2FwdHVyZSBraXQgdjMqKi4KCi0gICBUaGUgYW5hbHlzaXMgd2FzIHBlcmZvcm1lZCB3aXRoICoqYENlbGwgUmFuZ2VyIHYzYCoqLCBtYXBwaW5nIG9uIHRoZQogICAgKipgR1JDaDM4LTIwMjAtQWAqKiBtYW51ZmFjdHVyZXIgcmVmZXJlbmNlLgoKLSAgIFRoZSBkYXRhIGNvbnNpc3RzIGludG8gdGhlIHR5cGljYWwgMTB4IDMtZmlsZXMgc3RydWN0dXJlLCBob3N0ZWQgaW4KICAgIGEgWmVub2RvIHJlc3Bvc2l0b3J5IChJZCA6CiAgICBbYHIgemVuX2lkYF0oaHR0cHM6Ly96ZW5vZG8ub3JnL3JlY29yZHMvJTYwciUyMHplbl9pZCU2MCAiWmVub2RvIFByZXByb2MuMSIpe3RhcmdldD0iX2JsYW5rIn0pCgogICAgLSAgIGBiYXJjb2Rlcy50c3YuZ3pgIDogdGhlIGxpc3Qgb2YgcHV0YXRpdmUgZHJvcGxldCAqKmJhcmNvZGVzKioKICAgIC0gICBgZmVhdHVyZXMudHN2Lmd6YCA6IHRoZSBsaXN0IG9mICoqZmVhdHVyZXMqKiAoZ2VuZXMpCiAgICAtICAgYG1hdHJpeC5tdHguZ3pgIDogdGhlICoqZmVhdHVyZSoqIHggKipiYXJjb2RlKiogZXhwcmVzc2lvbiBjb3VudAogICAgICAgICoqbWF0cml4KioKCiMjIERvd25sb2FkIGRhdGEKCi0gICBUaGUgcmVsYXRpdmVseSBjb21wbGV4IGNvZGUgY2h1bmsgaGVyZWFmdGVyIGhhcyBiZWVuIHdyaXR0ZW4gdG8gOgoKICAgIC0gICBsZXQgeW91IGRvd25sb2FkIHRoZSBkYXRhIGZpbGVzIGZyb20gOgoKICAgICAgICAtICAgWmVub2RvIGJ5IGRlZmF1bHQKCiAgICAgICAgLSAgIGlmIGBiYWNrdXAgPSBUUlVFYCwgYSBsb2NhbCBiYWNrdXAgaW5zdGVhZCAoaW4gY2FzZSB3ZSBsb29zZSB0aGUKICAgICAgICAgICAgaW50ZXJuZXQgY29ubmVjdGlvbikKCiAgICAtICAgaWYgdGhlIHJlcXVlc3RlZCBkYXRhIGFyZSBhbHJlYWR5IGF2YWlsYWJsZSBsb2NhbGx5IDoKCiAgICAgICAgLSAgIGxvYWQgdGhpcyBmcm9tIHRoZSBsb2NhbCBzb3VyY2UKCiAgICAgICAgLSAgIGlmIGBmb3JjZSA9IFRSVUVgIGEgbmV3IGRvd25sb2FkIGlzIGZvcmNlZAoKYGBge3IgZGx6ZW4xMHgsIG1lc3NhZ2UgPSBUUlVFLCBlY2hvID0gVFJVRX0KIyBkbHplbjEweAoKIyMjIE5hbWVkIGZpbGVzICh3aWxsIGJlIHVzZWQgbGF0ZXIgb24gISkKbXR4X2ZpbGUgPC0gIm1hdHJpeC5tdHguZ3oiCmZlYXR1cmVzX2ZpbGUgPC0gImZlYXR1cmVzLnRzdi5neiIKYmFyY29kZXNfZmlsZSA8LSAiYmFyY29kZXMudHN2Lmd6IgpzdW1tYXJ5X2ZpbGUgPC0gInBibWNfMTBrX3YzX3N1bW1hcnkuaHRtbCIKCiMjIEZpbGVuYW1lKHMpIHRvIHJldHJpZXZlCnRvZ2V0X2ZpbGVzIDwtIGMobXR4X2ZpbGUsCiAgICAgICAgICAgICAgICAgZmVhdHVyZXNfZmlsZSwKICAgICAgICAgICAgICAgICBiYXJjb2Rlc19maWxlLAogICAgICAgICAgICAgICAgIHN1bW1hcnlfZmlsZSkKCiMjIEZvbGRlciB0byBzdG9yZSByZXRyaWV2ZWQgZmlsZXMKbG9jYWxfZm9sZGVyIDwtIGlucHV0X2RpcgoKIyMgVXNlIGxvY2FsIGJhY2t1cCA/CmJhY2t1cCA8LSBGQUxTRQppZihiYWNrdXApIG1lc3NhZ2UoIlVzaW5nIGxvY2FsIGJhY2t1cCAhIikKCiMjIEZvcmNlIGRvd25sb2FkID8KZm9yY2UgPC0gRkFMU0UKaWYoZm9yY2UpIG1lc3NhZ2UoIkZvcmNpbmcgKHJlKWRvd25sb2FkICEiKQoKIyMjIERlZmluZSByZW1vdGUgZm9sZGVyCnJlbW90ZV9mb2xkZXIgPC0gaWYgKGJhY2t1cCkgcGFzdGUoIi9zaGFyZWQvcHJvamVjdHMvIiwgc2Vzc2lvbmlkLCAiL2F0ZWxpZXJfc2NybmFzZXEvVEQvQkFDS1VQLzEwWC8iKSBlbHNlIHBhc3RlMCgiaHR0cHM6Ly96ZW5vZG8ub3JnL3JlY29yZHMvIiwgemVuX2lkLCAiL2ZpbGVzLyIpCgojIyMgUmVjb25zdHJ1Y3QgdGhlIGlucHV0IHBhdGhzCnJlbW90ZV9wYXRoIDwtIHBhc3RlMChyZW1vdGVfZm9sZGVyLCAiLyIsIHRvZ2V0X2ZpbGVzKQoKIyMjIFJlY29uc3RydWN0IHRoZSBvdXRwdXQgcGF0aHMKbG9jYWxfcGF0aCA8LSBwYXN0ZTAobG9jYWxfZm9sZGVyLCAiLyIsIHRvZ2V0X2ZpbGVzKQoKIyMgUmV0cmlldmUgZmlsZXMgKGlmIHRoZXkgZG9uJ3QgZXhpc3QpLCBpbiBsb29wCmZvciAodGcgaW4gc2VxX2Fsb25nKHRvZ2V0X2ZpbGVzKSkgewogICMjIElmIHRoZSBmaWxlIGRvZXMgbm90IGxvY2FsbHkgZXhpc3QKICBpZiAoIWZpbGUuZXhpc3RzKGxvY2FsX3BhdGhbdGddKSB8IGZvcmNlKSB7CiAgICAjIyBSZXRyaWV2ZSBkYXRhCiAgICBpZihiYWNrdXApIHsKICAgICAgZmlsZS5jb3B5KGZyb20gPSByZW1vdGVfcGF0aFt0Z10sCiAgICAgICAgICAgICAgICB0byA9IGxvY2FsX3BhdGhbdGddKQogICAgfSBlbHNlIHsKICAgICAgZG93bmxvYWQuZmlsZSh1cmwgPSByZW1vdGVfcGF0aFt0Z10sIAogICAgICAgICAgICAgICAgICAgIGRlc3RmaWxlID0gbG9jYWxfcGF0aFt0Z10pCiAgICB9CiAgICAjIyBDaGVjayBpZiBkb3dubG9hZGVkIGZpbGVzIGV4aXN0IGxvY2FsbHkKICAgIGlmKGZpbGUuZXhpc3RzKGxvY2FsX3BhdGhbdGddKSkgbWVzc2FnZSgiXHRPSyIpCiAgfSBlbHNlIG1lc3NhZ2UocGFzdGUwKHRvZ2V0X2ZpbGVzW3RnXSwgIiBhbHJlYWR5IGRvd25sb2FkZWQgISIpKQp9CgpgYGAKCiMjIExvYWQgaW50byBSCgotICAgKipMb2FkKiogdGhlICoqMTBYKiogKipkYXRhKiogZmlsZXMgaW50byBSIDoKCiAgICBgYGB7ciBxX3IxMFgsIGV2YWwgPSBGQUxTRSwgY2xhc3Muc291cmNlPSJxdWVzdGlvbiIsIGV2YWwgPSBGQUxTRX0KICAgICMgcV9yMTBYCiAgICAKICAgIEhvdyA/CiAgICAKICAgIGBgYAoKICAgIC0gICBBIGdvb2Qgc3RhcnRlcgoKICAgIGBgYHtyIGFfcjEwWDEsIGNsYXNzLnNvdXJjZSA9IGMoImZvbGQtaGlkZSIsICJhbnN3ZXIiKSwgZWNobyA9IEZBTFNFfQogICAgIyBhX3IxMFgxCiAgICAKICAgIGNhdCgnaHR0cHM6Ly9sbWRkZ3RmeS5uZXQvP3E9bG9hZCUyMDEweCUyMGRhdGElMjBpbiUyMFInKQogICAgCiAgICBgYGAKCiAgICAtICAgQmFjayB0byB5b3VyIGxvY2FsIGhlbHAgcmVzc291cmNlCgogICAgYGBge3IgYV9yMTBYMiwgY2xhc3Muc291cmNlID0gYygiZm9sZC1oaWRlIiwgImFuc3dlciIpLCBldmFsID0gRkFMU0V9CiAgICAjIGFfcjEwWDIKICAgIAogICAgIyMgUmVhZGluZyB0aGUgZnVuY3Rpb24gaGVscCBwYWdlCiAgICA/U2V1cmF0OjpSZWFkMTBYCiAgICAKICAgIGBgYAoKICAgIDxicj48YnI+CiAgICAKICAgIFdlIGNhbiBub3cgYXBwbHkgaXQgdG8gb3VyIGRhdGEgOiAKICAgIAogICAgYGBge3IgbG9hZDEweCwgY2xhc3MgPSAiZm9sZC1oaWRlIn0KICAgICMgbG9hZDEweAogICAgCiAgICAjIyBMb2FkaW5nIHVuZmlsdGVyZWQgMTBYIG1hdHJpeAogICAgc2NtYXQgPC0gU2V1cmF0OjpSZWFkMTBYKGRhdGEuZGlyID0gaW5wdXRfZGlyKQogICAgCiAgICBgYGAKCiAgICA8YnI+CgotICAgSG93ZXZlciwgd2UgZG9uJ3Qga25vdyB3aGF0IHRoZSBsb2FkZWQgb2JqZWN0IGFjdHVhbGx5IGxvb2tzIGxpa2UgaW4KICAgIFIuLi4KCi0gICAqKlF1ZXN0aW9ucyoqIDoKCiAgICAtICAgRG8geW91IGtub3cgYSB3YXkgdG8gLi4uCgogICAgYGBge3IgcV93aGF0aXMsIGNsYXNzLnNvdXJjZT0icXVlc3Rpb24iLCBldmFsID0gRkFMU0V9CiAgICAjIHFfd2hhdGlzCiAgICAKICAgIC4uLiBrbm93IHdoYXQgSVMgdGhlIHR5cGUgb2Ygb2JqZWN0IHRoYXQgd2FzIGNyZWF0ZWQgPwogICAgCiAgICBgYGAKCiAgICA8YnI+CgogICAgYGBge3IgYV93aGF0aXMsIGNsYXNzLnNvdXJjZSA9IGMoImZvbGQtaGlkZSIsICJhbnN3ZXIiKSwgZXZhbCA9IEZBTFNFfQogICAgIyBhX3doYXRpcwogICAgCiAgICAjIyBUbyBrbm93IHRoZSB0eXBlIG9mIGFuIFIgb2JqZWN0IDogbWV0aG9kczo6aXMoKQogICAgP21ldGhvZHM6OmlzCiAgICAKICAgIGBgYAoKICAgIGBgYHtyIGdvaXMsIGNsYXNzLnNvdXJjZSA9ICJmb2xkLWhpZGUifQogICAgIyBnb2lzCiAgICAKICAgICMjIEdldCB0aGUgdHlwZQogICAgbWV0aG9kczo6aXMoc2NtYXQpCiAgICAKICAgIGBgYAoKICAgIDxicj48YnI+CgogICAgYGBge3IgcV9zdHJ1YywgY2xhc3Muc291cmNlPSJxdWVzdGlvbiIsIGV2YWwgPSBGQUxTRX0KICAgICMgcV9zdHJ1YwogICAgCiAgICAuLi4gb2JzZXJ2ZSB0aGUgU1RSdWN0dXJlIG9mIHRoZSByZXN1bHRpbmcgb2JqZWN0ID8KICAgIAogICAgYGBgCgogICAgPGJyPgoKICAgIGBgYHtyIGFfc3RydWMsIGNsYXNzLnNvdXJjZSA9IGMoImZvbGQtaGlkZSIsICJhbnN3ZXIiKSwgZXZhbCA9IEZBTFNFfQogICAgIyBhX3N0cnVjCiAgICAKICAgICMjIFRvIGdldCBhIGJhc2ljIHN0cnVjdHVyZSBkZXNjcmlwdGlvbiA6IHV0aWxzOjpzdHIoKQogICAgP3V0aWxzOjpzdHIKICAgIAogICAgYGBgCgogICAgYGBge3IgZ29fc3RyLCBjbGFzcy5zb3VyY2UgPSAiZm9sZC1oaWRlIn0KICAgICMgZ29fc3RyCiAgICAKICAgICMjIEdldCB0aGUgYmFzaWMgc3RydWN0dXJlCiAgICB1dGlsczo6c3RyKHNjbWF0KQogICAgCiAgICBgYGAKICAgIAogICAgPGJyPgogICAgCiAgICBOb3QgdGhhdCBtZWFuaW5nZnVsLCB0aG91Z2guLi4gQXQgbGVhc3Qgd2Uga25vdyB3ZSBoYXZlIGEgbWF0cml4LCBzbyBvbmUgY2FuCiAgICBnZXQgaXRzIGRpbWVuc2lvbnMgIQogICAgCiAgICBgYGB7ciBzY21hdF9kaW19CiAgICAjIHNjbWF0X2RpbQogICAgCiAgICAjIyBEaW1lbnNpb25zIG9mIGEgPjFEIG9iamVjdAogICAgZGltKHNjbWF0KQogICAgCiAgICBgYGAKICAgIAogICAgPGJyPgoKIyMgRmVhdHVyZXMgY2xlYW51cAoKYFNldXJhdGAgZG9lcyBub3QgbGlrZSBhdCBhbGwgKip1bmRlcnNjb3JlcyoqIChgX2ApIGluIGZlYXR1cmUgbmFtZXMgLi4uCgpXZSB3aWxsIGhhdmUgdG8gOgoKLSAgICoqY29udHJvbCoqIGlmIG91ciBkYXRhc2V0IGNvbnRhaW5zIHNvbWUgPyBIZXJlLCB3ZSB3aWxsIHVzZSB0aGUgYGJhc2U6OmdyZXAoKWAgZnVuY3Rpb24gdGhhdCBsb29rcyBpbiBjaGFyYWN0ZXIgc3RyaW5ncyBgeGAgaWYgdGhleSBjb250YWluIGEgcHJvdmlkZWQgc2V0IG9mIGNoYXJhY3RlcnMgYXMgYSBgcGF0dGVybmAsIGFuZCByZXR1cm5zIHRoZWlyIHBvc2l0aW9uICh3aGVuIGB2YWx1ZSA9IEZBTFNFYCwgZGVmYXVsdCkgb3IgdmFsdWUgaXRzZWxmIChgdmFsdWUgPSBUUlVFYCkKCiAgICBgYGB7ciB1c19jaGVjazF9CiAgICAjIHVzX2NoZWNrMQogICAgCiAgICAjIyBDaGVjayBmZWF0dXJlIG5hbWVzIHdpdGggdW5kZXJzY29yZShzKQogICAgYmFzZTo6Z3JlcChwYXR0ZXJuID0gJ18nLCAKICAgICAgICAgeCA9IHJvd25hbWVzKHNjbWF0KSwgCiAgICAgICAgIHZhbHVlID0gVFJVRSkKICAgIAogICAgYGBgCgogICAgPGJyPgoKLSAgIGlmIHNvLCAqKmNsZWFuKiogaXQgISBIZXJlIHdlIHdpbGwgdXNlIHRoZSBgYmFzZTo6Z3N1YigpYCBmdW5jdGlvbiB0aGF0IGxvb2tzIGZvciBhIGBwYXR0ZXJuYCBpbiBhIGNoYXJhY3RlciBzdHJpbmcgYHhgIGFuZCBzdWJzdGl0dXRlcyBpdCB3aXRoIGEgYHJlcGxhY2VtZW50YCAgY2hhcmFjdGVyIHN0cmluZy4KCiAgICBgYGB7ciB1c19jbGVhbn0KICAgICMgdXNfY2xlYW4KICAgIAogICAgIyMgQ2xlYW4gZmVhdHVyZXMgOiByZXBsYWNpbmcgdW5kZXJzY29yZXMgYnkgZGFzaGVzCiAgICByb3duYW1lcyhzY21hdCkgPC0gYmFzZTo6Z3N1YigKICAgICAgcGF0dGVybiA9ICJfIiwgCiAgICAgIHJlcGxhY2VtZW50ID0gIi0iLCAKICAgICAgeCA9IHJvd25hbWVzKHNjbWF0KSkKICAgIAogICAgYGBgCgotICAgdGhlbiwgKipjb250cm9sKiogdGhlIGVmZmljaWVuY3kgb2Ygb3VyIGNsZWFuaW5nICh3aXRoIGEgbmV3IHVzZSBvZiBgYmFzZTo6Z3JlcGApCgogICAgYGBge3IgdXNfY2hlY2syLCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiAgICAjIHVzX2NoZWNrMgogICAgCiAgICAjIyBQb3N0LWNsZWFudXAgY29udHJvbAogICAgYmFzZTo6Z3JlcChwYXR0ZXJuID0gIl8iLCAKICAgICAgICAgICAgICAgeCA9IHJvd25hbWVzKHNjbWF0KSwgCiAgICAgICAgICAgICAgIHZhbHVlID0gVFJVRSkKICAgIAogICAgYGBgCgo8YnI+CgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyBFbXB0eSBkcm9wbGV0cwoKLSAgIFdoaWxlIGEgbGFyZ2UgcGFydCBvZiBvdXIgYmFyY29kZXMgY29udGFpbiAqKm5vIGNvdW50KiogYXQgYWxsLAogICAgYmFyY29kZXMgd2l0aCBhdCBsZWFzdCBvbmUgdG90YWwgY291bnQgKiptYXkgbm90KiogbWFuZGF0b3JpbHkKICAgIGNvcnJlc3BvbmQgdG8gKip0cnVlIGNlbGxzKiosIGR1ZSB0byBhbWJpZW50IFJOQS4KCi0gICBIZXJlIHdlIHdhbnQgdG8gKippZGVudGlmeSoqIGVtcHR5IGRyb3BsZXRzLCB0aGVuICoqcmVtb3ZlKiogdGhlbQogICAgZnJvbSB0aGUgbWF0cml4LgoKIyMgRGV0ZWN0aW9uCgotICAgVG8gZGV0ZWN0IHRoZSBtYXJnaW5zIG9mICJ0cnVlIiBjZWxscyB2ZXJzdXMgZW1wdHkgZHJvcGxldHMsIHdlIHdpbGwKICAgIHJlbHkgb24gdGhlIHNvLWNhbGxlZCAqKmBkb3VibGUta25lZXBsb3RgKiogOiBhIHBsb3Qgb2YgKipVTUkgY291bnRzIHBlciAKICAgIGJhcmNvZGUqKiwgaW4gZnVuY3Rpb24gb2YgdGhlaXIgaW5jcmVhc2luZ2x5ICoqb3JkZXJlZCByYW5rKiouCgo8YnI+CjxjZW50ZXI+IVtdKGVsYm93X2tuZWVfS05FRS5wbmcpPC9jZW50ZXI+Cjxicj4KCi0gICBUbyBkcmF3IHN1Y2ggcGxvdCwgd2UgbmVlZCB0byAqKnJhbmsgdGhlIGJhcmNvZGVzKiogdGhhbmtzIHRvIHRoZWlyIAogICAgdG90YWwgY291bnRzLgoKIyMjIFJhbmtpbmcgYmFyY29kZXMKCi0gICBXZSB3aWxsIHVzZSB0aGUKICAgIFsqKkRyb3BsZXRVdGlscyoqXShodHRwczovL2Jpb2NvbmR1Y3Rvci5vcmcvcGFja2FnZXMvcmVsZWFzZS9iaW9jL2h0bWwvRHJvcGxldFV0aWxzLmh0bWwpe3RhcmdldD0iX2JsYW5rIn0gcGFja2FnZQoKICAgIGBgYHtyIGhfYmFyY29kZVJhbmtzLCBldmFsID0gRkFMU0V9CiAgICAjIGhfYmFyY29kZVJhbmtzCiAgICAKICAgICMjIFJlYWQgdGhlIGZ1bmN0aW9uIGhlbHAgcGFnZQogICAgP0Ryb3BsZXRVdGlsczo6YmFyY29kZVJhbmtzKCkKICAgIAogICAgYGBgCgogICAgYGBge3IgYmNyYW5rc30KICAgICMgYmNyYW5rcwogICAgCiAgICAjIyBHZW5lcmF0ZSB0aGUgcmFuayBzdGF0aXN0aWNzCiAgICBiY19yYW5rIDwtIERyb3BsZXRVdGlsczo6YmFyY29kZVJhbmtzKHNjbWF0KQogICAgCiAgICBgYGAKCi0gICBDaGVjayB0aGUgcmFua2luZyBvZiBiYXJjb2RlcywgYWNjb3JkaW5nIHRvIHRoZWlyIGNvdW50cwoKICAgIGBgYHtyIGJjcmFua190b3AsIGNsYXNzLnNvdXJjZT0ibm90cnVuIiwgY2xhc3Mub3V0cHV0PSJub3RydW5vIiwgZm9sZC5vdXRwdXQgPSBGQUxTRX0KICAgICMgYmNyYW5rX3RvcAogICAgCiAgICAjIyBTaG93IGEgc3VtbWFyeSBvZiB0aGUgcmFua3MgKG9yZGVyZWQgZm9yIGRpc3BsYXkgY29udmVuaWVuY2UpCiAgICBwcmludChiY19yYW5rW29yZGVyKGJjX3JhbmskcmFuayksXSkKICAgIAogICAgYGBgCgo8YnI+CgojIyMgS25lZXBsb3QKCk5vdywgd2UgY2FuICoqZHJhdyoqIG91ciBrbmVlcGxvdCA6CgpgYGB7ciBrbmVlcGxvdDEsIGNsYXNzLnNvdXJjZT0ibm90cnVuIiwgY2xhc3Mub3V0cHV0PSJub3RydW5vIiwgZm9sZC5wbG90ID0gRkFMU0V9CiMga25lZXBsb3QxCgojIyBQbG90CmdyYXBoaWNzOjpwbG90KHggPSBiY19yYW5rJHJhbmssIAogICAgICAgICAgICAgICB5ID0gYmNfcmFuayR0b3RhbCArIDEsIAogICAgICAgICAgICAgICBsb2cgPSAieHkiLCAKICAgICAgICAgICAgICAgeGxhYiA9ICJCYXJjb2RlIHJhbmsiLCAKICAgICAgICAgICAgICAgeWxhYiA9ICJUb3RhbCBVTUlzIiwgCiAgICAgICAgICAgICAgIGNvbCA9ICJibGFjayIsIAogICAgICAgICAgICAgICBwY2ggPSAiLiIsIAogICAgICAgICAgICAgICBjZXggPSA1LCAKICAgICAgICAgICAgICAgbWFpbiA9ICJLbmVlcGxvdCIpCgpgYGAKCjxicj4KCi0gICAqKlF1ZXN0aW9uKiogOgoKICAgIGBgYHtyIHFfa25lZSwgY2xhc3Muc291cmNlPSJxdWVzdGlvbiIsIGV2YWwgPSBGQUxTRX0KICAgICMgcV9rbmVlCiAgICAKICAgIExvb2tpbmcgYXQgdGhlIGtuZWVwbG90LCBjYW4geW91IHJvdWdobHkgcHJlZGljdCBob3cgbWFueSBiYXJjb2RlcyB3aWxsIGJlIAogICAga2VwdCBhcyBjZWxscyAoaWUsIGFzIG5vbi1lbXB0eSBkcm9wbGV0cykgPwogICAgCiAgICBgYGAKCiAgICA8YnI+CgogICAgYGBge3IgYV9rbmVlLCBjbGFzcy5zb3VyY2UgPSBjKCJmb2xkLWhpZGUiLCAiYW5zd2VyIiksIGV2YWwgPSBGQUxTRX0KICAgICMgYV9rbmVlCiAgICAKICAgICMjIC4gVGhlICJjbGlmZiIgaW4gdGhlIGtuZWVwbG90IGlzIGF0IGEgCiAgICAjIyAgcmFuayB2YWx1ZSBhIGJpdCBiZWxvdyB+IDEsMDAwLgogICAgCiAgICBgYGAKClRoZW4sIHdlIGNhbiB1c2UgdGhlc2UgcmFua3MgdG8gKippZGVudGlmeSAidHJ1ZSIgY2VsbHMqKiA6CgpgYGB7ciBoX2VtcHR5RHJvcHMsIGNsYXNzLnNvdXJjZT0ibm90cnVuIiwgY2xhc3Mub3V0cHV0PSJub3RydW5vIiwgZXZhbCA9IEZBTFNFfQojIGhfZW1wdHlEcm9wcwoKIyMgUmVhZCB0aGUgZnVuY3Rpb24gaGVscCBwYWdlCj9Ecm9wbGV0VXRpbHM6OmVtcHR5RHJvcHMKCmBgYAoKVGhlIGZ1bmN0aW9uIHdlIHdpbGwgdXNlIHJlbGllcyBvbiByYW5kb20gbnVtYmVyIGdlbmVyYXRpb24sIHNvIHdlIHdpbGwgZml4IHRoZSBSTkcgc2VlZCwgZm9yIHJlcHJvZHVjaWJpbGl0eQoKYGBge3IgZW1wdHlkcm9wczF9CiMgZW1wdHlkcm9wczEKCiMjIFNldCB0aGUgUk5HIHNlZWQKc2V0LnNlZWQobXlfc2VlZCkKCiMjIElkZW50aWZ5IGVtcHR5IGRyb3BsZXRzCmJjX3JhbmsyIDwtIERyb3BsZXRVdGlsczo6ZW1wdHlEcm9wcyhzY21hdCkKCiMjIFdoYXRzIGlzIHRoZSBzdHJ1Y3R1cmUgb2YgdGhlIHJldHVybmVkIG91dHB1dCA/CnN0cihiY19yYW5rMikKYGBgCgpXZSBjYW4gc2hvdyBhIHN1bW1hcnkgb2YgdGhlIG9yZGVyZWQgYmFyY29kZXMgOgoKYGBge3IgZW1wdHlzbXJ5LCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiMgZW1wdHlzbXJ5CgojIyBTaG93IGEgc3VtbWFyeSBvZiB0aGUgcXVhbGlmaWVkIHJhbmtzIChvcmRlcmVkIGZvciBkaXNwbGF5IGNvbnZlbmllbmNlKQpwcmludChiY19yYW5rMltvcmRlcihiY19yYW5rMiRUb3RhbCwgZGVjcmVhc2luZyA9IFRSVUUpLF0pCgpgYGAKCjxicj4KCiMjIyBTZWxlY3Rpb24KCkJhcmNvZGVzIGNvbnNpZGVyZWQgYXMgZW1wdHkgZHJvcGxldHMgaGF2ZSAqKm5vIHAtdmFsdWUqKiAoYEZEUmAgPSBOQSkKCmBgYHtyIGJjcjJfbmF9CiMgYmNyMl9uYQoKIyMgQ3JlYXRpbmcgYSB2YWxpZGF0aW9uIGNvbHVtbiAoZW1wdHkgZHJvcGxldHMgaGF2ZSBOQSBhcyBwLXZhbHVlKQpiY19yYW5rMiRWQUxJRCA8LSAhaXMubmEoYmNfcmFuazIkRkRSKSAKCiMjIFF1YW50aWZ5ICJyZWFsIGNlbGxzIgp0YWJsZShiY19yYW5rMiRWQUxJRCkKCmBgYAoKPGJyPgoKV2UgY2FuICoqY29udHJvbCB0aGUgc3RyaW5nZW5jeSoqIGJ5IHNldHRpbmcgYSBGRFItYWRqdXN0ZWQgcC12YWx1ZQoqKnRocmVzaG9sZCoqCgpgYGB7ciB2YWxpZF9mZHJ9CiMgdmFsaWRfZmRyCgojIyBSZWNhbGwgdGhlIG1heCBwLXZhbHVlIHdlIHNldCBhcyBhIHZhcmlhYmxlIGVhcmx5IG9uCm1heF9wCgojIyBXZSd2ZSBzZXQgdGhlIG1heCBwLXZhbHVlIHRvIDFFLTAzCmJjX3JhbmsyJFZBTElEW2JjX3JhbmsyJEZEUiA+PSBtYXhfcF0gPC0gRkFMU0UKCiMjIENvbnRyb2wKdGFibGUoYmNfcmFuazIkVkFMSUQpCgpgYGAKCjxicj4KCk5vdywgd2UgY2FuICoqZGlzcGxheSoqIG91ciAqKnRydWUgY2VsbHMqKiBvbiB0aGUga25lZXBsb3QgOgoKYGBge3Iga25lZXBsb3QyLCBmaWcuYWxpZ249ImNlbnRlciJ9CiMga25lZXBsb3QyCgojIyBQbG90IHdpdGggaGlnaGxpZ2h0IG9mIHRoZSBzZWxlY3RlZCBiYXJjb2RlcyAocmVkKQpncmFwaGljczo6cGxvdCh4ID0gYmNfcmFuayRyYW5rLCAKICAgICAgICAgICAgICAgeSA9IGJjX3JhbmskdG90YWwrMSwgCiAgICAgICAgICAgICAgIGxvZyA9ICJ4eSIsIAogICAgICAgICAgICAgICB4bGFiID0gIkJhcmNvZGUgcmFuayIsIAogICAgICAgICAgICAgICB5bGFiID0gIlRvdGFsIFVNSXMiLCAKICAgICAgICAgICAgICAgY29sID0gaWZlbHNlKGJjX3JhbmsyJFZBTElELCAicmVkIiwgImJsYWNrIiksIAogICAgICAgICAgICAgICBwY2ggPSAiLiIsIAogICAgICAgICAgICAgICBjZXggPSBpZmVsc2UoYmNfcmFuazIkVkFMSUQsIDcsIDMpLCAKICAgICAgICAgICAgICAgbWFpbiA9IHBhc3RlMCgiS25lZXBsb3QgKCIsIAogICAgICAgICAgICAgICAgICAgICAgICAgICAgIGxlbmd0aCh3aGljaChiY19yYW5rMiRWQUxJRCkpLCAKICAgICAgICAgICAgICAgICAgICAgICAgICAgICAiIGJhcmNvZGVzIGtlcHQpIikpCgpgYGAKCjxicj4KCi0gICAqKlF1ZXN0aW9uKiogOgoKICAgIGBgYHtyIHFfa25lZXJlZCwgY2xhc3Muc291cmNlPSJxdWVzdGlvbiIsIGV2YWwgPSBGQUxTRX0KICAgICMgcV9rbmVlcmVkCiAgICAKICAgIElzIHRoZXJlIHNvbWV0aGluZyB1bmV4cGVjdGVkIGZvciB0aGUgcmV0YWluZWQgIlRSVUUiIGNlbGxzIChyZWQpIGluIHRoaXMga25lZXBsb3QgPwoKICAgIGBgYAoKICAgIDxicj4KCiAgICBgYGB7ciBhX2tuZWVyZWQsIGNsYXNzLnNvdXJjZSA9IGMoImZvbGQtaGlkZSIsICJhbnN3ZXIiKSwgZXZhbCA9IEZBTFNFfQogICAgIyBhX2tuZWVyZWQKICAgIAogICAgIyMgLiBZb3UgbWF5IG9ic2VydmUgdGhhdCB0aGUgc2VsZWN0aW9uIG9mICJ0cnVlIiBjZWxscyAocmVkKQogICAgIyMgICBkb2VzIG5vdCBzdHJpY3RseSBjb3JyZXNwb25kIHRvIGEgImN1dCIgaW4gdGhlIGtuZWVwbG90CiAgICAjIyAgIGRlc2NlbmRpbmcgY3VydmUuCiAgICAjIwogICAgIyMgLiBUaGlzIGlzIGJlY2F1c2UgZW1wdHlEcm9wcyBkb2VzIG5vdCBvbmx5ICJmb2xsb3ciIHRoZSBjdXJ2ZQogICAgIyMgICB0byBkZXRlcm1pbmUgYSB0aHJlc2hvbGQgY29ycmVzcG9uZGluZyB0byBhICJtaW5pbWFsIGNvdW50IiAKICAgICMjICAgdmFsdWUgdG8gY29uc2lkZXIgYSBiYXJjb2RlIGFzIGEgY2VsbCwgYnV0IGFsc28gcGVyZm9ybXMgYSAKICAgICMjICAgc3RhdGlzdGljYWwgYW5hbHlzaXMgZm9yIGVhY2ggYmFyY29kZSwgY29uc2lkZXJpbmcgaG93IG11Y2ggCiAgICAjIyAgIGl0cyBleHByZXNzaW9uIHByb2ZpbGUgcmVzc2VtYmxlcyB0aGUgb25lIG9mIG90aGVyIGNlbGxzLCAKICAgICMjICAgZXZlbiB3aXRoIGEgdmVyeSBkaWZmZXJlbnQgKGxvd2VyKSBnbG9iYWwgbGV2ZWwgb2YgZXhwcmVzc2lvbi4KICAgIAogICAgYGBgCgo8YnI+PGJyPgoKIyMgUmVtb3ZhbAoKV2UganVzdCBuZWVkIHRvIHJlc3RyaWN0IG91ciBjb3VudCBtYXRyaXggdG8gInRydWUiIGNlbGxzLgoKLSAgIEhvdyBtYW55IGJhcmNvZGVzIGRvIHdlIGhhdmUgYXQgdGhlIG1vbWVudCA/CgogICAgYGBge3Igc2NkaW0zLCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiAgICAjIHNjZGltMwogICAgCiAgICBkaW0oc2NtYXQpCiAgICAKICAgIGBgYAoKICAgIDxicj4KCi0gICBGaWx0ZXIgZW1wdHkgZHJvcGxldHMgb3V0IDoKCiAgICBgYGB7ciBlZGZpbHRlcn0KICAgICMgZWRmaWx0ZXIKICAgIAogICAgIyMgUmVzdHJpY3QgdG8gdmFsaWQgYmFyY29kZXMKICAgIHNjbWF0X2NlbGxzIDwtIHNjbWF0WywgYmNfcmFuazIkVkFMSURdCiAgICAKICAgICMjIENvbnRyb2wKICAgIGRpbShzY21hdF9jZWxscykKICAgIAogICAgYGBgCgo8YnI+CgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKIyBBbWJpZW50IFJOQQoKLSAgIFRoZSBjb3VudHMgbWVhc3VyZWQgaW4gdGhpcyBtYXRyaXggb2YgKG5vdyBjb25zaWRlcmVkICJ0cnVlIikgc2luZ2xlCiAgICBjZWxscyBkbyAqKm5vdCoqIHJlZmxlY3QgdGhlIG1lYXN1cmVtZW50IG9mIHRoZXNlIGNlbGxzCiAgICBleHByZXNzaW9uICoqb25seSoqCgotICAgRHVlIHRvIG90aGVyICoqb3Zlci1seXNlZCAvIGR5aW5nIGNlbGxzKiogaW4gdGhlIGVtdWxzaW9uLCB3ZSBhbHNvCiAgICBtZWFzdXJlZCBmcmFjdGlvbiBvZiAqKmFtYmllbnQgUk5BKiogdGhhdCBhZGRzIHVwIHRvIHRoZSAicmVhbCIKICAgIGV4cHJlc3Npb24KCi0gICBXZSBtYXkgKiplc3RpbWF0ZSoqIGl0IGJ5IG9ic2VydmluZyBjb3VudHMgZXhpc3RpbmcgYXQgYSAqKm1pbmltYWwKICAgIGxldmVsKiogYWNyb3NzIGNlbGxzIHRoYXQgZ3JlYXRseSBkaWZmZXIgaW4gdGhlaXIgZXhwcmVzc2lvbiBwcm9maWxlCiAgICAKLSAgIEp1c3QgaW4gY2FzZSwgZHVlIHRvIHRoZSBmdW4gdG9vbCBuYW1lIHdlIHdpbGwgdXNlLCBgU291cFhgLCBJJ2xsCiAgICB3cml0ZSBkb3duICJzb3VwIiBhcyBhbiBhbGlhcyBmb3IgImFtYmllbnQgUk5BIiwgYnV0IGJvdGggdGVybXMgaGF2ZQogICAgdGhlIHNhbWUgbWVhbmluZy4KCjxicj4KCiMjIEVzdGltYXRlIHRoZSAic291cCIKClsqKlNvdXBYKipdKGh0dHBzOi8vY3Jhbi5yLXByb2plY3Qub3JnL3dlYi9wYWNrYWdlcy9Tb3VwWC9yZWFkbWUvUkVBRE1FLmh0bWwpe3RhcmdldD0iX2JsYW5rIn0KY2FuIGVzdGltYXRlIHRoZSByYXRlIGFuZCBjb21wb3NpdGlvbiBvZiBhbWJpZW50IFJOQSwgaW4gc2V2ZXJhbCB3YXlzIDoKCi0gICAqKk1hbnVhbCBtb2RlKiogOiBwcm92aWRpbmcgYSBsaXN0IG9mIGZlYXR1cmVzIGV4cGVjdGVkIHRvIHJlZmxlY3QKICAgIHRoZSBzb3VwIChnZW5lcyBleHByZXNzZWQgYXQgaGlnaCBsZXZlbHMgaW4gYSBtYWpvcml0eSBvZiB0aGUKICAgIGV4cGVjdGVkIGNlbGwgdHlwZXMpCgotICAgKipBdXRvbWF0aWMgbW9kZXMqKiwgdXNpbmcgdGhlIGVtcHR5IGRyb3BsZXRzLSoqZmlsdGVyZWQqKiByYXcgY291bnQKICAgIG1hdHJpeCAuLi4KCiAgICAtICAgLi4uIGFsb25nICoqd2l0aCB0aGUgdW5maWx0ZXJlZCoqIG1hdHJpeCAoaWUsIGNvbnRhaW5pbmcgYWxsCiAgICAgICAgcXVhbnRpZmllZCBkcm9wbGV0cykuIFRoaXMgaXMgdGhlIG1vc3QgZWZmaWNpZW50IG1vZGUuCgogICAgLSAgIC4uLiAqKnNvbGVseSoqIDogdGhpcyBpcyBsZXNzIGVmZmljaWVudCwgYnV0IHVzZWZ1bCB3aGVuIHRoZQogICAgICAgIHVuZmlsdGVyZWQgbWF0cml4IGlzIG5vdCBhdmFpbGFibGUuCgpgYGB7ciBoX1NvdXBYX2F1dG8sIGV2YWwgPSBGQUxTRSwgY2xhc3Muc291cmNlPSJub3RydW4iLCBjbGFzcy5vdXRwdXQ9Im5vdHJ1bm8ifQojIGhfU291cFhfYXV0bwoKIyMgUmVhZGluZyB0aGUgU291cFggcGFja2FnZSBtYWluIGhlbHAgcGFnZQo/U291cFg6OlNvdXBYCgpgYGAKCkxldCdzIHJ1biBTb3VwWCwganVzdCBhc2tpbmcgZm9yIGFuIGVzdGltYXRpb24gb2YgdGhlICJzb3VwIiBwcm9wb3J0aW9uCjoKCmBgYHtyIHNvdXB4cmhvRX0KIyBzb3VweHJob0UKCiMjIFRJUCA6IEluIGNhc2Ugb2YgZXJyb3IgOgojIGluc3RhbGwucGFja2FnZXMoImlybGJhIiwgdHlwZT0ic291cmNlIiwgZm9yY2U9VFJVRSkKCiMjIFNvdXBYIGZvciBhbWJpZW50IGZyYWN0aW9uIGVzdGltYXRpb24Kc291cF9mcmFjIDwtIFNvdXBYX2F1dG8oCiAgc2NtYXRfZmlsdCA9IHNjbWF0X2NlbGxzLCAKICBzY21hdF9yYXcgPSBzY21hdCkKCmBgYAoKPGJyPgoKKioqTk9URSoqKiA6IFdlIGNhbiBvYnNlcnZlClsqKk1BTEFUMSoqXShodHRwczovL3d3dy5iaW9yeGl2Lm9yZy9jb250ZW50LzEwLjExMDEvMjAyNC4wNy4xNC42MDM0Njl2Mil7dGFyZ2V0PSJfYmxhbmsifSBhcyB0aGUgdG9wIGZpcnN0IGdlbmUgKGFzIHdlbGwgYXMgbXVsdGlwbGUgcmlib3Byb3RlaW4tY29kaW5nIGdlbmVzKS4KCjxicj4KCldoYXQncyB0aGUgZXN0aW1hdGVkIGZyYWN0aW9uIG9mIHNvdXAgaW4gb3VyIGNvdW50cyA/CgpgYGB7ciBzb3VweGZyYWN9CiMgc291cHhmcmFjCgpjYXQoIlNvdXAgZnJhY3Rpb24gOiAiLCBzb3VwX2ZyYWMpCgpgYGAKCjxicj4KCi0gICAqKlF1ZXN0aW9uKiogOgoKICAgIGBgYHtyIHFfc291cHJhdGUsIGNsYXNzLnNvdXJjZT0icXVlc3Rpb24iLCBldmFsID0gRkFMU0V9CiAgICAjIHFfc291cHJhdGUKICAgIAogICAgU2hvdWxkIHdlIHRyeSB0byByZW1vdmUgc3VjaCBhbiBhbW91bnQgb2YgYW1iaWVudCBSTkEgPwogICAgCiAgICBgYGAKCiAgICA8YnI+CgogICAgYGBge3IgYV9zb3VwcmF0ZSwgY2xhc3Muc291cmNlID0gYygiZm9sZC1oaWRlIiwgImFuc3dlciIpLCBldmFsID0gRkFMU0V9CiAgICAjIGFfc291cHJhdGUKICAgIAogICAgVGhlIGJlc3Qgd2F5IHRvIGtub3cgaXMgdG8gdHJ5IGl0IG91dCA9RAogICAgCiAgICBgYGAKCjxicj4KCiMjIFJlbW92ZSB0aGUgInNvdXAiCgpTb3VwWCBjYW4gKipyZW1vdmUqKiB0aGUgc291cCAoc3VidHJhY3QgY291bnRzIGZvciBpZGVudGlmaWVkIHNvdXAgZ2VuZXMKdG8gYWxsIGJhcmNvZGVzL2NlbGxzKSB1c2luZyBkaWZmZXJlbnQgbWV0aG9kcyA6CgotICAgKipBdXRvbWF0aWNhbGx5KiogKHNhZmVyKQoKLSAgIFNwZWNpZnlpbmcgYSAqKmZyYWN0aW9uIG9mIGNvdW50cyoqIHRvIHJlbW92ZS4KCiAgICAtICAgKlByb3MqIDoKCiAgICAgICAgLSAgIEJ5IHJlbW92aW5nIFggJSBvZiByZWFkcyA6CgogICAgICAgICAgICAtICAgaWYgc291cCB3YXMgZWZmZWN0aXZlbHkgcHJlc2VudCAoYXQgYSByYXRlIFw8PSBYKSwgdGhpcwogICAgICAgICAgICAgICAgd2lsbCAqKm1vc3RseSBhZmZlY3Qgc291cCoqIGdlbmVzIGZpcnN0LgoKICAgICAgICAgICAgLSAgIElmIG5vdCwgaXQgd2lsbCAqKmRlY3JlYXNlIGNvdW50cyoqIGdsb2JhbGx5LCB0aHVzIGhhdmUKICAgICAgICAgICAgICAgIG5lZ2xpZ2libGUgZWZmZWN0IGluIGRpZmZlcmVuY2VzIGJldHdlZW4gY2VsbHMgLyBjZWxsCiAgICAgICAgICAgICAgICB0eXBlcy4KCiAgICAtICAgKkNvbnMqIDoKCiAgICAgICAgLSAgICoqQ291bnRzIGFyZSoqIGFscmVhZHkgKipyYXJlKiogaW4gZHJvcGxldC1iYXNlZCBzaW5nbGUgY2VsbAogICAgICAgICAgICB0ZWNobm9sb2dpZXMuLi4gUmVtb3ZpbmcgdGhlbSAiYnkgZGVmYXVsdCIgbWFrZXMgdGhlbSAqKmV2ZW4KICAgICAgICAgICAgbW9yZSByYXJlKiouLi4KCiAgICAgICAgLSAgIElmIGFjdHVhbCBzb3VwICoqcmF0ZSBpcyBoaWdoZXIqKiB0aGFuIHRoZSBkZWZpbmVkIHJhdGUgWCB0bwogICAgICAgICAgICByZW1vdmUsIHRoaXMgbW9kZSB3aWxsIHJlbWFpbiAqKnVuZWZmZWN0aXZlKiouLi4KCmBgYHtyIHNvdXB4cm19CiMgc291cHhybQoKIyMgUnVuIFNvdXBYIGluICJyZW1vdmFsIiBtb2RlCnNjbWF0X3Vuc291cCA8LSBTb3VwWF9hdXRvKAogIHNjbWF0X2ZpbHQgPSBzY21hdF9jZWxscywgCiAgc2NtYXRfcmF3ID0gc2NtYXQsIAogIHJldHVybl9vYmplY3QgPSBUUlVFKQoKYGBgCgo8YnI+CgojIyBFZmZlY3Qgb2YgdGhlICJzb3VwIiByZW1vdmFsCgotICAgTm93LCB3ZSB3YW50IHRvICoqY2hhcmFjdGVyaXplIHRoZSBlZmZlY3QqKiBvZiB0aGUgcmVtb3ZhbCBvZiB0aGUKICAgIGFtYmllbnQgUk5BLCBieSBjb21wYXJpbmcgdGhlIHByZS0gYW5kIHBvc3QtIHNvdXBYIG1hdHJpY2VzLgogICAgCi0gICBBbiBlYXN5IHdheSBpcyB0byAqKmRlc2NyaWJlKiosIHRoZW4gKip2aXN1YWxpemUqKiBzdWNoIGFuIGVmZmVjdC4KCjxicj48YnI+CgojIyMgRGVzY3JpYmUgYW5kIHZpc3VhbGl6ZSB7LnRhYnNldCAudGFic2V0LWZhZGUgLnRhYnNldC1waWxsc30KCiMjIyMgQmVmb3JlIFNvdXBYIChERU1PKQoKV2UgY2FuIGRlc2NyaWJlIG91ciBjb3VudCBtYXRyaXggdGhhbmtzIHRvIHNvbWUgdXNlZnVsIG1ldHJpY3MgOgoKLSAgIFRvdGFsIGNvdW50cwoKYGBge3IgZGVzY0JTX3RvdGNvdW50LCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiMgZGVzY0JTX3RvdGNvdW50CgojIyBBbGwgY291bnRzCmJhc2U6OnN1bShzY21hdF9jZWxscykKCmBgYAoKLSAgIFNwYXJzaXR5IChmcmFjdGlvbiBvZiBgMGBzKQoKYGBge3IgZGVzY0JTX3NwYXJzaXR5LCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiMgZGVzY0JTX3NwYXJzaXR5CgojIyBDb21wdXRlIHNwYXJzaXR5IChudW1iZXIgb2YgbWF0cml4IHZhbHVlcyBhdCAwIC8gY3Jvc3MtcHJvZHVjdCBvZiBtYXRyaXggZGltKQpiYXNlOjpzdW0oc3BhcnNlTWF0cml4U3RhdHM6OmNvbENvdW50cygKICB4ID0gc2NtYXRfY2VsbHMsIHZhbHVlID0gMCkKICApIC8gYmFzZTo6cHJvZChkaW0oc2NtYXRfY2VsbHMpKQoKYGBgCgotICAgRGlzdHJpYnV0aW9uIG9mIGNvdW50cyBpbiBjZWxscwoKYGBge3IgZGVzY0JTX25Db3VudCwgY2xhc3Muc291cmNlPSJub3RydW4iLCBjbGFzcy5vdXRwdXQ9Im5vdHJ1bm8ifQojIGRlc2NCU19uQ291bnQKCiMjIFN1bXMgdXAgY291bnRzIGZvciBlYWNoIGNlbGwKYmFzZTo6c3VtbWFyeShzcGFyc2VNYXRyaXhTdGF0czo6Y29sU3VtczIoeCA9IHNjbWF0X2NlbGxzLCBuYS5ybSA9IFRSVUUpKQoKYGBgCgotICAgRGlzdHJpYnV0aW9uIG9mIGV4cHJlc3NlZCBnZW5lcyBpbiBjZWxscwoKYGBge3IgZGVzY0JTX25GZWF0dXJlLCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiMgZGVzY0JTX25GZWF0dXJlCgojIyBBbGwgZ2VuZXMgbWludXMgMC1jb3VudCBnZW5lcwpiYXNlOjpzdW1tYXJ5KG5yb3coc2NtYXRfY2VsbHMpIC0gc3BhcnNlTWF0cml4U3RhdHM6OmNvbENvdW50cygKICB4ID0gc2NtYXRfY2VsbHMsIAogIHZhbHVlID0gMCwgCiAgbmEucm0gPSBUUlVFKSkKCmBgYAoKIyMjIyBBZnRlciBTb3VwWAoKV2UgY2FuIGRlc2NyaWJlIG91ciBjb3VudCBtYXRyaXggdGhhbmtzIHRvIHNvbWUgdXNlZnVsIG1ldHJpY3MgOgoKLSAgIFRvdGFsIGNvdW50cwoKYGBge3IgZGVzY0FTX3RvdGNvdW50fQojIGRlc2NBU190b3Rjb3VudAoKIyMgQWxsIGNvdW50cwpiYXNlOjpzdW0oc2NtYXRfdW5zb3VwKQoKYGBgCgotICAgU3BhcnNpdHkgKGZyYWN0aW9uIG9mIGAwYHMpCgpgYGB7ciBkZXNjQVNfc3BhcnNpdHl9CiMgZGVzY0FTX3NwYXJzaXR5CgojIyBDb21wdXRlIHNwYXJzaXR5IChudW1iZXIgb2YgbWF0cml4IHZhbHVlcyBhdCAwIC8gY3Jvc3MtcHJvZHVjdCBvZiBtYXRyaXggZGltKQpiYXNlOjpzdW0oc3BhcnNlTWF0cml4U3RhdHM6OmNvbENvdW50cygKICB4ID0gc2NtYXRfdW5zb3VwLCB2YWx1ZSA9IDApCiAgKSAvIGJhc2U6OnByb2QoZGltKHNjbWF0X3Vuc291cCkpCgpgYGAKCi0gICBEaXN0cmlidXRpb24gb2YgY291bnRzIGluIGNlbGxzCgpgYGB7ciBkZXNjQVNfbkNvdW50fQojIGRlc2NBU19uQ291bnQKCiMjIFN1bXMgdXAgY291bnRzIGZvciBlYWNoIGNlbGwKYmFzZTo6c3VtbWFyeShzcGFyc2VNYXRyaXhTdGF0czo6Y29sU3VtczIoeCA9IHNjbWF0X3Vuc291cCwgbmEucm0gPSBUUlVFKSkKCmBgYAoKLSAgIERpc3RyaWJ1dGlvbiBvZiBleHByZXNzZWQgZ2VuZXMgaW4gY2VsbHMKCmBgYHtyIGRlc2NBU19uRmVhdHVyZX0KIyBkZXNjQVNfbkZlYXR1cmUKCiMjIENvbXB1dGUgc3BhcnNpdHkgKG51bWJlciBvZiBtYXRyaXggdmFsdWVzIGF0IDAgLyBjcm9zcy1wcm9kdWN0IG9mIG1hdHJpeCBkaW0pCmJhc2U6OnN1bW1hcnkobnJvdyhzY21hdF9jZWxscykgLSBzcGFyc2VNYXRyaXhTdGF0czo6Y29sQ291bnRzKAogIHggPSBzY21hdF91bnNvdXAsIAogIHZhbHVlID0gMCwgCiAgbmEucm0gPSBUUlVFKSkKCmBgYAoKIyMgIHsudW5udW1iZXJlZH0KCjxicj4KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgo8YnI+CgotICAgKipRdWVzdGlvbnMqKiA6CgogICAgYGBge3IgcV9kZXNjZGlmZiwgY2xhc3Muc291cmNlPSJxdWVzdGlvbiIsIGV2YWwgPSBGQUxTRX0KICAgICMgcV9kZXNjZGlmZgogICAgCiAgICBDYW4geW91IGRlc2NyaWJlIHRoZSBkaWZmZXJlbmNlIChlc3BlY2lhbGx5IGZvciBjb3VudHMgYW5kIHNwYXJzaXR5KSA/CiAgICAgIAogICAgV2FzIGl0IGV4cGVjdGVkID8KICAgIAogICAgYGBgCgogICAgPGJyPgoKICAgIGBgYHtyIGFfZGVzY2RpZmYsIGNsYXNzLnNvdXJjZSA9IGMoImZvbGQtaGlkZSIsICJhbnN3ZXIiKSwgZXZhbCA9IEZBTFNFfQogICAgIyBhX2Rlc2NkaWZmCiAgICAKICAgICMjIC4gVG90YWwgY291bnRzLCBtYXggdmFsdWUsIGNvdW50cyBwZXIgY2VsbCBhbmQgbnVtYmVyIG9mCiAgICAjIyAgIGV4cHJlc3NlZCBmZWF0dXJlcyBwZXIgY2VsbCwgYWxsIGhhdmUgZGVjcmVhc2VkLgogICAgIyMKICAgICMjIC4gVGhlIHNwYXJzaXR5IGxldmVsIGhhcyBzbGlnaHRseSBpbmNyZWFzZWQgKHJlbW92aW5nIHRoZSBzb3VwCiAgICAjIyAgIG9ibGl0ZXJhdGVkIGFsbCBjb3VudHMgZm9yIHNvbWUgZmVhdHVyZXMgaW4gc29tZSBjZWxscykuCiAgICAjIwogICAgIyMgLiBBbGwgb2YgdGhpcyBpcyBhYnNvbHV0ZWx5IGV4cGVjdGVkLgogICAgCiAgICBgYGAKCjxicj4KCldoYXQgYXJlIHRoZSBmZWF0dXJlcyBtb3N0IGFmZmVjdGVkIGJ5IHRoZSBzb3VwIHJlbW92YWwgPwoKYGBge3Igc291cGZlYXR1cmVzLCBjbGFzcy5zb3VyY2U9Im5vdHJ1biIsIGNsYXNzLm91dHB1dD0ibm90cnVubyJ9CiMgc291cGZlYXR1cmVzCgojIyBDb21wdXRlIHRoZSBmcmFjdGlvbiBtYXRyaXggKFBSRSAvIFBPU1QpCiMjICBOT1RFIDogd2UgbWF5IGhhdmUgMCBjb3VudHMgaW4gdGhlIGRpdmlkZXIsIAojIyAgc28gd2UgaW5jcmVtZW50IGJvdGggbWF0cml4IGJ5ICsxICEKY29udF9mcmFjIDwtIChzY21hdF9jZWxscyArIDEpIC8gKHNjbWF0X3Vuc291cCArIDEpCgojIyBGcmFjdGlvbiBvZiBjb3VudHMgZGVjcmVhc2VkIGJ5IHNvdXAgcmVtb3ZhbCwgcGVyIGZlYXR1cmUgKG9yZGVyZWQpCmZlYXRfZnJhYyA8LSBzb3J0KHNwYXJzZU1hdHJpeFN0YXRzOjpyb3dNZWFuczIoY29udF9mcmFjKSwgZGVjcmVhc2luZyA9IFRSVUUpCgojIyMgRGlzcGxheSBUb3AgNSBmZWF0dXJlcwp1dGlsczo6aGVhZChmZWF0X2ZyYWMsIGRlY3JlYXNpbmcgPSBUUlVFKQoKIyMgUmVtb3Zpbmcgb2JqZWN0cyB0byBmcmVlIHNvbWUgUkFNCnJtKGNvbnRfZnJhYywgZmVhdF9mcmFjKQoKYGBgCgpgYGB7ciBxc291cGZlYXQsIGNsYXNzLnNvdXJjZT0icXVlc3Rpb24iLCBldmFsID0gRkFMU0V9CiMgcXNvdXBmZWF0CgpEb2VzIHNvbWV0aGluZyBhYm91dCB0aGVzZSBnZW5lcyBzdHJpa2UgeW91ID8KCmBgYAoKSGludCA6ICoqZnJhY3Rpb24qKiBvciAqKmRpZmZlcmVuY2UqKiA/CgoKYGBge3Igc291cGZlYXR1cmVzZGlmZiwgY2xhc3Muc291cmNlPWMoImZvbGQtaGlkZSIsICJhbnN3ZXIiKSwgY2xhc3Mub3V0cHV0PSJhbnN3ZXJvIn0KIyBzb3VwZmVhdHVyZXNkaWZmCgojIyBDb21wdXRlIHRoZSBESUZGRVJFTkNFIG1hdHJpeCAoUFJFIC8gUE9TVCkKIyMgIE5PVEUgOiB3ZSBtYXkgaGF2ZSAwIGNvdW50cyBpbiB0aGUgZGl2aWRlciwgCiMjICBzbyB3ZSBpbmNyZW1lbnQgYm90aCBtYXRyaXggYnkgKzEgIQpjb250X2RpZiA8LSBzY21hdF9jZWxscyAtIHNjbWF0X3Vuc291cAoKIyMgRnJhY3Rpb24gb2YgY291bnRzIGRlY3JlYXNlZCBieSBzb3VwIHJlbW92YWwsIHBlciBmZWF0dXJlIChvcmRlcmVkKQpmZWF0X2RpZiA8LSBzb3J0KHNwYXJzZU1hdHJpeFN0YXRzOjpyb3dNZWFuczIoY29udF9kaWYpLCBkZWNyZWFzaW5nID0gVFJVRSkKCiMjIyBEaXNwbGF5IFRvcCA1IGZlYXR1cmVzCnV0aWxzOjpoZWFkKGZlYXRfZGlmLCBkZWNyZWFzaW5nID0gVFJVRSkKCiMjIFJlbW92aW5nIG9iamVjdHMgdG8gZnJlZSBzb21lIFJBTQpybShjb250X2RpZiwgZmVhdF9kaWYpCgpgYGAKCjxicj4KCiMjIyBQbG90IHRoZSAidG9wIDEiIGZlYXR1cmUgey50YWJzZXQgLnRhYnNldC1mYWRlIC50YWJzZXQtcGlsbHN9CgpXZSB3YW50IHRvIHZpc3VhbGl6ZSB0aGUgZXhwcmVzc2lvbiBvZiB0aGUgInRvcCIgc291cCBnZW5lLCBMWVosIGluIG91ciBkYXRhc2V0LgoKVG8gZG8gc28sIHdlIHdpbGwgZ2VuZXJhdGUgYSBwcm9qZWN0aW9uIG9mIHRoZSBkYXRhIGludG8gYSByZWR1Y2VkICgyLUQpIAptYXRoZW1hdGljYWwgc3BhY2UuIFRoaXMgaXMgZG9uZSBieSB1c2luZyBhIHNlcmllcyBvZiBSIGZ1bmN0aW9uIGZyb20gdGhlIApTZXVyYXQgcGFja2FnZSwgdGhhdCB3ZSB3b24ndCBkZXNjcmliZSBoZXJlIGFuZCBub3cgOiB1bmRlcnN0YW5kaW5nIHRoZXNlIApjb21tYW5kcyBiZWxvdyBpcyBub3QgdGhlIHNjb3BlIG9mIG91ciBjdXJyZW50IHNlc3Npb24sIGJ1dCBpdCB3aWxsIGJlIGR1cmluZyAKdGhlIG5leHQgZmV3IGRheXMuCgpBY3R1YWxseSwgdGhlc2UgY29tbWFuZHMgcGVyZm9ybSBhbGwgdGhlIHN0ZXBzIHdlJ2xsIHNlZSBpbiB0aGUgbmV4dCBmZXcgZGF5cywgCmluIGEgcm93LCBidXQgd2l0aCBfZGVmYXVsdCB2YWx1ZXNfICh0aHVzLCBhbG1vc3QgY2VydGFpbmx5IGluYWRlcXVhdGUgdG8gb3VyIApjdXJyZW50IGRhdGEpLgoKSW4gY29uc2VxdWVuY2UsIHBsZWFzZSBqdXN0IGNvbnNpZGVyIHRoZSBvdXRwdXQgYXMgYSB3YXkgdG8gcXVpY2tseSBhbmQgZGlydGlseSAKdmlzdWFsaXplIHRoZSBkYXRhIDogdGhlIGRlcGljdGVkIHRvcG9sb2d5IHNob3VsZCBiZSBjb25zaWRlcmVkIGFzIHZlcnkgcm91Z2ggCmFuZCBtb3N0bHkgaW1wZXJmZWN0LgoKIyMjIyBCZWZvcmUgU291cFgKCmBgYHtyIGJzb3VweCwgZmlnLmhlaWdodCA9IDYsIGZpZy53aWR0aCA9IDYsIGNsYXNzLnNvdXJjZT0ibm90cnVuIiwgY2xhc3Mub3V0cHV0PSJub3RydW5vIn0KIyBic291cHgKCiMjIFJvbGwgcHJvY2VzcywgdXNpbmcgYSB0ZW1wb3JhcnkgInNvYmpfcHJlc291cHgiIFNldXJhdCBvYmplY3QKc29ial9wcmVzb3VweCA8LSBTZXVyYXQ6OkNyZWF0ZVNldXJhdE9iamVjdChjb3VudHMgPSBzY21hdF9jZWxscywgcHJvamVjdCA9ICJQQk1DMTBLX0JlZm9yZVNvdXBYIikKc29ial9wcmVzb3VweCA8LSBTZXVyYXQ6Ok5vcm1hbGl6ZURhdGEob2JqZWN0ID0gc29ial9wcmVzb3VweCwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX3ByZXNvdXB4IDwtIFNldXJhdDo6U2NhbGVEYXRhKG9iamVjdCA9IHNvYmpfcHJlc291cHgsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9wcmVzb3VweCA8LSBTZXVyYXQ6OkZpbmRWYXJpYWJsZUZlYXR1cmVzKG9iamVjdCA9IHNvYmpfcHJlc291cHgsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9wcmVzb3VweCA8LSBTZXVyYXQ6OlJ1blBDQShvYmplY3QgPSBzb2JqX3ByZXNvdXB4LCBucGNzID0gMjEsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9wcmVzb3VweCA8LSBTZXVyYXQ6OlJ1blVNQVAob2JqZWN0ID0gc29ial9wcmVzb3VweCwgZGltcyA9IGMoMToyMCksIHZlcmJvc2UgPSBGQUxTRSkKCiMjIFRvcDEgZmVhdHVyZSBwbG90CmZwX3ByZSA8LSBTZXVyYXQ6OkZlYXR1cmVQbG90KG9iamVjdCA9IHNvYmpfcHJlc291cHgsIGZlYXR1cmVzID0gJ0xZWicpICsgU2V1cmF0OjpEYXJrVGhlbWUoKQpwcmludChmcF9wcmUpCgpgYGAKCjxicj4KCiMjIyMgQWZ0ZXIgU291cFgKCmBgYHtyIGFzb3VweCwgZmlnLmhlaWdodCA9IDYsIGZpZy53aWR0aCA9IDZ9CiMgYXNvdXB4CgojIyBSb2xsIHByb2Nlc3MsIHVzaW5nIGEgdGVtcG9yYXJ5ICJzb2JqX3Bvc3Rzb3VweCIgU2V1cmF0IG9iamVjdApzb2JqX3Bvc3Rzb3VweCA8LSBTZXVyYXQ6OkNyZWF0ZVNldXJhdE9iamVjdChjb3VudHMgPSBzY21hdF91bnNvdXAsIHByb2plY3QgPSAiUEJNQzEwS19BZnRlclNvdXBYIikKc29ial9wb3N0c291cHggPC0gU2V1cmF0OjpOb3JtYWxpemVEYXRhKG9iamVjdCA9IHNvYmpfcG9zdHNvdXB4LCB2ZXJib3NlID0gRkFMU0UpCnNvYmpfcG9zdHNvdXB4IDwtIFNldXJhdDo6U2NhbGVEYXRhKG9iamVjdCA9IHNvYmpfcG9zdHNvdXB4LCB2ZXJib3NlID0gRkFMU0UpCnNvYmpfcG9zdHNvdXB4IDwtIFNldXJhdDo6RmluZFZhcmlhYmxlRmVhdHVyZXMob2JqZWN0ID0gc29ial9wb3N0c291cHgsIHZlcmJvc2UgPSBGQUxTRSkKc29ial9wb3N0c291cHggPC0gU2V1cmF0OjpSdW5QQ0Eob2JqZWN0ID0gc29ial9wb3N0c291cHgsIG5wY3MgPSAyMSwgdmVyYm9zZSA9IEZBTFNFKQpzb2JqX3Bvc3Rzb3VweCA8LSBTZXVyYXQ6OlJ1blVNQVAob2JqZWN0ID0gc29ial9wb3N0c291cHgsIGRpbXMgPSBjKDE6MjApLCB2ZXJib3NlID0gRkFMU0UpCgojIyBUb3AxIGZlYXR1cmUgcGxvdApmcF9wb3N0IDwtIFNldXJhdDo6RmVhdHVyZVBsb3Qob2JqZWN0ID0gc29ial9wb3N0c291cHgsIGZlYXR1cmVzID0gJ0xZWicpICsgU2V1cmF0OjpEYXJrVGhlbWUoKQpwcmludChmcF9wb3N0KQoKYGBgCgo8YnI+CgojIyAgey51bm51bWJlcmVkfQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCk1lcmdpbmcgcGxvdHMgZm9yIGVhc2Ugb2YgdXNlIDoKCmBgYHtyIHNvdXB4X3VtYXBzLCBmaWcud2lkdGggPSAxMiwgZmlnLmhlaWdodCA9IDYsIGNsYXNzLnNvdXJjZT0iYmV5b25kIiwgY2xhc3Mub3V0cHV0PSJiZXlvbmRvIn0KIyBzb3VweF91bWFwcwoKIyMgVXNpbmcgdGhlIHBhdGNod29yayBwYWNrYWdlIHRvIG1lcmdlIHBsb3RzIChhbmQgZ2dwbG90MiB0byBhZGQgdGl0bGVzKQpwYXRjaHdvcms6OndyYXBfcGxvdHMoCiAgbGlzdCgKICAgIGZwX3ByZSAmIGdncGxvdDI6OmdndGl0bGUobGFiZWwgPSAiTFlaIGJlZm9yZSBTb3VwWCIpLCAKICAgIGZwX3Bvc3QgJiBnZ3Bsb3QyOjpnZ3RpdGxlKGxhYmVsID0gIkxZWiBhZnRlciBTb3VwWCIpKSwgCiAgbnJvdyA9IDEpCgpgYGAKCjxicj4KCi0gICAqKlF1ZXN0aW9uKioKCiAgICBgYGB7ciBxX3VtYXBzb3VweCwgY2xhc3Muc291cmNlPSJxdWVzdGlvbiIsIGV2YWwgPSBGQUxTRX0KICAgICMgcV91bWFwc291cHgKICAgIAogICAgQ2FuIHlvdSBjb21wYXJlIHRoZSB0d28gcGxvdHMgOgogICAgICAtIGZvciB0aGUgTFlaIGV4cHJlc3Npb24gYWNyb3NzIHRoZSAyLUQgc3BhY2UgPwogICAgICAtIGFib3V0IHRoZSBzcGFjZSB0b3BvbG9naWVzID8KICAgICAgLSBmb3IgdGhlIExZWiBleHByZXNzaW9uIGdsb2JhbGx5ID8KICAgICAgCiAgICBgYGAKCiAgICA8YnI+CgogICAgYGBge3IgYV91bWFwc291cHgsIGNsYXNzLnNvdXJjZSA9IGMoImZvbGQtaGlkZSIsICJhbnN3ZXIiKSwgZXZhbCA9IEZBTFNFfQogICAgIyBhX3VtYXBzb3VweAogICAgCiAgICAjIyAuIFRoZSBsZXZlbCBvZiBleHByZXNzaW9uIG1lYXN1cmVkIGluIGNlbGxzIHdpdGggc21hbGxlciBleHByZXNzaW9uCiAgICAjIyAgIGJlZm9yZSBTb3VwWCBoYXMgZ29uZSB0byBhbG1vc3Qgbm9uZSBhZnRlciA6IHRoZSBhbWJpZW50IGV4cHJlc3Npb24KICAgICMjICAgb2YgdGhpcyBnZW5lIGhhcyBlZmZpY2llbnRseSBiZWVuIHJlbW92ZWQuCiAgICAjIwogICAgIyMgLiBIb3BlZnVsbHksIHRoZSBjbHVzdGVyKHMpIHdpdGggaGlnaGVyIGV4cHJlc3Npb24gcmVtYWluKHMpIGhpZ2guCiAgICAjIwogICAgIyMgLiBUaGUgZXhwcmVzc2lvbiBzY2FsZSB1cHBlciBib3VuZCBhZnRlciBTb3VwWCBpcyBoaWdoZXIgdGhhbiBiZWZvcmUgPyEgCiAgICAjIyAgIEJ1dCB3ZSBSRU1PVkVEIHNvbWUgY291bnRzID8hIFdoaWxlIGJlaW5nIGNvdW50ZXItaW50dWl0aXZlLCB0aGlzIGlzCiAgICAjIyAgIGEgcG9zaXRpdmUgY29uc2VxdWVuY2Ugb2YgdGhlIGFtYmllbnQgcmVtb3ZhbCBvbiB0aGUgc2NhbGluZyBvZgogICAgIyMgICBiYXJjb2RlcyBleHByZXNzaW9uIChzZWUgbGF0ZXIpLiAKICAgICMjCiAgICAjIyAuIFRoZSByYW5nZSBvZiBib3RoIFggYW5kIFkgYXhlcyBoYXZlIGluY3JlYXNlZCwgZGVwaWN0aW5nIGEgYmV0dGVyCiAgICAjIyAgIHNlcGFyYXRpb24gb2YgLyBkaXN0YW5jZSBiZXR3ZWVuIHNvbWUgb2YgdGhlIGNsdXN0ZXJzLgogICAgIyMKICAgICMjIC4gVGhlIGdsb2JhbCB0b3BvbG9neSBvZiBjbHVzdGVycyByZW1haW5zIGNsb3NlIChidXQgbm90IGlkZW50aWNhbCkuCiAgICAjIyAKICAgICMjIC4gQ2x1c3RlcnMgY29tcGFjaXR5IGhhcyBpbmNyZWFzZWQuCiAgICAKICAgIGBgYAoKPGJyPgoKIyMgQ29uY2x1c2lvbgoKLSAgIFJlbW92aW5nIGFtYmllbnQgUk5BIGhhcyBhICoqcG9zaXRpdmUgZWZmZWN0KiogOgoKICAgIC0gICBvbiB0aGUgKipxdWFsaXR5IG9mIHRoZSBleHByZXNzaW9uKiogbGV2ZWwgbWVhc3VyZW1lbnRzCiAgICAKICAgIC0gICBvbiB0aGUgb2JzZXJ2ZWQgKip0b3BvbG9neSoqCiAgICAKICAgIC0gICAuLi4gZGVzcGl0ZSBhbiBlc3RpbWF0aW9uIG9mIG9ubHkgfjUlICEKCi0gICBBY3R1YWxseSwgYW1iaWVudCBSTkEgbGV2ZWwgYW5kIGNvbXBvc2l0aW9uIGlzIChvbmUgb2YpIHRoZSBtYWpvciBzb3VyY2UocykgCm9mIHRoZSAqKmJhdGNoIGVmZmVjdCBiaWFzKiogdGhhdCBtYXkgYWx0ZXIgdGhlIGludGVncmF0aW9uIG9mIGRpZmZlcmVudCAKc2FtcGxlcy4KCjxicj4KCi0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLQoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgo8YnI+PGJyPjxicj4KCioqQmV5b25kKiogOgoKMS4gIE9uZSBjYW4gY2hlY2sgaWYgdGhlICJ1bnNvdXBlZCIgKGllLCBwb3N0LVNvdXBYKSBtYXRyaXggcmV0YWluIHNvbWUgc291cCA/LiBIb3cgd291bGQgeW91IGRvIGl0ID8KCiAgICBgYGB7ciBiX3NvdXB4MiwgY2xhc3Muc291cmNlID0gYygiZm9sZC1oaWRlIiwgImJleW9uZCIpLCBldmFsID0gRkFMU0V9CiAgICAjIGJfc291cHgyCiAgICAKICAgICMjIFNvdXBYIGZvciBhbWJpZW50IGZyYWN0aW9uIGVzdGltYXRpb24KICAgIHNvdXAyX2ZyYWMgPC0gU291cFhfYXV0bygKICAgICAgc2NtYXRfZmlsdCA9IHNjbWF0X3Vuc291cCwgCiAgICAgIHNjbWF0X3JhdyA9IHNjbWF0KQogICAgCiAgICAjIyBSdW4gU291cFggaW4gInJlbW92YWwiIG1vZGUKICAgIHNjbWF0X3Vuc291cDIgPC0gU291cFhfYXV0bygKICAgICAgc2NtYXRfZmlsdCA9IHNjbWF0X3Vuc291cCwgCiAgICAgIHNjbWF0X3JhdyA9IHNjbWF0LCAKICAgICAgcmV0dXJuX29iamVjdCA9IFRSVUUpCiAgICAKICAgIGBgYAoKMi4gIFNhdmUgdGhlIGBzb2JqX3ByZXNvdXB4YCBSIG9iamVjdCBvbiBkaXNrLgoKICAgIGBgYHtyIGJfc2F2ZV9hbnN3ZXIxLCBjbGFzcy5zb3VyY2UgPSBjKCJmb2xkLWhpZGUiLCAiYmV5b25kIiksIGV2YWwgPSBGQUxTRX0KICAgICMgYl9zYXZlX2Fuc3dlcjEKICAgIAogICAgIyMgQSBzaW5nbGUgb2JqZWN0IDogdXNlIFJEUyBmb3JtYXQgd2l0aCBzYXZlUkRTKCkKICAgIHNhdmVSRFMob2JqZWN0ID0gc29ial9wcmVzb3VweCwgZmlsZSA9IHBhc3RlMChvdXRwdXRfZGlyLCAiL3NvYmpfcHJlc291cHguUkRTIikpCiAgICAKICAgIGBgYAoKMy4gIFNhdmUgYm90aCB0aGUgYHNvYmpfcHJlc291cHhgIGFuZCBgc29ial9wb3N0c291cHhgIG9iamVjdHMgaW4gYSBzaW5nbGUgYXJjaGl2ZQogICAgb24gZGlzay4KCiAgICBgYGB7ciBiX3NhdmVfYW5zd2VyMiwgY2xhc3Muc291cmNlID0gYygiZm9sZC1oaWRlIiwgImJleW9uZCIpLCBldmFsID0gRkFMU0V9CiAgICAjIGJfc2F2ZV9hbnN3ZXIyCiAgICAKICAgICMjIE11bHRpcGxlIG9iamVjdHMgOiB1c2UgUmRhdGEgZm9ybWF0IHdpdGggc2F2ZSgpCiAgICBzYXZlKGxpc3QgPSBjKHNvYmpfcHJlc291cHgsIHNvYmpfcG9zdHNvdXB4KSwgZmlsZSA9IHBhc3RlMChvdXRwdXRfZGlyLCAiL3NvYmplY3RzLlJEQSIpKQogICAgCiAgICBgYGAKCjQuICBTYW1lIGFzIDEuIGJ1dCB1c2UgdGhlIGBiemlwMmAgY29tcHJlc3Npb24gYWxnb3JpdGhtCgogICAgYGBge3IgYl9zYXZlX2Fuc3dlcjMsIGNsYXNzLnNvdXJjZSA9IGMoImZvbGQtaGlkZSIsICJiZXlvbmQiKSwgZXZhbCA9IEZBTFNFfQogICAgIyBiX3NhdmVfYW5zd2VyMwogICAgCiAgICAjIyBDb21wcmVzcyB3aXRoIGJ6aXAyCiAgICBzYXZlUkRTKG9iamVjdCA9IHNvYmosIGZpbGUgPSBwYXN0ZTAob3V0cHV0X2RpciwgIi9zb2JqLlJEUyIpLCBjb21wcmVzcyA9ICJiemlwMiIpCiAgICAKICAgIGBgYAoKPGJyPgoKLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tCgotLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KCjxicj48YnI+PGJyPgoKIyBSc2Vzc2lvbgoKRm9yIHJlcHJvZHVjaWJpbGl0eSBhbmQgY29udGV4dCwgaXQgaXMgcmVjb21tZW5kZWQgdG8gaW5jbHVkZSBpbiB5b3VyIFJNYXJrZG93biB0aGUgbGlzdCBvZiBsb2FkZWQgcGFja2FnZXMgYW5kIHRoZWlyIHZlcnNpb24uCgpgYGB7ciByc2Vzc2lvbiwgY2xhc3Muc291cmNlPSJub3RydW4iLCBjbGFzcy5vdXRwdXQ9Im5vdHJ1bm8ifQojIHJzZXNzaW9uCgp1dGlsczo6c2Vzc2lvbkluZm8oKQoKYGBgCg==