@@ -42,8 +42,8 @@ def findomega(stab_fh):
4242 pparams ['u0' ] = 1.0
4343 swparams = {}
4444 #swparams['collocation_class'] = collclass.CollGaussLobatto
45- swparams ['collocation_class' ] = collclass .CollGaussLegendre
46- # swparams['collocation_class'] = collclass.CollGaussRadau_Right
45+ # swparams['collocation_class'] = collclass.CollGaussLegendre
46+ swparams ['collocation_class' ] = collclass .CollGaussRadau_Right
4747 swparams ['num_nodes' ] = 3
4848 K = 5
4949 dirk_order = K
@@ -130,7 +130,7 @@ def findomega(stab_fh):
130130 fig = plt .figure ()
131131 plt .plot (k_vec , (U_speed + c_speed )+ np .zeros (np .size (k_vec )), '--' , color = 'k' , linewidth = 1.5 , label = 'Exact' )
132132 plt .plot (k_vec , phase [1 ,:], '-' , color = 'g' , linewidth = 1.5 , label = 'DIRK(' + str (dirkts .order )+ ')' )
133- plt .plot (k_vec , phase [2 ,:], '-+' , color = 'r' , linewidth = 1.5 , label = 'RK- IMEX(' + str (rkimex .order )+ ')' , markevery = 5 , mew = 1.0 )
133+ plt .plot (k_vec , phase [2 ,:], '-+' , color = 'r' , linewidth = 1.5 , label = 'IMEX(' + str (rkimex .order )+ ')' , markevery = 5 , mew = 1.0 )
134134 plt .plot (k_vec , phase [0 ,:], '-o' , color = 'b' , linewidth = 1.5 , label = 'SDC(' + str (K )+ ')' , markevery = 5 , markersize = fs / 2 )
135135 plt .xlabel ('Wave number' , fontsize = fs , labelpad = 0.25 )
136136 plt .ylabel ('Phase speed' , fontsize = fs , labelpad = 0.5 )
@@ -147,7 +147,7 @@ def findomega(stab_fh):
147147 fig = plt .figure ()
148148 plt .plot (k_vec , 1.0 + np .zeros (np .size (k_vec )), '--' , color = 'k' , linewidth = 1.5 , label = 'Exact' )
149149 plt .plot (k_vec , amp_factor [1 ,:], '-' , color = 'g' , linewidth = 1.5 , label = 'DIRK(' + str (dirkts .order )+ ')' )
150- plt .plot (k_vec , amp_factor [2 ,:], '-+' , color = 'r' , linewidth = 1.5 , label = 'RK- IMEX(' + str (rkimex .order )+ ')' , markevery = 5 , mew = 1.0 )
150+ plt .plot (k_vec , amp_factor [2 ,:], '-+' , color = 'r' , linewidth = 1.5 , label = 'IMEX(' + str (rkimex .order )+ ')' , markevery = 5 , mew = 1.0 )
151151 plt .plot (k_vec , amp_factor [0 ,:], '-o' , color = 'b' , linewidth = 1.5 , label = 'SDC(' + str (K )+ ')' , markevery = 5 , markersize = fs / 2 )
152152 plt .xlabel ('Wave number' , fontsize = fs , labelpad = 0.25 )
153153 plt .ylabel ('Amplification factor' , fontsize = fs , labelpad = 0.5 )
0 commit comments