Mercurial > repos > siwaa > redelac_stics_t

comparison redelacSticsTool.R @ 19:8848db582840 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://forgemia.inra.fr/redelac commit 9f38541b8f84bf16077819d47c96dbc69e3d1675"
author siwaa
date Thu, 09 Jul 2026 14:30:55 +0000
parents 8293d67cd1a9
children
comparison
equal deleted inserted replaced
18:8293d67cd1a9 19:8848db582840
135 135
136 afterGenTxtTime <- Sys.time() 136 afterGenTxtTime <- Sys.time()
137 durGenT <- 137 durGenT <-
138 lubridate::as.duration(lubridate::interval(beforeGenTxtTime, afterGenTxtTime)) 138 lubridate::as.duration(lubridate::interval(beforeGenTxtTime, afterGenTxtTime))
139 139
140 endTime <- Sys.time()
141
142 dur <-
143 lubridate::as.period(lubridate::as.duration(lubridate::interval(startTime, endTime)))
144 durGenXP <- lubridate::as.period(durGenX)
145 durGenTP <- lubridate::as.period(durGenT)
146
147 print("REDELACSticsInputGenerator Report")
148 print("=================================")
149 print(paste("number of usms: ", nrow(usms_param)))
150 print(paste("redelac optimization:", redelac))
151 print(paste("process in parallel:", parallel))
152 print(paste("nbslots:", nbSlots))
153
154 sprintf(
155 'overall duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
156 lubridate::day(dur),
157 lubridate::hour(dur),
158 lubridate::minute(dur),
159 lubridate::second(dur)
160 )
161
162 sprintf(
163 'inside gen_usms_xml duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
164 lubridate::day(durGenXP),
165 lubridate::hour(durGenXP),
166 lubridate::minute(durGenXP),
167 lubridate::second(durGenXP)
168 )
169
170 sprintf(
171 'inside gen_usms_xml2txt duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
172 lubridate::day(durGenTP),
173 lubridate::hour(durGenTP),
174 lubridate::minute(durGenTP),
175 lubridate::second(durGenTP)
176 )
177 print(
178 "==================================================================================="
179 )
180
181 # 140 #
182 # Simulation part 141 # Simulation part
183 # 142 #
184 143
185 beforeSimsTime <- Sys.time() 144 beforeSimsTime <- Sys.time()
249 cat(paste0("--> succession : ", successive_usms_tab_split[[j]][i, 1], " launched\n")) 208 cat(paste0("--> succession : ", successive_usms_tab_split[[j]][i, 1], " launched\n"))
250 209
251 return(stics_wrapper(model_options = sim_options, situation = list_usms)) 210 return(stics_wrapper(model_options = sim_options, situation = list_usms))
252 211
253 } 212 }
254 213
255 afterSimsTime <- Sys.time() 214 afterSimsTime <- Sys.time()
256 durSimsT <- 215 durSimsT <-
257 lubridate::as.duration(lubridate::interval(beforeSimsTime, afterSimsTime)) 216 lubridate::as.duration(lubridate::interval(beforeSimsTime, afterSimsTime))
258 durSimsTP <- lubridate::as.period(durSimsT) 217 durSimsTP <- lubridate::as.period(durSimsT)
259 218
260 sprintf(
261 'Simulations duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
262 lubridate::day(durSimsTP),
263 lubridate::hour(durSimsTP),
264 lubridate::minute(durSimsTP),
265 lubridate::second(durSimsTP)
266 )
267 print(
268 "==================================================================================="
269 )
270
271 # 219 #
272 # Output part 220 # Output part
273 # 221 #
274 222
275 beforeOutputsTime <- Sys.time() 223 beforeOutputsTime <- Sys.time()
276 224
277 #for (j in 1:length(casesUpcSol)) { 225 for (j in 1:length(casesUpcSol)) {
278 226
279 foreach( 227 successive_usms_tab_split[[j]]$Treatment <- NULL
228 successive_usms_tab_split[[j]]$Number_USM <- NULL
229
230 }
231
232 foreach(
280 j = 1:length(casesUpcSol), 233 j = 1:length(casesUpcSol),
281 .options.future = list(packages =c("foreach", "tidyr", "arrow") 234 .options.future = list(packages =c("foreach", "tidyr", "arrow"))
282 ) %dordopar% { 235 ) %dordopar% {
283 236
284 splited <- strsplit(casesUpcSol[j], split = "-")[[1]] 237 splited <- strsplit(casesUpcSol[j], split = "-")[[1]]
285 238
286 upc <- splited[1] 239 upc <- splited[1]
399 ".parquet" 352 ".parquet"
400 353
401 )) 354 ))
402 } 355 }
403 356
404 successive_usms_tab_split[[j]]$Treatment <- NULL 357 #successive_usms_tab_split[[j]]$Treatment <- NULL
405 successive_usms_tab_split[[j]]$Number_USM <- NULL 358 #successive_usms_tab_split[[j]]$Number_USM <- NULL
406 359
407 list_usms <- 360 list_usms <-
408 unlist(unname(as.list( 361 unlist(unname(as.list(
409 transpose(successive_usms_tab_split[[j]]) 362 transpose(successive_usms_tab_split[[j]])
410 ))) 363 )))
549 )) 502 ))
550 } 503 }
551 } 504 }
552 505
553 afterOutputsTime <- Sys.time() 506 afterOutputsTime <- Sys.time()
507
554 durOutT <- 508 durOutT <-
555 lubridate::as.duration(lubridate::interval(beforeOutputsTime, afterOutputsTime)) 509 lubridate::as.duration(lubridate::interval(beforeOutputsTime, afterOutputsTime))
510
511
512 durGenXP <- lubridate::as.period(durGenX)
513 durGenTP <- lubridate::as.period(durGenT)
556 durOutTP <- lubridate::as.period(durOutT) 514 durOutTP <- lubridate::as.period(durOutT)
515
516 endTime <- Sys.time()
517
518 dur <-
519 lubridate::as.period(lubridate::as.duration(lubridate::interval(startTime, endTime)))
520
521 print("REDELACSticsInputGenerator Report")
522 print("=================================")
523 print(paste("number of usms: ", nrow(usms_param)))
524 print(paste("redelac optimization:", redelac))
525 print(paste("process in parallel:", parallel))
526 print(paste("nbslots:", nbSlots))
527
528 sprintf(
529 'overall duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
530 lubridate::day(dur),
531 lubridate::hour(dur),
532 lubridate::minute(dur),
533 lubridate::second(dur)
534 )
535 print(
536 "==================================================================================="
537 )
538
539 sprintf(
540 'inside gen_usms_xml duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
541 lubridate::day(durGenXP),
542 lubridate::hour(durGenXP),
543 lubridate::minute(durGenXP),
544 lubridate::second(durGenXP)
545 )
546 print(
547 "==================================================================================="
548 )
549
550 sprintf(
551 'inside gen_usms_xml2txt duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
552 lubridate::day(durGenTP),
553 lubridate::hour(durGenTP),
554 lubridate::minute(durGenTP),
555 lubridate::second(durGenTP)
556 )
557 print(
558 "==================================================================================="
559 )
560
561 sprintf(
562 'Simulations duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
563 lubridate::day(durSimsTP),
564 lubridate::hour(durSimsTP),
565 lubridate::minute(durSimsTP),
566 lubridate::second(durSimsTP)
567 )
568 print(
569 "==================================================================================="
570 )
557 571
558 sprintf( 572 sprintf(
559 'outputs generating duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)', 573 'outputs generating duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
560 lubridate::day(durOutTP), 574 lubridate::day(durOutTP),
561 lubridate::hour(durOutTP), 575 lubridate::hour(durOutTP),