Last change
on this file since 5286 was 5286, checked in by Mathieu Morlighem, 15 years ago |
Fixed enormous bug in Misfitx and added c/modules/ThicknessAbsMisfitx needed for Balancedthicknesses"
|
File size:
921 bytes
|
Rev | Line | |
---|
[5286] | 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( double* pJ, Elements* elements,Nodes* nodes, Vertices* vertices, Loads* loads, Materials* materials,Parameters* parameters,bool process_units){
|
---|
| 13 |
|
---|
| 14 | /*Intermediary*/
|
---|
| 15 | int i;
|
---|
| 16 | Element* element=NULL;
|
---|
| 17 |
|
---|
| 18 | /*output: */
|
---|
| 19 | double J=0;
|
---|
| 20 | double J_sum;
|
---|
| 21 |
|
---|
| 22 | /*Compute Misfit: */
|
---|
| 23 | for (i=0;i<elements->Size();i++){
|
---|
| 24 | element=(Element*)elements->GetObjectByOffset(i);
|
---|
| 25 | J+=element->ThicknessAbsMisfit(process_units);
|
---|
| 26 | }
|
---|
| 27 |
|
---|
| 28 | /*Sum all J from all cpus of the cluster:*/
|
---|
| 29 | MPI_Reduce (&J,&J_sum,1,MPI_DOUBLE,MPI_SUM,0,MPI_COMM_WORLD );
|
---|
| 30 | MPI_Bcast(&J_sum,1,MPI_DOUBLE,0,MPI_COMM_WORLD);
|
---|
| 31 | J=J_sum;
|
---|
| 32 |
|
---|
| 33 | /*Assign output pointers: */
|
---|
| 34 | *pJ=J;
|
---|
| 35 | }
|
---|
Note:
See
TracBrowser
for help on using the repository browser.