example1.log

./example1.out ex.dat 

 "mfb.txt" No.00 -- focused beam --
medium density                       [kg/m^3] :            1000
speed of sound in the medium            [m/s] :            1500
frequency                                [Hz] :          100000
sound pressure amplitude                 [Pa] :     10+      0I
radius of focused transducer              [m] :           0.045
curvature radius of the transducer        [m] :           0.072
x-component of translation vector         [m] :               0
y-component of translation vector         [m] :               0
z-component of translation vector         [m] :               0
rotation parameter theta                [rad] :               0
rotation parameter psi                  [rad] :               0
sampling number for Gauss-Legendre quadrature :              64
-- additional information --
wavelength                                [m] :           0.015
energy passing through a plane orthogonal to the acoustic axis [W]
                                              : 1.4965495613794e-06

--- medium data ---
medium data file name      = medium.txt
domain id  1 rho0 [kg/m^3] =     970
             c0      [m/s] =    1004

--- mesh data ---
mesh data file name    = sphere_m1_05z.msh
node number            =      602
defined element number =     1200


sound pressure and particle velocity at r=( -0.03, 0.01, 0.002 )
type=0 setting ( 4 point Gauss-Legendre )
p  =  7.83664515499439e-02 +6.96197074165065e-02 I 
vx = -5.40468252235035e-08 -2.49067291806228e-08 I 
vy =  1.80152705159352e-08 +8.30153026037150e-09 I 
vz =  6.81407244242909e-09 +2.32283134469617e-08 I 
sound pressure and particle velocity at r=( -0.03, 0.01, 0.002 )
type=1 setting ( 9 or 7 point Gauss-Legendre )
p  =  7.83664726333253e-02 +6.96197058340226e-02 I 
vx = -5.40468313392011e-08 -2.49067002728027e-08 I 
vy =  1.80152761596773e-08 +8.30154469622858e-09 I 
vz =  6.81412629103293e-09 +2.32282955478716e-08 I 

Radiation force and torque
type=0 setting ( 4 point Gauss-Legendre ) 
F = ( 1.82208380178588e-29,-1.81530452838163e-29, 9.13571182459779e-15)
N = ( 2.24370838521113e-32, 6.23519819495368e-33,-2.14799493318978e-32), center of rotation ( 0, 0, 0)
type=1 setting ( 9 or 7 point Gauss-Legendre ) 
F = ( 1.80020523761764e-29,-1.81461119360165e-29, 9.14470213627012e-15)
N = ( 2.98519141380022e-32, 1.50222535662204e-32,-1.90179939635356e-32), center of rotation ( 0, 0, 0)