b2element_klt_compute_constraint.H

//------------------------------------------------------------------------
// b2element_klt_compute_constraint.H --
//
//     Compute constraint and derivatives between neighbouring elements,
//     for G1 compatibility.
//
// Copyright (c) 2016
//     SMR Engineering & Development SA
//     2502 Bienne, Switzerland
//
// All Rights Reserved.  Proprietary source code.  The contents of
// this file may not be disclosed to third parties, copied or
// duplicated in any form, in whole or in part, without the prior
// written permission of SMR.
//------------------------------------------------------------------------
 
#ifndef _B2ELEMENT_KLT_COMPUTE_CONSTRAINT_H_
#define _B2ELEMENT_KLT_COMPUTE_CONSTRAINT_H_
 
#include <cassert>
#include <cmath>
 
#include "elements/klt_shells/b2element_klt_util.H"
 
namespace b2000 { namespace klt {
 
struct ComputeConstraint {
    const size_t ndof = 12;
    DofScalar b_0;
    DofScalar b_1;
    DofScalar b_2;
    DofScalar aa;
    DofScalar bb;
    DofScalar ab;
    DofScalar q;  // Contains the constraint + derivatives
 
    ComputeConstraint(const bool first, const bool second)
        : b_0(ndof, first, second),
          b_1(ndof, first, second),
          b_2(ndof, first, second),
          aa(ndof, first, second),
          bb(ndof, first, second),
          ab(ndof, first, second),
          q(ndof, first, second) {}
 
    void compute(const double X[4][3], const double u[4][3], const bool first, const bool second);
};
 
}}  // namespace b2000::klt
 
#endif /* _B2ELEMENT_KLT_COMPUTE_CONSTRAINT_H_ */