- File givens_base.h
- Make sure create_givens() in givens.h is documented.
- File householder_base.h
- Better documentation for the Householder functions.
- Class o2scl::astep_gsl< vec_y_t, vec_dydx_t, vec_yerr_t, func_t >
Document what happens when the stepper function returns a non-zero value, as it's different now with the ode-initval2 function.
Document count, failed_steps, etc.
- Class o2scl::cli
- Long options cannot be one letter long, or else process_args() will fail, thus the class should throw if a long option with only one letter is given.
- Global o2scl::cli::process_args (std::string s, std::vector< cmd_line_arg > &ca, int debug=0)
- There's a typecast in this function to (char *) from (const char *) which needs reworking.
- Class o2scl::contour
Convert distances in find_next_point functions to be scaled by grid spacing.
Copy constructor
- Global o2scl::count_words (std::string str)
- Make consistent with split_string().
- Global o2scl::deriv_eqi< func_t, vec_t >::deriv_vector (size_t nv, double dx, const vec_t &y, vec_t &dydx)
- generalize to other values of npoints.
- Class o2scl::fit_nonlin< func_t, vec_t, mat_t >
Allow the user to specify the derivatives
Fix so that the user can specify automatic scaling of the fitting parameters, where the initial guess are used for scaling so that the fitting parameters are near unity.
- Class o2scl::hist
Check implementation of hist::extend_lhs
.
More testing.
- Global o2scl::inte_cauchy_cern< func_t >::integ_err (func_t &func, double a, double b, double &res, double &err)
- Fix converge error issue here.
- Global o2scl::inte_gauss_cern< func_t >::integ_err (func_t &func, double a, double b, double &res, double &err)
- Fix converge error issue here.
- Class o2scl::inte_qawf_gsl_sin< func_t >
- More documentation and examples for the qawf, qawo and qawc integrators.
- Class o2scl::mcarlo_vegas< func_t, vec_t, rng_t, rng_dist_t >
- Mode = importance only doesn't give the same answer as GSL yet.
- Global o2scl::mcarlo_vegas< func_t, vec_t, rng_t, rng_dist_t >::vegas_minteg_err (int stage, func_t &func, size_t ndim, const vec_t &xl, const vec_t &xu, double &res, double &err)
Should stage be passed by reference?
There was an update between gsl-1.12 and 1.15 which has not been implemented here yet.
- Class o2scl::mmin_bfgs2< func_t, vec_t, dfunc_t, auto_grad_t, def_auto_grad_t >
- While BFGS does well in the
ex_mmin
example with the initial guess of
it seems to converge more poorly for the spring function than the other minimizers with other initial guesses, and I think this will happen in the GSL versions too. I need to examine this more closely with some code designed to clearly show this.
- Class o2scl::mmin_fix_params< func_t, vec_t >
- This hasn't been converted yet with the new C++11 function types.
- Global o2scl::mroot_hybrids< func_t, vec_t, mat_t, jfunc_t >::iterate ()
- This now just returns zero during convergence errors, because it caused problems with tov_eos_ts, but this should be fixed so that it returns a non-zero value, and functions which call iterate() need to handle the return value appropriately.
- Class o2scl::ode_bv_multishoot< func_t, vec_t, alloc_vec_t, alloc_t, vec_int_t, mat_t >
- Improve documentation a little and create testing code
- Class o2scl::ode_control_gsl< vec_y_t, vec_dydx_t, vec_yerr_t >
- Double check that the improvements in the ode-initval2 routines are available here
- Global o2scl::ode_iv_solve< func_t, vec_t >::solve_final_value (double x0, double x1, double h, size_t n, vec_t &ystart, vec_t ¥d, func_t &derivs)
- Document if yend can be the same as ystart.
- Class o2scl::ode_rkf45_gsl< vec_y_t, vec_dydx_t, vec_yerr_t, func_t >
- Check this because it may not give exact dydt_out.
- Class o2scl::other_todos_and_bugs
- Make sure an error message is printed to ensure users do 'make o2scl-test' and not 'make check'
- Make sure cstdlib is included wherever exit() is used [5/22/11 - Most of these are taken care of now.]
- Class o2scl::root_de< func_t, dfunc_t >
- At the moment, the functions solve() and solve_bkt() are not implemented for derivative solvers.
- Class o2scl::search_vec_ext< vec_t >
- The constructor is too restrictive, as it actually calls the parent search_vec constructor and thus requires nn<2 instead of nn<1. Fix this.
- Class o2scl::smooth_gsl
- Needs a bit more error checking and more documentation.
- Global o2scl::split_string (std::string str, std::vector< std::string > &sv)
- More documentation
- Add user-specified delimiters?
- Add version which ignores quotes
- Use this function in acol
- Class o2scl::table< vec_t >
- Specify somewhere what kind of vector types can be used for the template parameter. ublas objects work, but what about armadillo and Eigen vectors? The main reason the default type is std::vector is because of HDF5 I/O.
- Global o2scl::table< vec_t >::function_find_row (std::string function) const
- There is a fixme entry in this function associated with error handling in function parser class.
- Global o2scl::table< vec_t >::sort_column (std::string scol)
- Use vector_sort() rather than qsort().
- Global o2scl::table< vec_t >::subtable (std::string list, size_t top, size_t bottom) const
- Modify to return a shared pointer.
- Class o2scl::tensor< vec_t, vec_size_t_int >
Specify somewhere what kind of vector types can be used for the template parameter. ublas objects work, but what about armadillo and Eigen vectors? The main reason the default type is std::vector is because of HDF5 I/O.
The o2scl::tensor::vector_slice() function should clearly work for uBlas vectors, and seems to work with std::vector objects also, but latter use has not been fully tested.
- Class o2scl::tensor_grid
Make this a template like the o2scl::tensor class?
It is possible for the user to create a tensor_grid object, upcast it to a tensor object, and then use tensor::resize() to resize the tensor, failing to resize the grid. This probably needs fixing.
- Global o2scl::vector_invert_enclosed_sum (double sum, size_t n, vec_t &x, vec2_t &y, double &lev, int verbose=0)
- This function may also require that all of the y-vector values have the same sign or are all positive. Check this.
- Class o2scl_acol::acol_manager
Rework so that the error handler isn't called prematurely. Ensure table::make_cols() and table::calc() don't call the error handler.
Ensure add() copies constants and interpolation type, and units.
- Global o2scl_hdf::hdf_input (hdf_file &hf, o2scl::table_units< vec_t > &t, std::string name)
- Removed default value for
name
for compiling at nersc
- Global o2scl_hdf::hdf_input (hdf_file &hf, o2scl::table< vec_t > &t, std::string name)
- Removed default value for
name
for compiling at nersc
- Global o2scl_linalg::SV_decomp (size_t M, size_t N, mat_t &A, mat2_t &V, vec_t &S, vec2_t &work)
- Test N=1 case, N=2 case, and non-square matrices.