00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #include "SundanceOut.hpp"
00032 #include "PlayaTabs.hpp"
00033 #include "SundanceFunctionalAssemblyKernel.hpp"
00034 #include "SundanceIntegralGroup.hpp"
00035 #include "Teuchos_Time.hpp"
00036 #include "Teuchos_TimeMonitor.hpp"
00037
00038 using namespace Sundance;
00039 using namespace Teuchos;
00040 using Playa::MPIComm;
00041 using Playa::MPIDataType;
00042 using Playa::MPIOp;
00043
00044 using std::setw;
00045 using std::endl;
00046
00047 FunctionalAssemblyKernel::FunctionalAssemblyKernel(
00048 const MPIComm& comm,
00049 double* value,
00050 int verb
00051 )
00052 : AssemblyKernelBase(verb),
00053 comm_(comm),
00054 value_(value),
00055 localValue_(0.0)
00056 {}
00057
00058 void FunctionalAssemblyKernel::postLoopFinalization()
00059 {
00060 Tabs tab;
00061 SUNDANCE_MSG3(verb(), tab << "reducing functional values across processors");
00062 SUNDANCE_MSG3(verb(), tab << "local value=" << localValue_);
00063 double globalVal = localValue_;
00064
00065 comm_.allReduce((void*) &localValue_, (void*) &globalVal, 1,
00066 MPIDataType::doubleType(), MPIOp::sumOp());
00067
00068 *value_ = globalVal;
00069
00070 SUNDANCE_MSG3(verb(), tab << "reduced value = " << *value_);
00071 }
00072
00073 void FunctionalAssemblyKernel::fill(
00074 bool isBC,
00075 const IntegralGroup& group,
00076 const RCP<Array<double> >& localValues)
00077 {
00078 Tabs tab;
00079 SUNDANCE_MSG2(verb(), tab << "adding local increment " << (*localValues)[0]
00080 << " to local value" << std::endl << tab << "Before: " << localValue_);
00081
00082 localValue_ += (*localValues)[0];
00083 SUNDANCE_MSG2(verb(), tab << "After: " << localValue_);
00084 }
00085
00086