k-Wave

1. subhajit
   

   Member
   Joined: Apr '20
   Posts: 3

   Hi,
   I was comparing the photoacoustic pressure waves computed by matalb and c++ binary(...-OMP.exe) under identical geometry and other conditions.The waveforms look almost similar except
   1) the amplitude of the one computed by Matlab is 1.341 times of the other one
   2)-OMP.exe computed waveform is little noisy.

   Grid: 128x128x128, Nt: 936, dt:1.818 ns, dx=dy=dz:10um. Medium: Heterogeneous, Max speed: 1540m/s

   I followed the example file in k-wave to prepare the hdf5 input file.

   It would be of great help if the forum members kindly extend suggestions/comment son the observed difference.
   Regards,

   - Subhajit

   Posted 3 years ago #

2. subhajit
   

   Member
   Joined: Apr '20
   Posts: 3

   Corrected as ...extend suggestions/comments on the...

   Posted 3 years ago #

3. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,640

   Hi Subhajit,

   We have a pretty large test suite to make sure the output from the MATLAB code and C++ code are the same, although there are a few options (e.g., medium.alpha_mode) that are not supported by the C++ code. Do you generate the input file using the 'SaveToDisk' option of kspaceFirstOrder3D?

   Brad

   Posted 3 years ago #

4. s148275
   

   Member
   Joined: Mar '20
   Posts: 15

   Hi Brad,

   Are you planning to include medium.alpha_mode and medium.alpha_sign in the C++ code?

   Best,
   Frank

   Posted 3 years ago #

5. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,640

   Hi Frank,

   It's not currently on the list - is it something that you need?

   Brad.

   Posted 3 years ago #

6. s148275
   

   Member
   Joined: Mar '20
   Posts: 15

   Hi Brad,

   Yes this would be really nice to have. The 3D simulation I am currently running takes much longer with the Maltlab solver. And since I am using medium.alpha_sign I cannot use the C++ code. I do not need the support for medium.alpha_mode.

   Best,
   Frank

   Posted 3 years ago #

7. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,640

   Ok - I'll generate an issue. I can't promise when it will be done though!

   Are you also using medium.alpha_filter?

   Posted 3 years ago #

8. s148275
   

   Member
   Joined: Mar '20
   Posts: 15

   No I am currently not using the alpha_filter.

   Posted 3 years ago #

RSS feed for this topic

Reply

You must log in to post.