source: issm/trunk/src/c/modules/DragCoefficientAbsGradientx/DragCoefficientAbsGradientx.cpp@ 13975

Last change on this file since 13975 was 13975, checked in by Mathieu Morlighem, 12 years ago

merged trunk-jpl and trunk for revision 13974

File size: 1.0 KB
Line 
1/*!\file DragCoefficientAbsGradientx
2 * \brief: compute misfit between observations and model
3 */
4
5#include "./DragCoefficientAbsGradientx.h"
6
7#include "../../shared/shared.h"
8#include "../../include/include.h"
9#include "../../toolkits/toolkits.h"
10#include "../../EnumDefinitions/EnumDefinitions.h"
11
12void DragCoefficientAbsGradientx( IssmDouble* pJ, Elements* elements,Nodes* nodes, Vertices* vertices, Loads* loads, Materials* materials,Parameters* parameters,bool process_units, int weight_index){
13
14 /*Intermediary*/
15 int i;
16 Element* element=NULL;
17
18 /*output: */
19 IssmDouble J=0;
20 IssmDouble J_sum;
21
22 /*Compute Misfit: */
23 for (i=0;i<elements->Size();i++){
24 element=dynamic_cast<Element*>(elements->GetObjectByOffset(i));
25 J+=element->DragCoefficientAbsGradient(process_units,weight_index);
26 }
27
28 /*Sum all J from all cpus of the cluster:*/
29 #ifdef _HAVE_MPI_
30 MPI_Reduce (&J,&J_sum,1,MPI_DOUBLE,MPI_SUM,0,IssmComm::GetComm() );
31 MPI_Bcast(&J_sum,1,MPI_DOUBLE,0,IssmComm::GetComm());
32 J=J_sum;
33 #endif
34
35 /*Assign output pointers: */
36 *pJ=J;
37}
Note: See TracBrowser for help on using the repository browser.