a_ocriteria_ex.cxx

Go to the documentation of this file.

/*
Copyright 2010-2011 Pierre SMARS (smars@yuntech.edu.tw)
This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 2 of the License, or
(at your option) any later version.
 
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.
 
You should have received a copy of the GNU General Public License
along with this program.  If not, see <http://www.gnu.org/licenses/>.
*/
#include "a_ocriteria_ex.h"
#include "a_structure.h"
#include "a_face_2d4.h"
 
#include "a_debug.h"
 
//---------------------------------------------------------------------------
double a_ocriteria_ex::f(const a_structure * s) const
{
        int nf = 0;
        double etot;
        for (int i=0; i<s->nb(); i++)
        {
                const a_block * b = dynamic_cast<const a_block *>(s->b(i));
                for (int j=0; j<b->ni(); j++)
                {
//CL(.)
                        a_face_2d4 * f = cf_2d4(b->i(j));
                        double ex = fabs(f->exrx());
                        double ey = fabs(f->exry());
                        if (ey>ex)
                                ex = ey;
                        etot += ex;
                        nf++;
                }
        }
//VE(etot)
//VE(nf)
//VE(etot/nf)
        return etot/nf;
}
//---------------------------------------------------------------------------