User:Antonio Pino
Public application template
A: An introduction
- Please describe yourself in three sentences, one of them regarding your current studies.
I am Antonio Pino Robles and I am currently an undergraduate student at the University of the Basque Country, in the Basque Country, Northern Spain, Atlantic French border. Thus I an a native Basque and Spanish speaker, I can express myself in French and Italian (after a year long stay in Piemonte), besides being fluent in English.
I have always enjoyed literature and math, which are things I like to do when I enjoy the heavy rain pouring in my home town.
- What's your overall background?
The undergraduate programme I am following leads me to a major in Electronic Engineering and a minor in Physics, that means that I have experience in subjects such as Control Theory, Signal Analysis on one side, and Quantum Physics and Thermodynamics on the other; all with a strong mathematical background.
- Why do you want to participate in the Google Summer of Code? What do you hope to gain by doing so?
I hope to learn how programming is done outside of the sand boxed problems I have done until now (at Uni) while contributing to free software. Besides, I have been following the FSF for a long time, and been a free software user for longer, so this can be a way of giving back while receiving.
- Why are you choosing Octave?
Last year, I took two courses the relied heavily on the matrix exponential (Control Theory and Non-linear dynamics) to solve systems (or approximate them) by
.
Needless to say, they tried to push Matlab to all the students, and I got by with Gnu Octave, for which I am thankful to the whole developer team.
C: Contact
- Please state the (unique and identical where possible) nick you use on IRC and any other communication channel related to Octave.
While at IRC my nickname will be AntonioPino, and I am readily available by email at data.script93 <at> gmail.com
- Which time zone do you live in? Will that change over GSoC duration?
In Spain we are at UTC+1 (DST UTC+2), although there have been talks to change to an UTC+0 (DST UTC+1) frame, I am not aware of those taking place soon, nor do I plan to travel outside UTC+2 this summer.
- Please state the timeframe (in UTC+0) when you feel most comfortable working during GSoC. Where are your time buffers?
I am really a morning person, so a eight to six (UTC+0) schedule will do for me.
E: Coding experience
This part is one of the more important ones in your application. You are allowed to be as verbose as you want, as long as you stay on topic ;-)
- Please describe your experience with C++, Octave or Matlab m-scripts, OpenGL and Qt.
I have been writing m-scripts for two years now, to solve specific engineering problems, 100-200 lines long on average. Apart from that, I have a basic level of C++ (not with the octave api) and I have not used OpenGL nor Qt before.
- Please describe your experience with other programming languages.
A big part of the programming done at Uni has been done in Python: from the basic algorithm courses to symbolic problem solving with sympy. It has also been the tool of choice for a little kludge here and there. On the other hand, I have used VHDL for building and testing a discrete signal filter. And c to do a little systems and microcontroller programming.
- Please describe your experience with being in a development team.
If selected this will be my first experience as a free software developer.
- Please describe the biggest project you have written code for and what you learned by doing so. Also describe your role in that project over time.
The biggest project I have worked on was the VHDL discrete filter, carried out along with two colleagues. The idea was to make use of all we learned about digital circuits so we implemented a RAM, adders, various types of registers, finite state machines, and we all went through the nightmare of 'looking for the bug' and creating testbenchs for vhdl. But more importantly that an innocent looking modification may break others' job.
By the way, knowing about the nuts and bolts of how are numbers stored in the bare metal gave me an appreciation for what are the differences between an idealized (almighty) function and the real poor man's version of it.
- Please state the commits and patches you already contributed to Octave.
No commits or patches done, although I have successfully built it from source and written a bare bones matrix cosine function cosm.m based on Higham and Smith.
F: Feeling fine
Please describe (in short) your experience with the following tools:
- IRC and mailing lists
In spite of not having much experience with relay chat, I have been following several mailing lists in the free software world (fsf, fsfe, octave).
- Mercurial or other source code management systems
Basic usage of git and hg only (just in order to clone-update-build).
- Mediawiki or other wiki software
Definitely not an expert, but comfortable editing with it.
- make, gcc, gdb or other development tools
My experience with these is limited, having used GCC, make and autotools to build free software packages.
- What will make you actively stay in our community after this GSoC is over?
If I stick to my discipline, I will be using numerics for a long time and octave will be one of my main tools, to which I hope to contribute.
O: Only out of interest
- Did you ever hear about Octave before?
I have been using octave for a long time now, in varying courses such as Linear Algebra, Control Theory, Quantum Physics.
- What was the first question concerning Octave you could not find an answer to rather quickly?
Finding a way to systematically test a functions was not very easy. I also missed a guide on the differences with Matlab (specially syntax), which should be crucial as Matlab is sadly more promoted than Octave in universities and people discover the latter when looking for an alternative to the former, just like me.
P: Prerequisites
I have access as root to two laptops running Debian Gnu/Linux and Ubuntu, with a reasonably working internet connection during the whole day.
S: Self-assessment
- Please describe how useful criticism looks from your point of view as committing student.
I am a person that prefers to listen than to talk, but let me put it this way: I would rather be told straight away that something isn't working than end up in a misunderstanding, even if they need to be a little rude with me.
- How autonomous are you when developing?
I often start with a proof of concept and 'see how it turns out' is usually my approach to a problem, I never start doing it blindly or without a clue, though.
Y: Your task
The project I intend to do is Improve logm, sqrtm, funm; its aim is to improve the existing implementations of Matrix Functions in Octave based on the algorithms developed by a team lead by Prof. Higham (project entitled Numerical Analysis of Matrix Functions, NAMF) at the University of Manchester. At this point in time, in Octave there are the following: expm makes use of Padé approximant, logm uses a Schur-Parlett algorithm, and sqrtm using a variant of the algorithm in A New sqrtm for MATLAB[1]. On the other hand, in Octave-Forge there are funm and trigonometric and hyperbolic matrix functions. For a general survey-introduction to matrix functions (or matrix computation in general) refer to Golub & Van Loan[2].
I believe this is of interest to Gnu Octave first, due to the goal of overall MATLAB compatibility and second, because more and more systems are being described by a matrix equation lately.
Upon completion Gnu Octave should have a working funm based on the Schur-Parlett algorithms by Higham et al., that calls to specific matrix functions if these have an instance of their own: expm, logm, sqrtm etc.
Update: Part of the work is already done by Prof. N.J. Higham and is available under a GPLv3+ license: The Matrix Function Toolbox[3] which is closely related to the book by the same author[4]. A toolbox for matrix computations[4] (The Matrix Computation Toolbox) is also provided by the same author, under the same license. Finally, a funm function is provided in the page of the NAMF project under GPLv3+. One might suggest that there is still room for improvement; because as Marco Caliari noted the toolboxes are from 2008. A review of the literature needs to be done in order to use more recent algorithms when writing the new functions.
TENTATIVE TIMELINE
weeks 1-2
The start should be soft for I am having the finals in this period. At this point the list of algorithms to be used must be completely defined; that is, a final review of the literature will be done. For that purpose, the access to my university's bibliographical will be used.
Work with the toolboxes starts here.
weeks 3-4
Milestone 0: the toolboxes are working and packaged.
weeks 4-7
funm
Milestone 1: general purpose funm based on a Schur-Parlett algorithm)
weeks 8-9
expm and logm
weeks 10-11
sqrtm
p-th roots (p = integer) algorithm [?]
Milestone 2: logm, sqrtm, funm
weeks 12
implement trigonometric functions (I submitted a snippet of cosm to the mailing list, updated)
weeks 13-14
Pencils down. Run tests on the Matrix Functions and write/review their documentation.
TENTATIVE TIMELINE(previous)
It appears to me that the natural task division for this project is to first improve and test funm (this should be ready by the mid-term), and then get to the rest, picking logm first. weeks 1-2
The start should be soft for I am having the finals in this period. At this point the list of algorithms to be used must be completely defined.
weeks 3-7
funm (Milestone 1: general purpose funm based on a Schur-Parlett algorithm)
weeks 8-9
expm and logm
weeks 10-11
sqrtm p-th roots (p = integer) algorithm [?] (Milestone 2: logm, sqrtm, funm)
weeks 12
implement trigonometric functions (I submitted a snippet of cosm to the mailing list)
weeks 13-14
Pencils down. Run tests on the Matrix Functions and write/review their documentation.
[1] N.J. Higham. A New sqrtm for MATLAB. Numerical Analysis Report No. 336, Manchester Centre for Computational Mathematics, Manchester, England, January 1999.
[2] G.H. Golub and C.F. Van Loan. Matrix Computations, 4th Edition. The Johns Hopkins University Press, Baltimore, USA, 2013.
[3] N. J. Higham. The Matrix Function Toolbox. http://www.ma.man.ac.uk/~higham/mftoolbox
[4] N. J. Higham, Functions of Matrices: Theory and Computation, Society for Industrial and Applied Mathematics, Philadelphia, USA, 2008.
[5] N. J. Higham. The Matrix Computation Toolbox. http://www.ma.man.ac.uk/~higham/mctoolbox
This section is being reworked in a sandbox, the other wiki. Will add it here when done. Note that the previous is still going though changes.
Z: submitted proposal
(note: personal contact info. has been removed and the info. below is has mostly been extracted form above)
Title: Improving Matrix Functions in GNU Octave
Organization: GNU Project
Abstract: Matrix Functions are widely used to describe dynamical systems, this project aims at implementing those based on algorithms developed at the University of Manchester by Dr. Philip Davies and Matthew Smith, led by Prof. Nick Higham. Upon completion Gnu Octave should have a reliable working funm based on the Schur-Parlett algorithms by Higham et al., that calls to specific matrix functions if these have an instance of their own.
Additional info: http://wiki.octave.org/User:Antonio_Pino
I am Antonio Pino Robles, and if selected I will be working on the project for about 6-8 hours a day (8:00-16:00 UTC+0), with the exception of the first two weeks when I will be taking my final exams. Therefore, mine will be a soft start. As for IRC, I may be contacted anytime during the day at AntonioPino in the #octave channel at freenode, via a mobile connection. That work plan shall not be altered, as I don't plan to travel this summer.
On the other hand, I am a fluent English speaker, hold a Cambridge CAE certificate and have taken several university courses taught entirely in the English language. Besides, I am an avid reader, with a passion for the English literature. What I mean is that I do not expect any problem communicating with the mentors or organizers, either in a written fashion or orally.
Benefit-Deliverables
The project's aim is to improve the existing implementations of Matrix Functions in Gnu Octave based on the algorithms developed by a team lead by Prof. Higham at the University of Manchester. At this point in time, in Gnu Octave there are the following: expm makes use of Padé approximant, logm uses a Schur-Parlett algorithm, and sqrtm using a variant of the algorithm in A New sqrtm for MATLAB[1]. On the other hand, in Octave-Forge there are funm and trigonometric and hyperbolic matrix functions.
I believe this is of interest to Gnu Octave first, due to the goal of overall MATLAB compatibility and second, because more and more systems are being described by a matrix equation lately.
funm should be ready by the mid-term. The rest will follow. The tentative time line can be read in the "Additional Info URL".
Qualification
Despite I consider myself an introvert person, I get most of my time when working with people, for example at study groups in university. I have been let down by harsh or ambiguous comments in the past, I still get an enormous benefit from properly filtered feedback, which I try to give when I think its necessary.
On one side, I hope to learn how programming is done outside of the sand boxed problems I have done until now (at Uni) while contributing to free software. Besides, I have been following the FSF for a long time, and been a free software user for longer, so this is in fact a way of giving back while receiving.
On the other, more and more dynamical systems can be approximated by a Matrix differential equation, and to have efficient Matrix Function implementations in Gnu Octave will help any free software user (including myself) approach these systems.
Finally, I have been writing m-scripts for two years now, to solve specific engineering problems, a strong python background. I order to fulfil it I will have to refresh my C++ knowledge by becoming familiar with the source tree, and get used to mercurial (hg), the version control system used to manage Gnu Octave.
[1] N.J. Higham. A New sqrtm for MATLAB. Numerical Analysis Report No. 336, Manchester Centre for Computational Mathematics, Manchester, England, January 1999.