void fill(bool isBC,
   const IntegralGroup& group,
   const RCP<Array<double> >& localValues) 
   {
     if (group.isOneForm())
     {
       vecKernel_->fill(isBC, group, localValues);
     }
     else if (group.isZeroForm())
     {
       funcKernel_->fill(isBC, group, localValues);
     }
     else
     {
       TEST_FOR_EXCEPT(group.isTwoForm());
     }
   }
void VectorAssemblyKernel::fill(
  bool isBC, 
  const IntegralGroup& group,
  const RCP<Array<double> >& localValues) 
{
  Tabs tab0;
  SUNDANCE_MSG1(verb(), tab0 << "in VectorAssemblyKernel::fill()");

  TEUCHOS_TEST_FOR_EXCEPT(!group.isOneForm());

  bool useCofacets = group.usesMaximalCofacets();

  if (group.isOneForm())
  {
    insertLocalVectorBatch(isBC, useCofacets, 
      group.testID(), group.testBlock(), group.mvIndices(),
      *localValues);
  }

  SUNDANCE_MSG1(verb(), tab0 << "done VectorAssemblyKernel::fill()");
}