search

HomerReid / scuff-em

A comprehensive and full-featured computational physics suite for boundary-element analysis of electromagnetic scattering, fluctuation-induced phenomena (Casimir forces and radiative heat transfer), nanophotonics, RF device engineering, electrostatics, and more. Includes a core library with C++ and python APIs as well as many command-line applications.
http://www.homerreid.com/scuff-em
GNU General Public License v2.0
125 stars 50 forks source link

[scuff-intergrate]error: Interp1D: grid must have 2 or more data points (aborting) #180

Open Zhao1236 opened 5 years ago

Zhao1236 commented 5 years ago

Hi, I'm calculating the radiative heat transfer between two objects using "scuff-em". After running "scuff-neq", I get the datafile "x.SIFlux.EMTPFT". Then I'm trying to specify the temperatures and do the integrating using scuff-integrate --Datafile x.SIFlux.EMTPFT --Temperature 0 300 --Temperature 1 800 --Temperature 2 300 but there is an error: Interp1D: grid must have 2 or more data points (aborting)

Maybe is there any commands I used wrong ?

Thanks in advance!
Zhao
HomerReid commented 5 years ago

This looks to me like your .EMTPFT file only has a single frequency point, so the integrator can't evaluate the integral over frequencies. You want to specify an --OmegaFile to scuff-neq that lists multiple sample frequencies (quadrature points for the numerical quadrature over frequency). The more points you give, the more accurately it will be able to estimate the frequency integral.

Note that the documentation for scuff-neq and scuff-integrate is a little incomplete at the moment, because I am in the midst of making major improvements to scuff-neq and other codes in the suite and haven't finished updating the documentation. Stay tuned for some revisions in the near future.

Feel free to ask more questions, or close the issue if satisfied.