add new function

NonLinearSolver/modelReduction/DeepMaterialNetworks.cpp

@@ -922,6 +922,21 @@ void BimaterialHomogenization::computeBimaterial(const fullMatrix<double>& C1, c
   };
 };
 
+void BimaterialHomogenization::compute(const std::vector<fullMatrix<double> >& Cin, const std::vector<double>& f, 
+                        const std::vector<SVector3>& normal,
+                        fullMatrix<double>& Ceff, 
+                        bool stiff, 
+                        std::vector<hyperFullMatrix>& DCeffDcin, 
+                        std::vector<fullMatrix<double> >& DCeffDf, 
+                        std::vector<std::vector<fullMatrix<double> > >&DCeffDnormal) const
+{
+  if (DCeffDnormal.size() != 1)
+  {
+    DCeffDnormal.resize(1);
+  }
+  compute(Cin,f,normal[0],Ceff,stiff,DCeffDcin,DCeffDf,DCeffDnormal[0]);
+};
+
 
 
 void BimaterialHomogenization::compute(const std::vector<fullMatrix<double> >& Cin, const std::vector<double>& f, const SVector3& normal,

NonLinearSolver/modelReduction/DeepMaterialNetworks.h

@@ -90,12 +90,21 @@ class Homogenization
   public:
     #ifndef SWIG
     virtual ~Homogenization(){}
+    virtual void compute(const std::vector<fullMatrix<double> >& Cin, const std::vector<double>& f, 
+                        const std::vector<SVector3>& normal,
+                        fullMatrix<double>& Ceff, 
+                        bool stiff, 
+                        std::vector<hyperFullMatrix>& DCeffDcin, 
+                        std::vector<fullMatrix<double> >& DCeffDf, 
+                        std::vector<std::vector<fullMatrix<double> > >&DCeffDnormal) const = 0; 
+                        
     virtual void compute(const std::vector<fullMatrix<double> >& Cin, const std::vector<double>& f, const SVector3& normal,
                         fullMatrix<double>& Ceff, 
                         bool stiff, 
                         std::vector<hyperFullMatrix>& DCeffDcin, 
                         std::vector<fullMatrix<double>>& DCeffDf, 
                         std::vector<fullMatrix<double>>&DCeffDnormal) const = 0;
+                        
     virtual void printInfos() const = 0;
     virtual bool withPlaneStrain() const  = 0;
     #endif //SWIG
@@ -113,6 +122,15 @@ class BimaterialHomogenization : public Homogenization
                         fullMatrix<double>& DCeffDf1, fullMatrix<double>& DCeffDf2 ) const;
     #ifndef SWIG
     virtual ~BimaterialHomogenization(){}
+    
+    virtual void compute(const std::vector<fullMatrix<double> >& Cin, const std::vector<double>& f, 
+                        const std::vector<SVector3>& normal,
+                        fullMatrix<double>& Ceff, 
+                        bool stiff, 
+                        std::vector<hyperFullMatrix>& DCeffDcin, 
+                        std::vector<fullMatrix<double> >& DCeffDf, 
+                        std::vector<std::vector<fullMatrix<double> > >&DCeffDnormal) const; 
+    
     virtual void compute(const std::vector<fullMatrix<double> >& Cin, const std::vector<double>& f, const SVector3& normal,
                         fullMatrix<double>& Ceff, 
                         bool stiff,