The VectorSetAction output functor class can be used with a GenericProjector to set projection values (which must be of type Val) as coefficients of the given NumericVector. More...

#include <generic_projector.h>

Public Member Functions
	VectorSetAction (NumericVector< Val > &target_vec)

void	insert (const FEMContext &c, unsigned int var_num, const DenseVector< Val > &Ue)

Private Attributes
NumericVector< Val > &	target_vector

Detailed Description

template<typename Val>
class libMesh::VectorSetAction< Val >

The VectorSetAction output functor class can be used with a GenericProjector to set projection values (which must be of type Val) as coefficients of the given NumericVector.

Author: Roy H. Stogner

Date: 2016

Definition at line 110 of file generic_projector.h.

Constructor & Destructor Documentation

◆ VectorSetAction()

template<typename Val >

libMesh::VectorSetAction< Val >::VectorSetAction ( NumericVector< Val > & target_vec )

Definition at line 116 of file generic_projector.h.

116 :

117 target_vector(target_vec) {}

libMesh::VectorSetAction::target_vector

NumericVector< Val > & target_vector

Definition: generic_projector.h:113

Member Function Documentation

◆ insert()

template<typename Val >

void libMesh::VectorSetAction< Val >::insert	(	const FEMContext &	c,
		unsigned int	var_num,
		const DenseVector< Val > &	Ue
	)

Definition at line 119 of file generic_projector.h.

References libMesh::NumericVector< T >::first_local_index(), libMesh::DiffContext::get_dof_indices(), libMesh::NumericVector< T >::last_local_index(), libMesh::NumericVector< T >::set(), libMesh::DenseVector< T >::size(), libMesh::Threads::spin_mtx, and libMesh::VectorSetAction< Val >::target_vector.

   {
     const numeric_index_type
       first = target_vector.first_local_index(),
       last  = target_vector.last_local_index();
 
     const std::vector<dof_id_type> & dof_indices =
       c.get_dof_indices(var_num);
 
     unsigned int size = Ue.size();
 
     libmesh_assert_equal_to(size, dof_indices.size());
 
     // Lock the new vector since it is shared among threads.
     {
       Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
 
       for (unsigned int i = 0; i != size; ++i)
         if ((dof_indices[i] >= first) && (dof_indices[i] <  last))
           target_vector.set(dof_indices[i], Ue(i));
     }
   }