pdepe

[Mark Wistrom]

Hi,

 

I’ve been working on a feature that has the same functionality as “pdepe” in Matlab. Basically this: http://wiki.octave.org/GSoC_Project_Ideas#Implement_solver_for_initial-boundary_value_problems_for_parabolic-elliptic_PDEs_in_1D . My implementation is based upon the paper by Skeel and Berzins, just as the Matlab function is. The paper: http://dl.acm.org/citation.cfm?id=78481

 

I’ve tested against the examples from Matlab and some (not all seem to work) of the examples in the Skeel and Berzins paper. The code constructs a function based on the geometry and whether or not there is a singularity in the domain that approximates the solution at discrete points in one dimension . An ODE/DAE uses the function to integrate the solution over time in a “Method of Lines” fashion. I’m using daspk.

 

At this point, I’m not too sure what to do. I need to finish the documentation, and make one more optimization. The work seems to be winding down, but I’m sure there are several small things needed to integrate what I have. I would like to check in/submit code review/get feedback etc. relatively soon so that I can finish this up. What is next?

 

Thanks,

Mark