Re: Level set methods fail on non-uniform grids under FiPy 3.1

2014-07-30 Thread Daniel Wheeler
On Tue, Jul 29, 2014 at 10:57 AM, Jason Furtney jkfurt...@gmail.com wrote: Daniel and Robert, Any thoughts on a way forward? Thanks Jason for putting your ideas down. I agree with everything you have suggested. I certainly think you/we should tackle non-uniform grids before looking into

Re: Level set methods fail on non-uniform grids under FiPy 3.1

2014-07-29 Thread Jason Furtney
Daniel and Robert, I looked more closely at adapting scikit-fmm for non-uniform (but still structured) grids. For the first-order method this would be fairly simple as the stencil and quadratic used in the point update are basically the same. It would be a matter of programming to create a new

Re: Level set methods fail on non-uniform grids under FiPy 3.1

2014-07-25 Thread Daniel Wheeler
On Tue, Jul 22, 2014 at 5:35 PM, Jason Furtney jkfurt...@gmail.com wrote: Also, the Sethian book has the formulation for the fast marching method on unstructured grids which we could implement. Jason, I would be very interested in having the implemented (it would help my research

Re: Level set methods fail on non-uniform grids under FiPy 3.1

2014-07-22 Thread Daniel Wheeler
On Tue, Jul 22, 2014 at 12:52 AM, Robert R. Howell rhow...@uwyo.edu wrote: I've been using FiPy 3.0 level set methods for an application which requires non-uniform grids. I recently tried upgrading to version 3.1 but now get a variety of error messages when I call calcDistanceFunction(). To

Re: Level set methods fail on non-uniform grids under FiPy 3.1

2014-07-22 Thread Robert R. Howell
Hi Daniel I just spent a little time looking through the C++ code in the heart of Scikit-fmm before I replied, and while it's not as overwhelming as I feared, modifying it to handle non-uniform grids is probably more than I'm going to be able to take on in the immediate future. On a more