Hi,
I was wondering if the modified Pennes bioheat equation (kWaveDiffusion) could be optimized in theory with compiled binaries in the same fashion as kspaceFirstOrder3D? If so, would it be feasible to use the C++ code for kspaceFirstOrder3D (e.g. KSpaceFirstOrder3DSolver.cpp) as a starting point for a C++-based kWaveDiffusion, or is the problem much more complicated than that?
Thank you!
Hi Vich,
Yes in principle these codes are closely related (e.g., both use an explicit time stepping scheme and derivatives calculated using FFTs). We have it in the development pipeline to write a C++ code, but haven't started on it yet, so it might be some time away!
Brad
I realize it must be really time-consuming to maintain and add new features to a project as big as k-Wave has become... I appreciate your hard work and that of everyone who has contributed to k-Wave: theoretical work, programming, documentation, physical validation, all of that is extremely valuable to us. Thank you!
Hi,
I was wondering if nowadays there is an available C++_based KWaveDiffusion.
Congratulations for this great toolbox! Thank you!
Hi Davide,
It's in the development pipeline, but we won't have anything ready for release for a while I'm afraid. In general, the stability limit for the diffusion code is large (due to the k-space correction), so if you're trying to speed things up, you could conduct a convergence test to see how big you can make your time steps and still get accurate results.
Brad
You must log in to post.