{"action":"created","issue":{"url":"https://api.github.com/repos/OPM/opm-simulators/issues/1480","repository_url":"https://api.github.com/repos/OPM/opm-simulators","labels_url":"https://api.github.com/repos/OPM/opm-simulators/issues/1480/labels{/name}","comments_url":"https://api.github.com/repos/OPM/opm-simulators/issues/1480/comments","events_url":"https://api.github.com/repos/OPM/opm-simulators/issues/1480/events","html_url":"https://github.com/OPM/opm-simulators/pull/1480","id":325182060,"number":1480,"title":"[WIP and discussion] using invDXDirect from MSWellHelpers to avoid inverting matrix D","user":{"login":"GitPaean","id":6847941,"avatar_url":"https://avatars0.githubusercontent.com/u/6847941?v=4","gravatar_id":"","url":"https://api.github.com/users/GitPaean","html_url":"https://github.com/GitPaean","followers_url":"https://api.github.com/users/GitPaean/followers","following_url":"https://api.github.com/users/GitPaean/following{/other_user}","gists_url":"https://api.github.com/users/GitPaean/gists{/gist_id}","starred_url":"https://api.github.com/users/GitPaean/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/GitPaean/subscriptions","organizations_url":"https://api.github.com/users/GitPaean/orgs","repos_url":"https://api.github.com/users/GitPaean/repos","events_url":"https://api.github.com/users/GitPaean/events{/privacy}","received_events_url":"https://api.github.com/users/GitPaean/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"assignees":[],"milestone":null,"comments":42,"created_at":"2018-05-22T08:04:46Z","updated_at":"2018-05-28T19:55:42Z","closed_at":null,"author_association":"MEMBER","pull_request":{"url":"https://api.github.com/repos/OPM/opm-simulators/pulls/1480","html_url":"https://github.com/OPM/opm-simulators/pull/1480","diff_url":"https://github.com/OPM/opm-simulators/pull/1480.diff","patch_url":"https://github.com/OPM/opm-simulators/pull/1480.patch"},"body":"for StandardWell. Basically, as @blattms pointed out for more than one time, the `invert()` in DUNE is not very stable. This PR using a direct solver to avoid inverting the matrix D for the StandardWell. \r\n\r\nIt is for testing and discussion purpose and it is not in a good form for going in.  When we decide what we want, I will revise and make it in a good form to go in. \r\n\r\nPlease do not merge unless requested and agreed. "},"comment":{"url":"https://api.github.com/repos/OPM/opm-simulators/issues/comments/392594501","html_url":"https://github.com/OPM/opm-simulators/pull/1480#issuecomment-392594501","issue_url":"https://api.github.com/repos/OPM/opm-simulators/issues/1480","id":392594501,"user":{"login":"GitPaean","id":6847941,"avatar_url":"https://avatars0.githubusercontent.com/u/6847941?v=4","gravatar_id":"","url":"https://api.github.com/users/GitPaean","html_url":"https://github.com/GitPaean","followers_url":"https://api.github.com/users/GitPaean/followers","following_url":"https://api.github.com/users/GitPaean/following{/other_user}","gists_url":"https://api.github.com/users/GitPaean/gists{/gist_id}","starred_url":"https://api.github.com/users/GitPaean/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/GitPaean/subscriptions","organizations_url":"https://api.github.com/users/GitPaean/orgs","repos_url":"https://api.github.com/users/GitPaean/repos","events_url":"https://api.github.com/users/GitPaean/events{/privacy}","received_events_url":"https://api.github.com/users/GitPaean/received_events","type":"User","site_admin":false},"created_at":"2018-05-28T19:55:42Z","updated_at":"2018-05-28T19:55:42Z","author_association":"MEMBER","body":"> Could you post the error message you got?\r\n\r\n```\r\nThread 1 \"flow\" received signal SIGSEGV, Segmentation fault.\r\n0x00007ffff7527909 in umfdi_solve () from /usr/lib/x86_64-linux-gnu/libumfpack.so.5.7.1\r\n\r\n(gdb) backtrace \r\n#0  0x00007ffff7527909 in umfdi_solve () from /usr/lib/x86_64-linux-gnu/libumfpack.so.5.7.1\r\n#1  0x00007ffff75375f2 in umfpack_di_solve () from /usr/lib/x86_64-linux-gnu/libumfpack.so.5.7.1\r\n#2  0x000000000078f8ef in Dune::UMFPackMethodChooser<double>::solve<int, int*, int*, double*, double*, double*, void*, double*, double*> () at /home/kaib/local/dune2.5.2/release2.5.2_O3_mpi/include/dune/istl/umfpack.hh:108\r\n#3  Dune::UMFPack<Dune::BCRSMatrix<Dune::FieldMatrix<double, 3, 3>, std::allocator<Dune::FieldMatrix<double, 3, 3> > > >::apply (this=0x1b081e0, x=..., b=..., res=...)\r\n    at /home/kaib/local/dune2.5.2/release2.5.2_O3_mpi/include/dune/istl/umfpack.hh:329\r\n#4  0x00000000008f4a2c in Opm::wellhelpers::invDXDirect<Dune::UMFPack<Dune::BCRSMatrix<Dune::FieldMatrix<double, 3, 3>, std::allocator<Dune::FieldMatrix<double, 3, 3> > > >, Dune::BlockVector<Dune::FieldVector<double, 3>, std::allocator<Dune::FieldVector<double, 3> > > > (linsolver=..., x=...) at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/WellHelpers.hpp:209\r\n#5  0x00000000009014ed in Opm::StandardWell<Ewoms::Properties::TTag::EclFlowProblem>::solveEqAndUpdateWellState (this=0x1b07eb0, well_state=...)\r\n    at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/StandardWell_impl.hpp:1618\r\n#6  0x000000000090c2b4 in Opm::BlackoilWellModel<Ewoms::Properties::TTag::EclFlowProblem>::solveWellEq (this=this@entry=0x7fffffffce00, dt=dt@entry=86400)\r\n    at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/BlackoilWellModel_impl.hpp:494\r\n#7  0x000000000090e153 in Opm::BlackoilWellModel<Ewoms::Properties::TTag::EclFlowProblem>::assemble (this=0x7fffffffce00, iterationIdx=iterationIdx@entry=0, dt=86400)\r\n    at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/BlackoilWellModel_impl.hpp:281\r\n#8  0x000000000092eb49 in Opm::BlackoilModelEbos<Ewoms::Properties::TTag::EclFlowProblem>::assemble (iterationIdx=0, timer=..., this=0x1b0abb0) at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/BlackoilModelEbos.hpp:376\r\n#9  Opm::BlackoilModelEbos<Ewoms::Properties::TTag::EclFlowProblem>::nonlinearIteration<Opm::NonlinearSolver<Opm::BlackoilModelEbos<Ewoms::Properties::TTag::EclFlowProblem> > > (nonlinear_solver=..., timer=..., iteration=0, \r\n    this=0x1b0abb0) at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/BlackoilModelEbos.hpp:248\r\n#10 Opm::NonlinearSolver<Opm::BlackoilModelEbos<Ewoms::Properties::TTag::EclFlowProblem> >::step (this=this@entry=0x1a6b1d0, timer=..., initial_reservoir_state=..., initial_well_state=..., reservoir_state=..., well_state=...)\r\n    at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/NonlinearSolver_impl.hpp:145\r\n#11 0x0000000000931d95 in Opm::NonlinearSolver<Opm::BlackoilModelEbos<Ewoms::Properties::TTag::EclFlowProblem> >::step (well_state=..., reservoir_state=..., timer=..., this=0x1a6b1d0)\r\n    at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/NonlinearSolver_impl.hpp:111\r\n#12 Opm::AdaptiveTimeStepping::stepImpl<Opm::NonlinearSolver<Opm::BlackoilModelEbos<Ewoms::Properties::TTag::EclFlowProblem> >, Opm::BlackoilState, Opm::WellStateFullyImplicitBlackoil, Opm::BlackoilOutputEbos<Ewoms::Properties::TTag::EclFlowProblem> > (this=this@entry=0x1b0a640, simulatorTimer=..., solver=..., state=..., well_state=..., event=true, outputWriter=0x1a6e6e0, fipnum=0x0)\r\n    at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/simulators/timestepping/AdaptiveTimeStepping_impl.hpp:236\r\n#13 0x000000000094ce43 in Opm::AdaptiveTimeStepping::step<Opm::NonlinearSolver<Opm::BlackoilModelEbos<Ewoms::Properties::TTag::EclFlowProblem> >, Opm::BlackoilState, Opm::WellStateFullyImplicitBlackoil, Opm::BlackoilOutputEbos<Ewoms::Properties::TTag::EclFlowProblem> > (fipnum=0x0, outputWriter=..., event=<optimized out>, well_state=..., state=..., solver=..., simulatorTimer=..., this=<optimized out>)\r\n    at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/simulators/timestepping/AdaptiveTimeStepping_impl.hpp:177\r\n#14 Opm::SimulatorFullyImplicitBlackoilEbos<Ewoms::Properties::TTag::EclFlowProblem>::run (this=0x1b09b00, timer=...) at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/SimulatorFullyImplicitBlackoilEbos.hpp:250\r\n#15 0x000000000094ee8b in Opm::FlowMainEbos<Ewoms::Properties::TTag::EclFlowProblem>::runSimulator (this=this@entry=0x7fffffffdc10) at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/FlowMainEbos.hpp:517\r\n#16 0x00000000009537be in Opm::FlowMainEbos<Ewoms::Properties::TTag::EclFlowProblem>::execute (this=this@entry=0x7fffffffdc10, argc=<optimized out>, argv=<optimized out>)\r\n    at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/autodiff/FlowMainEbos.hpp:115\r\n#17 0x0000000000786fcd in Opm::flowEbosBlackoilMain (argc=3, argv=0x7fffffffe348) at /home/kaib/OPM-PR-test/debug/opm-simulators/opm/simulators/flow_ebos_blackoil.cpp:59\r\n#18 0x0000000000747f47 in main (argc=3, argv=0x7fffffffe348) at /home/kaib/OPM-PR-test/debug/opm-simulators/examples/flow.cpp:200\r\n```"}}
{"actor":{"display_login":"GitPaean"}}