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:
1015 bytes
|
Line | |
---|
1 | /*!\file ThicknessAbsMisfitx
|
---|
2 | * \brief: compute misfit between observations and model
|
---|
3 | */
|
---|
4 |
|
---|
5 | #include "./ThicknessAbsMisfitx.h"
|
---|
6 |
|
---|
7 | #include "../../shared/shared.h"
|
---|
8 | #include "../../include/include.h"
|
---|
9 | #include "../../toolkits/toolkits.h"
|
---|
10 | #include "../../EnumDefinitions/EnumDefinitions.h"
|
---|
11 |
|
---|
12 | void ThicknessAbsMisfitx( 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->ThicknessAbsMisfit(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.