@@ -76,26 +76,14 @@ namespace pfasst
7676
7777
7878 template <typename scalar>
79- static void init_config_options (po::options_description& opts )
79+ static void init_opts ( )
8080 {
81- opts.add_options ()
82- (" num_particles" , po::value<size_t >(), " number of particles in the cloud" )
83- (" epsilon" , po::value<scalar>(), " Boris' epsilon" )
84- (" omega_e" , po::value<scalar>(), " E-field constant" )
85- (" omega_b" , po::value<scalar>(), " B-field constant" )
86- ;
81+ pfasst::config::options::add_option<size_t >(" Boris-SDC" , " num_particles" , " number of particles in the cloud" );
82+ pfasst::config::options::add_option<scalar>(" Boris-SDC" , " epsilon" , " Boris' epsilon" );
83+ pfasst::config::options::add_option<scalar>(" Boris-SDC" , " omega_e" , " E-field constant" );
84+ pfasst::config::options::add_option<scalar>(" Boris-SDC" , " omega_b" , " B-field constant" );
8785 }
8886
89- template <typename scalar>
90- static void enable_config_options (size_t index = -1 )
91- {
92- config::Options::get_instance ()
93- .register_init_function (" Boris-SDC" ,
94- std::function<void (po::options_description&)>(pfasst::examples::boris::init_config_options<scalar>),
95- index);
96- }
97-
98-
9987 template <
10088 typename scalar,
10189 typename time
@@ -318,7 +306,7 @@ namespace pfasst
318306 time dt = this ->get_controller ()->get_time_step ();
319307
320308 C omega_tilde = sqrt (-2.0 * epsilon) * omega_e;
321- q.positions ()[0 ][2 ] = (z0 * cos (omega_tilde * (scalar)(dt))
309+ q.positions ()[0 ][2 ] = (z0 * cos (omega_tilde * (scalar)(dt))
322310 + w0 / omega_tilde * sin (omega_tilde * (scalar)(dt))).real ();
323311
324312 C sqrt_in_omega = sqrt (pow (omega_b, 2 ) + 4.0 * epsilon * pow (omega_e, 2 ));
@@ -440,7 +428,7 @@ namespace pfasst
440428
441429 // building rules for Q_E and Q_I:
442430 // Q_E is striclty lower diagonal matrix with delta nodes of column index
443- // Q_I is lower diagonal matrix with first row and column all zero and delta nodes of
431+ // Q_I is lower diagonal matrix with first row and column all zero and delta nodes of
444432 // column index minus one
445433 Matrix<time> qe_mat = Matrix<time>(nnodes, nnodes);
446434 qe_mat.fill (time (0.0 ));
0 commit comments