Final Project

Evaluates the likelihood of flaws with the statistic BMI with regards to the largely uncontrollable and genetically defined diameter of one’s bones. 

Evaluate accuracy of statistic BMI with regards to the waist diameter, an area on the woman’s body that largely is not increased by muscle mass but is increased by an increase in fat. 

Evaluate accuracy of statistic BMI with regards to the a measurement on the woman’s body that does generally change dramatically by muscle mass in athletic women and less so by body fat.

For further study, we want to evaluate the consistency of the BMI statistic between men and women to determine if the flaws we find in BMI measurements in women carry over to men body measurement statistics in men as well.

A categorical vector, 1 if the respondent is male, 0 if female. It is also the variable we held constant by only using female data.

Differences in body fat due to gender change BMI ratings 

A numerical vector, respondent's biiliac diameter (pelvic breadth) in centimeters.

The pelvic breadth of a woman correlates to the weight of a woman   but   not necessarily the physical fitness or health

A numerical vector, respondent's waist girth in centimeters,    measured at the narrowest part of torso below the rib cage as average   of contracted and relaxed position.

The waist girth of a woman usually is a good measurement of             women’s weight and overall health given that there is not a large       opportunity to build up muscle mass there

A numerical vector, respondent's thigh girth in centimeters, measured below gluteal fold as the average of right and left girths.

The girth of a woman’s thigh is correlated to muscle mass and strength of quadriceps which is generally a healthy trait yet raises the BMI of the individual. 

This is one of the categorical variables that we made and tested in our hypothesis tests. A BMI greater than 25 is categorized as overweight, and under 25 is not overweight. 

We are trying to prove that inaccuracy of the measurement BMI for physically active people because of the inability of the scale to take into account fat and muscle mass. Given that muscle weighs more per area that fat, active individuals with larger amounts of muscle mass and likely to be given a higher BMI and be categorized into overweight and obese categories despite their physical fitness and general health. This is especially true for women in that physically active women have a greater tendency to have larger thigh muscles which would greatly increase their BMI calculation. 

wgt/((hgt)^2)

Underweight = <18.5
Normal weight = 18.5–24.9 
Overweight = 25–29.9 
Obesity = BMI of 30 or greater

bdims <- bdims %>%
  mutate(htmeters = hgt/100)

bdims <- bdims %>%
  mutate(bmi = wgt/(htmeters*htmeters))

bmif <- bdims %>%
  filter(sex == "f")

bmif <- bdims %>%
  mutate(bmi_over = bmi)

bmif <- bdims %>% 
  mutate(bmi_over = ifelse(bmi > 25, "overweight", "not overweight"))

qplot(x = bii.di, y = bmi, data = bmif)

sex= f; and relationship between pelvic girth and bmi

This data shows that there is a variability in the BMI of women with the same pelvic girth demonstrating that there may not be a strong correlation between the variables. There is a slight upward trend that we would expect when taking into consideration the outliers which shows that a woman who has a larger pelvic bone (due to genetics) despite being healthy, still has a higher likelihood of being classified as overweight.

qplot(x = wai.gi, y = bmi, data = bmif)

sex= f; and relationship between waist girth and bmi

As the waist girth increases so does the BMI which shows the trend that an increase in width of waist ( a typically fat prone area on women) does lead to an increase in BMI. 

qplot(x = thi.gi, y = bmi, dat = bmif)

sex = f; and relationship between thigh girth and bmi

This data shows a strong correlation between the variables with an increase in thigh girth also contributing to a higher BMI. Given the assumption that these women are healthy, thigh girth growth is mostly due to an increase in muscle mass which should not necessarily be correlated with being considered less healthy or having a higher BMI yet the data does show a correlation with an increase in BMI.

bmif %>%
  group_by(bmi_over) %>%
  summarise(mean_bii.di = mean(bii.di))

## # A tibble: 2 x 2
##   bmi_over       mean_bii.di
##   <chr>                <dbl>
## 1 not overweight        27.4
## 2 overweight            28.9

bmif %>%
  group_by(bmi_over) %>%
  summarise(mean_wai.gi = mean(wai.gi))

## # A tibble: 2 x 2
##   bmi_over       mean_wai.gi
##   <chr>                <dbl>
## 1 not overweight        72.2
## 2 overweight            88.3

bmif %>%
  group_by(bmi_over) %>%
  summarise(mean_thi.gi = mean(thi.gi))

## # A tibble: 2 x 2
##   bmi_over       mean_thi.gi
##   <chr>                <dbl>
## 1 not overweight        55.2
## 2 overweight            60.8

bmif %>%
  group_by(bmi_over) %>%
  summarise(conditions = n())

## # A tibble: 2 x 2
##   bmi_over       conditions
##   <chr>               <int>
## 1 not overweight        357
## 2 overweight            150

inference(y = bii.di, x = bmi_over, data = bmif, statistic = "mean", type = "ht", null = 0,       alternative = "twosided", method = "theoretical")

## Response variable: numerical
## Explanatory variable: categorical (2 levels) 
## n_not overweight = 357, y_bar_not overweight = 27.3832, s_not overweight = 2.1377
## n_overweight = 150, y_bar_overweight = 28.8933, s_overweight = 1.9989
## H0: mu_not overweight =  mu_overweight
## HA: mu_not overweight != mu_overweight
## t = -7.6043, df = 149
## p_value = < 0.0001

inference(y = wai.gi, x = bmi_over, data = bmif, statistic = "mean", type = "ht", null = 0,       alternative = "twosided", method = "theoretical")

## Response variable: numerical
## Explanatory variable: categorical (2 levels) 
## n_not overweight = 357, y_bar_not overweight = 72.2269, s_not overweight = 7.9122
## n_overweight = 150, y_bar_overweight = 88.2907, s_overweight = 8.9088
## H0: mu_not overweight =  mu_overweight
## HA: mu_not overweight != mu_overweight
## t = -19.1389, df = 149
## p_value = < 0.0001

inference(y = thi.gi, x = bmi_over, data = bmif, statistic = "mean", type = "ht", null = 0,       alternative = "twosided", method = "theoretical")

## Response variable: numerical
## Explanatory variable: categorical (2 levels) 
## n_not overweight = 357, y_bar_not overweight = 55.191, s_not overweight = 3.3025
## n_overweight = 150, y_bar_overweight = 60.818, s_overweight = 4.364
## H0: mu_not overweight =  mu_overweight
## HA: mu_not overweight != mu_overweight
## t = -14.1779, df = 149
## p_value = < 0.0001

ggplot(data = bmif, aes(y= bii.di, x = bmi_over)) +
  geom_boxplot()

ggplot(data = bmif, aes(y= thi.gi, x = bmi_over)) +
  geom_boxplot()

ggplot(data = bmif, aes(y= wai.gi, x = bmi_over)) +
  geom_boxplot()

data(bmif)

## Warning in data(bmif): data set 'bmif' not found

m_all <- lm(bmi ~ bia.di + bii.di + bit.di + che.di + elb.di + wri.di + kne.di + ank.di + sho.gi + che.gi + wai.gi + nav.gi + hip.gi + thi.gi + bic.gi + for.gi + kne.gi + cal.gi + wri.gi + che.de, data = bmif)
summary(m_all)

## 
## Call:
## lm(formula = bmi ~ bia.di + bii.di + bit.di + che.di + elb.di + 
##     wri.di + kne.di + ank.di + sho.gi + che.gi + wai.gi + nav.gi + 
##     hip.gi + thi.gi + bic.gi + for.gi + kne.gi + cal.gi + wri.gi + 
##     che.de, data = bmif)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.6624 -0.7392 -0.0393  0.7245  4.0551 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -5.6502651  1.2183026  -4.638 4.53e-06 ***
## bia.di      -0.1833430  0.0347547  -5.275 2.00e-07 ***
## bii.di      -0.0845945  0.0351749  -2.405 0.016547 *  
## bit.di      -0.0686560  0.0503838  -1.363 0.173621    
## che.di       0.0309861  0.0455375   0.680 0.496542    
## elb.di      -0.2625983  0.1015438  -2.586 0.009998 ** 
## wri.di      -0.0126985  0.1241563  -0.102 0.918578    
## kne.di       0.3677408  0.0754217   4.876 1.47e-06 ***
## ank.di      -0.2843622  0.0820444  -3.466 0.000575 ***
## sho.gi      -0.0003324  0.0172157  -0.019 0.984605    
## che.gi       0.0705312  0.0212135   3.325 0.000952 ***
## wai.gi       0.1136994  0.0136685   8.318 9.05e-16 ***
## nav.gi       0.0324945  0.0129724   2.505 0.012575 *  
## hip.gi       0.0745908  0.0263443   2.831 0.004827 ** 
## thi.gi       0.1348821  0.0282729   4.771 2.43e-06 ***
## bic.gi       0.0569426  0.0459191   1.240 0.215551    
## for.gi       0.2101902  0.0753335   2.790 0.005476 ** 
## kne.gi      -0.0797466  0.0414291  -1.925 0.054826 .  
## cal.gi       0.1953225  0.0352945   5.534 5.12e-08 ***
## wri.gi      -0.1965782  0.1117607  -1.759 0.079220 .  
## che.de      -0.0312200  0.0383934  -0.813 0.416524    
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.171 on 486 degrees of freedom
## Multiple R-squared:  0.8772, Adjusted R-squared:  0.8721 
## F-statistic: 173.5 on 20 and 486 DF,  p-value: < 2.2e-16

m_all[["coefficients"]]

##   (Intercept)        bia.di        bii.di        bit.di        che.di 
## -5.6502650557 -0.1833430389 -0.0845945281 -0.0686560331  0.0309861079 
##        elb.di        wri.di        kne.di        ank.di        sho.gi 
## -0.2625982845 -0.0126985221  0.3677407889 -0.2843621592 -0.0003323548 
##        che.gi        wai.gi        nav.gi        hip.gi        thi.gi 
##  0.0705311659  0.1136993843  0.0324944847  0.0745908106  0.1348821144 
##        bic.gi        for.gi        kne.gi        cal.gi        wri.gi 
##  0.0569425657  0.2101902262 -0.0797465708  0.1953224615 -0.1965781982 
##        che.de 
## -0.0312200081

 qplot(x = .fitted, y = .resid, data = m_all) +
     geom_hline(yintercept = 0, linetype = "dashed") +
     xlab("Fitted values") +
     ylab("Residuals")

qplot(sample = .resid, data = m_all, stat = "qq")

## Warning: `stat` is deprecated

bmif.data <- bmif %>% select(-c('wgt','hgt', 'htmeters', 'bmi_over'))
n.obs <- dim(bmif.data)[1]

train.index <- sample(1:n.obs,floor(.8*n.obs),replace=FALSE)
bmif.train <- bmif.data[train.index,]
bmif.test <- bmif.data[-train.index,] 

sel.data <- bmif.train
m0 <- lm(bmi ~ ., data=sel.data)
summary(m0)

## 
## Call:
## lm(formula = bmi ~ ., data = sel.data)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.7138 -0.7504  0.0154  0.7147  3.8772 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -8.840368   1.463351  -6.041 3.64e-09 ***
## bia.di      -0.138248   0.040310  -3.430  0.00067 ***
## bii.di      -0.085033   0.038978  -2.182  0.02975 *  
## bit.di      -0.107745   0.056698  -1.900  0.05815 .  
## che.de      -0.043746   0.043490  -1.006  0.31511    
## che.di       0.020517   0.051950   0.395  0.69311    
## elb.di      -0.164255   0.113702  -1.445  0.14939    
## wri.di      -0.153904   0.137758  -1.117  0.26461    
## kne.di       0.415292   0.084381   4.922 1.28e-06 ***
## ank.di      -0.147536   0.094164  -1.567  0.11799    
## sho.gi       0.008842   0.019440   0.455  0.64948    
## che.gi       0.062484   0.023941   2.610  0.00941 ** 
## wai.gi       0.146250   0.016955   8.626  < 2e-16 ***
## nav.gi       0.009049   0.015800   0.573  0.56720    
## hip.gi       0.070091   0.029785   2.353  0.01912 *  
## thi.gi       0.108844   0.033110   3.287  0.00111 ** 
## bic.gi       0.089441   0.050676   1.765  0.07837 .  
## for.gi       0.237720   0.084470   2.814  0.00514 ** 
## kne.gi      -0.080924   0.047155  -1.716  0.08695 .  
## cal.gi       0.176656   0.042109   4.195 3.39e-05 ***
## ank.gi      -0.027815   0.061941  -0.449  0.65365    
## wri.gi      -0.075842   0.130461  -0.581  0.56136    
## age          0.003859   0.007642   0.505  0.61387    
## sexm        -1.454909   0.317857  -4.577 6.39e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.164 on 381 degrees of freedom
## Multiple R-squared:  0.8855, Adjusted R-squared:  0.8786 
## F-statistic: 128.2 on 23 and 381 DF,  p-value: < 2.2e-16

single.step.backwards <- function(data,response){
resp.indx <- which(names(data)==response)
y <- data[,resp.indx]
X <- data[,-resp.indx]
n.pred <- dim(X)[2]
if(n.pred > 1){
  for(i in 1:n.pred){
      print(paste0("Variable ", names(X)[i]," removed: Adjusted R-squared =  ",
          round(summary(lm(y~.,data=as.data.frame(X[,-i])))$adj.r.squared,5)))
      }
}
else{
  print("Model only contains one variable.")
}
}

model.best <- lm(bmi ~ ., data=bmif.train)
summary(model.best)

## 
## Call:
## lm(formula = bmi ~ ., data = bmif.train)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -3.7138 -0.7504  0.0154  0.7147  3.8772 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -8.840368   1.463351  -6.041 3.64e-09 ***
## bia.di      -0.138248   0.040310  -3.430  0.00067 ***
## bii.di      -0.085033   0.038978  -2.182  0.02975 *  
## bit.di      -0.107745   0.056698  -1.900  0.05815 .  
## che.de      -0.043746   0.043490  -1.006  0.31511    
## che.di       0.020517   0.051950   0.395  0.69311    
## elb.di      -0.164255   0.113702  -1.445  0.14939    
## wri.di      -0.153904   0.137758  -1.117  0.26461    
## kne.di       0.415292   0.084381   4.922 1.28e-06 ***
## ank.di      -0.147536   0.094164  -1.567  0.11799    
## sho.gi       0.008842   0.019440   0.455  0.64948    
## che.gi       0.062484   0.023941   2.610  0.00941 ** 
## wai.gi       0.146250   0.016955   8.626  < 2e-16 ***
## nav.gi       0.009049   0.015800   0.573  0.56720    
## hip.gi       0.070091   0.029785   2.353  0.01912 *  
## thi.gi       0.108844   0.033110   3.287  0.00111 ** 
## bic.gi       0.089441   0.050676   1.765  0.07837 .  
## for.gi       0.237720   0.084470   2.814  0.00514 ** 
## kne.gi      -0.080924   0.047155  -1.716  0.08695 .  
## cal.gi       0.176656   0.042109   4.195 3.39e-05 ***
## ank.gi      -0.027815   0.061941  -0.449  0.65365    
## wri.gi      -0.075842   0.130461  -0.581  0.56136    
## age          0.003859   0.007642   0.505  0.61387    
## sexm        -1.454909   0.317857  -4.577 6.39e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 1.164 on 381 degrees of freedom
## Multiple R-squared:  0.8855, Adjusted R-squared:  0.8786 
## F-statistic: 128.2 on 23 and 381 DF,  p-value: < 2.2e-16

qplot(x = .fitted, y = .resid, data = model.best) +
         geom_hline(yintercept = 0, linetype = "dashed") +
         xlab("Fitted values") +
        ylab("Residuals")

qplot(sample = .resid, data = model.best, stat = "qq")

## Warning: `stat` is deprecated

predictions.test <- predict(model.best,bmif.test)
 y <- bmif.test$bmi
 mse <- mean((y-predictions.test)^2)
 (rmse <- sqrt(mse))

## [1] 1.140995

bmif.test[2,]

## # A tibble: 1 x 24
##   bia.di bii.di bit.di che.de che.di elb.di wri.di kne.di ank.di sho.gi
##    <dbl>  <dbl>  <dbl>  <dbl>  <dbl>  <dbl>  <dbl>  <dbl>  <dbl>  <dbl>
## 1   44.3   29.9     34   18.4   28.2   13.9   11.2   20.9     15   104.
## # ... with 14 more variables: che.gi <dbl>, wai.gi <dbl>, nav.gi <dbl>,
## #   hip.gi <dbl>, thi.gi <dbl>, bic.gi <dbl>, for.gi <dbl>, kne.gi <dbl>,
## #   cal.gi <dbl>, ank.gi <dbl>, wri.gi <dbl>, age <int>, sex <fct>,
## #   bmi <dbl>

predict(model.best,bmif.test[2,],interval='prediction')

##        fit      lwr      upr
## 1 21.38582 19.00634 23.76531

LS0tCnRpdGxlOiAiRmluYWwgUHJvamVjdCIKYXV0aG9yOiAiS2F0aGxlZW4gVmVybiIKZGF0ZTogIjQvMTUvMjAxOSIKb3V0cHV0OiBvaWxhYnM6OmxhYl9yZXBvcnQKLS0tCgpgYGB7ciBzZXR1cCwgaW5jbHVkZT1GQUxTRX0Ka25pdHI6Om9wdHNfY2h1bmskc2V0KGVjaG8gPSBUUlVFKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShkcGx5cikKbGlicmFyeShnZ3Bsb3QyKQpsaWJyYXJ5KG9pbGFicykKZGF0YShiZGltcykKc2V0LnNlZWQoOTk5KQpgYGAKCiMjIyBUZWFtd29yayBSZXBvcnQKCnwgVGVhbSBtZW1iZXIgbmFtZSB8IEF0dGVuZGVuY2UgfCBBdXRob3IgfCBDb250cmlidXRpb24gJSB8CnwgOiAuLi4uLi4uLi4uLi4uLi4uLi4uLi4uLi4uLi58IDogLi4uLi4uLi58IDouLi4uLi4uLnwgOi4uLi4uLi4uLi4uCnwgS2F0aGxlZW4gVmVybiAgICAgIHwgeWVzIHwgeWVzIHwgMzAlCnwgQW1hbmRhIEJlYWNoICAgICAgIHwgeWVzIHwgbm8gIHwgMjAlCnwgTGV4aSBDbGlmZm9yZCAgICAgIHwgeWVzIHwgbm8gIHwgMjAlCnwgTWljaGVsbGUgUm9kcmlndWV6IHwgeWVzIHwgeWVzIHwgMzAlCnwgVG90YWwgfCB8IHwgMTAwJSB8CgoqICogKgoKIyMjIEdpdmVuOiAKCkJvZHkgZ2lydGggbWVhc3VyZW1lbnRzIGFuZCBza2VsZXRhbCBkaWFtZXRlciBtZWFzdXJlbWVudHMsIGFzIHdlbGwgYXMgYWdlLCB3ZWlnaHQsIGhlaWdodCBhbmQgZ2VuZGVyLCBhcmUgZ2l2ZW4gZm9yIDUwNyBwaHlzaWNhbGx5IGFjdGl2ZSBpbmRpdmlkdWFscyAtIDI0NyBtZW4gYW5kIDI2MCB3b21lbi4KCiMjIyBRdWVzdGlvbnM6IAoKCklzIHRoZXJlIGEgcmVsYXRpb25zaGlwIGJldHdlZW4gc2tlbGV0YWwgZGlhbWV0ZXIgKGJpaWxpYWMpIGFuZCBCTUkgaW4gd29tZW4/IAoKICAgIEV2YWx1YXRlcyB0aGUgbGlrZWxpaG9vZCBvZiBmbGF3cyB3aXRoIHRoZSBzdGF0aXN0aWMgQk1JIHdpdGggcmVnYXJkcyB0byB0aGUgbGFyZ2VseSB1bmNvbnRyb2xsYWJsZSBhbmQgZ2VuZXRpY2FsbHkgZGVmaW5lZCBkaWFtZXRlciBvZiBvbmXigJlzIGJvbmVzLiAKCgpJcyB0aGVyZSBhIHJlbGF0aW9uc2hpcCBiZXR3ZWVuIGJvZHkgZ2lydGggKHdhaXN0KSBhbmQgQk1JIGluIHdvbWVuPwoKICAgIEV2YWx1YXRlIGFjY3VyYWN5IG9mIHN0YXRpc3RpYyBCTUkgd2l0aCByZWdhcmRzIHRvIHRoZSB3YWlzdCBkaWFtZXRlciwgYW4gYXJlYSBvbiB0aGUgd29tYW7igJlzIGJvZHkgdGhhdCBsYXJnZWx5IGlzIG5vdCBpbmNyZWFzZWQgYnkgbXVzY2xlIG1hc3MgYnV0IGlzIGluY3JlYXNlZCBieSBhbiBpbmNyZWFzZSBpbiBmYXQuIAoKSXMgdGhlcmUgYSByZWxhdGlvbnNoaXAgYmV0d2VlbiBib2R5IGdpcnRoICh0aGlnaCkgYW5kIEJNSSBpbiB3b21lbj8KCiAgICBFdmFsdWF0ZSBhY2N1cmFjeSBvZiBzdGF0aXN0aWMgQk1JIHdpdGggcmVnYXJkcyB0byB0aGUgYSBtZWFzdXJlbWVudCBvbiB0aGUgd29tYW7igJlzIGJvZHkgdGhhdCBkb2VzIGdlbmVyYWxseSBjaGFuZ2UgZHJhbWF0aWNhbGx5IGJ5IG11c2NsZSBtYXNzIGluIGF0aGxldGljIHdvbWVuIGFuZCBsZXNzIHNvIGJ5IGJvZHkgZmF0LgoKUGVyY2VudGFnZXMgb2YgbWVuIHYgd29tZW4gd2hvIGFyZSBjb25zaWRlcmVkIG92ZXJ3ZWlnaHQgZ2l2ZW4gdGhlaXIgQk1JIGRlc3BpdGUgYmVpbmcgY29uc2lkZXJlZCBwaHlzaWNhbGx5IGFjdGl2ZT8gIAoKICAgIEZvciBmdXJ0aGVyIHN0dWR5LCB3ZSB3YW50IHRvIGV2YWx1YXRlIHRoZSBjb25zaXN0ZW5jeSBvZiB0aGUgQk1JIHN0YXRpc3RpYyBiZXR3ZWVuIG1lbiBhbmQgd29tZW4gdG8gZGV0ZXJtaW5lIGlmIHRoZSBmbGF3cyB3ZSBmaW5kIGluIEJNSSBtZWFzdXJlbWVudHMgaW4gd29tZW4gY2Fycnkgb3ZlciB0byBtZW4gYm9keSBtZWFzdXJlbWVudCBzdGF0aXN0aWNzIGluIG1lbiBhcyB3ZWxsLgoKIyMjIERhdGEgYW5kIFZhcmlhYmxlcyBDb25zaWRlcmVkOgoKIyMjIyBzZXg6IAoKICAgIEEgY2F0ZWdvcmljYWwgdmVjdG9yLCAxIGlmIHRoZSByZXNwb25kZW50IGlzIG1hbGUsIDAgaWYgZmVtYWxlLiBJdCBpcyBhbHNvIHRoZSB2YXJpYWJsZSB3ZSBoZWxkIGNvbnN0YW50IGJ5IG9ubHkgdXNpbmcgZmVtYWxlIGRhdGEuCgogICAgRGlmZmVyZW5jZXMgaW4gYm9keSBmYXQgZHVlIHRvIGdlbmRlciBjaGFuZ2UgQk1JIHJhdGluZ3MgCgojIyMjIGJpaS5kaTogCgogICAgQSBudW1lcmljYWwgdmVjdG9yLCByZXNwb25kZW50J3MgYmlpbGlhYyBkaWFtZXRlciAocGVsdmljIGJyZWFkdGgpIGluIGNlbnRpbWV0ZXJzLgoKICAgIFRoZSBwZWx2aWMgYnJlYWR0aCBvZiBhIHdvbWFuIGNvcnJlbGF0ZXMgdG8gdGhlIHdlaWdodCBvZiBhIHdvbWFuICAgYnV0ICAgbm90IG5lY2Vzc2FyaWx5IHRoZSBwaHlzaWNhbCBmaXRuZXNzIG9yIGhlYWx0aAoKIyMjIyB3YWkuZ2k6IAogICAgCiAgICBBIG51bWVyaWNhbCB2ZWN0b3IsIHJlc3BvbmRlbnQncyB3YWlzdCBnaXJ0aCBpbiBjZW50aW1ldGVycywgICAgbWVhc3VyZWQgYXQgdGhlIG5hcnJvd2VzdCBwYXJ0IG9mIHRvcnNvIGJlbG93IHRoZSByaWIgY2FnZSBhcyBhdmVyYWdlICAgb2YgY29udHJhY3RlZCBhbmQgcmVsYXhlZCBwb3NpdGlvbi4KCiAgICBUaGUgd2Fpc3QgZ2lydGggb2YgYSB3b21hbiB1c3VhbGx5IGlzIGEgZ29vZCBtZWFzdXJlbWVudCBvZiAgICAgICAgICAgICB3b21lbuKAmXMgd2VpZ2h0IGFuZCBvdmVyYWxsIGhlYWx0aCBnaXZlbiB0aGF0IHRoZXJlIGlzIG5vdCBhIGxhcmdlICAgICAgIG9wcG9ydHVuaXR5IHRvIGJ1aWxkIHVwIG11c2NsZSBtYXNzIHRoZXJlCgojIyMjIHRoaS5naTogCiAgICAKICAgIEEgbnVtZXJpY2FsIHZlY3RvciwgcmVzcG9uZGVudCdzIHRoaWdoIGdpcnRoIGluIGNlbnRpbWV0ZXJzLCBtZWFzdXJlZCBiZWxvdyBnbHV0ZWFsIGZvbGQgYXMgdGhlIGF2ZXJhZ2Ugb2YgcmlnaHQgYW5kIGxlZnQgZ2lydGhzLgoKICAgIFRoZSBnaXJ0aCBvZiBhIHdvbWFu4oCZcyB0aGlnaCBpcyBjb3JyZWxhdGVkIHRvIG11c2NsZSBtYXNzIGFuZCBzdHJlbmd0aCBvZiBxdWFkcmljZXBzIHdoaWNoIGlzIGdlbmVyYWxseSBhIGhlYWx0aHkgdHJhaXQgeWV0IHJhaXNlcyB0aGUgQk1JIG9mIHRoZSBpbmRpdmlkdWFsLiAKCgojIyMjIGJtaV9vdmVyCgogICAgVGhpcyBpcyBvbmUgb2YgdGhlIGNhdGVnb3JpY2FsIHZhcmlhYmxlcyB0aGF0IHdlIG1hZGUgYW5kIHRlc3RlZCBpbiBvdXIgaHlwb3RoZXNpcyB0ZXN0cy4gQSBCTUkgZ3JlYXRlciB0aGFuIDI1IGlzIGNhdGVnb3JpemVkIGFzIG92ZXJ3ZWlnaHQsIGFuZCB1bmRlciAyNSBpcyBub3Qgb3ZlcndlaWdodC4gCgoKIyMjIE92ZXJhbGw6IAoKICAgIFdlIGFyZSB0cnlpbmcgdG8gcHJvdmUgdGhhdCBpbmFjY3VyYWN5IG9mIHRoZSBtZWFzdXJlbWVudCBCTUkgZm9yIHBoeXNpY2FsbHkgYWN0aXZlIHBlb3BsZSBiZWNhdXNlIG9mIHRoZSBpbmFiaWxpdHkgb2YgdGhlIHNjYWxlIHRvIHRha2UgaW50byBhY2NvdW50IGZhdCBhbmQgbXVzY2xlIG1hc3MuIEdpdmVuIHRoYXQgbXVzY2xlIHdlaWdocyBtb3JlIHBlciBhcmVhIHRoYXQgZmF0LCBhY3RpdmUgaW5kaXZpZHVhbHMgd2l0aCBsYXJnZXIgYW1vdW50cyBvZiBtdXNjbGUgbWFzcyBhbmQgbGlrZWx5IHRvIGJlIGdpdmVuIGEgaGlnaGVyIEJNSSBhbmQgYmUgY2F0ZWdvcml6ZWQgaW50byBvdmVyd2VpZ2h0IGFuZCBvYmVzZSBjYXRlZ29yaWVzIGRlc3BpdGUgdGhlaXIgcGh5c2ljYWwgZml0bmVzcyBhbmQgZ2VuZXJhbCBoZWFsdGguIFRoaXMgaXMgZXNwZWNpYWxseSB0cnVlIGZvciB3b21lbiBpbiB0aGF0IHBoeXNpY2FsbHkgYWN0aXZlIHdvbWVuIGhhdmUgYSBncmVhdGVyIHRlbmRlbmN5IHRvIGhhdmUgbGFyZ2VyIHRoaWdoIG11c2NsZXMgd2hpY2ggd291bGQgZ3JlYXRseSBpbmNyZWFzZSB0aGVpciBCTUkgY2FsY3VsYXRpb24uIAogICAgCiMjIyBJbXBvcnRhbnQgRGF0YS9DYWxjdWxhdGlvbnM6IAoKSW1wb3J0YW50IERhdGE6IENhbGN1bGF0aW5nIEJNSQoKQk1JIGVxdWF0aW9uOiAKCiAgICB3Z3QvKChoZ3QpXjIpCiAgICAKaGd0OiBBIG51bWVyaWNhbCB2ZWN0b3IsIHJlc3BvbmRlbnQncyBoZWlnaHQgaW4gY2VudGltZXRlcnMg4oaSIGNvbnZlcnRlZCB0byBtZXRlcnMKCndndDogQSBudW1lcmljYWwgdmVjdG9yLCByZXNwb25kZW50J3Mgd2VpZ2h0IGluIGtpbG9ncmFtcy4KCkJNSSBDYXRlZ29yaWVzOiAKCiAgICBVbmRlcndlaWdodCA9IDwxOC41CiAgICBOb3JtYWwgd2VpZ2h0ID0gMTguNeKAkzI0LjkgCiAgICBPdmVyd2VpZ2h0ID0gMjXigJMyOS45IAogICAgT2Jlc2l0eSA9IEJNSSBvZiAzMCBvciBncmVhdGVyCiAgICAKIyMjIEVEQSBTZXQtdXA6IAoKYGBge3J9CmJkaW1zIDwtIGJkaW1zICU+JQogIG11dGF0ZShodG1ldGVycyA9IGhndC8xMDApCgpiZGltcyA8LSBiZGltcyAlPiUKICBtdXRhdGUoYm1pID0gd2d0LyhodG1ldGVycypodG1ldGVycykpCgpibWlmIDwtIGJkaW1zICU+JQogIGZpbHRlcihzZXggPT0gImYiKQoKYm1pZiA8LSBiZGltcyAlPiUKICBtdXRhdGUoYm1pX292ZXIgPSBibWkpCgpibWlmIDwtIGJkaW1zICU+JSAKICBtdXRhdGUoYm1pX292ZXIgPSBpZmVsc2UoYm1pID4gMjUsICJvdmVyd2VpZ2h0IiwgIm5vdCBvdmVyd2VpZ2h0IikpCmBgYAoKIyMjIEJhY2tncm91bmQ6IAoKIyMjIyBHcmFwaCAxOgoKYGBge3J9CnFwbG90KHggPSBiaWkuZGksIHkgPSBibWksIGRhdGEgPSBibWlmKQpgYGAKaHR0cHM6Ly9sYWJzLWF6LTAyLm9pdC5kdWtlLmVkdTozMDYyMy9ncmFwaGljcy9wbG90LnBuZz93aWR0aD00NDMmaGVpZ2h0PTI5OSZyYW5kb21pemVyPS0xNDYwNDgzNzU4CgpWYXJzIFJlbGF0aW9uc2hpcDogCgogICAgc2V4PSBmOyBhbmQgcmVsYXRpb25zaGlwIGJldHdlZW4gcGVsdmljIGdpcnRoIGFuZCBibWkKCgpEZXNjcmlwdGlvbjogCgogICAgVGhpcyBkYXRhIHNob3dzIHRoYXQgdGhlcmUgaXMgYSB2YXJpYWJpbGl0eSBpbiB0aGUgQk1JIG9mIHdvbWVuIHdpdGggdGhlIHNhbWUgcGVsdmljIGdpcnRoIGRlbW9uc3RyYXRpbmcgdGhhdCB0aGVyZSBtYXkgbm90IGJlIGEgc3Ryb25nIGNvcnJlbGF0aW9uIGJldHdlZW4gdGhlIHZhcmlhYmxlcy4gVGhlcmUgaXMgYSBzbGlnaHQgdXB3YXJkIHRyZW5kIHRoYXQgd2Ugd291bGQgZXhwZWN0IHdoZW4gdGFraW5nIGludG8gY29uc2lkZXJhdGlvbiB0aGUgb3V0bGllcnMgd2hpY2ggc2hvd3MgdGhhdCBhIHdvbWFuIHdobyBoYXMgYSBsYXJnZXIgcGVsdmljIGJvbmUgKGR1ZSB0byBnZW5ldGljcykgZGVzcGl0ZSBiZWluZyBoZWFsdGh5LCBzdGlsbCBoYXMgYSBoaWdoZXIgbGlrZWxpaG9vZCBvZiBiZWluZyBjbGFzc2lmaWVkIGFzIG92ZXJ3ZWlnaHQuCgojIyMjIEdyYXBoIDI6CgpgYGB7cn0KcXBsb3QoeCA9IHdhaS5naSwgeSA9IGJtaSwgZGF0YSA9IGJtaWYpCmBgYApodHRwczovL2xhYnMtYXotMDIub2l0LmR1a2UuZWR1OjMwNjIzL2dyYXBoaWNzL3Bsb3QucG5nP3dpZHRoPTQ0NiZoZWlnaHQ9Mjk5JnJhbmRvbWl6ZXI9MTI3MjExMDk4OSAKClZhcnMgUmVsYXRpb25zaGlwOiAKCiAgICBzZXg9IGY7IGFuZCByZWxhdGlvbnNoaXAgYmV0d2VlbiB3YWlzdCBnaXJ0aCBhbmQgYm1pCgpEZXNjcmlwdGlvbjogCgogICAgQXMgdGhlIHdhaXN0IGdpcnRoIGluY3JlYXNlcyBzbyBkb2VzIHRoZSBCTUkgd2hpY2ggc2hvd3MgdGhlIHRyZW5kIHRoYXQgYW4gaW5jcmVhc2UgaW4gd2lkdGggb2Ygd2Fpc3QgKCBhIHR5cGljYWxseSBmYXQgcHJvbmUgYXJlYSBvbiB3b21lbikgZG9lcyBsZWFkIHRvIGFuIGluY3JlYXNlIGluIEJNSS4gCiMjIyMgR3JhcGggMzoKCmBgYHtyfQpxcGxvdCh4ID0gdGhpLmdpLCB5ID0gYm1pLCBkYXQgPSBibWlmKQpgYGAKaHR0cHM6Ly9sYWJzLWF6LTAyLm9pdC5kdWtlLmVkdTozMDYyMy9ncmFwaGljcy9wbG90LnBuZz93aWR0aD00NDMmaGVpZ2h0PTI5OSZyYW5kb21pemVyPTE4NDAwNDk4OTQKClZhcnMgUmVsYXRpb25zaGlwOiAKCiAgICBzZXggPSBmOyBhbmQgcmVsYXRpb25zaGlwIGJldHdlZW4gdGhpZ2ggZ2lydGggYW5kIGJtaQoKRGVzY3JpcHRpb246CgogICAgVGhpcyBkYXRhIHNob3dzIGEgc3Ryb25nIGNvcnJlbGF0aW9uIGJldHdlZW4gdGhlIHZhcmlhYmxlcyB3aXRoIGFuIGluY3JlYXNlIGluIHRoaWdoIGdpcnRoIGFsc28gY29udHJpYnV0aW5nIHRvIGEgaGlnaGVyIEJNSS4gR2l2ZW4gdGhlIGFzc3VtcHRpb24gdGhhdCB0aGVzZSB3b21lbiBhcmUgaGVhbHRoeSwgdGhpZ2ggZ2lydGggZ3Jvd3RoIGlzIG1vc3RseSBkdWUgdG8gYW4gaW5jcmVhc2UgaW4gbXVzY2xlIG1hc3Mgd2hpY2ggc2hvdWxkIG5vdCBuZWNlc3NhcmlseSBiZSBjb3JyZWxhdGVkIHdpdGggYmVpbmcgY29uc2lkZXJlZCBsZXNzIGhlYWx0aHkgb3IgaGF2aW5nIGEgaGlnaGVyIEJNSSB5ZXQgdGhlIGRhdGEgZG9lcyBzaG93IGEgY29ycmVsYXRpb24gd2l0aCBhbiBpbmNyZWFzZSBpbiBCTUkuCgoKIyMjIyBTdW1tYXJ5IFN0YXRzIAoKYGBge3J9CmJtaWYgJT4lCiAgZ3JvdXBfYnkoYm1pX292ZXIpICU+JQogIHN1bW1hcmlzZShtZWFuX2JpaS5kaSA9IG1lYW4oYmlpLmRpKSkKCmJtaWYgJT4lCiAgZ3JvdXBfYnkoYm1pX292ZXIpICU+JQogIHN1bW1hcmlzZShtZWFuX3dhaS5naSA9IG1lYW4od2FpLmdpKSkKCmJtaWYgJT4lCiAgZ3JvdXBfYnkoYm1pX292ZXIpICU+JQogIHN1bW1hcmlzZShtZWFuX3RoaS5naSA9IG1lYW4odGhpLmdpKSkKCmJtaWYgJT4lCiAgZ3JvdXBfYnkoYm1pX292ZXIpICU+JQogIHN1bW1hcmlzZShjb25kaXRpb25zID0gbigpKQpgYGAKCkV4cGxhbmF0aW9uOiAKCgojIyMgSHlwb3RoZXNpcyB0ZXN0cyAKIyMjIyBJbmZlcmVuY2UgZm9yIGJpaS5kaSBhbmQgYm1pIGFzIG92ZXJ3ZWlnaHQKYGBge3J9CmluZmVyZW5jZSh5ID0gYmlpLmRpLCB4ID0gYm1pX292ZXIsIGRhdGEgPSBibWlmLCBzdGF0aXN0aWMgPSAibWVhbiIsIHR5cGUgPSAiaHQiLCBudWxsID0gMCwgICAgICAgYWx0ZXJuYXRpdmUgPSAidHdvc2lkZWQiLCBtZXRob2QgPSAidGhlb3JldGljYWwiKQpgYGAKCiMjIyMgSW5mZXJlbmNlIGZvciB3YWkuZ2kgYW5kIGJtaSBhcyBvdmVyd2VpZ2h0CmBgYHtyfQppbmZlcmVuY2UoeSA9IHdhaS5naSwgeCA9IGJtaV9vdmVyLCBkYXRhID0gYm1pZiwgc3RhdGlzdGljID0gIm1lYW4iLCB0eXBlID0gImh0IiwgbnVsbCA9IDAsICAgICAgIGFsdGVybmF0aXZlID0gInR3b3NpZGVkIiwgbWV0aG9kID0gInRoZW9yZXRpY2FsIikKYGBgCgojIyMjIEluZmVyZW5jZSBmb3IgdGhpLmdpIGFuZCBibWkgYXMgb3ZlcndlaWdodApgYGB7cn0KaW5mZXJlbmNlKHkgPSB0aGkuZ2ksIHggPSBibWlfb3ZlciwgZGF0YSA9IGJtaWYsIHN0YXRpc3RpYyA9ICJtZWFuIiwgdHlwZSA9ICJodCIsIG51bGwgPSAwLCAgICAgICBhbHRlcm5hdGl2ZSA9ICJ0d29zaWRlZCIsIG1ldGhvZCA9ICJ0aGVvcmV0aWNhbCIpCmBgYAoKIyMjIEVEQQojIyMjIEJveHBsb3RzIApgYGB7cn0KZ2dwbG90KGRhdGEgPSBibWlmLCBhZXMoeT0gYmlpLmRpLCB4ID0gYm1pX292ZXIpKSArCiAgZ2VvbV9ib3hwbG90KCkKCmdncGxvdChkYXRhID0gYm1pZiwgYWVzKHk9IHRoaS5naSwgeCA9IGJtaV9vdmVyKSkgKwogIGdlb21fYm94cGxvdCgpCgpnZ3Bsb3QoZGF0YSA9IGJtaWYsIGFlcyh5PSB3YWkuZ2ksIHggPSBibWlfb3ZlcikpICsKICBnZW9tX2JveHBsb3QoKQpgYGAKCiMjIyBNdWx0aXBsZSBSZWdyZXNzaW9uIE1vZGVsCldoZW4gY3JlYXRpbmcgb3VyIE11bHRpcGxlIFJlZ3Jlc3Npb24gTW9kZWwsIHdlIHJlbW92ZWQgImhlaWdodCIgYW5kICJ3ZWlnaHQiIHZhcmlhYmxlcyBmcm9tIGNvbnNpZGVyYXRpb24gYmVjYXVzZSBpdCBhcHBlYXJlZCBjb3VudGVyaW50dWl0aXZlIHRvIGluY2x1ZGUgdGhlc2UgZmFjdG9ycy0td2hpY2ggYXJlIHVzZWQgdG8gY2FsY3VsYXRlIEJNSS0tZnJvbSBjb25zaWRlcmF0aW9uIG9mIHdoYXQgZmFjdG9ycyBpbmZ1ZW5jZSBCTUkuIEZyb20gdGhlcmUsIHdlIGZvdW5kIHRoYXQgYmFja3dhcmRzIG1vZGVsbGluZyBjYXVzZWQgdGhlIGFkanVzdGVkIFItc3F1YXJlZCB2YWx1ZSB0byBkZWNyZWFzZSB3aGVuZXZlciBhIHZhcmlhYmxlIHdhcyByZW1vdmVkLCBzbyB0aGUgYmVzdCBsaW5lYXIgcmVncmVzc2lvbiBtb2RlbCBpbmNsdWRlcyBhbGwgb2YgdGhlIG9yaWdpbmFsIHZhcmlhYmxlcy4gQ29uY2VwdHVhbGx5LCB0aGlzIG1ha2VzIHNlbnNlIGJjIEJNSSBpcyBhIHZlcnkgc2Vuc2l0aXZlIGNhbGN1bGF0aW9uIGFuZCBhbnkgYWRkaXRpb25hbCBkYXRhIG1ha2VzIGl0IG1vcmUgYWNjdXJhdGUKYGBge3J9CmRhdGEoYm1pZikKbV9hbGwgPC0gbG0oYm1pIH4gYmlhLmRpICsgYmlpLmRpICsgYml0LmRpICsgY2hlLmRpICsgZWxiLmRpICsgd3JpLmRpICsga25lLmRpICsgYW5rLmRpICsgc2hvLmdpICsgY2hlLmdpICsgd2FpLmdpICsgbmF2LmdpICsgaGlwLmdpICsgdGhpLmdpICsgYmljLmdpICsgZm9yLmdpICsga25lLmdpICsgY2FsLmdpICsgd3JpLmdpICsgY2hlLmRlLCBkYXRhID0gYm1pZikKc3VtbWFyeShtX2FsbCkKbV9hbGxbWyJjb2VmZmljaWVudHMiXV0KYGBgCgojIyMgTGluZWFyaXR5IFRlc3RzIGZvciBvdXIgTW9kZWw6CkFsbCBvZiB0aGUgY29uZGl0aW9ucyBmb3IgbGluZWFyIHJlZ3Jlc3Npb24gYXJlIG1ldDoKTGluZWFyaXR5OiBUaGUgcmVzaWR1YWxzIHBsb3RzIGRvZXMgbm90IHZhcnkgbXVjaCBhcm91bmQgdGhlIGxpbmUgYW5kIHRoZXJlZm9yZSB0aGUgY29ycmVsYXRpb24gbGluZSB0aGF0IHdlIGNyZWF0ZWQgaXMgbm90IGEgZ29vZCBwcmVkaWN0b3Igb2YgYWN0dWFsIEJNSSBhbmQgbV9mdWxsIGlzIGEgZ29vZCBwcmVkaWN0b3Igb2YgQk1JCk5lYXJseSBOb3JtYWwgUmVzaWR1YWxzOiBUaGUgbm9ybWFsIHByb2JhYmlsaXR5IHBsb3QgaXMgY2VudGVyZWQgYXJvdW5kIDAKQ29uc3RhbnQgVmFyaWFiaWxpdHk6IGJhc2VkIG9uIHRoZSByZXNpZHVhbHMgdnMuIGZpdHRlZCBwbG90LCB0aGUgY29uc3RhbnQgdmFyaWFiaWxpdHkgY29uZGl0aW9uIGFwcGVhcnMgdG8gYmUgbWV0IGJlY2F1c2UgYXMgeCB2YXJpZXMgdGhlIHZhcmlhYmlsaXR5IGRvZXMgbm90IHZhcnkgc3Vic3RhbnRpYWxseS4KRWFjaCB2YXJpYWJsZSBpcyBsaW5lYXJseSByZWxhdGVkIHRvIHRoZSBvdXRjb21lOiB0aGUgcmVzaWR1YWxzIGFyZSByYW5kb21seSBzY2F0dGVyZWQgYXJvdW5kIDAsIHNvIHRoaXMgY29uZGl0aW9uIGlzIG1ldApgYGB7cn0KIHFwbG90KHggPSAuZml0dGVkLCB5ID0gLnJlc2lkLCBkYXRhID0gbV9hbGwpICsKICAgICBnZW9tX2hsaW5lKHlpbnRlcmNlcHQgPSAwLCBsaW5ldHlwZSA9ICJkYXNoZWQiKSArCiAgICAgeGxhYigiRml0dGVkIHZhbHVlcyIpICsKICAgICB5bGFiKCJSZXNpZHVhbHMiKQpxcGxvdChzYW1wbGUgPSAucmVzaWQsIGRhdGEgPSBtX2FsbCwgc3RhdCA9ICJxcSIpCmBgYApodHRwczovL2xhYnMtYXotMDIub2l0LmR1a2UuZWR1OjMwNjIzL2dyYXBoaWNzL3Bsb3QucG5nP3dpZHRoPTQ0NiZoZWlnaHQ9Mjk5JnJhbmRvbWl6ZXI9MTc3MjI1NDA1MgoKaHR0cHM6Ly9sYWJzLWF6LTAyLm9pdC5kdWtlLmVkdTozMDYyMy9ncmFwaGljcy9wbG90LnBuZz93aWR0aD00NDYmaGVpZ2h0PTI5OSZyYW5kb21pemVyPS0yMTI4OTEzNTQ3IAoKIyMjIE1vZGVsIFNlbGVjdGlvbiBhbmQgUGVyZm9ybWFuY2UKUlN0dWRpbyBjcmVhdGVkIHRoZSBzYW1lIGJlc3QtZml0IG1vZGVsIGFzIHdlIGRpZDogb25lIHRoYXQgaW5jbHVkZXMgYWxsIHZhcmlhYmxlcyB0aGF0IGNvbnRyaWJ1dGUgdG8gQk1JIChleGNsdWRpbmcgd2VpZ2h0IGFuZCBoZWlnaHQpCmBgYHtyfQpibWlmLmRhdGEgPC0gYm1pZiAlPiUgc2VsZWN0KC1jKCd3Z3QnLCdoZ3QnLCAnaHRtZXRlcnMnLCAnYm1pX292ZXInKSkKbi5vYnMgPC0gZGltKGJtaWYuZGF0YSlbMV0KCnRyYWluLmluZGV4IDwtIHNhbXBsZSgxOm4ub2JzLGZsb29yKC44Km4ub2JzKSxyZXBsYWNlPUZBTFNFKQpibWlmLnRyYWluIDwtIGJtaWYuZGF0YVt0cmFpbi5pbmRleCxdCmJtaWYudGVzdCA8LSBibWlmLmRhdGFbLXRyYWluLmluZGV4LF0gCgpzZWwuZGF0YSA8LSBibWlmLnRyYWluCm0wIDwtIGxtKGJtaSB+IC4sIGRhdGE9c2VsLmRhdGEpCnN1bW1hcnkobTApCgpzaW5nbGUuc3RlcC5iYWNrd2FyZHMgPC0gZnVuY3Rpb24oZGF0YSxyZXNwb25zZSl7CnJlc3AuaW5keCA8LSB3aGljaChuYW1lcyhkYXRhKT09cmVzcG9uc2UpCnkgPC0gZGF0YVsscmVzcC5pbmR4XQpYIDwtIGRhdGFbLC1yZXNwLmluZHhdCm4ucHJlZCA8LSBkaW0oWClbMl0KaWYobi5wcmVkID4gMSl7CiAgZm9yKGkgaW4gMTpuLnByZWQpewogICAgICBwcmludChwYXN0ZTAoIlZhcmlhYmxlICIsIG5hbWVzKFgpW2ldLCIgcmVtb3ZlZDogQWRqdXN0ZWQgUi1zcXVhcmVkID0gICIsCiAgICAgICAgICByb3VuZChzdW1tYXJ5KGxtKHl+LixkYXRhPWFzLmRhdGEuZnJhbWUoWFssLWldKSkpJGFkai5yLnNxdWFyZWQsNSkpKQogICAgICB9Cn0KZWxzZXsKICBwcmludCgiTW9kZWwgb25seSBjb250YWlucyBvbmUgdmFyaWFibGUuIikKfQp9CmBgYAoKIyMjIFByZWRpY3Rpb24gVXNpbmcgQmVzdCBGaXQgTW9kZWwKT3VyIFJvb3QgTWVhbiBTcXVhcmUgRXJyb3IgKFJNU0UpIGlzIDEuMDkxNzkzLCB3aGljaCBtZWFucyB0aGF0IG91ciBtb2RlbCB3b3VsZCBjYWxjdWxhdGUgYSBuZXcgQk1JIHRoYXQgd291bGQgYmUgb2ZmIGJ5IDEuMTIwMDMyIG9uIGF2ZXJhZ2UuCk91ciBtb2RlbCBjYWxjdWxhdGVkIHRoZSBCTUkgZm9yIGEgcmFuZG9tIHNhbXBsZSBhcyAyNC45NzAyMiwgd2l0aCBhIGxvd2VyIGxpbWl0IG9mIDIyLjYzMzk3IGFuZCBhbiB1cHBlciBsaW1pdCBvZiAyNy4zMDY0Ni4gT3VyIG1vZGVsIG92ZXJlc3RpbWF0ZWQgdGhlIHJlYWwgQk1JIGJ5IDMuODE4NTEgYW5kIG91ciBwcmVkaWN0aW9uIGludGVydmFsIGRpZCBub3QgaW5jbHVkZSB0aGUgcmVhbCB2YWx1ZSBvZiAyMS4wOTE4LiAKYGBge3J9Cm1vZGVsLmJlc3QgPC0gbG0oYm1pIH4gLiwgZGF0YT1ibWlmLnRyYWluKQpzdW1tYXJ5KG1vZGVsLmJlc3QpCgpxcGxvdCh4ID0gLmZpdHRlZCwgeSA9IC5yZXNpZCwgZGF0YSA9IG1vZGVsLmJlc3QpICsKICAgICAgICAgZ2VvbV9obGluZSh5aW50ZXJjZXB0ID0gMCwgbGluZXR5cGUgPSAiZGFzaGVkIikgKwogICAgICAgICB4bGFiKCJGaXR0ZWQgdmFsdWVzIikgKwogICAgICAgIHlsYWIoIlJlc2lkdWFscyIpCnFwbG90KHNhbXBsZSA9IC5yZXNpZCwgZGF0YSA9IG1vZGVsLmJlc3QsIHN0YXQgPSAicXEiKQoKcHJlZGljdGlvbnMudGVzdCA8LSBwcmVkaWN0KG1vZGVsLmJlc3QsYm1pZi50ZXN0KQogeSA8LSBibWlmLnRlc3QkYm1pCiBtc2UgPC0gbWVhbigoeS1wcmVkaWN0aW9ucy50ZXN0KV4yKQogKHJtc2UgPC0gc3FydChtc2UpKQogCmJtaWYudGVzdFsyLF0KcHJlZGljdChtb2RlbC5iZXN0LGJtaWYudGVzdFsyLF0saW50ZXJ2YWw9J3ByZWRpY3Rpb24nKQpgYGAKUk1TRTogMS4wOTE3OTMKCkNvbmRpdGlvbnMgQ2hlY2tzOiAKaHR0cHM6Ly9sYWJzLWF6LTAyLm9pdC5kdWtlLmVkdTozMDYyMy9ncmFwaGljcy9wbG90LnBuZz93aWR0aD00NDMmaGVpZ2h0PTI5OSZyYW5kb21pemVyPTExNTE3Njg2NzYKaHR0cHM6Ly9sYWJzLWF6LTAyLm9pdC5kdWtlLmVkdTozMDYyMy9ncmFwaGljcy9wbG90LnBuZz93aWR0aD00NDMmaGVpZ2h0PTI5OSZyYW5kb21pemVyPS0xNzgxMjk5NDA1IAoKVGVzdDoKICBmaXQgICAgICBsd3IgICAgICB1cHIKMSAyNC45NzAyMiAyMi42MzM5NyAyNy4zMDY0Ngo+IAo=