SundanceCellVectorExpr.hpp
Go to the documentation of this file.
00001 /* @HEADER@ */
00002 // ************************************************************************
00003 // 
00004 //                             Sundance
00005 //                 Copyright 2011 Sandia Corporation
00006 // 
00007 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
00008 // the U.S. Government retains certain rights in this software.
00009 //
00010 // Redistribution and use in source and binary forms, with or without
00011 // modification, are permitted provided that the following conditions are
00012 // met:
00013 //
00014 // 1. Redistributions of source code must retain the above copyright
00015 // notice, this list of conditions and the following disclaimer.
00016 //
00017 // 2. Redistributions in binary form must reproduce the above copyright
00018 // notice, this list of conditions and the following disclaimer in the
00019 // documentation and/or other materials provided with the distribution.
00020 //
00021 // 3. Neither the name of the Corporation nor the names of the
00022 // contributors may be used to endorse or promote products derived from
00023 // this software without specific prior written permission.
00024 //
00025 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
00026 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
00027 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
00028 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
00029 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
00030 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
00031 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
00032 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
00033 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
00034 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
00035 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00036 //
00037 // Questions? Contact Kevin Long (kevin.long@ttu.edu)
00038 // 
00039 
00040 /* @HEADER@ */
00041 
00042 #ifndef SUNDANCE_CELLVECTOREXPR_H
00043 #define SUNDANCE_CELLVECTOREXPR_H
00044 
00045 #include "SundanceEvaluatorFactory.hpp"
00046 #include "SundanceCellVectorEvaluator.hpp"
00047 #include "SundanceExpr.hpp"
00048 
00049 namespace Sundance
00050 {
00051 using namespace Sundance;
00052 using namespace Sundance;
00053 
00054 /** */
00055 enum CellVectorExprType {CellNormalVector, CellTangentSpace};
00056 
00057 /**
00058  * Expression that returns a component of a cell-determined vector for
00059  * each cell on 
00060  * which it is evaluated. This makes sense only for boundary cells. 
00061  */
00062 class CellVectorExpr
00063   : public EvaluatableExpr,
00064     public GenericEvaluatorFactory<CellVectorExpr, CellVectorEvaluator>
00065 {
00066 public:
00067   /** */
00068   CellVectorExpr(
00069     int normalComponentIndex, 
00070     int dimension,
00071     const std::string& name
00072     );
00073 
00074   /** */
00075   CellVectorExpr(
00076     int tangentBasisIndex, 
00077     int tangentComponentIndex,
00078     int dimension,
00079     const std::string& name
00080     );
00081 
00082 
00083     
00084   /** */
00085   virtual ~CellVectorExpr() {;}
00086 
00087   /** */
00088   bool isTangent() const {return type_ == CellTangentSpace;}
00089 
00090   /** */
00091   bool isNormal() const {return !isTangent();}
00092 
00093   /** */
00094   int componentIndex() const {return componentIndex_;}
00095 
00096   /** */
00097   int basisMemberIndex() const {return basisMemberIndex_;}
00098 
00099   /** */
00100   int dimension() const {return dim_;}
00101 
00102   /** */
00103   virtual XMLObject toXML() const ;
00104 
00105   const std::string& name() const {return name_;}
00106 
00107   /** Write a simple text description suitable 
00108    * for output to a terminal */
00109   virtual std::ostream& toText(std::ostream& os, bool paren) const ;
00110     
00111 
00112   /** */
00113   virtual Set<MultipleDeriv> 
00114   internalFindW(int order, const EvalContext& context) const ;
00115 
00116   /** */
00117   virtual RCP<ExprBase> getRcp() {return rcp(this);}
00118 
00119   /** Ordering operator for use in transforming exprs to standard form */
00120   virtual bool lessThan(const ScalarExpr* other) const ;
00121 private:
00122   std::string name_;
00123   int dim_;
00124   CellVectorExprType type_;
00125   int basisMemberIndex_;
00126   int componentIndex_;
00127   
00128 };
00129 
00130 
00131 /** \relates Expr \relates CellVectorExpr */
00132 Expr CellNormalExpr(int dimension, const std::string& name);
00133 
00134 /** \relates Expr \relates CellVectorExpr */
00135 Expr CellTangentExpr(int dimension, const std::string& name);
00136 
00137 }
00138 
00139 #endif

Site Contact