|
40 | 40 | }, |
41 | 41 | { |
42 | 42 | "cell_type": "code", |
43 | | - "execution_count": 1, |
| 43 | + "execution_count": 2, |
44 | 44 | "metadata": {}, |
45 | 45 | "outputs": [], |
46 | 46 | "source": [ |
|
253 | 253 | }, |
254 | 254 | { |
255 | 255 | "cell_type": "code", |
256 | | - "execution_count": 1, |
| 256 | + "execution_count": 5, |
257 | 257 | "metadata": { |
258 | 258 | "scrolled": false |
259 | 259 | }, |
260 | 260 | "outputs": [ |
261 | 261 | { |
262 | | - "ename": "NameError", |
263 | | - "evalue": "name 'SPOD_streaming' is not defined", |
264 | | - "output_type": "error", |
265 | | - "traceback": [ |
266 | | - "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", |
267 | | - "\u001b[0;31mNameError\u001b[0m Traceback (most recent call last)", |
268 | | - "\u001b[0;32m<ipython-input-1-8189a6bbb6a8>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m 1\u001b[0m \u001b[0;31m# Perform SPOD analysis using low storage module\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 2\u001b[0;31m \u001b[0mSPOD_analysis\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mSPOD_streaming\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mX\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mX\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mparams\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mparams\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdata_handler\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mFalse\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvariables\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mvariables\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 3\u001b[0m \u001b[0mspod\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mSPOD_analysis\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfit\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n", |
269 | | - "\u001b[0;31mNameError\u001b[0m: name 'SPOD_streaming' is not defined" |
| 262 | + "name": "stdout", |
| 263 | + "output_type": "stream", |
| 264 | + "text": [ |
| 265 | + "\n", |
| 266 | + "SPOD parameters\n", |
| 267 | + "------------------------------------\n", |
| 268 | + "Spectrum type : one-sided (real-valued signal)\n", |
| 269 | + "No. of snaphots per block : 64\n", |
| 270 | + "Block overlap : 32\n", |
| 271 | + "No. of blocks : 30\n", |
| 272 | + "Windowing fct. (time) : hamming\n", |
| 273 | + "Weighting fct. (space) : user-specified\n", |
| 274 | + "Mean : blockwise\n", |
| 275 | + "Time-step : 0.2\n", |
| 276 | + "Number of Frequencies : 33\n", |
| 277 | + "------------------------------------\n", |
| 278 | + "\n", |
| 279 | + "DATA MATRIX DIMENSIONS\n", |
| 280 | + "------------------------------------\n", |
| 281 | + "Problem size : 0.013113021850585958 GB. (double)\n", |
| 282 | + "data matrix dimensions: (1, 20, 88, 1)\n", |
| 283 | + "Make sure that first column of data matrix is time and last column is number of variables. \n", |
| 284 | + "First column dimension: 1 must correspond to number of time snapshots.\n", |
| 285 | + "Last column dimension: 1 must correspond to number of variables.\n", |
| 286 | + "------------------------------------\n", |
| 287 | + " \n", |
| 288 | + "Calculating temporal DFT (streaming)\n", |
| 289 | + "------------------------------------\n", |
| 290 | + "--> Updating left singular vectors Time 63 / block 1\n" |
| 291 | + ] |
| 292 | + }, |
| 293 | + { |
| 294 | + "name": "stderr", |
| 295 | + "output_type": "stream", |
| 296 | + "text": [ |
| 297 | + "/usr/local/lib/python3.8/site-packages/pyspod/spod_streaming.py:203: RuntimeWarning: divide by zero encountered in true_divide\n", |
| 298 | + " mse_prev[block_i,:,:] = (np.abs(S_hat_prev**2 - self._eigs**2)**2) / (S_hat_prev**2)\n", |
| 299 | + "/usr/local/lib/python3.8/site-packages/pyspod/spod_streaming.py:203: RuntimeWarning: invalid value encountered in true_divide\n", |
| 300 | + " mse_prev[block_i,:,:] = (np.abs(S_hat_prev**2 - self._eigs**2)**2) / (S_hat_prev**2)\n" |
| 301 | + ] |
| 302 | + }, |
| 303 | + { |
| 304 | + "name": "stdout", |
| 305 | + "output_type": "stream", |
| 306 | + "text": [ |
| 307 | + "--> Updating left singular vectors Time 95 / block 2\n", |
| 308 | + "--> Updating left singular vectors Time 128 / block 3\n", |
| 309 | + "--> Updating left singular vectors Time 160 / block 4\n", |
| 310 | + "--> Updating left singular vectors Time 193 / block 5\n", |
| 311 | + "--> Updating left singular vectors Time 225 / block 6\n", |
| 312 | + "--> Updating left singular vectors Time 258 / block 7\n", |
| 313 | + "--> Updating left singular vectors Time 290 / block 8\n", |
| 314 | + "--> Updating left singular vectors Time 323 / block 9\n", |
| 315 | + "--> Updating left singular vectors Time 355 / block 10\n", |
| 316 | + "--> Updating left singular vectors Time 388 / block 11\n", |
| 317 | + "--> Updating left singular vectors Time 420 / block 12\n", |
| 318 | + "--> Updating left singular vectors Time 453 / block 13\n", |
| 319 | + "--> Updating left singular vectors Time 485 / block 14\n", |
| 320 | + "--> Updating left singular vectors Time 518 / block 15\n", |
| 321 | + "--> Updating left singular vectors Time 550 / block 16\n", |
| 322 | + "--> Updating left singular vectors Time 583 / block 17\n", |
| 323 | + "--> Updating left singular vectors Time 615 / block 18\n", |
| 324 | + "--> Updating left singular vectors Time 648 / block 19\n", |
| 325 | + "--> Updating left singular vectors Time 680 / block 20\n", |
| 326 | + "--> Updating left singular vectors Time 713 / block 21\n", |
| 327 | + "--> Updating left singular vectors Time 745 / block 22\n", |
| 328 | + "--> Updating left singular vectors Time 778 / block 23\n", |
| 329 | + "--> Updating left singular vectors Time 810 / block 24\n", |
| 330 | + "--> Updating left singular vectors Time 843 / block 25\n", |
| 331 | + "--> Updating left singular vectors Time 875 / block 26\n", |
| 332 | + "--> Updating left singular vectors Time 908 / block 27\n", |
| 333 | + "--> Updating left singular vectors Time 940 / block 28\n", |
| 334 | + "--> Updating left singular vectors Time 973 / block 29\n", |
| 335 | + "--> Data stream ended.\n", |
| 336 | + "Elapsed time: 3.8960890769958496 s.\n" |
270 | 337 | ] |
271 | 338 | } |
272 | 339 | ], |
|
297 | 364 | }, |
298 | 365 | { |
299 | 366 | "cell_type": "code", |
300 | | - "execution_count": 7, |
| 367 | + "execution_count": 6, |
301 | 368 | "metadata": {}, |
302 | | - "outputs": [], |
| 369 | + "outputs": [ |
| 370 | + { |
| 371 | + "name": "stdout", |
| 372 | + "output_type": "stream", |
| 373 | + "text": [ |
| 374 | + "- RAM required for loading all modes ~ 0.0001573562622070315 GB\n", |
| 375 | + "- Available RAM memory ~ 14.751468658447289 GB\n" |
| 376 | + ] |
| 377 | + } |
| 378 | + ], |
303 | 379 | "source": [ |
304 | 380 | "# Show results\n", |
305 | 381 | "T_approx = 12.5 # approximate period\n", |
|
316 | 392 | }, |
317 | 393 | { |
318 | 394 | "cell_type": "code", |
319 | | - "execution_count": 8, |
| 395 | + "execution_count": 7, |
320 | 396 | "metadata": { |
321 | 397 | "scrolled": true |
322 | 398 | }, |
|
347 | 423 | }, |
348 | 424 | { |
349 | 425 | "cell_type": "code", |
350 | | - "execution_count": 9, |
| 426 | + "execution_count": 8, |
351 | 427 | "metadata": {}, |
352 | 428 | "outputs": [ |
353 | 429 | { |
|
370 | 446 | }, |
371 | 447 | { |
372 | 448 | "cell_type": "code", |
373 | | - "execution_count": 10, |
| 449 | + "execution_count": 9, |
374 | 450 | "metadata": {}, |
375 | 451 | "outputs": [ |
376 | 452 | { |
|
0 commit comments