Commit 9bbc0fc 1 parent 7ef9714 commit 9bbc0fc Copy full SHA for 9bbc0fc
File tree 1 file changed +3
-1
lines changed
opm/simulators/flow/python
1 file changed +3
-1
lines changed Original file line number Diff line number Diff line change @@ -115,6 +115,7 @@ getFluidStateVariable(const std::string &name) const
115
115
116
116
Model &model = this ->ebos_simulator_ ->model ();
117
117
auto size = model.numGridDof ();
118
+ std::cout << " model size=" << size << std::endl;
118
119
std::vector<double > array (size);
119
120
const auto & grid_view = this ->ebos_simulator_ ->vanguard ().gridView ();
120
121
/* NOTE: grid_view.size(0) should give the same value as
@@ -129,7 +130,7 @@ getFluidStateVariable(const std::string &name) const
129
130
for (; elem_itr != elem_end_itr; ++elem_itr) {
130
131
std::cout << " i=" << i << std::endl;
131
132
const Element& elem = *elem_itr;
132
- // elem_ctx.updatePrimaryStencil(elem);
133
+ elem_ctx.updatePrimaryStencil (elem);
133
134
elem_ctx.updatePrimaryIntensiveQuantities (/* timeIdx=*/ 0 );
134
135
for (unsigned dof_idx = 0 ; dof_idx < elem_ctx.numPrimaryDof (/* timeIdx=*/ 0 ); ++dof_idx) {
135
136
std::cout << " j=" << dof_idx << std::endl;
@@ -138,6 +139,7 @@ getFluidStateVariable(const std::string &name) const
138
139
unsigned global_dof_idx = elem_ctx.globalSpaceIndex (dof_idx, /* timeIdx=*/ 0 );
139
140
array[global_dof_idx] = getVariableValue_ (fs, var_type, name);
140
141
}
142
+ i+=1 ;
141
143
}
142
144
return array;
143
145
}
You can’t perform that action at this time.
0 commit comments