~vcs-imports/escript-finley/trunk

 

#include "Util.h"

#include "Assemble.h"

#ifdef _OPENMP

#include <omp.h>

#endif

            requireWrite(out);

            if  (out_data_type == FINLEY_NODES) {

               coupler=Paso_Coupler_alloc(nodes->degreesOfFreedomConnector,numComps);

                    /* It is not immediately clear whether coupler can be trusted with constant data so I'll assume RW */

                    /* Also, it holds pointers so it might not be safe to use on lazy data anyway?*/

                    requireWrite(in);

               Paso_Coupler_free(coupler);

            } else if  (out_data_type == FINLEY_REDUCED_NODES) {

               coupler=Paso_Coupler_alloc(nodes->degreesOfFreedomConnector,numComps);

                    requireWrite(in);   /* See comment above about coupler and const */

                    Paso_Coupler_startCollect(coupler,getDataRW(in));

                    recv_buffer=Paso_Coupler_finishCollect(coupler);

             Finley_setError(TYPE_ERROR,"Finley_Assemble_CopyNodalData: cannot copy from reduced degrees of freedom to nodes.");

           } else if (out_data_type == FINLEY_REDUCED_NODES) {

               coupler=Paso_Coupler_alloc(nodes->reducedDegreesOfFreedomConnector,numComps);

                    upperBound=Paso_Distribution_getMyNumComponents(nodes->reducedDegreesOfFreedomDistribution);

                    requireWrite(in);                   /* See comment about coupler and const */

                    Paso_Coupler_startCollect(coupler,getDataRW(in));