Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fit to Wilhoit for H2 atom does not converge #128

Closed
shamelmerchant opened this issue Jul 23, 2013 · 2 comments
Closed

Fit to Wilhoit for H2 atom does not converge #128

shamelmerchant opened this issue Jul 23, 2013 · 2 comments

Comments

@shamelmerchant
Copy link
Contributor

I get this error when I try to compute the thermo for H2 atom in cantherm ....

Generating Wilhoit thermo model for H2...
c:\Python27\lib\site-packages\scipy\optimize\optimize.py:988: RuntimeWarning: in
valid value encountered in absolute
q = abs(q)
c:\Python27\lib\site-packages\scipy\optimize\optimize.py:993: RuntimeWarning: in
valid value encountered in absolute
if ( (abs(p) < abs(0.5_q_r)) and (p > q*(a-xf)) and \

The line where it is failing is
scipy.optimize.fminbound(self.__residual, 300.0, 3000.0, args=(Tdata, Cpdata, Cp0, CpInf, H298, S298))

Currently I dont know what is the best way to fix this.

@rwest
Copy link
Member

rwest commented Jul 23, 2013

Are p q and r moments of inertia by any chance? H2 is linear, and I guess that abs(0) is invalid.
My two cents, without looking at any code.

@rwest
Copy link
Member

rwest commented Aug 30, 2013

Closed by 2230912, I think.
(@shamelmercant, please reopen if I'm wrong)

@rwest rwest closed this as completed Aug 30, 2013
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants