/*!\file ProcessParamsx
 * \brief: process parameters using partitioning vector. 
 * Go through all parameters in the 'parameters' dataset. For each parameter that holds a doublevec (ie, a double* vector synchronized across 
 * the MPI ring of a cluster), partition the vector so that the new node partitioning decided in ModelProcessor is applied. Otherwise, 
 * parameters coming directly from Matlab would be serially partitioned, which means could not be recognized by each individual node or element. 
 * The partition needs to be the parallel partitionting.
 */

#include "./ProcessParamsx.h"

#undef __FUNCT__ 
#define __FUNCT__ "ProcessParamsx"

#include "../shared/shared.h"
#include "../include/macros.h"
#include "../toolkits/toolkits.h"
#include "../DataSet/DataSet.h"
#include "../EnumDefinitions/EnumDefinitions.h"

void ProcessParamsx( DataSet* parameters, Vec  part){

	
	double* partition=NULL;
	Param*  param=NULL;

	/*diverse: */
	int     numberofnodes;
	int i,j,k;
	int ndof;

	/*processing of parameters: */
	double* parameter=NULL;
	double* new_parameter=NULL;

	/*Some parameters needed: */
	parameters->FindParam(&numberofnodes,"numberofnodes");

	/*serialize partition vector: */
	if(part)VecToMPISerial(&partition,part);

	for(i=0;i<parameters->Size();i++){

		Param* param=(Param*)parameters->GetObjectByOffset(i);

		if (param->GetType()==DOUBLEVEC){

			ndof=param->GetNdof();

			if (ndof!=0){ /*ok, we are dealing with a parameter that needs to be repartitioned, for ndof degrees of freedom: */

				new_parameter=(double*)xcalloc(numberofnodes*ndof,sizeof(double));
				param->GetParameterValue(&parameter);

				for(j=0;j<ndof;j++){

					for(k=0;k<numberofnodes;k++){
					
						new_parameter[(int)(ndof*partition[k]+j)]=parameter[ndof*k+j];

					}

				}
			
				
				/*Now, replace old parameter with new parameter: */
				param->SetDoubleVec(new_parameter,ndof*numberofnodes,ndof);

				/*Free ressources: */
				xfree((void**)&new_parameter);
				xfree((void**)&parameter);

			}
		}
	}
	xfree((void**)&partition);

}