+ // direction
+ std::vector<double> direction;
+ direction.resize(dim*dim);
+ for(int i=0; i<dim; i++)
+ {
+ for(int j=0; j<dim; j++)
+ {
+ if (i == j)
+ direction[i*dim+j] = 1;
+ else
+ direction[i*dim+j] = 0;
+ }
+ }
+ if (args_info.transformMatrix_given) {
+ if (args_info.transformMatrix_given != dim*dim) {
+ std::cerr << "ERROR : please give the same number of values for --transfomMatrix and --spacing." << std::endl;
+ exit(-1);
+ }
+ for(int i=0; i<dim; i++)
+ for(int j=0; j<dim; j++)
+ direction[i*dim+j] = args_info.transformMatrix_arg[i*dim+j];
+ }
+