"Fossies" - the Fresh Open Source Software Archive  

Source code changes of the file "Kernel/Gauss_Prism.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.

Gauss_Prism.cpp  (getdp-3.4.0-source.tgz)	:	Gauss_Prism.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.
#include "Gauss_Prism.h"		#include "Gauss_Prism.h"
#include "Message.h"		#include "Message.h"
/* Gauss integration over a prism */		/* Gauss integration over a prism */
void Gauss_Prism(int Nbr_Points, int Num,		void Gauss_Prism(int Nbr_Points, int Num, double *u, double *v, double *w,
double *u, double *v, double *w, double *wght)		double *wght)
{		{
switch (Nbr_Points) {		switch(Nbr_Points) {
		case 6:
case 6 :		*u = upri6[Num];
*u = upri6 [Num] ; *v = vpri6 [Num] ; *w = wpri6 [Num] ;		*v = vpri6[Num];
*wght = ppri6 [Num] ; break ;		*w = wpri6[Num];
		*wght = ppri6[Num];
		break;
case 9 :		case 9:
*u = upri9 [Num] ; *v = vpri9 [Num] ; *w = wpri9 [Num] ;		*u = upri9[Num];
*wght = ppri9 [Num] ; break ;		*v = vpri9[Num];
		*w = wpri9[Num];
		*wght = ppri9[Num];
		break;
case 21 :		case 21:
*u = upri21[Num] ; *v = vpri21[Num] ; *w = wpri21[Num] ;		*u = upri21[Num];
*wght = ppri21[Num] ; break ;		*v = vpri21[Num];
		*w = wpri21[Num];
		*wght = ppri21[Num];
		break;
case 42 :		case 42:
*u = upri42[Num] ; *v = vpri42[Num] ; *w = wpri42[Num] ;		*u = upri42[Num];
*wght = ppri42[Num] ; break ;		*v = vpri42[Num];
		*w = wpri42[Num];
		*wght = ppri42[Num];
		break;
default :		default:
Message::Error("Wrong number of Gauss points for Prism: "		Message::Error("Wrong number of Gauss points for Prism: "
"valid choices: 6, 9, 21, 42");		"valid choices: 6, 9, 21, 42");
break;		break;
}		}
}		}
End of changes. 8 change blocks.
19 lines changed or deleted		29 lines changed or added

---