GSoC - Community Bonding

“Make it work, make it right, make it efficient, make it fast.”

My project " Developing an optimization Library for Libroadrunner " got selected in GSoC 2019 under NRNB, here is my proposal.
 

About the Project


LibRoadrunner is a high-performance SBML based simulator that uses LLVM to generate very efficient runtime code. This enables libRoadrunner to simulate models on par with compiled C/C++ code. By combining libroadrunner with standard optimization algorithms it is possible to use libroadrunner to fit models to data. At present this is done by writing code to link the standard Python optimizers available via scipy with libroadrunner. Although this works, it is inefficient and for large models, it is not practical. In this project, we would like to develop a C/C++ based optimization library that can be used directly by libroadrunner without having to go via Python. This would enable us to provide high-performance optimization capabilities.

Community Bonding Period

 The Community bonding period starts from 6th May to 27th May. I have my last exam at 8th May so I started from 8th May. My mentor Herbert Sir arranged a meeting and discussed about the project and potential problems.

 By meeting experienced people and having an intense discussion about our project was a great experience for me. I have assigned a task to build a plugin which add two number and that must need not to be changed if Libroadrunner gets updated. This plugin will be called from Libroadrunner. Will keep on updating new results as the project goes on.

Comments

Post a Comment

Popular posts from this blog

Monte Carlo Bootstrap

Image
 The main work of Monte Carlo is to " Get some confidence   estimates".   The Monte Carlo plugin is used to obtain estimates of the confidence limits for a model’s parameters. This is in the context where experimental data exists and a parameter minimization method, such as Levenberg-Marquardt or Nelder-Mead has already been used in order to find a parameter minimum. The Monte Carlo algorithm is used subsequently at this minimum and will give an estimate of parameter confidence limits corresponding to the variance in the original experimental data. The plugin has properties such as the size of the Monte Carlo population, minimization algorithm to use (e.g. Nelder-Mead or Levenberg-Marquardt), and on output, confidence limits for each involved parameter. With some past updates, Monte Carlo is not working mainly due to the modifications in rrplugins and change of python version from 2.x to 3.x. I have modified the wrapper to let it work. Majorly the problem was with the passing
Read more

Update old plugins to Work in the current system

Image
 I have updated all the plugins ( https://github.com/sys-bio/rrplugins/tree/master/plugins/released ) working in the old system. Mostly they are not working because of the wrapper. When we pass a python string to C through CPython we need to it the format ctypes.c_char_p(String.encode('utf-8')) where String is the "string" we need to pass. Currently, the optimizer is not working if we need to fit the data by skipping any species. It is the first bad design I have seen in rrplugins. Now I am trying to port all the plugins from the old system to the new system.  One of the optimization example where I have to fit rate K1 and K2 based on species S1, S2, and S3 is:  Model: antimony_str = """ J1: S1 -> S2; k1*S1 J2: S2 -> S3; k2*S2    S1 = 1; S2 = 0; S3 = 0;  k1 = 0; k2 = 0;  """ The exact fitted value is k1 =0.45 and k2=0.15.  My output comes out to be  ChiSquare =  1.6501935363049038 Reduced ChiSquare =  0.017555250386222382 k1
Read more

End of GSoC 2.0

    My second innings of Google Summer of Code comes to an end. I have learned a lot of things. There is still some work left which I will complete after my exams. The link to my summer internship report is https://gist.github.com/debashish05/e1ec288a81d77e0fe61cbf336d147109 . I have learned a lot of things from my mentors and develop an inclination towards research. Works Done during this period: 1) RRPlugins have been integrated with Libroadrunner. So we don't require a separate build for rrplugins. Earlier this task was achieved in my last year project. But the problem is the change in roadrunner will force the plugins to need to rebuild, which is not the case now. 2) Build time will reduce since the plugins and roadrunner are not interdependent. Updated all the plugins (https://github.com/sys-bio/rrplugins/tree/master/plugins/released), in the current system. Earlier these plugins were not working due to outdated wrapper and some minor bugs. 3) Updated Python wrapper for rrplug
Read more