Commit 15962e03 authored by Backofen, Rainer's avatar Backofen, Rainer
Browse files

definition change in usage of GridWriter, details ask rainer

parent 8c0fc620
...@@ -38,9 +38,9 @@ namespace AMDiS { ...@@ -38,9 +38,9 @@ namespace AMDiS {
* Produces a grid-based output of a dof-vector or vector of dof-vectors * Produces a grid-based output of a dof-vector or vector of dof-vectors
* \param p array of world coordinates. * \param p array of world coordinates.
* - p[0] defines origin of grid-system. * - p[0] defines origin of grid-system.
* - p[1] defines first basis vector by p[1] - p[0]. * - p[1] defines first basis vector by p[1] relatively to p[0].
* - p[2] defines second basis vector by p[2] - p[0] (dim >= 2). * - p[2] defines second basis vector by p[2] relatively to p[0]. (dim >= 2).
* - p[3] defines third basis vector by p[3] - p[0] (dim == 3). * - p[3] defines third basis vector by p[3] relatively to p[0]. (dim == 3).
* \param numPoints number of grid points for the different dirctions * \param numPoints number of grid points for the different dirctions
* - numPoints[0] number of grid points in the first direction * - numPoints[0] number of grid points in the first direction
* - numPoints[1] number of grid points in the second direction * - numPoints[1] number of grid points in the second direction
......
...@@ -48,7 +48,7 @@ namespace AMDiS { ...@@ -48,7 +48,7 @@ namespace AMDiS {
// get basis of grid // get basis of grid
for (int i = 0; i < dim; i++) { for (int i = 0; i < dim; i++) {
for (int j = 0; j < dim; j++) { for (int j = 0; j < dim; j++) {
basis[i][j] = p[i + 1][j] - p[0][j]; basis[i][j] = p[i + 1][j] ;
lengthBasis[i] += basis[i][j] * basis[i][j]; lengthBasis[i] += basis[i][j] * basis[i][j];
} }
lengthBasis[i] = sqrt(lengthBasis[i]); lengthBasis[i] = sqrt(lengthBasis[i]);
...@@ -160,7 +160,7 @@ namespace AMDiS { ...@@ -160,7 +160,7 @@ namespace AMDiS {
// get basis of grid // get basis of grid
for (int i = 0; i < dim; i++) { for (int i = 0; i < dim; i++) {
for (int j = 0; j < dim; j++) { for (int j = 0; j < dim; j++) {
basis[i][j] = p[i + 1][j] - p[0][j]; basis[i][j] = p[i + 1][j] ;
lengthBasis[i] += basis[i][j] * basis[i][j]; lengthBasis[i] += basis[i][j] * basis[i][j];
} }
lengthBasis[i] = sqrt(lengthBasis[i]); lengthBasis[i] = sqrt(lengthBasis[i]);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment