SundanceVectorAssemblyKernel.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_VECTORASSEMBLYKERNEL_H
00043 #define SUNDANCE_VECTORASSEMBLYKERNEL_H
00044 
00045 #include "SundanceDefs.hpp"
00046 #include "SundanceVectorFillingAssemblyKernel.hpp"
00047 
00048 namespace Sundance
00049 {
00050 using namespace Teuchos;
00051 
00052 
00053 /**
00054  * VectorAssemblyKernel builds load vectors. It derives from VectorFillingAssemblyKernel,
00055  * which implements the functions shared between all kernels that build vectors,
00056  * such as FunctionalAndGradientAssemblyKernel.
00057  */
00058 class VectorAssemblyKernel : public VectorFillingAssemblyKernel
00059 {
00060 public:
00061   
00062   /**
00063    * Ctor takes several arguments:
00064    * \param dofMap is an array of DOFMap ptrs, one for each block 
00065    *
00066    * \param isBCIndex is an array of ptrs to arrays of ints (bools). The value 
00067    * (*isBCIndex[b])[d] indicates whether dof #d in block #b is or is not 
00068    * an essential BC dof. 
00069    *
00070    * \param lowestLocalIndex stores the lowest locally-owned DOF index for each 
00071    * block 
00072    * 
00073    * \param b multivector to be filled
00074    *
00075    * \param partitionBC whether dirichlet BCs are stored in a separate block
00076    *
00077    * \param verb verbosity level
00078    */
00079   VectorAssemblyKernel(
00080     const Array<RCP<DOFMapBase> >& dofMap,
00081     const Array<RCP<Array<int> > >& isBCIndex,
00082     const Array<int>& lowestLocalIndex,
00083     Array<Vector<double> >& b,
00084     bool partitionBCs,
00085     int verb
00086     );
00087 
00088   /** */
00089   virtual ~VectorAssemblyKernel(){;}
00090 
00091   /** */
00092   virtual void prepareForWorkSet(
00093     const Array<Set<int> >& requiredTests,
00094     const Array<Set<int> >& requiredUnks,
00095     RCP<StdFwkEvalMediator> mediator) ;
00096 
00097   /** */
00098   virtual void fill(bool isBC,
00099     const IntegralGroup& group,
00100     const RCP<Array<double> >& localValues) ;   
00101 
00102 };
00103 
00104 }
00105 
00106 
00107 
00108 #endif

Site Contact