Go to the documentation of this file.00001 #ifndef OPENTISSUE_DYNAMICS_FEM_FEM_CLEAR_STIFFNESS_ASSEMBLY_H
00002 #define OPENTISSUE_DYNAMICS_FEM_FEM_CLEAR_STIFFNESS_ASSEMBLY_H
00003
00004
00005
00006
00007
00008
00009
00010 #include <OpenTissue/configuration.h>
00011
00012 namespace OpenTissue
00013 {
00014 namespace fem
00015 {
00016 namespace detail
00017 {
00023 template < typename node_iterator >
00024 inline void clear_stiffness_assembly(node_iterator begin, node_iterator end)
00025 {
00026 for (node_iterator node = begin; node != end; ++node)
00027 {
00028 typedef typename node_iterator::value_type::matrix_iterator matrix_iterator;
00029 node->m_f0.clear();
00030 for (matrix_iterator Kij = node->Kbegin() ; Kij != node->Kend(); ++Kij )
00031 Kij->second.clear();
00032 }
00033 }
00034
00035 }
00036 }
00037 }
00038
00039
00040 #endif