private Surface CalculateSurface(PointKsiEta point1, PointKsiEta point2, double detJ) { Surface result = new Surface(); result.shapeFuncPc1 = CalculateShapeFunctions(point1); result.shapeFuncPc2 = CalculateShapeFunctions(point2); result.pc1 = result.shapeFuncPc1.Multiply((-alpha) * ambientTemperature); result.pc2 = result.shapeFuncPc2.Multiply((-alpha) * ambientTemperature); result.sum = result.pc1.Add(result.pc2); result.sum = result.sum.Multiply(detJ); return(result); }
private Surface CalculateSurface(PointKsiEta point1, PointKsiEta point2, double detJ) { Surface result = new Surface(); result.shapeFuncPc1 = CalculateShapeFunctions(point1); result.shapeFuncPc2 = CalculateShapeFunctions(point2); result.pc1 = result.shapeFuncPc1.Transpose().Multiply(result.shapeFuncPc1); result.pc1 = result.pc1.Multiply(alpha); result.pc2 = result.shapeFuncPc2.Transpose().Multiply(result.shapeFuncPc2); result.pc2 = result.pc2.Multiply(alpha); result.sum = result.pc1.Add(result.pc2); result.sum = result.sum.Multiply(detJ); return(result); }
private Matrix <double> CalculateShapeFunctions(PointKsiEta integralPoint) { return(Matrix <double> .Build.DenseOfRowArrays(new double[] { 0.25 * (1 - integralPoint.ksi) * (1 - integralPoint.eta), 0.25 * (1 + integralPoint.ksi) * (1 - integralPoint.eta), 0.25 * (1 + integralPoint.ksi) * (1 + integralPoint.eta), 0.25 * (1 - integralPoint.ksi) * (1 + integralPoint.eta) })); }