Coverage for pySDC / projects / Second_orderSDC / harmonic_oscillator_run_stability.py: 0%

1from pySDC.projects.Second_orderSDC.harmonic_oscillator_params import get_default_harmonic_oscillator_description 

2from pySDC.projects.Second_orderSDC.stability_simulation import StabilityImplementation 

3 

4if __name__ == '__main__': 

5 """ 

6 To implement Stability region for the Harmonic Oscillator problem 

7 Run for 

8 SDC stability region: model_stab.run_SDC_stability() 

9 Picard stability region: model_stab.run_Picard_stability() 

10 Runge-Kutta-Nzström stability region: model_run_RKN_stability() 

11 To implement spectral radius of iteration matrix 

12 Run: 

13 Iteration matrix of SDC method: model_stab.run_Ksdc() 

14 Iteration matrix of Picard method: model_stab.run_Kpicard() 

15 

16 """ 

17 # This code checks if the "data" folder exists or not. 

18 exec(open("check_data_folder.py").read()) 

19 # Execute the stability analysis for the damped harmonic oscillator 

20 description = get_default_harmonic_oscillator_description() 

21 # ============================================================================= 

22 # maxiter can be changed here manually. By default, maxiter is 50 

23 description['step_params']['maxiter'] = 50 

24 # ============================================================================= 

25 

26 model_stab = StabilityImplementation(description, kappa_max=18, mu_max=18, Num_iter=(200, 200)) 

27 

28 model_stab.run_SDC_stability() 

29 model_stab.run_Picard_stability() 

30 model_stab.run_RKN_stability() 

31 model_stab.run_Ksdc() 

32 # model_stab.run_Kpicard()