mirror of
https://github.com/OPM/opm-simulators.git
synced 2024-12-28 02:00:59 -06:00
commit
778a7e2343
@ -46,4 +46,4 @@ void printFlowTrailer(int nprocs, int nthreads,
|
||||
|
||||
} // namespace Opm
|
||||
|
||||
#endif // OPM_FLOW_BANNERS_EBOS_HEADER_INCLUDED
|
||||
#endif // OPM_FLOW_BANNERS_HEADER_INCLUDED
|
||||
|
@ -210,7 +210,7 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
|
||||
|
||||
void initialize()
|
||||
{
|
||||
OPM_TIMEBLOCK(IstlSolverEbos);
|
||||
OPM_TIMEBLOCK(IstlSolver);
|
||||
|
||||
if (parameters_[0].linsolver_ == "hybrid") {
|
||||
// Experimental hybrid configuration.
|
||||
@ -320,7 +320,7 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
|
||||
|
||||
// update matrix entries for solvers.
|
||||
if (firstcall) {
|
||||
// ebos will not change the matrix object. Hence simply store a pointer
|
||||
// model will not change the matrix object. Hence simply store a pointer
|
||||
// to the original one with a deleter that does nothing.
|
||||
// Outch! We need to be able to scale the linear system! Hence const_cast
|
||||
matrix_ = const_cast<Matrix*>(&M);
|
||||
@ -349,7 +349,7 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
|
||||
|
||||
void prepare(const Matrix& M, Vector& b)
|
||||
{
|
||||
OPM_TIMEBLOCK(istlSolverEbosPrepare);
|
||||
OPM_TIMEBLOCK(istlSolverPrepare);
|
||||
|
||||
initPrepare(M,b);
|
||||
|
||||
@ -382,7 +382,7 @@ std::unique_ptr<Matrix> blockJacobiAdjacency(const Grid& grid,
|
||||
|
||||
bool solve(Vector& x)
|
||||
{
|
||||
OPM_TIMEBLOCK(istlSolverEbosSolve);
|
||||
OPM_TIMEBLOCK(istlSolverSolve);
|
||||
++solveCount_;
|
||||
// Write linear system if asked for.
|
||||
const int verbosity = prm_[activeSolverNum_].get("verbosity", 0);
|
||||
|
@ -198,7 +198,7 @@ public:
|
||||
#if HAVE_OPENCL
|
||||
// update matrix entries for solvers.
|
||||
if (firstcall && bdaBridge_) {
|
||||
// ebos will not change the matrix object. Hence simply store a pointer
|
||||
// model will not change the matrix object. Hence simply store a pointer
|
||||
// to the original one with a deleter that does nothing.
|
||||
// Outch! We need to be able to scale the linear system! Hence const_cast
|
||||
// setup sparsity pattern for jacobi matrix for preconditioner (only used for openclSolver)
|
||||
|
@ -39,7 +39,7 @@ namespace Opm
|
||||
//=====================================================================
|
||||
// Implementation for ISTL-matrix based operators
|
||||
// Note: the classes WellModelMatrixAdapter and
|
||||
// WellModelGhostLastMatrixAdapter were moved from ISTLSolverEbos.hpp
|
||||
// WellModelGhostLastMatrixAdapter were moved from ISTLSolver.hpp
|
||||
// and subsequently modified.
|
||||
//=====================================================================
|
||||
|
||||
|
@ -27,8 +27,7 @@
|
||||
<h3>Solvers and simulators</h3>
|
||||
|
||||
The opm-simulators module is the home for the Flow reservoir
|
||||
simulator, the Ebos experimental simulator and library functionality
|
||||
on which they both depend.
|
||||
simulators and library functionality on which they depend.
|
||||
|
||||
*/
|
||||
|
||||
|
@ -88,7 +88,7 @@ namespace Opm
|
||||
|
||||
// It makes no sense adding time points. Therefore, do not
|
||||
// overwrite the value of global_time which gets set in
|
||||
// NonlinearSolverEbos.hpp by the line:
|
||||
// NonlinearSolver.hpp by the line:
|
||||
// report.global_time = timer.simulationTimeElapsed();
|
||||
}
|
||||
|
||||
|
Loading…
Reference in New Issue
Block a user