EpetraExt Development
Loading...
Searching...
No Matches
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
45
46namespace GLpApp {
47
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
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
81Teuchos::RCP<AdvDiffReactOptModel>
83 const Teuchos::RCP<const Epetra_Comm> &comm
84 ,std::ostream *out
85 ) const
86{
87 return Teuchos::rcp(
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
Teuchos::RCP< AdvDiffReactOptModel > createModel(const Teuchos::RCP< const Epetra_Comm > &comm, std::ostream *out=NULL) const
void setupCLP(Teuchos::CommandLineProcessor *clp)
PDE-constrained inverse problem based on a 2D discretization of a diffusion/reaction system.