Skip to content
Snippets Groups Projects
Commit cbcf66c7 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

add a comment around a test copied as-is from the old Gradient plugin
parent 24a4bbd9
No related branches found
No related tags found
No related merge requests found
// $Id: Lambda2.cpp,v 1.4 2004-12-27 19:18:12 geuzaine Exp $ // $Id: Lambda2.cpp,v 1.5 2004-12-27 20:39:14 geuzaine Exp $
// //
// Copyright (C) 1997-2004 C. Geuzaine, J.-F. Remacle // Copyright (C) 1997-2004 C. Geuzaine, J.-F. Remacle
// //
...@@ -104,6 +104,10 @@ static void eigen(List_T *inList, int inNb, ...@@ -104,6 +104,10 @@ static void eigen(List_T *inList, int inNb,
int nb = List_Nbr(inList) / inNb; int nb = List_Nbr(inList) / inNb;
for(int i = 0; i < List_Nbr(inList); i += nb) { for(int i = 0; i < List_Nbr(inList); i += nb) {
// FIXME: there was this test in the old Plugin(Gradient)...
// double *yy = (double *)List_Pointer_Fast(inList, i + nbNod);
// if(yy[0] > 0){
// copy node coordinates // copy node coordinates
for(int j = 0; j < 3 * nbNod; j++) for(int j = 0; j < 3 * nbNod; j++)
List_Add(outList, List_Pointer_Fast(inList, i + j)); List_Add(outList, List_Pointer_Fast(inList, i + j));
...@@ -133,8 +137,7 @@ static void eigen(List_T *inList, int inNb, ...@@ -133,8 +137,7 @@ static void eigen(List_T *inList, int inNb,
// val contains the velocities: compute the gradient tensor // val contains the velocities: compute the gradient tensor
// from them // from them
const int MAX_NOD = 4; const int MAX_NOD = 4;
const int MAX_COMP= 3; double val[3][MAX_NOD];
double val[MAX_COMP][MAX_NOD];
for(int k = 0; k < nbNod; k++){ for(int k = 0; k < nbNod; k++){
double *v = (double *)List_Pointer_Fast(inList, i + 3 * nbNod + double *v = (double *)List_Pointer_Fast(inList, i + 3 * nbNod +
nbNod * nbComp * j + nbComp * k); nbNod * nbComp * j + nbComp * k);
...@@ -143,8 +146,8 @@ static void eigen(List_T *inList, int inNb, ...@@ -143,8 +146,8 @@ static void eigen(List_T *inList, int inNb,
} }
} }
// compute gradient of shape functions // compute gradient of shape functions
double GradPhi_x[4][3]; double GradPhi_x[MAX_NOD][3];
double GradPhi_ksi[4][3]; double GradPhi_ksi[MAX_NOD][3];
double dx_dksi[3][3]; double dx_dksi[3][3];
double dksi_dx[3][3]; double dksi_dx[3][3];
double det; double det;
...@@ -220,6 +223,10 @@ static void eigen(List_T *inList, int inNb, ...@@ -220,6 +223,10 @@ static void eigen(List_T *inList, int inNb,
} }
(*outNb)++; (*outNb)++;
// FIXME: end of the yy[0]>0 test in the old Plugin(Gradient)...
// }
} }
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment