Added inizialization of ok, which correct error for solves without wells. The only thing which is check is if the well object have incompressible form.
This commit is contained in:
@@ -585,7 +585,7 @@ assemble_incompressible(struct UnstructuredGrid *G ,
|
|||||||
int res_is_neumann, wells_are_rate;
|
int res_is_neumann, wells_are_rate;
|
||||||
|
|
||||||
double s;
|
double s;
|
||||||
|
*ok=1;
|
||||||
csrmatrix_zero( h->A);
|
csrmatrix_zero( h->A);
|
||||||
vector_zero (h->A->m, h->b);
|
vector_zero (h->A->m, h->b);
|
||||||
|
|
||||||
@@ -694,7 +694,6 @@ ifs_tpfa_assemble(struct UnstructuredGrid *G ,
|
|||||||
/* ---------------------------------------------------------------------- */
|
/* ---------------------------------------------------------------------- */
|
||||||
{
|
{
|
||||||
int system_singular, ok;
|
int system_singular, ok;
|
||||||
|
|
||||||
assemble_incompressible(G, F, trans, gpress, h, &system_singular, &ok);
|
assemble_incompressible(G, F, trans, gpress, h, &system_singular, &ok);
|
||||||
|
|
||||||
if (ok && system_singular) {
|
if (ok && system_singular) {
|
||||||
|
|||||||
Reference in New Issue
Block a user