EpetraExt Package Browser (Single Doxygen Collection)
Development
Loading...
Searching...
No Matches
example
model_evaluator
GLpApp
GLpApp_AdvDiffReactOptModelCreator.cpp
Go to the documentation of this file.
1
/*
2
//@HEADER
3
// ***********************************************************************
4
//
5
// EpetraExt: Epetra Extended - Linear Algebra Services Package
6
// Copyright (2011) Sandia Corporation
7
//
8
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
9
// the U.S. Government retains certain rights in this software.
10
//
11
// Redistribution and use in source and binary forms, with or without
12
// modification, are permitted provided that the following conditions are
13
// met:
14
//
15
// 1. Redistributions of source code must retain the above copyright
16
// notice, this list of conditions and the following disclaimer.
17
//
18
// 2. Redistributions in binary form must reproduce the above copyright
19
// notice, this list of conditions and the following disclaimer in the
20
// documentation and/or other materials provided with the distribution.
21
//
22
// 3. Neither the name of the Corporation nor the names of the
23
// contributors may be used to endorse or promote products derived from
24
// this software without specific prior written permission.
25
//
26
// THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
27
// EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
28
// IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
29
// PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
30
// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
31
// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
32
// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
33
// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
34
// LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
35
// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
36
// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
37
//
38
// Questions? Contact Michael A. Heroux (maherou@sandia.gov)
39
//
40
// ***********************************************************************
41
//@HEADER
42
*/
43
44
#include "
GLpApp_AdvDiffReactOptModelCreator.hpp
"
45
46
namespace
GLpApp
{
47
48
AdvDiffReactOptModelCreator::AdvDiffReactOptModelCreator
()
49
:len_x_(1.0)
50
,len_y_(1.0)
51
,local_nx_(3)
52
,local_ny_(4)
53
,geomFileBase_(
""
)
54
,np_(1)
55
,normalizeBasis_(false)
56
,beta_(0.0)
57
,reactionRate_(1.0)
58
,x0_(0.0)
59
,p0_(1.0)
60
,supportDerivatives_(true)
61
{}
62
63
void
AdvDiffReactOptModelCreator::setupCLP
(
64
Teuchos::CommandLineProcessor *clp
65
)
66
{
67
clp->setOption(
"len-x"
, &
len_x_
,
"Mesh dimension in the x direction (Overridden by --geom-file-base)."
);
68
clp->setOption(
"len-y"
, &
len_y_
,
"Mesh dimension in the y direction (Overridden by --geom-file-base)."
);
69
clp->setOption(
"local-nx"
, &
local_nx_
,
"Number of local discretization segments in the x direction (Overridden by --geom-file-base)."
);
70
clp->setOption(
"local-ny"
, &
local_ny_
,
"Number of local discretization segments in the y direction (Overridden by --geom-file-base)."
);
71
clp->setOption(
"geom-file-base"
, &
geomFileBase_
,
"Base name of geometry file to read the mesh from."
);
72
clp->setOption(
"np"
, &
np_
,
"The number of optimization parameters p (If < 0 then all of boundary is used)"
);
73
clp->setOption(
"normalize-basis"
,
"no-normalize-basis"
, &
normalizeBasis_
,
"Normalize the basis for the parameters p or not."
);
74
clp->setOption(
"beta"
, &
beta_
,
"Regularization."
);
75
clp->setOption(
"reaction-rate"
, &
reactionRate_
,
"The rate of the reaction"
);
76
clp->setOption(
"x0"
, &
x0_
,
"Initial guess for the state."
);
77
clp->setOption(
"p0"
, &
p0_
,
"Initial guess or nonminal value for optimization parameters."
);
78
clp->setOption(
"support-derivatives"
,
"no-support-derivatives"
,&
supportDerivatives_
,
"Support derivatives or not."
);
79
}
80
81
Teuchos::RCP<AdvDiffReactOptModel>
82
AdvDiffReactOptModelCreator::createModel
(
83
const
Teuchos::RCP<const Epetra_Comm> &comm
84
,std::ostream *out
85
)
const
86
{
87
return
Teuchos::rcp(
88
new
GLpApp::AdvDiffReactOptModel
(
89
comm,
beta_
,
len_x_
,
len_y_
,
local_nx_
,
local_ny_
,
geomFileBase_
.c_str()
90
,
np_
,
x0_
,
p0_
,
reactionRate_
,
normalizeBasis_
,
supportDerivatives_
91
)
92
);
93
}
94
95
}
// namespace GLpApp
GLpApp_AdvDiffReactOptModelCreator.hpp
GLpApp::AdvDiffReactOptModelCreator::len_x_
double len_x_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:70
GLpApp::AdvDiffReactOptModelCreator::createModel
Teuchos::RCP< AdvDiffReactOptModel > createModel(const Teuchos::RCP< const Epetra_Comm > &comm, std::ostream *out=NULL) const
Definition
GLpApp_AdvDiffReactOptModelCreator.cpp:82
GLpApp::AdvDiffReactOptModelCreator::setupCLP
void setupCLP(Teuchos::CommandLineProcessor *clp)
Definition
GLpApp_AdvDiffReactOptModelCreator.cpp:63
GLpApp::AdvDiffReactOptModelCreator::local_ny_
int local_ny_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:73
GLpApp::AdvDiffReactOptModelCreator::x0_
double x0_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:79
GLpApp::AdvDiffReactOptModelCreator::reactionRate_
double reactionRate_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:77
GLpApp::AdvDiffReactOptModelCreator::beta_
double beta_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:78
GLpApp::AdvDiffReactOptModelCreator::normalizeBasis_
bool normalizeBasis_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:76
GLpApp::AdvDiffReactOptModelCreator::AdvDiffReactOptModelCreator
AdvDiffReactOptModelCreator()
Definition
GLpApp_AdvDiffReactOptModelCreator.cpp:48
GLpApp::AdvDiffReactOptModelCreator::geomFileBase_
std::string geomFileBase_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:74
GLpApp::AdvDiffReactOptModelCreator::p0_
double p0_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:80
GLpApp::AdvDiffReactOptModelCreator::supportDerivatives_
bool supportDerivatives_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:81
GLpApp::AdvDiffReactOptModelCreator::local_nx_
int local_nx_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:72
GLpApp::AdvDiffReactOptModelCreator::np_
int np_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:75
GLpApp::AdvDiffReactOptModelCreator::len_y_
double len_y_
Definition
GLpApp_AdvDiffReactOptModelCreator.hpp:71
GLpApp::AdvDiffReactOptModel
PDE-constrained inverse problem based on a 2D discretization of a diffusion/reaction system.
Definition
GLpApp_AdvDiffReactOptModel.hpp:165
GLpApp
Definition
GLpApp_AdvDiffReactOptModel.cpp:93
Generated by
1.10.0