max/iphreeqc
1
0
Fork 0
mirror of https://git.gfz-potsdam.de/naaice/iphreeqc.git synced 2025-12-15 16:18:22 +01:00
Code Issues Packages Projects Releases Wiki Activity
iphreeqc/unit/TestInterface.cpp
Scott R Charlton 2c4422593b # IPQ_DLL_EXPORT int GetDumpOn(int id); 
#	IPQ_DLL_EXPORT int         GetDumpFileOn(int id);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetDumpOn(int id, int dump_on);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetDumpFileOn(int id, int dump_on);
s/etDumpOn/etDumpFileOn/g
s/etdumpon/etdumpfileon/g
s/ETDUMPON/ETDUMPFILEON/g
#	IPQ_DLL_EXPORT int         GetErrorOn(int id);
#	IPQ_DLL_EXPORT int         GetErrorFileOn(int id);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetErrorOn(int id, int error_on);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetErrorFileOn(int id, int error_on);
s/etErrorOn/etErrorFileOn/g
s/eterroron/eterrorfileon/g
s/ETERRORON/ETERRORFILEON/g
#	IPQ_DLL_EXPORT int         GetLogOn(int id);
#	IPQ_DLL_EXPORT int         GetLogFileOn(int id);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetLogOn(int id, int log_on);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetLogFileOn(int id, int log_on);
s/etLogOn/etLogFileOn/g
s/etlogon/etlogfileon/g
s/ETLOGON/ETLOGFILEON/g
#	IPQ_DLL_EXPORT int         GetOutputOn(int id);
#	IPQ_DLL_EXPORT int         GetOutputFileOn(int id);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetOutputOn(int id, int output_on);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetOutputFileOn(int id, int output_on);
s/etOutputOn/etOutputFileOn/g
s/etoutputon/etoutputfileon/g
s/ETOUTPUTON/ETOUTPUTFILEON/g
#	IPQ_DLL_EXPORT int         GetSelectedOutputOn(int id);
#	IPQ_DLL_EXPORT int         GetSelectedOutputFileOn(int id);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetSelectedOutputOn(int id, int sel_on);
#	IPQ_DLL_EXPORT IPQ_RESULT  SetSelectedOutputFileOn(int id, int sel_on);
s/etSelectedOutputOn/etSelectedOutputFileOn/g
s/etselectedoutputon/etselectedoutputfileon/g
s/ETSELECTEDOUTPUTON/ETSELECTEDOUTPUTFILEON/g
#	IPQ_DLL_EXPORT const char* GetDumpLine(int id, int n);
#	IPQ_DLL_EXPORT const char* GetDumpStringLine(int id, int n);
#	IPQ_DLL_EXPORT int         GetDumpLineCount(int id);
#	IPQ_DLL_EXPORT int         GetDumpStringLineCount(int id);
s/GetDumpLine/GetDumpStringLine/g
s/getdumpline/getdumpstringline/g
s/GETDUMPLINE/GETDUMPSTRINGLINE/g
#	IPQ_DLL_EXPORT const char* GetErrorLine(int id, int n);
#	IPQ_DLL_EXPORT const char* GetErrorStringLine(int id, int n);
#	IPQ_DLL_EXPORT int         GetErrorLineCount(int id);
#	IPQ_DLL_EXPORT int         GetErrorStringLineCount(int id);
s/GetErrorLine/GetErrorStringLine/g
s/geterrorline/geterrorstringline/g
s/GETERRORLINE/GETERRORSTRINGLINE/g
#	IPQ_DLL_EXPORT const char* GetWarningLine(int id, int n);
#	IPQ_DLL_EXPORT const char* GetWarningStringLine(int id, int n);
#	IPQ_DLL_EXPORT int         GetWarningLineCount(int id);
#	IPQ_DLL_EXPORT int         GetWarningStringLineCount(int id);
s/GetWarningLine/GetWarningStringLine/g
s/getwarningline/getwarningstringline/g
s/GETWARNINGLINE/GETWARNINGSTRINGLINE/g
#	IPQ_DLL_EXPORT void        OutputError(int id);
#	IPQ_DLL_EXPORT void        OutputErrorString(int id);
s/OutputError/OutputErrorString/g
s/outputerror/outputerrorstring/g
s/OUTPUTERROR/OUTPUTERRORSTRING/g
#	IPQ_DLL_EXPORT void        OutputLines(int id);
#	IPQ_DLL_EXPORT void        OutputAccumulatedLines(int id);
s/OutputLines/OutputAccumulatedLines/g
s/outputlines/outputaccumulatedlines/g
s/OUTPUTLINES/OUTPUTACCUMULATEDLINES/g
#	IPQ_DLL_EXPORT void        OutputWarning(int id);
#	IPQ_DLL_EXPORT void        OutputWarningString(int id);
s/OutputWarning/OutputWarningString/g
s/outputwarning/outputwarningstring/g
s/OUTPUTWARNING/OUTPUTWARNINGSTRING/g


git-svn-id: svn://136.177.114.72/svn_GW/IPhreeqc/trunk@4411 1feff8c3-07ed-0310-ac33-dd36852eb9cd
2010-05-18 01:59:33 +00:00

1803 lines
60 KiB
C++
Raw Blame History

#include "TestInterface.h"
#include "IPhreeqc.hpp"
#undef true
#undef false
#include "IPhreeqc.h"
#include "CVar.hxx"
#include <iostream>
#include <cmath>
#include <cfloat>
#if defined(_WIN32) || defined(__CYGWIN32__)
// DeleteFile defined in <windows.h>
#else
int DeleteFile(const char* szPathName);
#endif
bool FileExists(const char *szPathName);
VRESULT SOLUTION(double C, double Ca, double Na);
VRESULT EQUILIBRIUM_PHASES(const char* phase, double si, double amount);
VRESULT USER_PUNCH(const char* element, int max);
VRESULT SELECTED_OUTPUT(void);
VRESULT DUMP(void);
void TestOnOff(const char* FILENAME, int output_on, int error_on, int log_on, int selected_output_on, int dump_on);
TestInterface::TestInterface()
{
}
TestInterface::~TestInterface()
{
}
void TestInterface::TestLoadDatabase()
{
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
#if defined(_WIN32)
CPPUNIT_ASSERT_EQUAL(0, ::_fcloseall());
#endif
}
void TestInterface::TestLoadDatabaseString()
{
const char ex15_dat[] =
"SOLUTION_MASTER_SPECIES\n"
"C CO2 2.0 61.0173 12.0111\n"
"Cl Cl- 0.0 Cl 35.453\n"
"Co Co+2 0.0 58.93 58.93 \n"
"E e- 0.0 0.0 0.0\n"
"H H+ -1. 1.008 1.008\n"
"H(0) H2 0.0 1.008\n"
"H(1) H+ -1. 1.008\n"
"N NH4+ 0.0 14.0067 14.0067\n"
"Na Na+ 0.0 Na 22.9898\n"
"Nta Nta-3 3.0 1. 1.\n"
"O H2O 0.0 16.00 16.00\n"
"O(-2) H2O 0.0 18.016\n"
"O(0) O2 0.0 16.00\n"
"SOLUTION_SPECIES\n"
"2H2O = O2 + 4H+ + 4e- \n"
" log_k -86.08; -gamma 1e7 0.0\n"
"2 H+ + 2 e- = H2\n"
" log_k -3.15; -gamma 1e7 0.0\n"
"H+ = H+\n"
" log_k 0.0; -gamma 1e7 0.0\n"
"e- = e-\n"
" log_k 0.0; -gamma 1e7 0.0\n"
"H2O = H2O\n"
" log_k 0.0; -gamma 1e7 0.0\n"
"CO2 = CO2\n"
" log_k 0.0; -gamma 1e7 0.0\n"
"Na+ = Na+\n"
" log_k 0.0; -gamma 1e7 0.0\n"
"Cl- = Cl-\n"
" log_k 0.0; -gamma 1e7 0.0\n"
"Co+2 = Co+2\n"
" log_k 0.0; -gamma 1e7 0.0\n"
"NH4+ = NH4+\n"
" log_k 0.0; -gamma 1e7 0.0\n"
"Nta-3 = Nta-3\n"
" log_k 0.0; -gamma 1e7 0.0\n"
"Nta-3 + 3H+ = H3Nta\n"
" log_k 14.9; -gamma 1e7 0.0\n"
"Nta-3 + 2H+ = H2Nta-\n"
" log_k 13.3; -gamma 1e7 0.0\n"
"Nta-3 + H+ = HNta-2\n"
" log_k 10.3; -gamma 1e7 0.0\n"
"Nta-3 + Co+2 = CoNta-\n"
" log_k 11.7; -gamma 1e7 0.0\n"
"2 Nta-3 + Co+2 = CoNta2-4\n"
" log_k 14.5; -gamma 1e7 0.0\n"
"Nta-3 + Co+2 + H2O = CoOHNta-2 + H+\n"
" log_k 0.5; -gamma 1e7 0.0\n"
"Co+2 + H2O = CoOH+ + H+\n"
" log_k -9.7; -gamma 1e7 0.0\n"
"Co+2 + 2H2O = Co(OH)2 + 2H+\n"
" log_k -22.9; -gamma 1e7 0.0\n"
"Co+2 + 3H2O = Co(OH)3- + 3H+\n"
" log_k -31.5; -gamma 1e7 0.0\n"
"CO2 + H2O = HCO3- + H+\n"
" log_k -6.35; -gamma 1e7 0.0\n"
"CO2 + H2O = CO3-2 + 2H+\n"
" log_k -16.68; -gamma 1e7 0.0\n"
"NH4+ = NH3 + H+\n"
" log_k -9.3; -gamma 1e7 0.0\n"
"H2O = OH- + H+\n"
" log_k -14.0; -gamma 1e7 0.0\n"
"END\n";
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabaseString(ex15_dat));
}
void TestInterface::TestLoadDatabaseMissingFile()
{
CPPUNIT_ASSERT_EQUAL(false, ::FileExists("missing.file"));
CPPUNIT_ASSERT_EQUAL(1, ::LoadDatabase("missing.file"));
CPPUNIT_ASSERT_EQUAL(1, ::LoadDatabase("missing.file"));
CPPUNIT_ASSERT_EQUAL(1, ::LoadDatabase("missing.file"));
const char expected[] =
"ERROR: LoadDatabase: Unable to open:\"missing.file\".\n"
"Stopping.\n";
const char* err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL( std::string(expected), std::string(err) );
}
void TestInterface::TestLoadDatabaseWithErrors()
{
for (int i = 0; i < 5; ++i)
{
CPPUNIT_ASSERT_EQUAL(true, ::FileExists("missing_e.dat"));
CPPUNIT_ASSERT_EQUAL(6, ::LoadDatabase("missing_e.dat"));
static const char *expected =
"ERROR: Could not reduce equation to primary master species, CH4.\n"
"ERROR: Could not reduce equation to primary master species, Cu+.\n"
"ERROR: Could not reduce equation to primary master species, Fe+3.\n"
"ERROR: Could not reduce equation to primary master species, H2.\n"
"ERROR: Could not reduce equation to primary master species, Mn+3.\n"
"ERROR: Could not reduce equation to primary master species, NH4+.\n"
"ERROR: Could not reduce equation to primary master species, N2.\n"
"ERROR: Could not reduce equation to primary master species, NO2-.\n"
"ERROR: Could not reduce equation to primary master species, O2.\n"
"ERROR: Could not reduce equation to primary master species, HS-.\n"
"ERROR: Could not reduce equation to secondary master species, e-.\n"
"ERROR: Non-master species in secondary reaction, e-.\n"
"ERROR: No master species for element e.\n"
"ERROR: Could not find primary master species for e.\n"
"ERROR: No master species for element e.\n"
"ERROR: Could not reduce equation to secondary master species, Hausmannite.\n"
"ERROR: Could not reduce equation to secondary master species, Manganite.\n"
"ERROR: Could not reduce equation to secondary master species, Pyrite.\n"
"ERROR: Could not reduce equation to secondary master species, Pyrolusite.\n"
"ERROR: Could not reduce equation to secondary master species, Sulfur.\n"
"ERROR: e-, primary master species for E-, not defined.\n"
"ERROR: Calculations terminating due to input errors.\n"
"Stopping.\n";
const char* err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL( std::string(expected), std::string(err) );
}
}
void TestInterface::TestRun()
{
int files_on = 0;
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine("solution 12"));
::SetOutputFileOn(files_on);
::SetErrorFileOn(files_on);
::SetLogFileOn(files_on);
::SetSelectedOutputFileOn(files_on);
::SetDumpFileOn(files_on);
CPPUNIT_ASSERT_EQUAL(0, ::RunAccumulated());
}
void TestInterface::TestRunWithErrors()
{
const char dump_file[] = "error.inp";
// remove dump file if it exists
//
if (::FileExists(dump_file))
{
CPPUNIT_ASSERT(::DeleteFile(dump_file));
}
int files_on = 0;
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine("SOLUTION 1"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine(" pH 7"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine(" Na 1"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine("PHASES"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine(" Fix_H+"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine(" H+ = H+"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine(" log_k 0"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine("EQUILIBRIUM_PHASES"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine(" Fix_H+ -10 HCl 10"));
CPPUNIT_ASSERT_EQUAL(VR_OK, ::AccumulateLine("END"));
::SetOutputFileOn(files_on);
::SetErrorFileOn(files_on);
::SetLogFileOn(files_on);
::SetSelectedOutputFileOn(files_on);
::SetDumpFileOn(files_on);
CPPUNIT_ASSERT_EQUAL(1, ::RunAccumulated());
const char expected[] =
"ERROR: Numerical method failed on all combinations of convergence parameters\n"
"Stopping.\n";
const char* err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL( std::string(expected), std::string(err) );
CPPUNIT_ASSERT_EQUAL(true, ::FileExists(dump_file) );
CPPUNIT_ASSERT(::DeleteFile(dump_file));
}
void TestInterface::TestRunFile()
{
const char dump_file[] = "error.inp";
// remove dump file if it exists
//
if (::FileExists(dump_file))
{
CPPUNIT_ASSERT(::DeleteFile(dump_file));
}
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL(1, ::RunFile("conv_fail.in"));
const char expected[] =
"ERROR: Numerical method failed on all combinations of convergence parameters\n"
"Stopping.\n";
const char* err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL(std::string(expected), std::string(err));
CPPUNIT_ASSERT_EQUAL(true, ::FileExists(dump_file));
CPPUNIT_ASSERT(::DeleteFile(dump_file));
}
void TestInterface::TestRunString(void)
{
const char input[] =
"TITLE Example 1.--Add uranium and speciate seawater.\n"
"SOLUTION 1 SEAWATER FROM NORDSTROM ET AL. (1979)\n"
" units ppm\n"
" pH 8.22\n"
" pe 8.451\n"
" density 1.023\n"
" temp 25.0\n"
" redox O(0)/O(-2)\n"
" Ca 412.3\n"
" Mg 1291.8\n"
" Na 10768.0\n"
" K 399.1\n"
" Fe 0.002\n"
" Mn 0.0002 pe\n"
" Si 4.28\n"
" Cl 19353.0\n"
" Alkalinity 141.682 as HCO3\n"
" S(6) 2712.0\n"
" N(5) 0.29 gfw 62.0\n"
" N(-3) 0.03 as NH4\n"
" U 3.3 ppb N(5)/N(-3)\n"
" O(0) 1.0 O2(g) -0.7\n"
"SOLUTION_MASTER_SPECIES\n"
" U U+4 0.0 238.0290 238.0290\n"
" U(4) U+4 0.0 238.0290\n"
" U(5) UO2+ 0.0 238.0290\n"
" U(6) UO2+2 0.0 238.0290\n"
"SOLUTION_SPECIES\n"
" #primary master species for U\n"
" #is also secondary master species for U(4)\n"
" U+4 = U+4\n"
" log_k 0.0\n"
" U+4 + 4 H2O = U(OH)4 + 4 H+\n"
" log_k -8.538\n"
" delta_h 24.760 kcal\n"
" U+4 + 5 H2O = U(OH)5- + 5 H+\n"
" log_k -13.147\n"
" delta_h 27.580 kcal\n"
" #secondary master species for U(5)\n"
" U+4 + 2 H2O = UO2+ + 4 H+ + e-\n"
" log_k -6.432\n"
" delta_h 31.130 kcal\n"
" #secondary master species for U(6)\n"
" U+4 + 2 H2O = UO2+2 + 4 H+ + 2 e-\n"
" log_k -9.217\n"
" delta_h 34.430 kcal\n"
" UO2+2 + H2O = UO2OH+ + H+\n"
" log_k -5.782\n"
" delta_h 11.015 kcal\n"
" 2UO2+2 + 2H2O = (UO2)2(OH)2+2 + 2H+\n"
" log_k -5.626\n"
" delta_h -36.04 kcal\n"
" 3UO2+2 + 5H2O = (UO2)3(OH)5+ + 5H+\n"
" log_k -15.641\n"
" delta_h -44.27 kcal\n"
" UO2+2 + CO3-2 = UO2CO3\n"
" log_k 10.064\n"
" delta_h 0.84 kcal\n"
" UO2+2 + 2CO3-2 = UO2(CO3)2-2\n"
" log_k 16.977\n"
" delta_h 3.48 kcal\n"
" UO2+2 + 3CO3-2 = UO2(CO3)3-4\n"
" log_k 21.397\n"
" delta_h -8.78 kcal\n"
"PHASES\n"
" Uraninite\n"
" UO2 + 4 H+ = U+4 + 2 H2O\n"
" log_k -3.490\n"
" delta_h -18.630 kcal\n"
"END\n"
"\n";
if (::FileExists("phreeqc.out"))
{
CPPUNIT_ASSERT(::DeleteFile("phreeqc.out"));
}
CPPUNIT_ASSERT_EQUAL(false, ::FileExists("phreeqc.out"));
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
::SetOutputFileOn(1);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL(false, ::FileExists("phreeqc.out"));
CPPUNIT_ASSERT_EQUAL(0, ::RunString(input));
CPPUNIT_ASSERT_EQUAL(true, ::FileExists("phreeqc.out"));
if (::FileExists("phreeqc.out"))
{
CPPUNIT_ASSERT(::DeleteFile("phreeqc.out"));
}
}
void TestInterface::TestGetSelectedOutputRowCount()
{
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("llnl.dat"));
int max = 6;
CPPUNIT_ASSERT_EQUAL(VR_OK, SOLUTION(1.0, 1.0, 1.0));
CPPUNIT_ASSERT_EQUAL(VR_OK, EQUILIBRIUM_PHASES("calcite", 0.0, 0.010));
CPPUNIT_ASSERT_EQUAL(VR_OK, USER_PUNCH("Ca", max));
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL(0, ::RunAccumulated());
CPPUNIT_ASSERT_EQUAL(3, ::GetSelectedOutputRowCount()); // rows + header
}
void TestInterface::TestGetSelectedOutputValue()
{
int col;
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("llnl.dat"));
int max = 6;
CPPUNIT_ASSERT_EQUAL(VR_OK, SOLUTION(1.0, 1.0, 1.0));
CPPUNIT_ASSERT_EQUAL(VR_OK, EQUILIBRIUM_PHASES("calcite", 0.0, 0.010));
CPPUNIT_ASSERT_EQUAL(VR_OK, USER_PUNCH("Ca", max));
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL(0, ::RunAccumulated());
/*
EXPECTED selected.out:
sim state soln dist_x time step pH pe C Ca Na m_CO3-2 m_CaOH+ m_NaCO3- la_CO3-2 la_CaOH+ la_NaCO3- Calcite d_Calcite si_CO2(g) si_Siderite pressure total mol volume g_CO2(g) g_N2(g) k_Albite dk_Albite k_Pyrite dk_Pyrite s_CaSO4 s_SrSO4 1.name 1.type 1.moles 2.name 2.type 2.moles 3.name 3.type 3.moles 4.name 4.type 4.moles 5.name 5.type 5.moles 6.name 6.type 6.moles
1 i_soln 1 -99 -99 -99 7 4 1.0000e-003 1.0000e-003 1.0000e-003 4.2975e-007 1.1819e-009 1.1881e-009 -6.4686e+000 -8.9530e+000 -8.9507e+000 0.0000e+000 0.0000e+000 -2.2870 -999.9990 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 Ca+2 aq 9.9178e-004 CaHCO3+ aq 7.5980e-006 CaCO3 aq 6.2155e-007 CaOH+ aq 1.1819e-009
1 react 1 -99 0 1 7.86135 10.18 1.1556e-003 1.1556e-003 1.0000e-003 4.2718e-006 9.7385e-009 1.1620e-008 -5.4781e+000 -8.0388e+000 -7.9621e+000 9.8444e-003 -1.5555e-004 -3.0192 -999.9990 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 0.0000e+000 calcite equi 9.8444e-003 Ca+2 aq 1.1371e-003 CaHCO3+ aq 1.1598e-005 CaCO3 aq 6.8668e-006 CaOH+ aq 9.7385e-009
*/
CVar v;
CPPUNIT_ASSERT_EQUAL(VR_INVALIDROW, ::GetSelectedOutputValue(-1, 0, &v));
CPPUNIT_ASSERT_EQUAL(TT_ERROR, v.type);
CPPUNIT_ASSERT_EQUAL(VR_INVALIDROW, v.vresult);
CPPUNIT_ASSERT_EQUAL(VR_INVALIDROW, ::GetSelectedOutputValue(::GetSelectedOutputRowCount(), 0, &v));
CPPUNIT_ASSERT_EQUAL(TT_ERROR, v.type);
CPPUNIT_ASSERT_EQUAL(VR_INVALIDROW, v.vresult);
CPPUNIT_ASSERT_EQUAL(VR_INVALIDCOL, ::GetSelectedOutputValue(0, -1, &v));
CPPUNIT_ASSERT_EQUAL(TT_ERROR, v.type);
CPPUNIT_ASSERT_EQUAL(VR_INVALIDCOL, v.vresult);
CPPUNIT_ASSERT_EQUAL(VR_INVALIDCOL, ::GetSelectedOutputValue(0, ::GetSelectedOutputColumnCount(), &v));
CPPUNIT_ASSERT_EQUAL(TT_ERROR, v.type);
CPPUNIT_ASSERT_EQUAL(VR_INVALIDCOL, v.vresult);
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, 0, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("sim"), std::string(v.sVal));
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, 1, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("state"), std::string(v.sVal));
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, 2, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("soln"), std::string(v.sVal));
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, 3, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("dist_x"), std::string(v.sVal));
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, 4, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("time"), std::string(v.sVal));
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, 5, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("step"), std::string(v.sVal));
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, 6, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("pH"), std::string(v.sVal));
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, 7, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("pe"), std::string(v.sVal));
col = 7;
// -totals C Ca Na
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("C(mol/kgw)"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("Ca(mol/kgw)"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("Na(mol/kgw)"), std::string(v.sVal));
// -molalities CO3-2 CaOH+ NaCO3-
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("m_CO3-2(mol/kgw)"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("m_CaOH+(mol/kgw)"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("m_NaCO3-(mol/kgw)"), std::string(v.sVal));
// -activities CO3-2 CaOH+ NaCO3-
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("la_CO3-2"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("la_CaOH+"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("la_NaCO3-"), std::string(v.sVal));
// -equilibrium_phases Calcite
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("Calcite"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("d_Calcite"), std::string(v.sVal));
// -saturation_indices CO2(g) Siderite
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("si_CO2(g)"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("si_Siderite"), std::string(v.sVal));
// -gases CO2(g) N2(g)
// pressure "total mol" volume g_CO2(g) g_N2(g)
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("pressure"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("total mol"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("volume"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("g_CO2(g)"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("g_N2(g)"), std::string(v.sVal));
// -kinetic_reactants Albite Pyrite
// k_Albite dk_Albite k_Pyrite dk_Pyrite
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("k_Albite"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("dk_Albite"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("k_Pyrite"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("dk_Pyrite"), std::string(v.sVal));
// -solid_solutions CaSO4 SrSO4
// s_CaSO4 s_SrSO4
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("s_CaSO4"), std::string(v.sVal));
++col;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("s_SrSO4"), std::string(v.sVal));
for (int i = 0; i < max; ++i)
{
std::ostringstream oss1, oss2, oss3;
// 1.name
//
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col + 1 + (i*3) , &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
oss1 << i + 1 << ".name";
CPPUNIT_ASSERT_EQUAL( oss1.str(), std::string(v.sVal));
// 1.type
//
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col + 2 + (i*3), &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
oss2 << i + 1 << ".type";
CPPUNIT_ASSERT_EQUAL( oss2.str(), std::string(v.sVal));
// 1.moles
//
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, col + 3 + (i*3), &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
oss3 << i + 1 << ".moles";
CPPUNIT_ASSERT_EQUAL( oss3.str(), std::string(v.sVal));
}
// sim
//
col = 0;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_LONG );
CPPUNIT_ASSERT_EQUAL( 1L, v.lVal );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_LONG );
CPPUNIT_ASSERT_EQUAL( 1L, v.lVal );
// state
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("i_soln"), std::string(v.sVal) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("react"), std::string(v.sVal) );
// soln
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_LONG );
CPPUNIT_ASSERT_EQUAL( 1L, v.lVal );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_LONG );
CPPUNIT_ASSERT_EQUAL( 1L, v.lVal );
// dist_x -- sometimes as double sometimes as long (depends on state)
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_LONG );
CPPUNIT_ASSERT_EQUAL( -99L, v.lVal );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_LONG );
CPPUNIT_ASSERT_EQUAL( -99L, v.lVal );
// time -- sometimes as double sometimes as long (depends on state)
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_LONG );
CPPUNIT_ASSERT_EQUAL( -99L, v.lVal );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.0, v.dVal, ::pow(10., -DBL_DIG) );
// step
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_LONG );
CPPUNIT_ASSERT( v.lVal == -99 );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_LONG );
CPPUNIT_ASSERT( v.lVal == 1 );
// pH
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 7.0, v.dVal, ::pow(10., -DBL_DIG) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 7.861354, v.dVal, ::pow(10., -6) );
// pe
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 4.0, v.dVal, ::pow(10., -DBL_DIG) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 10.4, v.dVal, ::pow(10., -1) );
//
// -totals C Ca Na
//
// C(mol/kgw)
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0000e-003, v.dVal, ::pow(10., -DBL_DIG) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.1556e-003, v.dVal, ::pow(10., -7) );
// Ca(mol/kgw)
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0000e-003, v.dVal, ::pow(10., -DBL_DIG) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.1556e-003, v.dVal, ::pow(10., -7) );
// Na(mol/kgw)
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0000e-003, v.dVal, ::pow(10., -DBL_DIG) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.0000e-003, v.dVal, ::pow(10., -7) );
// -molalities CO3-2 CaOH+ NaCO3-
col += 3;
// -activities CO3-2 CaOH+ NaCO3-
col += 3;
// -equilibrium_phases Calcite
col += 2;
// -saturation_indices CO2(g) Siderite
col += 2;
// -gases CO2(g) N2(g)
col += 5;
// -kinetic_reactants Albite Pyrite
// k_Albite dk_Albite k_Pyrite dk_Pyrite
col += 4;
// -solid_solutions CaSO4 SrSO4
// s_CaSO4 s_SrSO4
col += 2;
// 1.name
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("Ca+2"), std::string(v.sVal));
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("calcite"), std::string(v.sVal));
// 1.type
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("aq"), std::string(v.sVal));
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("equi"), std::string(v.sVal));
// 1.moles
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 9.9177923E-04, v.dVal, ::pow(10., -11) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 9.8444477E-03, v.dVal, ::pow(10., -10) );
// 2.name
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT_EQUAL( TT_STRING, v.type );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("CaHCO3+"), std::string(v.sVal));
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("Ca+2"), std::string(v.sVal));
// 2.type
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("aq"), std::string(v.sVal));
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("aq"), std::string(v.sVal));
// 2.moles
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 7.5980e-006, v.dVal, ::pow(10., -10) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.1371e-003, v.dVal, ::pow(10., -7) );
// 3.name
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT_EQUAL( TT_STRING, v.type );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("CaCO3"), std::string(v.sVal));
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("CaHCO3+"), std::string(v.sVal));
// 3.type
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("aq"), std::string(v.sVal));
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("aq"), std::string(v.sVal));
// 3.moles
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 6.2155e-007, v.dVal, ::pow(10., -11) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.1598e-005, v.dVal, ::pow(10., -9) );
// 4.name
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT_EQUAL( TT_STRING, v.type );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("CaOH+"), std::string(v.sVal));
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("CaCO3"), std::string(v.sVal));
// 4.type
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("aq"), std::string(v.sVal));
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("aq"), std::string(v.sVal));
// 4.moles
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 1.1819e-009, v.dVal, ::pow(10., -13) );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 6.8668e-006, v.dVal, ::pow(10., -10) );
// 5.name
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT_EQUAL( TT_EMPTY, v.type );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("CaOH+"), std::string(v.sVal));
// 5.type
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_EMPTY );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string("aq"), std::string(v.sVal));
// 5.moles
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_EMPTY );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_DOUBLE );
CPPUNIT_ASSERT_DOUBLES_EQUAL( 9.7385e-009, v.dVal, ::pow(10., -13) );
// 6.name
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT_EQUAL( TT_EMPTY, v.type );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT_EQUAL( TT_EMPTY, v.type );
// 6.type
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_EMPTY );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT_EQUAL( TT_EMPTY, v.type );
// 6.moles
//
++col;
// i_soln
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, col, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_EMPTY );
// react
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, col, &v) == VR_OK );
CPPUNIT_ASSERT_EQUAL( TT_EMPTY, v.type );
}
void TestInterface::TestGetSelectedOutputColumnCount()
{
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("llnl.dat"));
CPPUNIT_ASSERT_EQUAL( 0, ::GetSelectedOutputColumnCount() );
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
CPPUNIT_ASSERT_EQUAL( VR_OK, EQUILIBRIUM_PHASES("calcite", 1.0, 1.0) );
CPPUNIT_ASSERT_EQUAL( VR_OK, USER_PUNCH("Ca", 10) );
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( 62, ::GetSelectedOutputColumnCount() );
}
void TestInterface::TestAddError()
{
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
// make sure initialized to empty
//
const char* err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(err) );
// make sure initialized to empty
//
const char *expected = "TESTING AddError\n";
CPPUNIT_ASSERT_EQUAL(1u, ::AddError(expected));
// check 1
//
err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL( std::string(expected), std::string(err) );
// check increment
//
const char *expected2 = "XXXXXX\n";
CPPUNIT_ASSERT_EQUAL(2u, ::AddError(expected2));
// check concatenation
//
err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL( std::string(expected) + std::string(expected2), std::string(err) );
// clear errors
//
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
// make sure back to empty
//
err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(err) );
}
void TestInterface::TestAccumulateLine()
{
// TODO
}
void TestInterface::TestOutputLastError()
{
// TODO
}
void TestInterface::TestRunWithCallback()
{
// TODO
}
void TestInterface::TestRunNoDatabaseLoaded()
{
UnLoadDatabase();
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 1, ::RunAccumulated() );
const char expected[] =
"ERROR: RunAccumulated: No database is loaded\n"
"Stopping.\n";
const char* err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL( std::string(expected), std::string(err) );
}
void TestInterface::TestRunFileNoDatabaseLoaded()
{
UnLoadDatabase();
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 1, ::RunFile("dummy") );
const char expected[] =
"ERROR: RunFile: No database is loaded\n"
"Stopping.\n";
const char* err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL( std::string(expected), std::string(err) );
}
void TestInterface::TestCase1()
{
// Case 1 (see do_run)
// pr.punch == TRUE
// punch.new_def == FALSE
// output_isopen(OUTPUT_PUNCH) == FALSE
// selected_output_on == TRUE
// remove punch file if it exists
if (::FileExists("selected.out"))
{
CPPUNIT_ASSERT(::DeleteFile("selected.out"));
}
CPPUNIT_ASSERT_EQUAL( false, ::FileExists("selected.out") );
// clear all flags
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// COMMENT: {3/4/2010 6:28:53 PM} CPPUNIT_ASSERT_EQUAL( FALSE, punch.in);
// COMMENT: {3/4/2010 6:28:53 PM} CPPUNIT_ASSERT_EQUAL( TRUE, pr.punch);
CPPUNIT_ASSERT_EQUAL( FALSE, IPhreeqc::LibraryInstance()->punch.in);
CPPUNIT_ASSERT_EQUAL( TRUE, IPhreeqc::LibraryInstance()->pr.punch);
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
CPPUNIT_ASSERT_EQUAL( VR_OK, USER_PUNCH("Ca", 10) );
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(1);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( true, ::FileExists("selected.out") );
CPPUNIT_ASSERT_EQUAL( 62, ::GetSelectedOutputColumnCount() );
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(1);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( true, ::FileExists("selected.out") );
CPPUNIT_ASSERT_EQUAL( 62, ::GetSelectedOutputColumnCount() );
}
void TestInterface::TestCase2()
{
// Case 2 (see do_run)
// pr.punch == TRUE
// punch.new_def == TRUE
// output_isopen(OUTPUT_PUNCH) == FALSE
// selected_output_on == TRUE
// remove punch files if they exists
//
if (::FileExists("selected.out"))
{
::DeleteFile("selected.out");
}
if (::FileExists("case2.punch"))
{
::DeleteFile("case2.punch");
}
CPPUNIT_ASSERT_EQUAL( false, ::FileExists("selected.out") );
CPPUNIT_ASSERT_EQUAL( false, ::FileExists("case2.punch") );
// clear all flags
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// COMMENT: {3/4/2010 6:29:02 PM} CPPUNIT_ASSERT_EQUAL( FALSE, punch.in);
// COMMENT: {3/4/2010 6:29:02 PM} CPPUNIT_ASSERT_EQUAL( TRUE, pr.punch);
CPPUNIT_ASSERT_EQUAL( FALSE, IPhreeqc::LibraryInstance()->punch.in);
CPPUNIT_ASSERT_EQUAL( TRUE, IPhreeqc::LibraryInstance()->pr.punch);
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
CPPUNIT_ASSERT_EQUAL( VR_OK, USER_PUNCH("Ca", 10) );
CPPUNIT_ASSERT_EQUAL( VR_OK, ::AccumulateLine("-file case2.punch") ); // force have_punch_name to TRUE (see read_selected_ouput)
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(1);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( false, ::FileExists("selected.out") );
CPPUNIT_ASSERT_EQUAL( true, ::FileExists("case2.punch") );
CPPUNIT_ASSERT_EQUAL( 62, ::GetSelectedOutputColumnCount() );
// remove punch files if they exists
//
if (::FileExists("selected.out"))
{
::DeleteFile("selected.out");
}
if (::FileExists("case2.punch"))
{
::DeleteFile("case2.punch");
}
CPPUNIT_ASSERT_EQUAL( false, ::FileExists("selected.out") );
CPPUNIT_ASSERT_EQUAL( false, ::FileExists("case2.punch") );
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
CPPUNIT_ASSERT_EQUAL( VR_OK, USER_PUNCH("Ca", 10) );
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(1);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( false, ::FileExists("selected.out") );
CPPUNIT_ASSERT_EQUAL( true, ::FileExists("case2.punch") );
CPPUNIT_ASSERT_EQUAL( 62, ::GetSelectedOutputColumnCount() );
if (::FileExists("case2.punch"))
{
::DeleteFile("case2.punch");
}
CPPUNIT_ASSERT_EQUAL( false, ::FileExists("case2.punch") );
}
void TestInterface::TestPrintSelectedOutputFalse()
{
// remove punch files if they exists
//
if (::FileExists("selected.out"))
{
::DeleteFile("selected.out");
}
CPPUNIT_ASSERT_EQUAL( false, ::FileExists("selected.out") );
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// add selected output block
CPPUNIT_ASSERT_EQUAL( VR_OK, SELECTED_OUTPUT() );
// turn off selected output
CPPUNIT_ASSERT_EQUAL( VR_OK, ::AccumulateLine("PRINT; -selected_output false \n") );
// run
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(1);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( 0, ::GetSelectedOutputColumnCount() );
CPPUNIT_ASSERT_EQUAL( 0, ::GetSelectedOutputRowCount() );
// reset pr.punch to TRUE
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// add selected output block
CPPUNIT_ASSERT_EQUAL( VR_OK, SELECTED_OUTPUT() );
// run
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(1);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( 11, ::GetSelectedOutputColumnCount() );
CPPUNIT_ASSERT_EQUAL( 2, ::GetSelectedOutputRowCount() );
}
void TestInterface::TestOutputOnOff()
{
int onoff[5];
onoff[0] = 1; // output_on
onoff[1] = 0; // error_on
onoff[2] = 0; // log_on
onoff[3] = 0; // selected_output_on
onoff[4] = 0; // dump_on
TestOnOff("phreeqc.out", onoff[0], onoff[1], onoff[2], onoff[3], onoff[4]);
}
void TestInterface::TestErrorOnOff()
{
int onoff[5];
onoff[0] = 0; // output_on
onoff[1] = 1; // error_on
onoff[2] = 0; // log_on
onoff[3] = 0; // selected_output_on
onoff[4] = 0; // dump_on
TestOnOff("phreeqc.err", onoff[0], onoff[1], onoff[2], onoff[3], onoff[4]);
}
void TestInterface::TestLogOnOff()
{
int onoff[5];
onoff[0] = 0; // output_on
onoff[1] = 0; // error_on
onoff[2] = 1; // log_on
onoff[3] = 0; // selected_output_on
onoff[4] = 0; // dump_on
TestOnOff("phreeqc.log", onoff[0], onoff[1], onoff[2], onoff[3], onoff[4]);
}
void TestInterface::TestDumpOn()
{
int onoff[5];
onoff[0] = 0; // output_on
onoff[1] = 0; // error_on
onoff[2] = 0; // log_on
onoff[3] = 0; // selected_output_on
onoff[4] = 1; // dump_on
TestOnOff("dump.out", onoff[0], onoff[1], onoff[2], onoff[3], onoff[4]);
}
void TestInterface::TestSelOutOnOff()
{
int onoff[5];
onoff[0] = 0; // output_on
onoff[1] = 0; // error_on
onoff[2] = 0; // log_on
onoff[3] = 1; // selected_output_on
onoff[4] = 0; // dump_on
TestOnOff("selected.out", onoff[0], onoff[1], onoff[2], onoff[3], onoff[4]);
}
VRESULT
SOLUTION(double C, double Ca, double Na)
{
std::ostringstream oss;
oss << "SOLUTION 1\n";
oss << "C " << C << "\n";
oss << "Ca " << Ca << "\n";
oss << "Na " << Na << "\n";
return ::AccumulateLine(oss.str().c_str());
}
VRESULT
EQUILIBRIUM_PHASES(const char* phase, double si, double amount)
{
std::ostringstream oss;
oss << "EQUILIBRIUM_PHASES\n";
oss << phase << " " << si << " " << amount << "\n";
return ::AccumulateLine(oss.str().c_str());
}
VRESULT
USER_PUNCH(const char* element, int max)
{
std::ostringstream oss;
oss << "USER_PUNCH\n";
oss << "-head ";
for (int i = 1; i <= max; ++i)
{
oss << i << ".name " << i << ".type " << i << ".moles ";
}
oss << "\n";
oss << "-start" << "\n";
oss << "10 n = sys(\"" << element << "\"" << ", count, names$, types$, moles)" << "\n";
oss << "20 n = " << max << "\n";
oss << "30 if count < " << max << " then n = count" << "\n";
oss << "40 for i = 1 to count" << "\n";
oss << "50 PUNCH names$(i), types$(i), moles(i)" << "\n";
oss << "60 next i" << "\n";
oss << "70 list" << "\n";
oss << "-end" << "\n";
oss << "SELECTED_OUTPUT" << "\n";
oss << "-totals C Ca Na" << "\n";
oss << "-molalities CO3-2 CaOH+ NaCO3-" << "\n";
oss << "-activities CO3-2 CaOH+ NaCO3-" << "\n";
oss << "-equilibrium_phases Calcite" << "\n";
oss << "-saturation_indices CO2(g) Siderite" << "\n";
oss << "-gases CO2(g) N2(g)" << "\n";
oss << "-kinetic_reactants Albite Pyrite" << "\n";
oss << "-solid_solutions CaSO4 SrSO4" << "\n";
return ::AccumulateLine(oss.str().c_str());
}
VRESULT
SELECTED_OUTPUT(void)
{
std::ostringstream oss;
oss << "SELECTED_OUTPUT" << "\n";
oss << "-file selected.out" << "\n";
oss << "-totals C Ca Na" << "\n";
return ::AccumulateLine(oss.str().c_str());
}
VRESULT
DUMP(void)
{
std::ostringstream oss;
oss << "DUMP" << "\n";
oss << "-solution 1" << "\n";
return ::AccumulateLine(oss.str().c_str());
}
void TestOnOff(const char* FILENAME, int output_on, int error_on, int log_on, int selected_output_on, int dump_on)
{
//const char *FILENAME = "phreeqc.out";
// remove FILENAME if it exists
//
if (::FileExists(FILENAME))
{
::DeleteFile(FILENAME);
}
CPPUNIT_ASSERT_EQUAL( false, ::FileExists(FILENAME) );
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// add selected output block
CPPUNIT_ASSERT_EQUAL( VR_OK, SELECTED_OUTPUT() );
// add dump block
CPPUNIT_ASSERT_EQUAL( VR_OK, DUMP() );
// run all off
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( false, ::FileExists(FILENAME) );
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// add selected output block
CPPUNIT_ASSERT_EQUAL( VR_OK, SELECTED_OUTPUT() );
// add dump block
CPPUNIT_ASSERT_EQUAL( VR_OK, DUMP() );
// run
::SetOutputFileOn(output_on);
::SetErrorFileOn(error_on);
::SetLogFileOn(log_on);
::SetSelectedOutputFileOn(selected_output_on);
::SetDumpFileOn(dump_on);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( true, ::FileExists(FILENAME) );
CPPUNIT_ASSERT( ::DeleteFile(FILENAME) );
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// add selected output block
CPPUNIT_ASSERT_EQUAL( VR_OK, SELECTED_OUTPUT() );
// add dump block
CPPUNIT_ASSERT_EQUAL( VR_OK, DUMP() );
// run
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( false, ::FileExists(FILENAME) );
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// add selected output block
CPPUNIT_ASSERT_EQUAL( VR_OK, SELECTED_OUTPUT() );
// add dump block
CPPUNIT_ASSERT_EQUAL( VR_OK, DUMP() );
// run
::SetOutputFileOn(output_on);
::SetErrorFileOn(error_on);
::SetLogFileOn(log_on);
::SetSelectedOutputFileOn(selected_output_on);
::SetDumpFileOn(dump_on);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( true, ::FileExists(FILENAME) );
CPPUNIT_ASSERT( ::DeleteFile(FILENAME) );
}
void
TestInterface::TestLongHeadings()
{
char long_header[] = "this_is_a_long_header_0123456789012345678901234567890123456789";
char long_value[] = "this_is_a_long_value_01234567890123456789012345678901234567890";
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
std::ostringstream oss;
oss << "SOLUTION" << "\n";
oss << "SELECTED_OUTPUT" << "\n";
oss << "-reset false" << "\n";
oss << "USER_PUNCH" << "\n";
oss << "-head " << long_header << "\n";
oss << "-start" << "\n";
oss << "10 PUNCH \"" << long_value << "\"\n";
oss << "-end" << "\n";
CPPUNIT_ASSERT_EQUAL( VR_OK, ::AccumulateLine(oss.str().c_str()) );
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL(2, ::GetSelectedOutputRowCount());
CPPUNIT_ASSERT_EQUAL(1, ::GetSelectedOutputColumnCount());
CVar v;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(0, 0, &v) == VR_OK );
CPPUNIT_ASSERT( v.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string(long_header), std::string(v.sVal));
CVar v1;
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, 0, &v1) == VR_OK );
CPPUNIT_ASSERT( v1.type == TT_STRING );
CPPUNIT_ASSERT_EQUAL( std::string(long_value), std::string(v1.sVal));
CPPUNIT_ASSERT( ::GetSelectedOutputValue(1, 1, &v1) == VR_INVALIDCOL );
CPPUNIT_ASSERT( ::GetSelectedOutputValue(2, 0, &v1) == VR_INVALIDROW );
}
void
TestInterface::TestDatabaseKeyword()
{
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
CPPUNIT_ASSERT_EQUAL(1, ::RunFile("dump"));
const char *expected =
"ERROR: Gas not found in PHASES data base, Amm(g).\n"
"ERROR: Calculations terminating due to input errors.\n"
"Stopping.\n";
const char* err = ::GetLastErrorString();
CPPUNIT_ASSERT_EQUAL(std::string(expected), std::string(err));
}
void
TestInterface::TestDumpString()
{
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// add dump block
CPPUNIT_ASSERT_EQUAL( VR_OK, DUMP() );
// run
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
::SetDumpStringOn(1);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
const char *expected =
#if defined(_MSC_VER)
"SOLUTION_RAW 1 \n"
" -temp 25\n"
" -total_h 111.0132593403\n"
" -total_o 55.509043478605\n"
" -cb 0.0021723831003929\n"
" -totals\n"
" C(4) 0.0010000000000376\n"
" Ca 0.001000000004331\n"
" H(0) 1.4149476909313e-025\n"
" Na 0.001\n"
" -Isotopes\n"
" -pH 7\n"
" -pe 4\n"
" -mu 0.0028961089894362\n"
" -ah2o 0.99994915105857\n"
" -mass_water 1\n"
" -total_alkalinity 0.00082761690826911\n"
" -activities\n"
" C(-4) -67.370522674574\n"
" C(4) -6.4415889265024\n"
" Ca -3.1040445240857\n"
" E -4\n"
" H(0) -25.15\n"
" Na -3.0255625287599\n"
" O(0) -42.080044167952\n"
" -gammas\n";
#endif
#if defined(__GNUC__)
"SOLUTION_RAW 1 \n"
" -temp 25\n"
" -total_h 111.0132593403\n"
" -total_o 55.509043478605\n"
" -cb 0.0021723831003928\n"
" -totals\n"
" C(4) 0.0010000000000376\n"
" Ca 0.001000000004331\n"
" H(0) 1.4149476909313e-25\n"
" Na 0.001\n"
" -Isotopes\n"
" -pH 7\n"
" -pe 4\n"
" -mu 0.0028961089894362\n"
" -ah2o 0.99994915105857\n"
" -mass_water 1\n"
" -total_alkalinity 0.00082761690826912\n"
" -activities\n"
" C(-4) -67.370522674574\n"
" C(4) -6.4415889265024\n"
" Ca -3.1040445240857\n"
" E -4\n"
" H(0) -25.15\n"
" Na -3.0255625287599\n"
" O(0) -42.080044167952\n"
" -gammas\n";
#endif
const char* dump_str = ::GetDumpString();
CPPUNIT_ASSERT_EQUAL(std::string(expected), std::string(dump_str));
}
void
TestInterface::TestGetDumpStringLineCount()
{
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// add dump block
CPPUNIT_ASSERT_EQUAL( VR_OK, DUMP() );
// run
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
::SetDumpStringOn(1);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( 26, ::GetDumpStringLineCount() );
}
void
TestInterface::TestGetDumpStringLine()
{
CPPUNIT_ASSERT_EQUAL(0, ::LoadDatabase("phreeqc.dat"));
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// add dump block
CPPUNIT_ASSERT_EQUAL( VR_OK, DUMP() );
// run
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
::SetDumpStringOn(1);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( 26, ::GetDumpStringLineCount() );
int line = 0;
#if defined(_MSC_VER)
CPPUNIT_ASSERT_EQUAL( std::string("SOLUTION_RAW 1 "), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -temp 25"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -total_h 111.0132593403"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -total_o 55.509043478605"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -cb 0.0021723831003929"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -totals"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" C(4) 0.0010000000000376"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" Ca 0.001000000004331"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" H(0) 1.4149476909313e-025"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" Na 0.001"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -Isotopes"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -pH 7"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -pe 4"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -mu 0.0028961089894362"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -ah2o 0.99994915105857"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -mass_water 1"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -total_alkalinity 0.00082761690826911"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -activities"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" C(-4) -67.370522674574"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" C(4) -6.4415889265024"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" Ca -3.1040445240857"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" E -4"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" H(0) -25.15"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" Na -3.0255625287599"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" O(0) -42.080044167952"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -gammas"), std::string(::GetDumpStringLine(line++)) );
#endif
#if defined(__GNUC__)
CPPUNIT_ASSERT_EQUAL( std::string("SOLUTION_RAW 1 "), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -temp 25"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -total_h 111.0132593403"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -total_o 55.509043478605"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -cb 0.0021723831003928"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -totals"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" C(4) 0.0010000000000376"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" Ca 0.001000000004331"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" H(0) 1.4149476909313e-25"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" Na 0.001"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -Isotopes"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -pH 7"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -pe 4"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -mu 0.0028961089894362"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -ah2o 0.99994915105857"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -mass_water 1"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -total_alkalinity 0.00082761690826912"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -activities"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" C(-4) -67.370522674574"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" C(4) -6.4415889265024"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" Ca -3.1040445240857"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" E -4"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" H(0) -25.15"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" Na -3.0255625287599"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" O(0) -42.080044167952"), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(" -gammas"), std::string(::GetDumpStringLine(line++)) );
#endif
// remaining lines should be empty
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetDumpStringLine(line++)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetDumpStringLine(line++)) );
// negative lines should be empty
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetDumpStringLine(-1)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetDumpStringLine(-2)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetDumpStringLine(-3)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetDumpStringLine(-4)) );
}
void
TestInterface::TestGetComponentCount(void)
{
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// run
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
::SetDumpStringOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( 3, ::GetComponentCount() );
}
void
TestInterface::TestGetComponent(void)
{
CPPUNIT_ASSERT_EQUAL( 0, ::LoadDatabase("phreeqc.dat") );
// add solution block
CPPUNIT_ASSERT_EQUAL( VR_OK, SOLUTION(1.0, 1.0, 1.0) );
// run
::SetOutputFileOn(0);
::SetErrorFileOn(0);
::SetLogFileOn(0);
::SetSelectedOutputFileOn(0);
::SetDumpFileOn(0);
::SetDumpStringOn(0);
CPPUNIT_ASSERT_EQUAL( 0, ::RunAccumulated() );
CPPUNIT_ASSERT_EQUAL( 3, ::GetComponentCount() );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetComponent(-2)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetComponent(-1)) );
CPPUNIT_ASSERT_EQUAL( std::string("C"), std::string(::GetComponent(0)) );
CPPUNIT_ASSERT_EQUAL( std::string("Ca"), std::string(::GetComponent(1)) );
CPPUNIT_ASSERT_EQUAL( std::string("Na"), std::string(::GetComponent(2)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetComponent(3)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetComponent(4)) );
CPPUNIT_ASSERT_EQUAL( std::string(""), std::string(::GetComponent(5)) );
}
Reference in New Issue View Git Blame Copy Permalink