Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 1 | # Copyright 2017 Open Networking Foundation (ONF) |
| 2 | # |
| 3 | # Please refer questions to either the onos test mailing list at <onos-test@onosproject.org>, |
| 4 | # the System Testing Plans and Results wiki page at <https://wiki.onosproject.org/x/voMg>, |
| 5 | # or the System Testing Guide page at <https://wiki.onosproject.org/x/WYQg> |
| 6 | # |
| 7 | # TestON is free software: you can redistribute it and/or modify |
| 8 | # it under the terms of the GNU General Public License as published by |
| 9 | # the Free Software Foundation, either version 2 of the License, or |
| 10 | # (at your option) any later version. |
| 11 | # |
| 12 | # TestON is distributed in the hope that it will be useful, |
| 13 | # but WITHOUT ANY WARRANTY; without even the implied warranty of |
| 14 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
| 15 | # GNU General Public License for more details. |
| 16 | # |
| 17 | # You should have received a copy of the GNU General Public License |
| 18 | # along with TestON. If not, see <http://www.gnu.org/licenses/>. |
| 19 | # |
| 20 | # If you have any questions, or if you don't understand R, |
Jeremy Ronquillo | b626884 | 2017-10-03 13:02:58 -0700 | [diff] [blame] | 21 | # please contact Jeremy Ronquillo: j_ronquillo@u.pacific.edu |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 22 | |
| 23 | # ********************************************************** |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 24 | # STEP 1: Data management. |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 25 | # ********************************************************** |
| 26 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 27 | print( "**********************************************************" ) |
| 28 | print( "STEP 1: Data management." ) |
| 29 | print( "**********************************************************" ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 30 | |
Jeremy Ronquillo | b626884 | 2017-10-03 13:02:58 -0700 | [diff] [blame] | 31 | # Command line arguments are read. |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 32 | print( "Reading commmand-line args." ) |
| 33 | args <- commandArgs( trailingOnly=TRUE ) |
| 34 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 35 | # ---------------- |
| 36 | # Import Libraries |
| 37 | # ---------------- |
| 38 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 39 | print( "Importing libraries." ) |
| 40 | library( ggplot2 ) |
| 41 | library( reshape2 ) |
| 42 | library( RPostgreSQL ) # For databases |
| 43 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 44 | # ------------------- |
| 45 | # Check CLI Arguments |
| 46 | # ------------------- |
| 47 | |
| 48 | print( "Verifying CLI args." ) |
| 49 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 50 | if ( is.na( args[ 9 ] ) ){ |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 51 | |
| 52 | print( paste( "Usage: Rscript SCPFflowTp1g.R", |
| 53 | "<has-flow-obj>", |
| 54 | "<database-host>", |
| 55 | "<database-port>", |
| 56 | "<database-user-id>", |
| 57 | "<database-password>", |
| 58 | "<test-name>", |
| 59 | "<branch-name>", |
| 60 | "<has-neighbors>", |
| 61 | "<directory-to-save-graphs>", |
| 62 | sep=" " ) ) |
| 63 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 64 | q() # basically exit(), but in R |
| 65 | } |
| 66 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 67 | # ----------------- |
| 68 | # Create File Names |
| 69 | # ----------------- |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 70 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 71 | print( "Creating filenames and title of graph." ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 72 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 73 | chartTitle <- "Flow Throughput Test" |
| 74 | fileNeighborsModifier <- "no" |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 75 | commandNeighborModifier <- "" |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 76 | fileFlowObjModifier <- "" |
| 77 | sqlFlowObjModifier <- "" |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 78 | if ( args[ 1 ] == 'y' ){ |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 79 | fileFlowObjModifier <- "_flowObj" |
| 80 | sqlFlowObjModifier <- "_fobj" |
| 81 | chartTitle <- paste( chartTitle, " with Flow Objectives", sep="" ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 82 | } |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 83 | |
| 84 | chartTitle <- paste( chartTitle, "\nNeighbors =", sep="" ) |
| 85 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 86 | if ( args[ 8 ] == 'y' ){ |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 87 | fileNeighborsModifier <- "all" |
Jeremy Ronquillo | 2d2649d | 2017-09-14 12:53:06 -0700 | [diff] [blame] | 88 | commandNeighborModifier <- "scale=1 OR NOT " |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 89 | chartTitle <- paste( chartTitle, "Cluster Size - 1" ) |
| 90 | } else { |
| 91 | chartTitle <- paste( chartTitle, "0" ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 92 | } |
| 93 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 94 | errBarOutputFile <- paste( args[ 9 ], |
| 95 | args[ 6 ], |
| 96 | "_", |
| 97 | args[ 7 ], |
| 98 | "_", |
| 99 | fileNeighborsModifier, |
| 100 | "-neighbors", |
| 101 | fileFlowObjModifier, |
| 102 | "_graph.jpg", |
| 103 | sep="" ) |
| 104 | # ------------------ |
| 105 | # SQL Initialization |
| 106 | # ------------------ |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 107 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 108 | print( "Initializing SQL" ) |
| 109 | |
| 110 | con <- dbConnect( dbDriver( "PostgreSQL" ), |
| 111 | dbname = "onostest", |
| 112 | host = args[ 2 ], |
| 113 | port = strtoi( args[ 3 ] ), |
| 114 | user = args[ 4 ], |
| 115 | password = args[ 5 ] ) |
| 116 | |
| 117 | # --------------------------- |
| 118 | # Flow Throughput SQL Command |
| 119 | # --------------------------- |
| 120 | |
| 121 | print( "Generating Flow Throughput SQL command." ) |
| 122 | |
| 123 | command <- paste( "SELECT scale, avg( avg ), avg( std ) FROM flow_tp", |
| 124 | sqlFlowObjModifier, |
| 125 | "_tests WHERE (", |
| 126 | commandNeighborModifier, |
| 127 | "neighbors = 0 ) AND branch = '", |
| 128 | args[ 7 ], |
| 129 | "' AND date IN ( SELECT max( date ) FROM flow_tp", |
| 130 | sqlFlowObjModifier, |
| 131 | "_tests WHERE branch='", |
| 132 | args[ 7 ], |
| 133 | "' ) GROUP BY scale ORDER BY scale", |
| 134 | sep="" ) |
| 135 | |
| 136 | print( "Sending SQL command:" ) |
| 137 | print( command ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 138 | |
| 139 | fileData <- dbGetQuery( con, command ) |
| 140 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 141 | # ********************************************************** |
| 142 | # STEP 2: Organize data. |
| 143 | # ********************************************************** |
| 144 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 145 | print( "**********************************************************" ) |
| 146 | print( "STEP 2: Organize Data." ) |
| 147 | print( "**********************************************************" ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 148 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 149 | # ------------ |
| 150 | # Data Sorting |
| 151 | # ------------ |
| 152 | |
| 153 | print( "Sorting data for Flow Throughput." ) |
| 154 | |
| 155 | colnames( fileData ) <- c( "scale", |
| 156 | "avg", |
| 157 | "std" ) |
| 158 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 159 | avgs <- c( fileData[ 'avg' ] ) |
| 160 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 161 | |
| 162 | # ---------------------------- |
| 163 | # Flow TP Construct Data Frame |
| 164 | # ---------------------------- |
| 165 | |
| 166 | print( "Constructing Flow TP data frame." ) |
| 167 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 168 | dataFrame <- melt( avgs ) # This is where reshape2 comes in. Avgs list is converted to data frame |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 169 | dataFrame$scale <- fileData$scale # Add node scaling to the data frame. |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 170 | dataFrame$std <- fileData$std |
| 171 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 172 | colnames( dataFrame ) <- c( "throughput", |
| 173 | "type", |
| 174 | "scale", |
| 175 | "std" ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 176 | |
Jeremy Ronquillo | 2d2649d | 2017-09-14 12:53:06 -0700 | [diff] [blame] | 177 | dataFrame <- na.omit( dataFrame ) # Omit any data that doesn't exist |
| 178 | |
| 179 | print( "Data Frame Results:" ) |
| 180 | print( dataFrame ) |
| 181 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 182 | # ********************************************************** |
| 183 | # STEP 3: Generate graphs. |
| 184 | # ********************************************************** |
| 185 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 186 | print( "**********************************************************" ) |
| 187 | print( "STEP 3: Generate Graph." ) |
| 188 | print( "**********************************************************" ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 189 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 190 | # ------------------ |
| 191 | # Generate Main Plot |
| 192 | # ------------------ |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 193 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 194 | print( "Generating main plot." ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 195 | # Create the primary plot here. |
| 196 | # ggplot contains the following arguments: |
| 197 | # - data: the data frame that the graph will be based off of |
| 198 | # - aes: the asthetics of the graph which require: |
| 199 | # - x: x-axis values (usually node scaling) |
| 200 | # - y: y-axis values (usually time in milliseconds) |
| 201 | # - fill: the category of the colored side-by-side bars (usually type) |
Jeremy Ronquillo | 2d2649d | 2017-09-14 12:53:06 -0700 | [diff] [blame] | 202 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 203 | mainPlot <- ggplot( data = dataFrame, aes( x = scale, |
| 204 | y = throughput, |
| 205 | ymin = throughput, |
| 206 | ymax = throughput + std, |
| 207 | fill = type ) ) |
| 208 | # ------------------------------ |
| 209 | # Fundamental Variables Assigned |
| 210 | # ------------------------------ |
Jeremy Ronquillo | 2d2649d | 2017-09-14 12:53:06 -0700 | [diff] [blame] | 211 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 212 | print( "Generating fundamental graph data." ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 213 | |
| 214 | # Formatting the plot |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 215 | theme_set( theme_grey( base_size = 22 ) ) # set the default text size of the graph. |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 216 | width <- 0.7 # Width of the bars. |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 217 | xScaleConfig <- scale_x_continuous( breaks = dataFrame$scale, |
| 218 | label = dataFrame$scale ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 219 | xLabel <- xlab( "Scale" ) |
Jeremy Ronquillo | 1bdaae5 | 2017-09-22 11:39:48 -0700 | [diff] [blame] | 220 | yLabel <- ylab( "Throughput (,000 Flows/sec)" ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 221 | fillLabel <- labs( fill="Type" ) |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 222 | imageWidth <- 15 |
| 223 | imageHeight <- 10 |
| 224 | imageDPI <- 200 |
| 225 | errorBarColor <- rgb( 140, 140, 140, maxColorValue=255 ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 226 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 227 | theme <- theme( plot.title = element_text( hjust = 0.5, |
| 228 | size = 32, |
| 229 | face = 'bold' ) ) |
| 230 | title <- ggtitle( chartTitle ) |
Jeremy Ronquillo | 4363d09 | 2017-10-13 13:28:47 -0700 | [diff] [blame] | 231 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 232 | # Store plot configurations as 1 variable |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 233 | fundamentalGraphData <- mainPlot + |
| 234 | xScaleConfig + |
| 235 | xLabel + |
| 236 | yLabel + |
| 237 | fillLabel + |
| 238 | theme + |
| 239 | title |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 240 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 241 | # --------------------------- |
| 242 | # Generating Bar Graph Format |
| 243 | # --------------------------- |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 244 | |
| 245 | # Create the stacked bar graph with error bars. |
| 246 | # geom_bar contains: |
| 247 | # - stat: data formatting (usually "identity") |
| 248 | # - width: the width of the bar types (declared above) |
| 249 | # geom_errorbar contains similar arguments as geom_bar. |
| 250 | print( "Generating bar graph with error bars." ) |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 251 | barGraphFormat <- geom_bar( stat = "identity", |
| 252 | width = width, |
| 253 | fill = "#FFAA3C" ) |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 254 | |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 255 | errorBarFormat <- geom_errorbar( width = width, |
| 256 | position = position_dodge(), |
| 257 | color = errorBarColor ) |
| 258 | |
| 259 | values <- geom_text( aes( x = dataFrame$scale, |
| 260 | y = dataFrame$throughput + 0.03 * max( dataFrame$throughput ), |
| 261 | label = format( dataFrame$throughput, |
| 262 | digits=3, |
| 263 | big.mark = ",", |
| 264 | scientific = FALSE ) ), |
| 265 | size = 7.0, |
| 266 | fontface = "bold" ) |
| 267 | |
| 268 | result <- fundamentalGraphData + |
| 269 | barGraphFormat + |
| 270 | errorBarFormat + |
| 271 | values |
| 272 | |
| 273 | # ----------------------- |
| 274 | # Exporting Graph to File |
| 275 | # ----------------------- |
| 276 | |
Jeremy Ronquillo | 6df8781 | 2017-08-28 16:17:36 +0000 | [diff] [blame] | 277 | print( paste( "Saving bar chart with error bars to", errBarOutputFile ) ) |
Jeremy Ronquillo | 7673f80 | 2017-10-30 09:42:44 -0700 | [diff] [blame^] | 278 | |
| 279 | ggsave( errBarOutputFile, |
| 280 | width = imageWidth, |
| 281 | height = imageHeight, |
| 282 | dpi = imageDPI ) |
| 283 | |
| 284 | print( paste( "[SUCCESS] Successfully wrote bar chart with error bars out to", errBarOutputFile ) ) |