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

Output all field values for subsamples

Browse Source
This commit is contained in:
Kari B. Skjerve committed 2011-11-23 15:28:07 +01:00
1 parent e5e18a8b76
commit fc992da21f
1 file changed
+37 -18
Download Patch File Download Diff File Expand all files Collapse all files
dune/common/CornerpointChopper.hpp
+37 -18
View File
@@ -335,32 +335,51 @@ namespace Dune
const std::string& keyword)
{
if (field.empty()) return;
os << keyword << '\n';
int sz = field.size();
T last = std::numeric_limits<T>::max();
int repeats = 0;
for (int i = 0; i < sz; ++i) {
T val = field[i];
if (val == last) {
++repeats;
} else {
if (repeats == 1) {
os << last << '\n';
} else if (repeats > 1) {
os << repeats << '*' << last << '\n';
}
last = val;
repeats = 1;
//int num_new_zcorn = new_ZCORN_.size();
assert(sz%20 == 0);
int num_full_rows=sz/20;
int num_extra_entries=sz-num_full_rows*20;
for (int i = 0; i < num_full_rows; ++i) {
for (int j = 0; j < 20; ++j) {
os << " " << field[8*i + j];
}
os << '\n';
}
if (repeats == 1) {
os << last << '\n';
} else if (repeats > 1) {
os << repeats << '*' << last << '\n';
for (int i = 0; i < num_extra_entries; ++i) {
os << " " << field[num_full_rows*20+i];
}
os << "/\n\n";
}
// os << keyword << '\n';
// int sz = field.size();
// T last = std::numeric_limits<T>::max();
// int repeats = 0;
// for (int i = 0; i < sz; ++i) {
// T val = field[i];
// if (val == last) {
// ++repeats;
// } else {
// if (repeats == 1) {
// os << last << '\n';
// } else if (repeats > 1) {
// os << repeats << '*' << last << '\n';
// }
// last = val;
// repeats = 1;
// }
// }
// if (repeats == 1) {
// os << last << '\n';
// } else if (repeats > 1) {
// os << repeats << '*' << last << '\n';
// }
// os << "/\n\n";
// }
template <typename T>
void filterField(const std::vector<T>& field,