OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Gives the same jacobi and residual for geocell for fixed as r 8522 in MATLAB modulo different transmissiblity calculations for pcflux.

Browse Source
This commit is contained in:
Ingeborg Ligaarden 2011-11-11 10:21:24 +01:00
parent 17895db0ff
commit 298280129c
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/ImplicitTransport.hpp
View File

@ -114,13 +114,17 @@ namespace Opm {
linsolve.solve(sys_.matrix(),
sys_.vector().residual(),
sys_.vector().writableIncrement());
Dune::writeMatrixToMatlab(sys_.matrix(),"matrix_matlab");
std::ofstream myfile2;
myfile2.open ("jacobi.txt");
Dune::printSparseMatrix(myfile2, sys_.matrix(), "", "", 26, 15);
myfile2.close();
Dune::writeMatrixToMatlab(sys_.matrix(),"matrix_matlab"); // write jacobi
std::ofstream myfile;
myfile.open ("residual.txt");
Dune::printvector(myfile,sys_.vector().residual(),"redidual","");
Dune::printvector(myfile,sys_.vector().residual(),"redidual","", 1, 10, 15); //write residual
myfile.close();
VNeg<vector_type>::negate(sys_.vector().writableIncrement());
//model_.finishIteration(state, g, sys_.vector());