-
Notifications
You must be signed in to change notification settings - Fork 62
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Exact representation of contour #87
Comments
Hi David, thanks for getting in touch with this. I need some time to digest this and do some testing to better understand what is going on. What you describe seems reasonable, let me have a look and comment back here soon. Thanks again |
Made a fork that makes it work: https://github.com/hamdav/scikit-fmm. I just did the algorithm part and renamed the previous initializeFrozen -> intializeFrozen2. So no option on which algorithm to use yet. I know python and I know c++, but I've never interfaced them before so I don't know how to do that. I'm sure it's not too hard though :) However, running the tests makes some of them fail, not sure if this is just because the distances changed or because they got worse... |
Thanks for this, is there a quick way I can see a diff of initializeFrozen and initializeFrozen2? I will be able to do a deeper dive on this soon. |
I don't think you would get much out of it. The conceptual diff is: First for loop (where i is marked frozen if phi[i]==0) is the same, the rest is different. There are of course similarities, both loop over the far values, both loop over the dimensions. The previous version also loop over the directions, though I skipped that loop and only go in one direction, so I only loop over each edge once. But I think the similarities are few enough that it's easier to just read it on it's own. I'd be happy to explain something in more detail if you want. |
Hey again! I am trying to compile skfmm on Windows, and I can't seem to get it to work. When I run
The output from the build is
|
It looks like the build completed. You can get that error (ModuleNotFoundError: No module named 'skfmm.cfmm') if you try to run python and import skfmm from the scikit-fmm source folder. In this case, just change to a different folder and it should work. I think in this case Python searches the local directory first. If you want to do some experiments locally, I think you can do |
Ah, that's a silly mistake by me, I was indeed stood in the source folder. Thanks :) And good to know about |
Hi!
When computing the distances to the 0 contour using this package, the returned array does not have exactly the same 0 contour. I am using it to keep track of an evolving interface, and so for me it is imperative that the contour is not moved when reinitializing. Trying to fix this might make the distances ever so slightly worse, however, they are not exact to begin with so I think that this is a reasonable tradeoff. It would also increase the complexity of the call to initializeFrozen, but probably not very significantly. Would you be open to a pull request for this? Maybe adding an option on which algorithm to use so as to not change the defaults? I have some ideas on how to compute the initialFrozen values...
The text was updated successfully, but these errors were encountered: