Mercurial > repos > siwaa > redelac_stics_t

comparison redelacSticsTool.R @ 18:8293d67cd1a9 draft

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 b972382172a3923486a799d46efe363d21e8ea1d"
author siwaa
date Thu, 09 Jul 2026 07:09:47 +0000
parents fa4e7fd1d3f7
children 8848db582840
comparison
equal deleted inserted replaced
17:fa4e7fd1d3f7 18:8293d67cd1a9
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 afterGenTxtTime <- Sys.time()
141 durGenT <-
142 lubridate::as.duration(lubridate::interval(beforeGenTxtTime, afterGenTxtTime))
143
144 endTime <- Sys.time() 140 endTime <- Sys.time()
145 141
146 dur <- 142 dur <-
147 lubridate::as.period(lubridate::as.duration(lubridate::interval(startTime, endTime))) 143 lubridate::as.period(lubridate::as.duration(lubridate::interval(startTime, endTime)))
148 durGenXP <- lubridate::as.period(durGenX) 144 durGenXP <- lubridate::as.period(durGenX)
180 ) 176 )
181 print( 177 print(
182 "===================================================================================" 178 "==================================================================================="
183 ) 179 )
184 180
185
186 # 181 #
187 # Simulation part 182 # Simulation part
188 # 183 #
184
185 beforeSimsTime <- Sys.time()
189 186
190 # List of successive USMs 187 # List of successive USMs
191 successive_usms_file <- ("successionPlan.csv") 188 successive_usms_file <- ("successionPlan.csv")
192 successive_usms_tab <- 189 successive_usms_tab <-
193 read.csv(file = file.path(workspace, successive_usms_file), 190 read.csv(file = file.path(workspace, successive_usms_file),
252 cat(paste0("--> succession : ", successive_usms_tab_split[[j]][i, 1], " launched\n")) 249 cat(paste0("--> succession : ", successive_usms_tab_split[[j]][i, 1], " launched\n"))
253 250
254 return(stics_wrapper(model_options = sim_options, situation = list_usms)) 251 return(stics_wrapper(model_options = sim_options, situation = list_usms))
255 252
256 } 253 }
257 254
258 for (j in 1:length(casesUpcSol)) { 255 afterSimsTime <- Sys.time()
256 durSimsT <-
257 lubridate::as.duration(lubridate::interval(beforeSimsTime, afterSimsTime))
258 durSimsTP <- lubridate::as.period(durSimsT)
259
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 #
272 # Output part
273 #
274
275 beforeOutputsTime <- Sys.time()
276
277 #for (j in 1:length(casesUpcSol)) {
278
279 foreach(
280 j = 1:length(casesUpcSol),
281 .options.future = list(packages =c("foreach", "tidyr", "arrow")
282 ) %dordopar% {
283
259 splited <- strsplit(casesUpcSol[j], split = "-")[[1]] 284 splited <- strsplit(casesUpcSol[j], split = "-")[[1]]
260 285
261 upc <- splited[1] 286 upc <- splited[1]
262 287
263 sol <- splited[2] 288 sol <- splited[2]
522 ".parquet" 547 ".parquet"
523 548
524 )) 549 ))
525 } 550 }
526 } 551 }
552
553 afterOutputsTime <- Sys.time()
554 durOutT <-
555 lubridate::as.duration(lubridate::interval(beforeOutputsTime, afterOutputsTime))
556 durOutTP <- lubridate::as.period(durOutT)
557
558 sprintf(
559 'outputs generating duration: %g jour(s) %g heure(s) %g minute(s) %g seconde(s)',
560 lubridate::day(durOutTP),
561 lubridate::hour(durOutTP),
562 lubridate::minute(durOutTP),
563 lubridate::second(durOutTP)
564 )
565 print(
566 "==================================================================================="
567 )