OilfieldToolsNet/opm-upscaling
4
0
Fork 0
mirror of https://github.com/OPM/opm-upscaling.git synced 2025-02-25 18:45:23 -06:00
Code Issues Packages Projects Releases Wiki Activity

Add perm-upscaling on the fly for subsamples. cpchop segfaults for the moment

Browse Source
This commit is contained in:
Håvard Berland
2010-10-26 13:29:46 +02:00
parent 8197d97cfe
commit 74a6465be1
1 changed files with 15 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

22
dune/upscaling/test/cpchop.cpp
View File

@@ -90,18 +90,26 @@ int main(int argc, char** argv)
Dune::SinglePhaseUpscaler upscaler;
upscaler.init(subparser, Dune::SinglePhaseUpscaler::Fixed, 1e-9);
//Dune::SinglePhaseUpscaler::permtensor_t upscaled_K = upscaler.upscaleSinglePhase();
//std::cout << "upscaling done, doing porosity " << std::endl;
porosities.push_back(upscaler.upscalePorosity());
//upscaled_K *= (1.0/(Dune::prefix::milli*Dune::unit::darcy));
//std::cout << sample << '\t' << porosity << std::endl; // << '\t' << upscaled_K(0,0) << std::endl;
}
Dune::SinglePhaseUpscaler::permtensor_t upscaled_K = upscaler.upscaleSinglePhase();
upscaled_K *= (1.0/(Dune::prefix::milli*Dune::unit::darcy));
porosities.push_back(upscaler.upscalePorosity());
permxs.push_back(upscaled_K(0,0));
permys.push_back(upscaled_K(1,1));
permzs.push_back(upscaled_K(2,2));
}
// Output results
// (TODO: Redirect to user-supplied filename)
for (int sample = 1; sample <= subsamples; ++sample) {
std::cout << sample << '\t' << porosities[sample-1] << std::endl;
std::cout << sample << '\t' << porosities[sample-1] <<
permxs[sample-1] << '\t' <<
permys[sample-1] << '\t' <<
permzs[sample-1] << '\t' <<
std::endl;
}
}