some information output

for debugging puropose.

opm/autodiff/BlackoilModelBase_impl.hpp

@@ -870,6 +870,13 @@ namespace detail {
         // 4. Store the results
         well_perforation_densities_ = Eigen::Map<const V>(cd.data(), nperf);
         well_perforation_pressure_diffs_ = Eigen::Map<const V>(cdp.data(), nperf);
+#if 0
+        std::cout << "well_perforation_densities_ " << std::endl;
+        std::cout << well_perforation_densities_ << std::endl;
+
+        std::cout << "well_perforation_pressure_diffs_ " << std::endl;
+        std::cout << well_perforation_pressure_diffs_ << std::endl;
+#endif
     }
 
 
@@ -2191,6 +2198,15 @@ namespace detail {
             const V dbhp_limited = sign(dbhp) * dbhp.abs().min(bhp_old.abs()*dpmaxrel);
             const V bhp = bhp_old - dbhp_limited;
             std::copy(&bhp[0], &bhp[0] + bhp.size(), well_state.bhp().begin());
+#if 0
+            std::cout << " output wr in updateWellState " << std::endl;
+            std::cout << wr << std::endl;
+            std::cin.ignore();
+
+            std::cout << " output bhp is updateWellState " << std::endl;
+            std::cout << bhp << std::endl;
+            std::cin.ignore();
+#endif
 
             //Get gravity for THP hydrostatic correction
             const double gravity = detail::getGravity(geo_.gravity(), UgGridHelpers::dimensions(grid_));

opm/autodiff/WellStateFullyImplicitBlackoil.hpp

@@ -32,6 +32,7 @@
 #include <map>
 #include <algorithm>
 #include <array>
+#include <iostream>
 
 namespace Opm
 {
@@ -167,6 +168,42 @@ namespace Opm
                     }
                 }
             }
+#if 0
+            // Debugging output.
+            std::cout << " output all the well state informations after initialization " << std::endl;
+            const int nperf_total = nperf;
+
+            std::cout << " number of phase : " << np << " nubmer of perforations " << nperf_total << std::endl;
+
+            std::cout << " bhps : " << std::endl;
+            for (int i = 0; i < nw; ++i) {
+                std::cout << bhp()[i] << std::endl;
+            }
+
+            std::cout << " thps : " << std::endl;
+
+            for (int i = 0; i < nw; ++i) {
+                std::cout << thp()[i] << std::endl;
+            }
+
+            std::cout << " well rates " << std::endl;
+            for (int i = 0; i < nw; ++i) {
+                std::cout << i;
+                for (int p = 0; p < np; ++p) {
+                    std::cout << " " << wellRates()[np * i + p];
+                }
+                std::cout << std::endl;
+            }
+
+            std::cout << " perf pressures and pref phase rates : " << std::endl;
+            for (int i = 0; i < nperf_total; ++i) {
+                std::cout << i << " " << perfPress()[i];
+                for (int p = 0; p < np; ++p) {
+                    std::cout << " " << perfPhaseRates()[np * i + p];
+                }
+                std::cout << std::endl;
+            }
+#endif
         }
 
         /// One rate per phase and well connection.

opm/autodiff/WellStateMultiSegment.hpp

@@ -428,7 +428,6 @@ namespace Opm
                 }
             }
             std::cout << " output the well state right after intialization is DONE! " << std::endl;
-            std::cin.ignore();
 #endif
         }