source: issm/trunk/src/c/modules/Dux/Dux.cpp@ 4091

/*!\file Dux
 * \brief: diff between observed and modeled velocity
 */

#include "./Dux.h"

#include "../../shared/shared.h"
#include "../../include/include.h"
#include "../../toolkits/toolkits.h"
#include "../../EnumDefinitions/EnumDefinitions.h"
#include "../SurfaceAreax/SurfaceAreax.h"

void Dux( Vec* pdu_g, DataSet* elements,DataSet* nodes, DataSet* vertices, DataSet* loads, DataSet* materials, Parameters* parameters){

        /*Intermediary*/
        int i;
        int gsize;
        int found;
        double fit=-1;
        double S;
        Element* element=NULL;

        int analysis_type;

        /*output: */
        Vec du_g=NULL;

        /*retrive parameters: */
        parameters->FindParam(&analysis_type,AnalysisTypeEnum);

        /*First, get elements and loads configured: */
        elements->  Configure(elements,loads, nodes, vertices, materials,parameters);
        nodes->     Configure(elements,loads, nodes, vertices, materials,parameters);
        parameters->Configure(elements,loads, nodes, vertices, materials,parameters);

        /*Compute surface area: */
        SurfaceAreax(&S,elements,nodes,vertices, loads,materials,parameters);

        /*add surface area to elements :*/
        for(i=0;i<elements->Size();i++){
                Element* element=(Element*)elements->GetObjectByOffset(i);
                element->InputUpdateFromVector(&S,SurfaceAreaEnum,ConstantEnum);
        }
        
        /*Get size of matrix: */
        gsize=nodes->NumberOfDofs(analysis_type);

        /*Allocate du_g: */
        du_g=NewVec(gsize);

        /*Compute velocity difference: */
        for (i=0;i<elements->Size();i++){
                element=(Element*)elements->GetObjectByOffset(i);
                element->Du(du_g);
        }

        /*Assemble vector: */
        VecAssemblyBegin(du_g);
        VecAssemblyEnd(du_g);

        /*Assign output pointers: */
        *pdu_g=du_g;
        
}