vcs_internal.h

Go to the documentation of this file.

/**
 *  @file vcs_internal.h Internal declarations for the VCSnonideal package
 */
 
// This file is part of Cantera. See License.txt in the top-level directory or
// at https://cantera.org/license.txt for license and copyright information.
 
#ifndef _VCS_INTERNAL_H
#define _VCS_INTERNAL_H
 
#include "cantera/base/global.h"
namespace Cantera
{
//! define this %Cantera function to replace printf
/*!
 * We can replace this with printf easily
 */
#define plogf writelogf
 
//! Class to keep track of time and iterations
/*!
 * class keeps all of the counters together.
 */
class VCS_COUNTERS
{
public:
    //! Total number of iterations in the main loop
    //! of vcs_TP() to solve for thermo equilibrium
    int T_Its;
 
    //! Current number of iterations in the main loop
    //! of vcs_TP() to solve for thermo equilibrium
    int Its;
 
    //! Total number of optimizations of the components basis set done
    int T_Basis_Opts;
 
    //! number of optimizations of the components basis set done
    int Basis_Opts;
 
    //! Current number of times the initial thermo equilibrium estimator has
    //! been called
    int T_Calls_Inest;
 
    //! Current number of calls to vcs_TP
    int T_Calls_vcs_TP;
};
 
//! determine the l2 norm of a vector of doubles
/*!
 * @param vec vector of doubles
 * @returns   the l2 norm of the vector
 */
double vcs_l2norm(const vector<double>& vec);
 
//! Returns a const char string representing the type of the species given by
//! the first argument
/*!
 * @param speciesStatus  Species status integer representing the type
 *                       of the species.
 * @param length         Maximum length of the string to be returned.
 *                       Shorter values will yield abbreviated strings.
 *                       Defaults to a value of 100.
 */
const char* vcs_speciesType_string(int speciesStatus, int length = 100);
 
//! Simple routine to check whether two doubles are equal up to roundoff error
/*!
 * Currently it's set to check for 10 digits of relative accuracy.
 *
 * @param d1 first double
 * @param d2 second double
 *
 * @returns true if the doubles are "equal" and false otherwise
 */
bool vcs_doubleEqual(double d1, double d2);
}
 
#endif