Exception of type: OPTION_INVALID in file "IpAlgBuilder.cpp" at line 271:
 Exception message: Selected linear solver MA27 not available.
Tried to obtain MA27 from shared library "libhsl.so", but the following error occured:
libhsl.so: cannot open shared object file: No such file or directory

EXIT: Invalid option encountered.

Update: Should terminate if any unknown options are set

• Ipopt
• Ceres
• AMICI (is now included via git subrepo)

e.g. https://github.com/robotology/icub-main/blob/master/conf/FindIPOPT.cmake https://github.com/cdcseacave/openMVS/blob/master/build/Modules/FindCERES.cmake

So far only happening for ipopt

Make as generic as possible

add to ParallelOptimizationProblem?

should be usable like ./optimizationOptions bla.h5 -s ceres/maxiter 1

Currently disabled because __node / STL problem.

Add signal handler to cancel all tasks and flush output.

E.g.

• finalize column names
• rules for observable names
• rules for scaling parameter names in model and data files (_offset_, scaling, sigma`)
• specification of timepoints, inf, units, ...
• document hdf5 output
• provide example data and sample measurements

• ensure sigmas provide in data file are used
• compute analytically

In intermediate function, evaluate model on test set and see if prediction likelihood improves

-> need to add validation set to dataprovider

Should be usable for both batch and minibatch optimizers

IpOpt may call Eval_F and Eval_Grad_F with identical parameters. Save previous results and check for new_x to avoid recomputation.

see amiciUserDataSerialization
same for ExpData, so worker can simply run simulations based on received messages.

Will simplify #10

of model, suitesparse, ipopt, ceres, ...

To allow for negative scaling parameters.

Should only be necessary for offset parameters?!

Specify amici options in hdf5 file

For IpOpt and Ceres

For Ipopt, mind difference between ApplicationReturnStatus and SolverReturn.

Need to reorganize training data

class MinibatchDataProvider : public MultiConditionDataProvider

Unlike MultiConditionDataProvider, need one instance per optimization to keep track of batches

Add class Optimizer, class LocalOptimizer : public Optimizer, class MiniBatchOptimizer : public LocalOptimizer
Add Optimizer::optimize(OptimizationProblem*)
-> need to refactor Ceres and Ipopt wrappers

Currently all dependencies are scanned as well.

• hdf5 import (simulation<->optimization parameter mapping, scaling factors)
• remove "genotypespecific" parameter, replace by parameter map

Will make it easier to minibatch methods later on

add factory method to OptimizationOptions

Subclass OptimizationOptions for each Optimizer to account for unique options?

Find problem and reenable mock support.

• https://worhp.de/ (need to check licensing conditions)
• https://www.math.uni-bremen.de/zetem/cms/media.php/279/2546042wassel.pdf
• FFSQP/CFSQP
• https://esa.github.io/pagmo2/docs/algorithm_list.html

Need to run simulation in separate thread on worker.

This will also facilitate sending cancel requests from master.

See /fixedParameters/referenceConditions

Don't wait for all jobs to finish.

Because execute_process is not run before every build.

TODO: Add dummy target as in https://stackoverflow.com/questions/13920072/how-to-always-run-command-when-building-regardless-of-any-dependency

Use h5copy for now, e.g.:

h5copy -i amici/examples/steadystate/data.h5 -s / -o steadystate_2017-06-28_rank00000.h5 -d /inputData

NOTE: https://support.hdfgroup.org/HDF5/doc/RM/RM_H5O.html#Object-Copy

http://people.sc.fsu.edu/~jburkardt/f_src/toms611/toms611.html

https://github.com/pyushkevich/cmrep/blob/master/extras/toms611/toms611.c
entry functions:
sumsl_
/* *** minimize general unconstrained objective function using *** /
/ *** analytic gradient and hessian approx. from secant update *** */

smsno_
/* *** minimize general unconstrained objective function using /
/ *** finite-difference gradients and secant hessian approximations. */

Ensure proper naming in result files

To be able to use same starting points for different optimizers

HDF5 file: measurements/y(sigma): condition x t x y
-> should enable chunked storage / compression
-> need to label timepoints: /measurements/[attr]timepoints double; [infinity] for steady-state data

Currently only integer values are used, which is hard to read and error-prone. Need string-to-enum. Available for ceres in types.cc.

To avoid recomputation in case of crash or termination due to exceeding wall-time limits.

EDIT: warm-start probably not possibly with all optimizers

To speed up tests

• libceres-dev
• coinor-libipopt-dev
• libcpputest-dev

Check versions

HDF5 generation script still looks for _genotypespecific_ in parameter names. get rid of that.