Hi every one.
I wrote FMM code in C++ by myself but for large number of grids, I think it take too much time. Would you please tell me know that typically how long it should take for running the efficient code in a mesh that contains 1 million grids.(100*100*100).
Thanks a lot for sharing your code! One thing though, quite basic stuff I am afraid, but hopefully you can help: I cannot compile rk4, I get the error below.
Any help will be more than welcomed!
Building with 'Xcode with Clang'.
Error using mex
Undefined symbols for architecture x86_64:
"_interpgrad3d", referenced from:
_RK4STEP_3D in rk4.o
ld: symbol(s) not found for architecture x86_64
clang: error: linker command failed with exit code 1 (use -v to see invocation)
Error in compile_c_files (line 12)
Thank you for sharing a seemingly great toolbox. But like some users, I also failed to compile mex files. Although I don't intend to use "perform_front_propagation_anisotropic", it calls an error:
Unable to complete successfully.
Error in compile_mex (line 7)
leaving perform_front_propagation_mesh, which I wish to use, not compiled.
So I just commented the line 7 in "compile_mex.m", which reads "% mex mex/anisotropic-fm//perform_front_propagation_anisotropic.cpp".
Then compile_mex.m returns a bunch of warnings but it ends without errors. and the functions seem to work, I guess. I wish the author update/modify the compile_mex.m to bypass such errors, and hope this helps others having this error, although I'm not sure if this introduces some insane behaviors or not.
06 Nov 2014
Accurate Fast Marching
Multistencils second order Fast Marching 2D and 3D including rk4 shortest path and skeletonize
Seems to have error for constant speed field tracking.
DistanceMap= msfm(SpeedImage, SourcePoint);
figure, imshow(DistanceMap,[0 3400])
hold on, plot(ShortestLine(:,2),ShortestLine(:,1),'r')