http://www.k-wave.org/forum/topic/comparing-my-own-fdtd-code-with-k-wave Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 00:41:30 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/comparing-my-own-fdtd-code-with-k-wave#post-8085 Wed, 17 Mar 2021 02:31:57 +0000 rainbowwhale 8085@http://www.k-wave.org/forum/ <p>k-Wave use staggered grid in both space and time.<br /> May be that causes those little delay.<br /> For time varying source, did you compared 'dirichlet' mode? </p> http://www.k-wave.org/forum/topic/comparing-my-own-fdtd-code-with-k-wave#post-8081 Tue, 09 Mar 2021 16:12:50 +0000 ZekuiJia 8081@http://www.k-wave.org/forum/ <p>In addition， the comparison when using initial pressure distribution is similar, the difference is more like a numerical error. However, when I use time varying pressure source, the result of kwave is even largerer, usually dozens of times of my own result. So I think there must be some kind of coefficient missing for my result. Any idea what could that be?</p> <p>Thanks</p> <p>Zekui </p> http://www.k-wave.org/forum/topic/comparing-my-own-fdtd-code-with-k-wave#post-8080 Mon, 08 Mar 2021 13:53:02 +0000 ZekuiJia 8080@http://www.k-wave.org/forum/ <p>Hi,</p> <p>Thanks for sharing this amazing toolbox.</p> <p>Recently I am writing my own 2D FDTD code. It's still in an early stage, so the code is very simple. I tested a simple case with my code and k-wave. I used the same initial pressure distribution for both codes. I controlled the simulation time so that the wave won't reach the boundary since my code doesn't support PML right now. The two results are very close but not exactly the same. The magnitude of the k-wave result tends to be larger and the whole result seems to be a little delayed comparing to my result. The grid and time spacing I chose should be fine for this simulation.</p> <p>In my opinion, even my code is very simple, they should handle these simple cases. I guess I am wrong. So I want to ask where my problem might be. </p> <p><a href="https://imgtu.com/i/6l1Inx"><img src="https://s3.ax1x.com/2021/03/08/6l1Inx.png" alt="6l1Inx.png" /></a></p> <p>Thanks</p> <p>Zekui </p>