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) }));
}