Skip to content

Commit

Permalink
-vxworks includes
Browse files Browse the repository at this point in the history
git-svn-id: https://openmodelica.org/svn/OpenModelica/trunk@22566 f25d12d1-65f4-0310-ae8a-bbce733d8d8e
  • Loading branch information
RuedKamp committed Oct 2, 2014
1 parent 0a00782 commit 1f1710e
Show file tree
Hide file tree
Showing 2 changed files with 8 additions and 2 deletions.
8 changes: 6 additions & 2 deletions SimulationRuntime/cpp/Solver/Newton/Newton.cpp
Expand Up @@ -12,6 +12,7 @@ Newton::Newton(IAlgLoop* algLoop, INonLinSolverSettings* settings)
, _yHelp (NULL)
, _f (NULL)
, _fHelp (NULL)
,_iHelp (NULL)
, _jac (NULL)
, _dimSys (0)
, _firstCall (true)
Expand All @@ -25,6 +26,7 @@ Newton::~Newton()
if(_yHelp) delete [] _yHelp;
if(_f) delete [] _f;
if(_fHelp) delete [] _fHelp;
if(_iHelp) delete [] _iHelp;
if(_jac) delete [] _jac;
}

Expand Down Expand Up @@ -53,12 +55,14 @@ void Newton::initialize()
if(_f) delete [] _f;
if(_yHelp) delete [] _yHelp;
if(_fHelp) delete [] _fHelp;
if(_iHelp) delete [] _iHelp;
if(_jac) delete [] _jac;

_y = new double[_dimSys];
_f = new double[_dimSys];
_yHelp = new double[_dimSys];
_fHelp = new double[_dimSys];
_iHelp = new long int[_dimSys];
_jac = new double[_dimSys*_dimSys];

_algLoop->getReal(_y);
Expand Down Expand Up @@ -129,7 +133,7 @@ void Newton::solve()
{
//calcFunction(_yHelp,_fHelp);
_algLoop->getSystemMatrix(_jac);
dgesv_(&_dimSys,&dimRHS,_jac,&_dimSys,_fHelp,_f,&_dimSys,&irtrn);
dgesv_(&_dimSys,&dimRHS,_jac,&_dimSys,_iHelp,_f,&_dimSys,&irtrn);
memcpy(_y,_f,_dimSys*sizeof(double));
_algLoop->setReal(_y);
_iterationStatus = DONE;
Expand All @@ -141,7 +145,7 @@ void Newton::solve()
calcJacobian();
}
// Solve linear System
dgesv_(&_dimSys,&dimRHS,_jac,&_dimSys,_fHelp,_f,&_dimSys,&irtrn);
dgesv_(&_dimSys,&dimRHS,_jac,&_dimSys,_iHelp,_f,&_dimSys,&irtrn);

if(irtrn!=0)
{
Expand Down
2 changes: 2 additions & 0 deletions SimulationRuntime/cpp/Solver/Newton/Newton.h
Expand Up @@ -89,4 +89,6 @@ class Newton : public IAlgLoopSolver
*_yHelp, ///< Temp - Auxillary variables
*_fHelp, ///< Temp - Auxillary variables
*_jac; ///< Temp - Jacobian

long int *_iHelp;
};

0 comments on commit 1f1710e

Please sign in to comment.