doublereal MineralEQ3::LookupGe(const std::string& elemName) { #ifdef OLDWAY int num = sizeof(geDataTable) / sizeof(struct GeData); string s3 = elemName.substr(0,3); for (int i = 0; i < num; i++) { //if (!std::strncmp(elemName.c_str(), aWTable[i].name, 3)) { if (s3 == geDataTable[i].name) { return (geDataTable[i].GeValue); } } throw CanteraError("LookupGe", "element " + s + " not found"); return -1.0; #else int iE = elementIndex(elemName); if (iE < 0) { throw CanteraError("PDSS_HKFT::LookupGe", "element " + elemName + " not found"); } doublereal geValue = entropyElement298(iE); if (geValue == ENTROPY298_UNKNOWN) { throw CanteraError("PDSS_HKFT::LookupGe", "element " + elemName + " doesn not have a supplied entropy298"); } geValue *= (-298.15); return geValue; #endif }
doublereal MineralEQ3::LookupGe(const std::string& elemName) { size_t iE = elementIndex(elemName); if (iE == npos) { throw CanteraError("PDSS_HKFT::LookupGe", "element " + elemName + " not found"); } doublereal geValue = entropyElement298(iE); if (geValue == ENTROPY298_UNKNOWN) { throw CanteraError("PDSS_HKFT::LookupGe", "element " + elemName + " does not have a supplied entropy298"); } geValue *= (-298.15); return geValue; }