YAC 3.18.0
Yet Another Coupler
Loading...
Searching...
No Matches
test_events.c
Go to the documentation of this file.
1// Copyright (c) 2024 The YAC Authors
2//
3// SPDX-License-Identifier: BSD-3-Clause
4
5#include <stdlib.h>
6#include <stdio.h>
7#include <string.h>
8#include "tests.h"
9#include "event.h"
10
16int main (void) {
17
18 // basic setup
19 initCalendar(PROLEPTIC_GREGORIAN);
20
21 {
22 struct {
23 enum yac_time_unit_type time_unit;
24 char const * value;
25 } tests[] =
26 {{.time_unit = C_MILLISECOND,
27 .value = "3600000"},
28 {.time_unit = C_SECOND,
29 .value = "3600"},
30 {.time_unit = C_MINUTE,
31 .value = "60"},
32 {.time_unit = C_HOUR,
33 .value = "1"},
34 {.time_unit = C_ISO_FORMAT,
35 .value = "PT01H"}};
36 enum {NUM_TESTS = sizeof(tests) / sizeof(tests[0])};
37 for (size_t i = 0; i < NUM_TESTS; ++i)
38 if (strcmp(
39 yac_time_to_ISO(
40 tests[i].value, tests[i].time_unit), "PT01H"))
41 PUT_ERR("ERROR in yac_time_to_ISO");
42 }
43
44 {
45 struct event * event = yac_event_new ();
46
47 char const * model_time_step = "5000"; // 5s per time step
48 char const * coupling_time_step = "15000"; // 15s per coupling time step
49
50 enum yac_time_unit_type time_unit = C_MILLISECOND;
51 int lag = 0;
52 enum yac_reduction_type time_operation = TIME_AVERAGE;
53
54 // 60s runtime
55 char const * start_date = "1850-01-01T00:00:00";
56 char const * stop_date = "1850-01-01T00:01:00";
57
58 char const * model_time_step_iso =
59 strdup(yac_time_to_ISO(model_time_step, time_unit));
60 char const * coupling_time_step_iso =
61 strdup(yac_time_to_ISO(coupling_time_step, time_unit));
62 yac_event_add(
63 event, model_time_step_iso, coupling_time_step_iso,
64 lag, time_operation, start_date, stop_date);
65
66 char time_step_buffer[MAX_TIMEDELTA_STR_LEN];
67 if (strcmp(
68 yac_get_event_coupling_timestep(event, time_step_buffer),
69 coupling_time_step_iso))
70 PUT_ERR("ERROR in yac_get_event_coupling_timestep");
71 if (strcmp(
72 yac_get_event_model_timestep(event, time_step_buffer),
73 model_time_step_iso))
74 PUT_ERR("ERROR in yac_get_event_model_timestep");
75 if (yac_get_event_lag(event) != lag)
76 PUT_ERR("ERROR in yac_get_event_lag");
77
78 free((void*)coupling_time_step_iso);
79 free((void*)model_time_step_iso);
80
81 enum {NUM_TIME_STEPS = 15};
82
83 enum yac_action_type ref_action[NUM_TIME_STEPS] =
84 {COUPLING, REDUCTION, REDUCTION,
85 COUPLING, REDUCTION, REDUCTION,
86 COUPLING, REDUCTION, REDUCTION,
87 COUPLING, REDUCTION, REDUCTION,
88 RESTART, OUT_OF_BOUND, OUT_OF_BOUND};
89
90 for (int t = 0; t < NUM_TIME_STEPS; ++t) {
91
92 if (yac_event_check(event) != ref_action[t]) PUT_ERR("wrong action\n");
93 yac_event_update(event);
94 }
95
96 yac_event_delete(event);
97 }
98
99 {
100 struct event * event = yac_event_new ();
101
102 char const * model_time_step = "5000"; // 5s per time step
103 char const * coupling_time_step = "15000"; // 15s per coupling time step
104
105 enum yac_time_unit_type time_unit = C_MILLISECOND;
106 int lag = 0;
107 enum yac_reduction_type time_operation = TIME_NONE;
108
109 // 60s runtime
110 char const * start_date = "1850-01-01T00:00:00";
111 char const * stop_date = "1850-01-01T00:01:00";
112
113 char const * model_time_step_iso =
114 strdup(yac_time_to_ISO(model_time_step, time_unit));
115 char const * coupling_time_step_iso =
116 strdup(yac_time_to_ISO(coupling_time_step, time_unit));
117 yac_event_add(
118 event, model_time_step_iso, coupling_time_step_iso,
119 lag, time_operation, start_date, stop_date);
120
121 char time_step_buffer[MAX_TIMEDELTA_STR_LEN];
122 if (strcmp(
123 yac_get_event_coupling_timestep(event, time_step_buffer),
124 coupling_time_step_iso))
125 PUT_ERR("ERROR in yac_get_event_coupling_timestep");
126 if (strcmp(
127 yac_get_event_model_timestep(event, time_step_buffer),
128 model_time_step_iso))
129 PUT_ERR("ERROR in yac_get_event_model_timestep");
130 if (yac_get_event_lag(event) != lag)
131 PUT_ERR("ERROR in yac_get_event_lag");
132
133 free((void*)coupling_time_step_iso);
134 free((void*)model_time_step_iso);
135
136 enum {NUM_TIME_STEPS = 15};
137
138 enum yac_action_type ref_action[NUM_TIME_STEPS] =
139 {COUPLING, NONE, NONE,
140 COUPLING, NONE, NONE,
141 COUPLING, NONE, NONE,
142 COUPLING, NONE, NONE,
143 RESTART, OUT_OF_BOUND, OUT_OF_BOUND};
144
145 for (int t = 0; t < NUM_TIME_STEPS; ++t) {
146
147 if (yac_event_check(event) != ref_action[t]) PUT_ERR("wrong action\n");
148 yac_event_update(event);
149 }
150
151 yac_event_delete(event);
152 }
153
154 {
155 struct event * event = yac_event_new ();
156
157 char * model_time_step = "PT30M";
158 char * coupling_time_step = "PT60M";
159
160 enum yac_time_unit_type time_unit = C_ISO_FORMAT;
161 int lag = 0;
162 enum yac_reduction_type time_operation = TIME_AVERAGE;
163
164 // 60s runtime
165 char const * start_date = "1850-01-01T00:00:00";
166 char const * stop_date = "1850-01-01T06:00:00";
167
168 char const * model_time_step_iso =
169 strdup(yac_time_to_ISO(model_time_step, time_unit));
170 char const * coupling_time_step_iso =
171 strdup(yac_time_to_ISO(coupling_time_step, time_unit));
172 yac_event_add(
173 event, model_time_step_iso, coupling_time_step_iso,
174 lag, time_operation, start_date, stop_date);
175
176 char time_step_buffer[MAX_TIMEDELTA_STR_LEN];
177 if (strcmp(
178 yac_get_event_coupling_timestep(event, time_step_buffer),
179 coupling_time_step_iso))
180 PUT_ERR("ERROR in yac_get_event_coupling_timestep");
181 if (strcmp(
182 yac_get_event_model_timestep(event, time_step_buffer),
183 model_time_step_iso))
184 PUT_ERR("ERROR in yac_get_event_model_timestep");
185 if (yac_get_event_lag(event) != lag)
186 PUT_ERR("ERROR in yac_get_event_lag");
187
188 free((void*)coupling_time_step_iso);
189 free((void*)model_time_step_iso);
190
191 enum {NUM_TIME_STEPS = 16};
192
193 enum yac_action_type ref_action[NUM_TIME_STEPS] =
194 {COUPLING, REDUCTION, COUPLING, REDUCTION,
195 COUPLING, REDUCTION, COUPLING, REDUCTION,
196 COUPLING, REDUCTION, COUPLING, REDUCTION,
197 RESTART, OUT_OF_BOUND, OUT_OF_BOUND, OUT_OF_BOUND};
198
199 for (int t = 0; t < NUM_TIME_STEPS; ++t) {
200
201 if (yac_event_check(event) != ref_action[t]) PUT_ERR("wrong action\n");
202 yac_event_update(event);
203 }
204
205 yac_event_delete(event);
206 }
207
208 {
209 char const * model_time_step = "5000"; // 5s per time step
210 char const * coupling_time_step = "15000"; // 15s per coupling time step
211
212 enum yac_time_unit_type time_unit = C_MILLISECOND;
213 enum yac_reduction_type time_operation = TIME_AVERAGE;
214
215 // 60s runtime
216 char const * start_date = "1850-01-01T00:00:00";
217 char const * stop_date = "1850-01-01T00:01:00";
218
219 enum {MIN_LAG = -6, MAX_LAG = 3};
220 enum {NUM_LAGS = MAX_LAG - MIN_LAG + 1};
221 enum {NUM_TIME_STEPS = 15};
222
223 enum yac_action_type ref_action[NUM_TIME_STEPS - MIN_LAG] =
224 {COUPLING, REDUCTION, REDUCTION,
225 COUPLING, REDUCTION, REDUCTION,
226 COUPLING, REDUCTION, REDUCTION,
227 COUPLING, REDUCTION, REDUCTION,
228 COUPLING, REDUCTION, REDUCTION,
229 COUPLING, REDUCTION, REDUCTION,
230 RESTART, OUT_OF_BOUND, OUT_OF_BOUND};
231
232 for (int lag = MIN_LAG; lag <= MAX_LAG; ++lag) {
233
234 struct event * event = yac_event_new ();
235
236 char const * model_time_step_iso =
237 strdup(yac_time_to_ISO(model_time_step, time_unit));
238 char const * coupling_time_step_iso =
239 strdup(yac_time_to_ISO(coupling_time_step, time_unit));
240 yac_event_add(
241 event, model_time_step_iso, coupling_time_step_iso,
242 lag, time_operation, start_date, stop_date);
243
244 char time_step_buffer[MAX_TIMEDELTA_STR_LEN];
245 if (strcmp(
246 yac_get_event_coupling_timestep(event, time_step_buffer),
247 coupling_time_step_iso))
248 PUT_ERR("ERROR in yac_get_event_coupling_timestep");
249 if (strcmp(
250 yac_get_event_model_timestep(event, time_step_buffer),
251 model_time_step_iso))
252 PUT_ERR("ERROR in yac_get_event_model_timestep");
253 if (yac_get_event_lag(event) != lag)
254 PUT_ERR("ERROR in yac_get_event_lag");
255
256 free((void*)coupling_time_step_iso);
257 free((void*)model_time_step_iso);
258
259 for (int t = lag; t < NUM_TIME_STEPS; ++t) {
260
261 enum yac_action_type action = yac_event_check(event);
262 if (action != ref_action[t - MIN_LAG])
263 PUT_ERR("wrong action\n");
264 yac_event_update(event);
265 }
266
267 yac_event_delete(event);
268 }
269 }
270
271 freeCalendar();
272
273 return TEST_EXIT_CODE;
274}
275
276
yac_time_to_ISO
char const * yac_time_to_ISO(char const *time, enum yac_time_unit_type time_unit)
Definition event.c:336
COUPLING
@ COUPLING
Definition config_yaml.c:113
yac_reduction_type
yac_reduction_type
Definition couple_config.h:22
TIME_NONE
@ TIME_NONE
Definition couple_config.h:23
TIME_AVERAGE
@ TIME_AVERAGE
Definition couple_config.h:25
yac_time_unit_type
yac_time_unit_type
Definition couple_config.h:10
C_SECOND
@ C_SECOND
Definition couple_config.h:12
C_HOUR
@ C_HOUR
Definition couple_config.h:14
C_MILLISECOND
@ C_MILLISECOND
Definition couple_config.h:11
C_MINUTE
@ C_MINUTE
Definition couple_config.h:13
C_ISO_FORMAT
@ C_ISO_FORMAT
Definition couple_config.h:18
yac_get_event_model_timestep
char * yac_get_event_model_timestep(struct event *event, char *timedelta_str)
Definition event.c:311
yac_get_event_coupling_timestep
char * yac_get_event_coupling_timestep(struct event *event, char *timedelta_str)
Definition event.c:298
yac_get_event_lag
int yac_get_event_lag(struct event *event)
Definition event.c:284
yac_event_delete
void yac_event_delete(struct event *event)
Definition event.c:265
yac_event_new
struct event * yac_event_new()
Definition event.c:36
yac_event_add
void yac_event_add(struct event *event, char const *delta_model_time, char const *delta_coupling_time, int lag, enum yac_reduction_type time_operation, const char *startdate, const char *stopdate)
Definition event.c:55
yac_event_check
enum yac_action_type yac_event_check(struct event *event)
Definition event.c:182
yac_event_update
void yac_event_update(struct event *event)
Definition event.c:170
yac_action_type
yac_action_type
Definition event.h:13
REDUCTION
@ REDUCTION
Definition event.h:15
RESTART
@ RESTART
Definition event.h:17
OUT_OF_BOUND
@ OUT_OF_BOUND
Definition event.h:22
NONE
@ NONE
Definition event.h:14
value
struct @23::@24 value
test_coupling1.t
int t
Definition test_coupling1.py:42
test_raw_exchange.i
int i
Definition test_raw_exchange.py:89
test_yac_plot_coupling_config.tests
list tests
Definition test_yac_plot_coupling_config.py:281
event
Definition event.c:18
event::time_unit
enum yac_time_unit_type time_unit
Definition event.c:26
event::time_operation
enum yac_reduction_type time_operation
Definition event.c:24
event::action
enum yac_action_type action
Definition event.c:25
event::model_time_step
struct _timedelta * model_time_step
Definition event.c:20
event::lag
int lag
Definition event.c:23
event::coupling_time_step
struct _timedelta * coupling_time_step
Definition event.c:21
NUM_TESTS
#define NUM_TESTS
TEST_EXIT_CODE
#define TEST_EXIT_CODE
Definition tests.h:15
PUT_ERR
#define PUT_ERR(string)
Definition tests.h:10