"Fossies" - the Fresh Open Source Software Archive  

Source code changes of the file "Kernel/Operation_ChangeOfCoordinates.cpp" between
getdp-3.4.0-source.tgz and getdp-3.5.0-source.tgz

About: GetDP is a general finite element solver using mixed elements to discretize de Rham-type complexes in one, two and three dimensions.

[ To the main getdp source changes report ]

Operation_ChangeOfCoordinates.cpp  (getdp-3.4.0-source.tgz):Operation_ChangeOfCoordinates.cpp  (getdp-3.5.0-source.tgz)
// GetDP - Copyright (C) 1997-2021 P. Dular and C. Geuzaine, University of Liege // GetDP - Copyright (C) 1997-2022 P. Dular and C. Geuzaine, University of Liege
// //
// See the LICENSE.txt file for license information. Please report all // See the LICENSE.txt file for license information. Please report all
// issues on https://gitlab.onelab.info/getdp/getdp/issues. // issues on https://gitlab.onelab.info/getdp/getdp/issues.
// //
// Contributor(s): // Contributor(s):
// Johan Gyselinck // Johan Gyselinck
// Ruth Sabariego // Ruth Sabariego
// //
#include <stdio.h> #include <stdio.h>
skipping to change at line 24 skipping to change at line 24
#include "ProData.h" #include "ProData.h"
#include "DofData.h" #include "DofData.h"
#include "GeoData.h" #include "GeoData.h"
#include "Pos_Search.h" #include "Pos_Search.h"
#include "SolvingOperations.h" #include "SolvingOperations.h"
#include "Cal_Quantity.h" #include "Cal_Quantity.h"
#include "ExtendedGroup.h" #include "ExtendedGroup.h"
#include "BF.h" #include "BF.h"
#include "Message.h" #include "Message.h"
extern struct Problem Problem_S ; extern struct Problem Problem_S;
extern struct CurrentData Current ; extern struct CurrentData Current;
/* ------------------------------------------------------------------------ */ /* ------------------------------------------------------------------------ */
/* O p e r a t i o n _ C h a n g e O f C o o r d i n a t e s */ /* O p e r a t i o n _ C h a n g e O f C o o r d i n a t e s */
/* ------------------------------------------------------------------------ */ /* ------------------------------------------------------------------------ */
void Operation_ChangeOfCoordinates(struct Resolution * Resolution_P, void Operation_ChangeOfCoordinates(struct Resolution *Resolution_P,
struct Operation * Operation_P, struct Operation *Operation_P,
struct DofData * DofData_P0, struct DofData *DofData_P0,
struct GeoData * GeoData_P0) struct GeoData *GeoData_P0)
{ {
int i, Nbr_Node, Num_Node, ThisIsTheNode ; int i, Nbr_Node, Num_Node, ThisIsTheNode;
double x=0., y=0., z=0.; double x = 0., y = 0., z = 0.;
struct Value Value, Value1, Value2 ; struct Value Value, Value1, Value2;
struct Group * Group_P ; struct Group *Group_P;
// Note: Current.{u,v,w} is not defined, so we cannot interpolate expressions // Note: Current.{u,v,w} is not defined, so we cannot interpolate expressions
// in the reference element. We thus set Current.Element=0 and rely on // in the reference element. We thus set Current.Element=0 and rely on
// Current.{x,y,z}. // Current.{x,y,z}.
struct Element *old = Current.Element; struct Element *old = Current.Element;
Current.Element = 0; Current.Element = 0;
Group_P = (struct Group *) Group_P = (struct Group *)List_Pointer(
List_Pointer(Problem_S.Group, Problem_S.Group, Operation_P->Case.ChangeOfCoordinates.GroupIndex);
Operation_P->Case.ChangeOfCoordinates.GroupIndex) ; if(!Group_P->ExtendedList) Generate_ExtendedGroup(Group_P);
if (!Group_P->ExtendedList) Generate_ExtendedGroup(Group_P) ; if(Group_P->FunctionType != NODESOF)
if (Group_P->FunctionType != NODESOF) Message::Error(
Message::Error("ChangeOfCoordinates: Group must be of NodesOf function type" "ChangeOfCoordinates: Group must be of NodesOf function type");
) ;
Nbr_Node = List_Nbr(Group_P->ExtendedList);
Nbr_Node = List_Nbr(Group_P->ExtendedList) ;
for(i = 0; i < Nbr_Node; i++) {
for (i=0 ; i < Nbr_Node ; i++) { List_Read(Group_P->ExtendedList, i, &Num_Node);
List_Read(Group_P->ExtendedList, i, &Num_Node) ;
if(!Group_P->InitialSuppList ||
if (!Group_P->InitialSuppList || !List_Search(Group_P->ExtendedSuppList, &Num_Node, fcmp_int)) {
! List_Search(Group_P->ExtendedSuppList, &Num_Node, fcmp_int)) { Geo_GetNodesCoordinates(1, &Num_Node, &Current.x, &Current.y, &Current.z);
Geo_GetNodesCoordinates(1, &Num_Node, &Current.x, &Current.y, &Current.z) if(Operation_P->Case.ChangeOfCoordinates.ExpressionIndex2 >= 0 &&
; Num_Node == Operation_P->Case.ChangeOfCoordinates.NumNode) {
x = Current.x;
if (Operation_P->Case.ChangeOfCoordinates.ExpressionIndex2 >= 0 && y = Current.y;
Num_Node == Operation_P->Case.ChangeOfCoordinates.NumNode) { z = Current.z;
x = Current.x ; Get_ValueOfExpressionByIndex(
y = Current.y ; Operation_P->Case.ChangeOfCoordinates.ExpressionIndex2, NULL, 0., 0.,
z = Current.z ; 0., &Value1);
Get_ValueOfExpressionByIndex ThisIsTheNode = 1;
(Operation_P->Case.ChangeOfCoordinates.ExpressionIndex2,
NULL, 0., 0., 0., &Value1) ;
ThisIsTheNode = 1 ;
} else ThisIsTheNode = 0 ;
Get_ValueOfExpressionByIndex
(Operation_P->Case.ChangeOfCoordinates.ExpressionIndex,
NULL, 0., 0., 0., &Value) ;
if (ThisIsTheNode) {
Current.x = Value.Val[0] ;
Current.y = Value.Val[1] ;
Current.z = Value.Val[2] ;
Get_ValueOfExpressionByIndex
(Operation_P->Case.ChangeOfCoordinates.ExpressionIndex2,
NULL, 0., 0., 0., &Value2) ;
Message::Debug("before x %e y %e z %e ||| after x %e y %e z %e\n",
x, y, z, Value.Val[0], Value.Val[1], Value.Val[2]);
Message::Debug(" before %e after %e", Value1.Val[0], Value2.Val[0]) ;
} }
else
ThisIsTheNode = 0;
Geo_SetNodesCoordinates(1, &Num_Node, Get_ValueOfExpressionByIndex(
&Value.Val[0], &Value.Val[1], &Value.Val[2]) ; Operation_P->Case.ChangeOfCoordinates.ExpressionIndex, NULL, 0., 0., 0.,
} &Value);
if(ThisIsTheNode) {
Current.x = Value.Val[0];
Current.y = Value.Val[1];
Current.z = Value.Val[2];
Get_ValueOfExpressionByIndex(
Operation_P->Case.ChangeOfCoordinates.ExpressionIndex2, NULL, 0., 0.,
0., &Value2);
Message::Debug("before x %e y %e z %e ||| after x %e y %e z %e\n",
x, y, z, Value.Val[0], Value.Val[1], Value.Val[2]);
Message::Debug(" before %e after %e", Value1.Val[0], Value2.Val[0]);
}
Geo_SetNodesCoordinates(1, &Num_Node, &Value.Val[0], &Value.Val[1],
&Value.Val[2]);
}
} }
Current.Element = old; Current.Element = old;
Free_SearchGrid(&Current.GeoData->Grid); Free_SearchGrid(&Current.GeoData->Grid);
} }
/* ------------------------------------------------------------------------ */ /* ------------------------------------------------------------------------ */
/* O p e r a t i o n _ D e f o r m e M e s h */ /* O p e r a t i o n _ D e f o r m e M e s h */
/* ------------------------------------------------------------------------ */ /* ------------------------------------------------------------------------ */
void Operation_DeformMesh(struct Resolution * Resolution_P, void Operation_DeformMesh(struct Resolution *Resolution_P,
struct Operation * Operation_P, struct Operation *Operation_P,
struct DofData * DofData_P0, struct DofData *DofData_P0,
struct GeoData * GeoData_P0) struct GeoData *GeoData_P0)
{ {
int Num_Node, NumBF_X=-1, NumBF_Y=-1, NumBF_Z=-1 ; int Num_Node, NumBF_X = -1, NumBF_Y = -1, NumBF_Z = -1;
double Value, un_x = 0., un_y = 0., un_z = 0. ; double Value, un_x = 0., un_y = 0., un_z = 0.;
struct Group *Group_P = (Operation_P->Case.DeformMesh.GroupIndex >=0) ? struct Group *Group_P =
(struct Group *)List_Pointer(Problem_S.Group, Operation_P->Case.DeformMesh.G (Operation_P->Case.DeformMesh.GroupIndex >= 0) ?
roupIndex) (struct Group *)List_Pointer(Problem_S.Group,
: NULL; Operation_P->Case.DeformMesh.GroupIndex) :
NULL;
if (Group_P && Group_P->FunctionType != NODESOF) if(Group_P && Group_P->FunctionType != NODESOF)
Message::Error("DeformMesh: Group must be of NodesOf function type") ; Message::Error("DeformMesh: Group must be of NodesOf function type");
struct DefineSystem *DS = (struct DefineSystem*)List_Pointer struct DefineSystem *DS = (struct DefineSystem *)List_Pointer(
(Resolution_P->DefineSystem, Operation_P->DefineSystemIndex) ; Resolution_P->DefineSystem, Operation_P->DefineSystemIndex);
if(List_Nbr(DS->FormulationIndex) > 1) if(List_Nbr(DS->FormulationIndex) > 1)
Message::Error("DeformMesh: Only one formulation must be associated to the s Message::Error(
ystem %s", "DeformMesh: Only one formulation must be associated to the system %s",
DS->Name) ; DS->Name);
struct Formulation *FO = (struct Formulation *) List_Pointer struct Formulation *FO = (struct Formulation *)List_Pointer(
(Problem_S.Formulation, *((int *)List_Pointer(DS->FormulationIndex, 0))) ; Problem_S.Formulation, *((int *)List_Pointer(DS->FormulationIndex, 0)));
struct DofData *DofData_P = DofData_P0 + Operation_P->DefineSystemIndex ; struct DofData *DofData_P = DofData_P0 + Operation_P->DefineSystemIndex;
int i; int i;
struct DefineQuantity *DQ_P = 0, *DQ2_P = 0, *DQ3_P = 0; struct DefineQuantity *DQ_P = 0, *DQ2_P = 0, *DQ3_P = 0;
if((i = List_ISearchSeq(FO->DefineQuantity, Operation_P->Case.DeformMesh.Quant if((i = List_ISearchSeq(FO->DefineQuantity,
ity, Operation_P->Case.DeformMesh.Quantity,
fcmp_DefineQuantity_Name)) < 0) fcmp_DefineQuantity_Name)) < 0)
Message::Error("Unknown Quantity '%s' in Formulation %s", Message::Error("Unknown Quantity '%s' in Formulation %s",
Operation_P->Case.DeformMesh.Quantity, FO->Name) ; Operation_P->Case.DeformMesh.Quantity, FO->Name);
else else
DQ_P = (struct DefineQuantity *) List_Pointer(FO->DefineQuantity, i) ; DQ_P = (struct DefineQuantity *)List_Pointer(FO->DefineQuantity, i);
if(Operation_P->Case.DeformMesh.Quantity2 && if(Operation_P->Case.DeformMesh.Quantity2 &&
(i = List_ISearchSeq(FO->DefineQuantity, Operation_P->Case.DeformMesh.Quant (i = List_ISearchSeq(FO->DefineQuantity,
ity2, Operation_P->Case.DeformMesh.Quantity2,
fcmp_DefineQuantity_Name)) < 0) fcmp_DefineQuantity_Name)) < 0)
Message::Error("Unknown Quantity '%s' in Formulation %s", Message::Error("Unknown Quantity '%s' in Formulation %s",
Operation_P->Case.DeformMesh.Quantity2, FO->Name) ; Operation_P->Case.DeformMesh.Quantity2, FO->Name);
else else
DQ2_P = (struct DefineQuantity *) List_Pointer(FO->DefineQuantity, i) ; DQ2_P = (struct DefineQuantity *)List_Pointer(FO->DefineQuantity, i);
if(Operation_P->Case.DeformMesh.Quantity3 && if(Operation_P->Case.DeformMesh.Quantity3 &&
(i = List_ISearchSeq(FO->DefineQuantity, Operation_P->Case.DeformMesh.Quant (i = List_ISearchSeq(FO->DefineQuantity,
ity3, Operation_P->Case.DeformMesh.Quantity3,
fcmp_DefineQuantity_Name)) < 0) fcmp_DefineQuantity_Name)) < 0)
Message::Error("Unknown Quantity '%s' in Formulation %s", Message::Error("Unknown Quantity '%s' in Formulation %s",
Operation_P->Case.DeformMesh.Quantity3, FO->Name) ; Operation_P->Case.DeformMesh.Quantity3, FO->Name);
else else
DQ3_P = (struct DefineQuantity *) List_Pointer(FO->DefineQuantity, i) ; DQ3_P = (struct DefineQuantity *)List_Pointer(FO->DefineQuantity, i);
struct FunctionSpace *FunctionSpace_P = 0, *FunctionSpace2_P = 0, *FunctionSpa struct FunctionSpace *FunctionSpace_P = 0, *FunctionSpace2_P = 0,
ce3_P = 0; *FunctionSpace3_P = 0;
if(DQ_P) if(DQ_P)
FunctionSpace_P = (struct FunctionSpace*)List_Pointer FunctionSpace_P = (struct FunctionSpace *)List_Pointer(
(Problem_S.FunctionSpace, DQ_P->FunctionSpaceIndex) ; Problem_S.FunctionSpace, DQ_P->FunctionSpaceIndex);
if(DQ2_P) if(DQ2_P)
FunctionSpace2_P = (struct FunctionSpace*)List_Pointer FunctionSpace2_P = (struct FunctionSpace *)List_Pointer(
(Problem_S.FunctionSpace, DQ2_P->FunctionSpaceIndex) ; Problem_S.FunctionSpace, DQ2_P->FunctionSpaceIndex);
if(DQ3_P) if(DQ3_P)
FunctionSpace3_P = (struct FunctionSpace*)List_Pointer FunctionSpace3_P = (struct FunctionSpace *)List_Pointer(
(Problem_S.FunctionSpace, DQ3_P->FunctionSpaceIndex) ; Problem_S.FunctionSpace, DQ3_P->FunctionSpaceIndex);
for(i = 0 ; i < List_Nbr(FunctionSpace_P->BasisFunction); i++){ for(i = 0; i < List_Nbr(FunctionSpace_P->BasisFunction); i++) {
if((void(*)())((struct BasisFunction*)List_Pointer(FunctionSpace_P->BasisFun if((void (*)())((struct BasisFunction *)List_Pointer(
ction, FunctionSpace_P->BasisFunction, i))
i))->Function == (void(*) ->Function == (void (*)())BF_NodeX)
())BF_NodeX) NumBF_X = ((struct BasisFunction *)List_Pointer(
NumBF_X = ((struct BasisFunction*)List_Pointer(FunctionSpace_P->BasisFunct FunctionSpace_P->BasisFunction, i))
ion, i))->Num; ->Num;
if((void(*)())((struct BasisFunction*)List_Pointer(FunctionSpace_P->BasisFun if((void (*)())((struct BasisFunction *)List_Pointer(
ction, FunctionSpace_P->BasisFunction, i))
i))->Function == (void(*) ->Function == (void (*)())BF_NodeY)
())BF_NodeY) NumBF_Y = ((struct BasisFunction *)List_Pointer(
NumBF_Y = ((struct BasisFunction*)List_Pointer(FunctionSpace_P->BasisFunct FunctionSpace_P->BasisFunction, i))
ion, i))->Num; ->Num;
if((void(*)())((struct BasisFunction*)List_Pointer(FunctionSpace_P->BasisFun if((void (*)())((struct BasisFunction *)List_Pointer(
ction, FunctionSpace_P->BasisFunction, i))
i))->Function == (void(*) ->Function == (void (*)())BF_NodeZ)
())BF_NodeZ) NumBF_Z = ((struct BasisFunction *)List_Pointer(
NumBF_Z = ((struct BasisFunction*)List_Pointer(FunctionSpace_P->BasisFunct FunctionSpace_P->BasisFunction, i))
ion, i))->Num; ->Num;
if((void(*)())((struct BasisFunction*)List_Pointer(FunctionSpace_P->BasisFun if((void (*)())((struct BasisFunction *)List_Pointer(
ction, FunctionSpace_P->BasisFunction, i))
i))->Function == (void(*) ->Function == (void (*)())BF_Node)
())BF_Node) NumBF_X = ((struct BasisFunction *)List_Pointer(
NumBF_X = ((struct BasisFunction*)List_Pointer(FunctionSpace_P->BasisFunct FunctionSpace_P->BasisFunction, i))
ion, i))->Num; ->Num;
} }
if(FunctionSpace2_P){ if(FunctionSpace2_P) {
for(i = 0 ; i < List_Nbr(FunctionSpace2_P->BasisFunction); i++){ for(i = 0; i < List_Nbr(FunctionSpace2_P->BasisFunction); i++) {
if((void(*)())((struct BasisFunction*)List_Pointer(FunctionSpace2_P->Basis if((void (*)())((struct BasisFunction *)List_Pointer(
Function, FunctionSpace2_P->BasisFunction, i))
i))->Function == (void( ->Function == (void (*)())BF_Node)
*)())BF_Node) NumBF_Y = ((struct BasisFunction *)List_Pointer(
NumBF_Y = ((struct BasisFunction*)List_Pointer(FunctionSpace2_P->BasisFu FunctionSpace2_P->BasisFunction, i))
nction, i))->Num; ->Num;
} }
} }
if(FunctionSpace3_P){ if(FunctionSpace3_P) {
for(i = 0 ; i < List_Nbr(FunctionSpace3_P->BasisFunction); i++){ for(i = 0; i < List_Nbr(FunctionSpace3_P->BasisFunction); i++) {
if((void(*)())((struct BasisFunction*)List_Pointer(FunctionSpace3_P->Basis if((void (*)())((struct BasisFunction *)List_Pointer(
Function, FunctionSpace3_P->BasisFunction, i))
i))->Function == (void ->Function == (void (*)())BF_Node)
(*)())BF_Node) NumBF_Z = ((struct BasisFunction *)List_Pointer(
NumBF_Z = ((struct BasisFunction*)List_Pointer(FunctionSpace3_P->BasisFu FunctionSpace3_P->BasisFunction, i))
nction, i))->Num; ->Num;
} }
} }
for(i = 0 ; i < List_Nbr(FunctionSpace_P->DofData->DofList) ; i++){ for(i = 0; i < List_Nbr(FunctionSpace_P->DofData->DofList); i++) {
if (((struct Dof*)List_Pointer(FunctionSpace_P->DofData->DofList, i))->NumTy if(((struct Dof *)List_Pointer(FunctionSpace_P->DofData->DofList, i))
pe == NumBF_X || ->NumType == NumBF_X ||
((struct Dof*)List_Pointer(FunctionSpace_P->DofData->DofList, i))->NumTy ((struct Dof *)List_Pointer(FunctionSpace_P->DofData->DofList, i))
pe == NumBF_Y || ->NumType == NumBF_Y ||
((struct Dof*)List_Pointer(FunctionSpace_P->DofData->DofList, i))->NumTy ((struct Dof *)List_Pointer(FunctionSpace_P->DofData->DofList, i))
pe == NumBF_Z){ ->NumType == NumBF_Z) {
Num_Node =
Num_Node = ((struct Dof*)List_Pointer(FunctionSpace_P->DofData->DofList, i ((struct Dof *)List_Pointer(FunctionSpace_P->DofData->DofList, i))
))->Entity ; ->Entity;
if (!Group_P || Check_IsEntityInExtendedGroup(Group_P, Num_Node, 0)) { if(!Group_P || Check_IsEntityInExtendedGroup(Group_P, Num_Node, 0)) {
Dof_GetRealDofValue Dof_GetRealDofValue(
(FunctionSpace_P->DofData, FunctionSpace_P->DofData,
((struct Dof*)List_Pointer(FunctionSpace_P->DofData->DofList, i)) , & ((struct Dof *)List_Pointer(FunctionSpace_P->DofData->DofList, i)),
Value) ; &Value);
/* Reference mesh */ /* Reference mesh */
Geo_SetCurrentGeoData(Current.GeoData = Geo_SetCurrentGeoData(Current.GeoData =
GeoData_P0 + Operation_P->Case.DeformMesh.GeoDataI GeoData_P0 +
ndex) ; Operation_P->Case.DeformMesh.GeoDataIndex);
Geo_GetNodesCoordinates(1, &Num_Node, &un_x, &un_y, &un_z) ; Geo_GetNodesCoordinates(1, &Num_Node, &un_x, &un_y, &un_z);
/* Mesh associated to the electromechanical system */ /* Mesh associated to the electromechanical system */
if(GeoData_P0 + DofData_P->GeoDataIndex != if(GeoData_P0 + DofData_P->GeoDataIndex !=
GeoData_P0 + Operation_P->Case.DeformMesh.GeoDataIndex) GeoData_P0 + Operation_P->Case.DeformMesh.GeoDataIndex)
Geo_SetCurrentGeoData(Current.GeoData = GeoData_P0 + DofData_P->GeoDat Geo_SetCurrentGeoData(Current.GeoData =
aIndex) ; GeoData_P0 + DofData_P->GeoDataIndex);
if (((struct Dof*)List_Pointer(FunctionSpace_P->DofData->DofList, i))->N if(((struct Dof *)List_Pointer(FunctionSpace_P->DofData->DofList, i))
umType == NumBF_X){ ->NumType == NumBF_X) {
un_x += Operation_P->Case.DeformMesh.Factor * Value ; un_x += Operation_P->Case.DeformMesh.Factor * Value;
Geo_SetNodesCoordinatesX(1, &Num_Node, &un_x) ; Geo_SetNodesCoordinatesX(1, &Num_Node, &un_x);
} }
if (((struct Dof*)List_Pointer(FunctionSpace_P->DofData->DofList, i))->N if(((struct Dof *)List_Pointer(FunctionSpace_P->DofData->DofList, i))
umType == NumBF_Y){ ->NumType == NumBF_Y) {
un_y += Operation_P->Case.DeformMesh.Factor * Value ; un_y += Operation_P->Case.DeformMesh.Factor * Value;
Geo_SetNodesCoordinatesY(1, &Num_Node, &un_y) ; Geo_SetNodesCoordinatesY(1, &Num_Node, &un_y);
} }
if (((struct Dof*)List_Pointer(FunctionSpace_P->DofData->DofList, i))->N if(((struct Dof *)List_Pointer(FunctionSpace_P->DofData->DofList, i))
umType == NumBF_Z){ ->NumType == NumBF_Z) {
un_z += Operation_P->Case.DeformMesh.Factor * Value ; un_z += Operation_P->Case.DeformMesh.Factor * Value;
Geo_SetNodesCoordinatesZ(1, &Num_Node, &un_z) ; Geo_SetNodesCoordinatesZ(1, &Num_Node, &un_z);
} }
} }
} }
} }
Current.GeoData = GeoData_P0 + Operation_P->Case.DeformMesh.GeoDataIndex; Current.GeoData = GeoData_P0 + Operation_P->Case.DeformMesh.GeoDataIndex;
Free_SearchGrid(&Current.GeoData->Grid); Free_SearchGrid(&Current.GeoData->Grid);
} }
 End of changes. 38 change blocks. 
178 lines changed or deleted 177 lines changed or added
To the C++ Programs section of the getdp source changes report or the top of this page
Mainly generated by pkgdiff using rfcdiff
Home  |  About  |  Features  |  All  |  Newest  |  Dox  |  Diffs  |  RSS Feeds  |  Screenshots  |  Comments  |  Imprint  |  Privacy  |  HTTP(S)