// Translate this file with
//
// g++ -O3 --std=c++11 assignment-2019.c -o assignment
//
// Run it with
//
// ./assignment
//
// There should be a result.pvd file that you can open with Paraview.
// Sometimes, Paraview requires to select the representation "Point Gaussian"
// to see something meaningful.
//
// (C) 2019 Tobias Weinzierl

#include <fstream>
#include <sstream>
#include <iostream>
#include <string>
#include <math.h>
#include <limits>
#include <stdalign.h>
// #pragma simd

double t       = 0;
double tFinal = 0;
double tPlot    = 0;
double tPlotDelta = 0;

int NumberOfBodies = 0;

/**
* Pointer to pointers. Each pointer in turn points to three coordinates, i.e.
* each pointer represents one molecule/particle/body.
*/
double** x;

/**
* Equivalent to x storing the velocities.
*/
double** v;

/**
* Global time step size used.
*/
double timeStepSize = 0.0001;

/**
* Maximum velocity of all particles.
*/
double maxV;

/**
* Minimum distance between two elements.
*/
double minDx;

/**
* Set up scenario from the command line.
*
* This operation is not to be changed in the assignment.
*/
void setUp(int argc, char** argv) {
NumberOfBodies = (argc-2) / 6;

x = new double*[NumberOfBodies];
v = new double*[NumberOfBodies];

int readArgument = 1;

tPlotDelta = std::stof(argv[readArgument]); readArgument++;
tFinal    = std::stof(argv[readArgument]); readArgument++;

for (int i=0; i<NumberOfBodies; i++) {
x[i] = new double[3];
v[i] = new double[3];

x[i][0] = std::stof(argv[readArgument]); readArgument++;
x[i][1] = std::stof(argv[readArgument]); readArgument++;
x[i][2] = std::stof(argv[readArgument]); readArgument++;

v[i][0] = std::stof(argv[readArgument]); readArgument++;
v[i][1] = std::stof(argv[readArgument]); readArgument++;
v[i][2] = std::stof(argv[readArgument]); readArgument++;
}

std::cout << "created setup with " << NumberOfBodies << " bodies" << std::endl;

if (tPlotDelta<=0.0) {
std::cout << "plotting switched off" << std::endl;
tPlot = tFinal + 1.0;
}
else {
std::cout << "plot initial setup plus every " << tPlotDelta << " time units" << std::endl;
tPlot = 0.0;
}
}

std::ofstream videoFile;

/**
* This operation is not to be changed in the assignment.
*/
void openParaviewVideoFile() {
videoFile.open( "result.pvd" );
videoFile << "<?xml version=\"1.0\"?>" << std::endl
         << "<VTKFile type=\"Collection\" version=\"0.1\" byte_order=\"LittleEndian\" compressor=\"vtkZLibDataCompressor\">" << std::endl
         << "<Collection>";
}

/**
* This operation is not to be changed in the assignment.
*/
void closeParaviewVideoFile() {
videoFile << "</Collection>"
         << "</VTKFile>" << std::endl;
}

/**
* The file format is documented at http://www.vtk.org/wp-content/uploads/2015/04/file-formats.pdf
*
* This operation is not to be changed in the assignment.
*/
void printParaviewSnapshot() {
static int counter = -1;
counter++;
std::stringstream filename;
filename << "result-" << counter << ".vtp";
std::ofstream out( filename.str().c_str() );
out << "<VTKFile type=\"PolyData\" >" << std::endl
   << "<PolyData>" << std::endl
   << " <Piece NumberOfPoints=\"" << NumberOfBodies << "\">" << std::endl
   << " <Points>" << std::endl
   << " <DataArray type=\"Float32\" NumberOfComponents=\"3\" format=\"ascii\">";

for (int i=0; i<NumberOfBodies; i++) {
out << x[i][0]
   << " "
   << x[i][1]
   << " "
   << x[i][2]
   << " ";
}

out << " </DataArray>" << std::endl
   << " </Points>" << std::endl
   << " </Piece>" << std::endl
   << "</PolyData>" << std::endl
   << "</VTKFile>" << std::endl;

videoFile << "<DataSet timestep=\"" << counter << "\" group=\"\" part=\"0\" file=\"" << filename.str() << "\"/>" << std::endl;
}

Superhero fight Part 1: Powers.h â€¢ This file has 4 sections: 1...