k-Wave

1. SergioGwirc
   

   Member
   Joined: Aug '13
   Posts: 3

   Dear Sir

   I intalled k-wave and trying to use with GPUmat and NVIDIA CUDA GTX 660.
   When I run the example "example_us_bmode_linear_transducer" after the GPUstart it gives an error and I not understand why. Could you help me?
   The command windows is as follows:

   Starting GPU
   - GPUmat version: 0.280
   - Required CUDA version: 5.0
   There is 1 device supporting CUDA
   CUDA Driver Version: 5.0
   CUDA Runtime Version: 5.0

   Device 0: "GeForce GTX 660"
   CUDA Capability Major revision number: 3
   CUDA Capability Minor revision number: 0
   Total amount of global memory: 2147483648 bytes
   - CUDA compute capability 3.0
   ...done
   - Loading module EXAMPLES_CODEOPT
   - Loading module EXAMPLES_NUMERICS
   -> numerics30.cubin
   - Loading module NUMERICS
   -> numerics30.cubin
   - Loading module RAND
   >> example_us_bmode_linear_transducer

   k-Wave Transducer Properties
   transducer position: [1 22 42]
   transducer width: 11.8519mm (64 grid points)
   number of elements: 32
   number of active elements: 32 (elements 1 to 32)
   element width: 370.3704um (2 grid points)
   element spacing (kerf): 0ym (0 grid points)
   element pitch: 370.3704um (2 grid points)
   element length: 4.4444mm (24 grid points)
   sound speed: 1540m/s
   focus distance: 20mm
   elevation focus distance: 19mm
   steering angle: 0 degrees
   steering angle max: auto degrees
   Computing scan line 1 of 96
   Running k-Wave simulation...
   start time: 28-Nov-2013 14:08:11
   reference sound speed: 1600m/s
   prepending transducer.input_signal with 17 leading zeros
   appending transducer.input_signal with 17 trailing zeros
   dt: 36.075ns, t_end: 57.1429us, time steps: 1585
   input grid size: 216 by 108 by 108 grid points (40 by 20 by 20mm)
   maximum supported frequency: 3.78MHz
   expanding computational grid...
   computational grid size: 256 by 128 by 128 grid points
   casting variables to kWaveGPUsingle type...
   precomputation completed in 4.5368s
   starting time loop...
   Undefined function 'ifftn' for input arguments of type 'GPUsingle'.

   Error in kspaceFirstOrder3D (line 822)
   ux_sgx = bsxfun(@times, pml_x_sgx, ...

   Error in example_us_bmode_linear_transducer (line 240)
   sensor_data = kspaceFirstOrder3D(kgrid, medium, transducer, transducer, input_args{:});

   Posted 10 years ago #

2. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi SergioGwirc

   If using the GPUmat toolbox with k-Wave, the functions fftn, ifftn, and bsxfun from the other/GPUmat" folder within the k-Wave Toolbox should be renamed to have .m extensions and placed in the /@GPUtype" folder within the GPUmat toolbox directory.

   Brad.

   Posted 10 years ago #

3. SergioGwirc
   

   Member
   Joined: Aug '13
   Posts: 3

   Hi Brad

   Thank you for your answer.

   I renamed the functions and placed in /@GPUtype folder when I installed the GPUmat, but the errors remain. I followed a bit the program and the problem seems to be with bsxfun in kspaceFirstOrder3D (line 822). What I don't understand is that the program stops when bsxfun is used for the third time in line 822.

   Sergio

   Posted 10 years ago #

4. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   Hi Sergio,

   If you cast a 3D matrix to the GPU from the command line, do the fftn and ifftn functions work? i.e., something like:

   x = GPUsingle(rand(64, 64, 64));
   y = fftn(x);

   Brad.

   Posted 10 years ago #

5. SergioGwirc
   

   Member
   Joined: Aug '13
   Posts: 3

   Hi Brad

   It does not work.
   That means that I made a mistake when I have renamed functions fftn, ifftn, and bsxfun and placed in the folder /@GPUtype ?

   Sergio

   Posted 10 years ago #

6. Bradley Treeby
   

   Developer
   Joined: Mar '10
   Posts: 1,643

   It would seem that way. The functions need to be renamed from .txt to .m and placed in the root @GPUtype folder in the GPUmat toolbox. The GPUmat toolbox also needs to be added to the MATLAB path. I hope you manage to get it to work.

   Brad.

   Posted 10 years ago #

RSS feed for this topic

Reply

You must log in to post.