|
1070 | 1070 | "outputs": [], |
1071 | 1071 | "source": [ |
1072 | 1072 | "nsnaps = 30\n", |
1073 | | - "factor = round(nt/nsnaps)\n", |
| 1073 | + "factor = int(np.ceil(nt/nsnaps))\n", |
1074 | 1074 | "\n", |
1075 | 1075 | "time_subsampled = ConditionalDimension('t_sub', parent=model.grid.time_dim, factor=factor)\n", |
1076 | 1076 | "\n", |
1077 | | - "u_save = TimeFunction(name='usave', grid=snapshotdomain, time_order=2, space_order=2,\n", |
| 1077 | + "u_save = TimeFunction(name='usave', grid=snapshotdomain, time_order=0, space_order=2,\n", |
1078 | 1078 | " save=nsnaps, time_dim=time_subsampled)\n", |
1079 | 1079 | "\n", |
1080 | 1080 | "# \"Normal\" snapshotting for comparison\n", |
1081 | | - "u_save_grid = TimeFunction(name='usavegrid', grid=model.grid, time_order=2, space_order=2,\n", |
| 1081 | + "u_save_grid = TimeFunction(name='usavegrid', grid=model.grid, time_order=0, space_order=2,\n", |
1082 | 1082 | " save=nsnaps, time_dim=time_subsampled)" |
1083 | 1083 | ] |
1084 | 1084 | }, |
|
2499 | 2499 | "metadata": {}, |
2500 | 2500 | "outputs": [], |
2501 | 2501 | "source": [ |
2502 | | - "assert np.isclose(np.linalg.norm(u_save.data), 53.853184, atol=0, rtol=1e-4)" |
| 2502 | + "assert np.isclose(np.linalg.norm(u_save.data), 54.026115, atol=0, rtol=1e-4)\n", |
| 2503 | + "assert np.isclose(np.linalg.norm(u_save_grid.data), 77.89853, atol=0, rtol=1e-4)" |
2503 | 2504 | ] |
2504 | 2505 | }, |
2505 | 2506 | { |
|
0 commit comments