// Membrane - A linear solver for membrane problem using FEM
// Copyright (C) 2010-2026 Eric Bechet
//
// See the LICENSE file for license information and contributions.
// Please report all bugs and problems to <bechet@cadxfem.org>.
// Description: stress for thin bodies formulation

#include "membraneStress.h"

// put here otherwise multiple definition Why
void stressData(genTensor4<double> &H, genScalar<genTensor2<double> >&Bvals, genScalar< genTensor2<double> > &tmpvals)
{
  tmpvals() = H*Bvals(); // OK as B is symmetric H(i,j,k,l)*Bvals(l,k) ??
}

void stressData(genTensor4<double> &H, genVector<genTensor2<double> > &Bvals, genVector<genTensor2<double> > &tmpvals)
{
  tmpvals.resize(Bvals.size(),false);
  for(int i=0;i<Bvals.size();i++)
    tmpvals(i) = H*(Bvals(i)); // OK as B is symmetric H(i,j,k,l)*Bvals(l,k) ??
}