cannot get the power transfer between two SiO2Spheres

Hi, Homer

The results are all 0.000000e+00 when I run the example about the power transfer between two SiO2Spheres, http://homerreid.github.io/buff-em-documentation/examples/SiO2Spheres/ The results are same for methods of --OPFT, --EMTPFT, and --MomentPFT . RunScript

#!/bin/bash

ARGS=""
ARGS="${ARGS} --geometry SiO2Spheres_677.buffgeo"
ARGS="${ARGS} --OmegaFile OmegaFile"
ARGS="${ARGS} --EMTPFT"
ARGS="${ARGS} --PAbs"

buff-neq ${ARGS}

results

# buff-neq run on jia-Precision-Tower-7810 (12/09/16::21:48:36)
# data file columns: 
# 1 omega 
# 2 transform tag
# 3 (sourceObject,destObject) 
# 4 absorbed power spectral density
# 5 radiated power spectral density
# 6 x-force flux spectral density
# 7 y-force flux spectral density
# 8 z-force flux spectral density
# 9 x-torque flux spectral density
# 10 y-torque flux spectral density
# 11 z-torque flux spectral density
DEFAULT 1.000000e-02 11 -1.670797e-09 1.670797e-09 2.965621e-11 -1.909764e-11 9.200753e-08 -1.000000e-02 2.561126e-10 -1.000000e-04 
DEFAULT 1.000000e-02 12 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 
DEFAULT 1.000000e-02 21 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 
DEFAULT 1.000000e-02 22 -1.670758e-09 1.670758e-09 1.992007e-14 -5.163125e-14 2.739386e-14 -2.334779e-12 5.712766e-13 1.269466e-12 
DEFAULT 1.258925e-02 11 -4.201814e-09 4.201814e-09 2.948099e-11 -1.836088e-11 9.224337e-08 -1.258925e-02 2.572150e-10 -1.584892e-04 
DEFAULT 1.258925e-02 12 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00 0.000000e+00

I hve looked through the codes, but i still do not know the reasons. Is it a problem about SCUFF_PFT_EMT(https://github.com/HomerReid/buff-em/blob/master/src/applications/buff-neq/CreateBNEQData.cc#L177)? Or anything else?

Thank you.

Zhihao

HomerReid / buff-em

cannot get the power transfer between two SiO2Spheres #5