get number of subfields in a SubdomainField

src/field/SubdomainField.cpp

@@ -62,6 +62,12 @@ namespace gmshddm
       return T_Form;
     }
 
+    template< class T_Scalar, gmshfem::field::Form T_Form >
+    unsigned long long SubdomainField< T_Scalar, T_Form >::size() const
+    {
+      return _fields.size();
+    }
+
     template< class T_Scalar, gmshfem::field::Form T_Form >
     const gmshfem::field::Field< T_Scalar, T_Form > &SubdomainField< T_Scalar, T_Form >::operator()(const unsigned int i) const
     {

src/field/SubdomainField.h

@@ -38,6 +38,7 @@ namespace gmshddm
       SubdomainField &operator=(const SubdomainField< T_Scalar, T_Form > &other);
 
       gmshfem::field::Form form() const;
+      unsigned long long size() const;
 
       const gmshfem::field::Field< T_Scalar, T_Form > &operator()(const unsigned int i) const;
       gmshfem::field::Field< T_Scalar, T_Form > &operator()(const unsigned int i);