1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123
|
#include "FEDataStructures.h"
#include <assert.h>
#include <iostream>
#include <mpi.h>
Grid::Grid()
{
this->NumPoints[0] = this->NumPoints[1] = this->NumPoints[2] = 0;
this->Spacing[0] = this->Spacing[1] = this->Spacing[2] = 0;
}
void Grid::Initialize(const unsigned int numPoints[3], const double spacing[3])
{
if (numPoints[0] == 0 || numPoints[1] == 0 || numPoints[2] == 0)
{
std::cerr << "Must have a non-zero amount of points in each direction.\n";
}
for (int i = 0; i < 3; i++)
{
this->NumPoints[i] = numPoints[i];
this->Spacing[i] = spacing[i];
}
int mpiRank = 0, mpiSize = 1;
MPI_Comm_rank(MPI_COMM_WORLD, &mpiRank);
MPI_Comm_size(MPI_COMM_WORLD, &mpiSize);
this->Extent[0] = mpiRank * numPoints[0] / mpiSize;
this->Extent[1] = (mpiRank + 1) * numPoints[0] / mpiSize;
if (mpiSize != mpiRank + 1)
{
this->Extent[1]++;
}
this->Extent[2] = this->Extent[4] = 0;
this->Extent[3] = numPoints[1];
this->Extent[5] = numPoints[2];
}
unsigned int Grid::GetNumberOfLocalPoints()
{
return (this->Extent[1] - this->Extent[0] + 1) * (this->Extent[3] - this->Extent[2] + 1) *
(this->Extent[5] - this->Extent[4] + 1);
}
unsigned int Grid::GetNumberOfLocalCells()
{
return (this->Extent[1] - this->Extent[0]) * (this->Extent[3] - this->Extent[2]) *
(this->Extent[5] - this->Extent[4]);
}
void Grid::GetLocalPoint(unsigned int pointId, double* point)
{
unsigned int logicalX = pointId % (this->Extent[1] - this->Extent[0] + 1);
assert(logicalX <= this->Extent[1]);
point[0] = this->Spacing[0] * logicalX;
unsigned int logicalY =
pointId % ((this->Extent[1] - this->Extent[0] + 1) * (this->Extent[3] - this->Extent[2] + 1));
logicalY /= this->Extent[1] - this->Extent[0] + 1;
assert(logicalY <= this->Extent[3]);
point[1] = this->Spacing[1] * logicalY;
unsigned int logicalZ =
pointId / ((this->Extent[1] - this->Extent[0] + 1) * (this->Extent[3] - this->Extent[2] + 1));
assert(logicalZ <= this->Extent[5]);
point[2] = this->Spacing[2] * logicalZ;
}
unsigned int* Grid::GetNumPoints()
{
return this->NumPoints;
}
unsigned int* Grid::GetExtent()
{
return this->Extent;
}
double* Grid::GetSpacing()
{
return this->Spacing;
}
Attributes::Attributes()
{
this->GridPtr = NULL;
}
void Attributes::Initialize(Grid* grid)
{
this->GridPtr = grid;
}
void Attributes::UpdateFields(double time)
{
unsigned int numPoints = this->GridPtr->GetNumberOfLocalPoints();
this->Velocity.resize(numPoints * 3);
for (unsigned int pt = 0; pt < numPoints; pt++)
{
double coord[3];
this->GridPtr->GetLocalPoint(pt, coord);
this->Velocity[pt] = coord[1] * time;
}
std::fill(this->Velocity.begin() + numPoints, this->Velocity.end(), 0.);
unsigned int numCells = this->GridPtr->GetNumberOfLocalCells();
this->Pressure.resize(numCells);
std::fill(this->Pressure.begin(), this->Pressure.end(), 1.);
}
double* Attributes::GetVelocityArray()
{
if (this->Velocity.empty())
{
return NULL;
}
return &this->Velocity[0];
}
float* Attributes::GetPressureArray()
{
if (this->Pressure.empty())
{
return NULL;
}
return &this->Pressure[0];
}
|