An error occurred while loading the file. Please try again.
-
Thibault Hallouin authoredThis is required for the bindings, because the 2D views on 4D tensors must be on the correct type (i.e. pytensor/rtensor/xtensor). Although, casting of the views to become xtensor could have been used instead. In doing so, there is no source file to be compiled anymore, so this becomes a header-only library. But to keep the "public" headers separate from the implementation, a subdirectory "detail/" is added.
ee8e2ef0
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
import unittest
from extreme_data.meteo_france_data.adamont_data.adamont.adamont_snowfall import AdamontSnowfall
from extreme_data.meteo_france_data.adamont_data.adamont_scenario import AdamontScenario
from extreme_data.meteo_france_data.scm_models_data.safran.safran import SafranSnowfall1Day
from projects.altitude_spatial_model.altitudes_fit.altitudes_studies import AltitudesStudies
from spatio_temporal_dataset.slicer.split import Split
class TestAltitudesStudies(unittest.TestCase):
def setUp(self) -> None:
super().setUp()
altitudes = [900, 1200]
study_class = SafranSnowfall1Day
self.studies = AltitudesStudies(study_class, altitudes, year_min=1959, year_max=1962)
class TestVisualization(TestAltitudesStudies):
def test_plot_maxima_time_series(self):
self.studies.plot_maxima_time_series(massif_names=['Vercors'], show=False)
class TestSpatioTemporalCoordinates(TestAltitudesStudies):
def test_temporal_split(self):
s_split_temporal = self.studies.random_s_split_temporal(train_split_ratio=0.75)
coordinates = self.studies.spatio_temporal_coordinates(s_split_temporal=s_split_temporal)
self.assertEqual(coordinates.coordinates_values(split=Split.train_temporal).shape, (6, 2))
self.assertEqual(coordinates.coordinates_values(split=Split.test_temporal).shape, (2, 2))
def test_spatial_split(self):
s_split_spatial = self.studies.random_s_split_spatial(train_split_ratio=0.5)
coordinates = self.studies.spatio_temporal_coordinates(s_split_spatial=s_split_spatial)
self.assertEqual(coordinates.coordinates_values(split=Split.train_spatial).shape, (4, 2))
self.assertEqual(coordinates.coordinates_values(split=Split.test_spatial).shape, (4, 2))
def test_spatio_temporal_split(self):
s_split_spatial = self.studies.random_s_split_spatial(train_split_ratio=0.5)
s_split_temporal = self.studies.random_s_split_temporal(train_split_ratio=0.75)
coordinates = self.studies.spatio_temporal_coordinates(s_split_spatial=s_split_spatial,
s_split_temporal=s_split_temporal)
self.assertEqual(coordinates.coordinates_values(split=Split.train_spatiotemporal).shape, (3, 2))
self.assertEqual(coordinates.coordinates_values(split=Split.test_spatiotemporal_spatial).shape, (3, 2))
self.assertEqual(coordinates.coordinates_values(split=Split.test_spatiotemporal_temporal).shape, (1, 2))
self.assertEqual(coordinates.coordinates_values(split=Split.test_spatiotemporal).shape, (1, 2))
class TestSpatioTemporalDataset(TestAltitudesStudies):
def setUp(self) -> None:
super().setUp()
self.massif_name = "Vercors"
def test_temporal_split(self):
s_split_temporal = self.studies.random_s_split_temporal(train_split_ratio=0.75)
dataset = self.studies.spatio_temporal_dataset(massif_name=self.massif_name,
s_split_temporal=s_split_temporal)
self.assertEqual(len(dataset.maxima_gev(split=Split.train_temporal)), 6)
self.assertEqual(len(dataset.maxima_gev(split=Split.test_temporal)), 2)
def test_spatial_split(self):
s_split_spatial = self.studies.random_s_split_spatial(train_split_ratio=0.5)
dataset = self.studies.spatio_temporal_dataset(massif_name=self.massif_name,
s_split_spatial=s_split_spatial)
self.assertEqual(len(dataset.maxima_gev(split=Split.train_spatial)), 4)
self.assertEqual(len(dataset.maxima_gev(split=Split.test_spatial)), 4)
def test_spatio_temporal_split(self):
s_split_spatial = self.studies.random_s_split_spatial(train_split_ratio=0.5)
s_split_temporal = self.studies.random_s_split_temporal(train_split_ratio=0.75)
dataset = self.studies.spatio_temporal_dataset(massif_name=self.massif_name,
s_split_spatial=s_split_spatial,
s_split_temporal=s_split_temporal)
self.assertEqual(len(dataset.maxima_gev(split=Split.train_spatiotemporal)), 3)
self.assertEqual(len(dataset.maxima_gev(split=Split.test_spatiotemporal)), 1)
self.assertEqual(len(dataset.maxima_gev(split=Split.test_spatiotemporal_temporal)), 1)
self.assertEqual(len(dataset.maxima_gev(split=Split.test_spatiotemporal_spatial)), 3)
class TestSpatioTemporalDatasetForClimateModels(unittest.TestCase):
def setUp(self) -> None:
super().setUp()
altitudes = [900, 1200]
study_class = AdamontSnowfall
self.studies = AltitudesStudies(study_class, altitudes,
year_min=2009, year_max=2012,
scenario=AdamontScenario.rcp85)
self.massif_name = "Vercors"
def test_dataset(self):
dataset = self.studies.spatio_temporal_dataset(self.massif_name)
self.assertEqual(len(dataset.coordinates.df_coordinate_climate_model.columns), 3)
if __name__ == '__main__':
unittest.main()