Skip to content

Create src #1

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
daniel4993 wants to merge 10 commits into
base: master
Choose a base branch
from
Open

Create src #1

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
0adf042
Create src
daniel4993 May 4, 2017
e517099
Add files via upload
daniel4993 May 4, 2017
9e8aaf4
warnings in the resolution
daniel4993 May 4, 2017
d0af3e3
Delete marginalpdf.jpg
daniel4993 May 4, 2017
871e019
Delete density.jpg
daniel4993 May 4, 2017
e612944
Delete drag_distribution.jpg
daniel4993 May 4, 2017
035d13d
Add files via upload
daniel4993 May 4, 2017
fe0e8c4
Delete file_sip_raw_mh1.jpg
daniel4993 May 4, 2017
3354a45
Delete fixed_fluid_density_no_fixed_infer_drag_radius_scenario.C
daniel4993 May 4, 2017
024bc57
Delete fixed_fluid_density_fixed_droplet_radius_position_infer_drag_s…
daniel4993 May 4, 2017
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Binary file removed figures/density.jpg
View file
Open in desktop
Binary file not shown.
Binary file removed figures/drag_distribution.jpg
View file
Open in desktop
Binary file not shown.
Binary file added figures/pdf_droplet_density.jpg
View file
Open in desktop
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file removed file_sip_raw_mh1.jpg
View file
Open in desktop
Binary file not shown.
117 changes: 0 additions & 117 deletions fixed_fluid_density_no_fixed_infer_drag_radius_scenario.C
View file
Open in desktop

This file was deleted.

Binary file removed marginalpdf.jpg
View file
Open in desktop
Binary file not shown.
7 changes: 4 additions & 3 deletions ...let_radius_position_infer_drag_scenario.C → ...let_radius_position_infer_drag_scenario.C
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -52,7 +52,8 @@ namespace DIP
const QUESO::BaseEnvironment & queso_env,
const DropletData<VectorType,MatrixType> & data )
: FixedFluidDensityFixedDropletRadiusPositionScenarioBase<VectorType,MatrixType>(input,data),
_n_params(1+this->_data.n_data())
//_n_params(1+this->_data.n_data())
_n_params(2)
//_n_marg_params(1)
{
// Setup parameter space, domain, and RV
Expand Down Expand Up @@ -90,7 +91,7 @@ namespace DIP
queso_assert_equal_to(model_output.sizeGlobal(), this->_data.n_data());

const double rho_d = this->_rho_d_nom*domain_vector[0];
//const double _cd = this->_cd*domain_vector[0];
const double cd = this->_cd*domain_vector[1];


for( unsigned int i = 0; i < this->_data.n_data(); i++ )
Expand All @@ -102,7 +103,7 @@ namespace DIP
model_output[i] = this->_model.position(this->_rho_f,
rho_d,
r_d,
this->_cd,
cd,
t);
}
}
Expand Down
28 changes: 28 additions & 0 deletions model/model_input.in
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,28 @@
# Input file for droplet inverse problem

[Data]
conversion_datasets = './data/slomo_1471360151_mmPixel.csv'
position_datasets = './data/slomo_1471360151_rise.csv'
frame_rate = '60'

# Number of pixels "fuzz" in a single position measurement (at the bubble boundary)
n_pixels_boundary = '10'
[]

[Droplet]
fluid_density = '997' #Water, [kg/m^3]

drag_coeff = '0.9'
drag_coeff_min = '0.8' #Normalized based on nominal radius determined by data
drag_coeff_max = '2.0' #Normalized based on nominal radius determined by data

droplet_density = '961' #Published value of caster oil, [kg/m^3]
droplet_density_min = '0.8'
droplet_density_max = '1.2'
[]

[InverseProblem]
scenario = 'fixed_rhof_fixed_rd_position_infer_cd'
solver_type = 'metropolis-hastings'
guess = '1.0 1.0'
[]
1 change: 1 addition & 0 deletions model/src
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@

21 changes: 21 additions & 0 deletions model/warnings
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,21 @@
WARNING In Algorithm<V,M>::alpha(x,y), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0: y.logTarget() = -nan, x.values() = 1.01259 0.905306 , y.values() = 1.03972 0.800748
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 740451, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.01295 0.898824 , [inputSize - 1]->values() = 1.03972 0.800748
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 740451, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.01259 0.905306 , [inputSize - 1]->values() = 1.03972 0.800748
WARNING In Algorithm<V,M>::alpha(x,y), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0: y.logTarget() = -nan, x.values() = 1.01147 0.926868 , y.values() = 1.03761 0.801581
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 741230, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.01434 0.884699 , [inputSize - 1]->values() = 1.03761 0.801581
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 741230, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.01147 0.926868 , [inputSize - 1]->values() = 1.03761 0.801581
WARNING In Algorithm<V,M>::alpha(x,y), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0: y.logTarget() = -nan, x.values() = 1.00747 1.06607 , y.values() = 1.03871 0.801129
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 742854, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.01126 0.945578 , [inputSize - 1]->values() = 1.03871 0.801129
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 742854, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.00747 1.06607 , [inputSize - 1]->values() = 1.03871 0.801129
WARNING In Algorithm<V,M>::alpha(x,y), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0: y.logTarget() = -nan, x.values() = 1.00125 1.30614 , y.values() = 1.03764 0.802089
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 743001, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 0.999488 1.44239 , [inputSize - 1]->values() = 1.03764 0.802089
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 743001, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.00125 1.30614 , [inputSize - 1]->values() = 1.03764 0.802089
WARNING In Algorithm<V,M>::alpha(x,y), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0: y.logTarget() = -nan, x.values() = 1.01384 0.820688 , y.values() = 1.04628 0.800015
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 745458, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.0163 0.809217 , [inputSize - 1]->values() = 1.04628 0.800015
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 745458, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.01384 0.820688 , [inputSize - 1]->values() = 1.04628 0.800015
WARNING In Algorithm<V,M>::alpha(x,y), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0: y.logTarget() = -nan, x.values() = 1.01325 0.850944 , y.values() = 1.04717 0.800092
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 745484, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.01224 0.846561 , [inputSize - 1]->values() = 1.04717 0.800092
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 745484, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.01325 0.850944 , [inputSize - 1]->values() = 1.04717 0.800092
WARNING In Algorithm<V,M>::alpha(x,y), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0: y.logTarget() = -nan, x.values() = 1.00906 1.0143 , y.values() = 1.04379 0.800475
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 746404, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.00747 1.07847 , [inputSize - 1]->values() = 1.04379 0.800475
WARNING In MetropolisHastingsSG<P_V,P_M>::alpha(vec), worldRank 0, fullRank 0, subEnvironment 0, subRank 0, inter0Rank 0, positionId = 746404, stageId = 1: inputSize = 2, inputPositionsData[inputSize - 1]->logTarget() = -nan, [0]->values() = 1.00906 1.0143 , [inputSize - 1]->values() = 1.04379 0.800475