Code covered by the BSD License
by Qun HAN
07 Jun 2009
Faster version of the MATLAB function bvp5c
| Watch this File
The inefficient sparse matrix operation in colloc_JacODE_region is optimized according to the method intruduced in Loren's blog (http://blogs.mathworks.com/loren/2007/03/01/creating-sparse-finite-element-matrices-in-matlab/). More than 10x efficiency improvement can be expected for large BVP systems.
The mex file SPARSE2.MEXW32 is compiled from SuiteSparse(http://www.cise.ufl.edu/research/sparse/SuiteSparse/). It is about 2x faster than SPARSE. If you encounters problem with it, just comment it out in bvp5c2 and use the MATLAB function sparse.
USAGE: Unzip to MATLABROOT\toolbox\matlab\funfun
For a benchmark see: http://junziyang.spaces.live.com/blog/cns!909609E321D6E37A!498.entry
Contact us at files@mathworks.com