research/cantera
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Replaced incorrect manual memory management with std::vector

Browse source
This commit is contained in:
Ray Speth 2012-02-17 20:28:58 +00:00
parent e2ac2afad5
commit 708b78f68b
2 changed files with 6 additions and 12 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
src/transport/LiquidTranInteraction.cpp
View file

@ -821,8 +821,8 @@ void LTI_StokesEinstein::getMatrixTransProp(DenseMatrix& mat, doublereal* specie
int nsp = m_thermo->nSpecies();
doublereal temp = m_thermo->temperature();
double* viscSpec = new double(nsp);
double* radiusSpec = new double(nsp);
array_fp viscSpec(nsp);
array_fp radiusSpec(nsp);
for (int k = 0; k < nsp; k++) {
viscSpec[k] = m_viscosity[k]->getSpeciesTransProp() ;
@ -834,8 +834,6 @@ void LTI_StokesEinstein::getMatrixTransProp(DenseMatrix& mat, doublereal* specie
for (int j = 0; j < nsp; j++) {
mat(i,j) = (6.0 * Pi * radiusSpec[i] * viscSpec[j]) / GasConstant / temp;
}
delete radiusSpec;
delete viscSpec;
}
doublereal LTI_MoleFracs_ExpT::getMixTransProp(doublereal* speciesValues, doublereal* speciesWeight)

12
src/transport/LiquidTransport.cpp
View file

@ -923,10 +923,10 @@ doublereal LiquidTransport::getElectricConduct()
set_Grad_X(&gradX[0]);
set_Grad_V(&gradV[0]);
doublereal* fluxes = new doublereal(m_nsp * m_nDim);
vector_fp fluxes(m_nsp * m_nDim);
doublereal current;
getSpeciesFluxesExt(m_nDim, fluxes);
getSpeciesFluxesExt(m_nDim, &fluxes[0]);
//sum over species charges, fluxes, Faraday to get current
// Since we want the scalar conductivity, we need only consider one-dim
@ -938,9 +938,7 @@ doublereal LiquidTransport::getElectricConduct()
//divide by unit potential gradient
current /= - gradV[i];
}
delete fluxes;
return current;
}
// Compute the electric current density in A/m^2
@ -977,9 +975,9 @@ void LiquidTransport::getElectricCurrent(int ndim,
set_Grad_X(grad_X);
set_Grad_V(grad_V);
doublereal* fluxes = new doublereal(m_nsp * m_nDim);
vector_fp fluxes(m_nsp * m_nDim);
getSpeciesFluxesExt(ldf, fluxes);
getSpeciesFluxesExt(ldf, &fluxes[0]);
//sum over species charges, fluxes, Faraday to get current
for (size_t i = 0; i < m_nDim; i++) {
@ -989,8 +987,6 @@ void LiquidTransport::getElectricCurrent(int ndim,
}
//divide by unit potential gradient
}
delete fluxes;
}
// Get the species diffusive velocities wrt to