OpenCASCADE Gauss Integration

eryar@163.com

Abstract. Numerical integration is the approximate computation of an integral using numerical techniques. The numerical computation of an integral is sometimes called quadrature. The most straightforward numerical integration technique uses the Newton-Cotes formulas(also called quadrature formulas), which approximate a function tabulated sequence of regularly spaced intervals by various degree polynomials. If the functions are known analytically instead of being tabulated at equally spaced intervals, the best numerical method of integrations is called Gauss Integration(Gaussian quadrature). By picking the abscissas at which to evaluate the function, Gaussian quadrature produces the most accurate approximations possible. In OpenCASCADE math package it implement the Gauss-Legendre integration. So I will focus on the usage of the class in OpenCASCADE.

Key Words. OpenCASCADE, Gauss Integration, Gauss-Legendre, Numerical Analysis

1. Introduction

在科学和工程计算问题中,经常要计算一些定积分或微分,它们的精确值无法算出或计算量太大,只能用数值的方法给出具有指定误差限的近似值。最直观的数值积分方法有Newton-Cotes,其将积分区间等分之,并取分点为积分节点。这种做法虽然简化了计算,但却降低了所得公式的代数精度。

Gauss型求积公式是一种高精度的数值积分公式。在求积节点数相同的情况下,即计算工作量相近的情况下,利用Gauss型求积公式往往可以获得准确程序较高的积分结果,只是它在不等距的无理数上计算被积函数。

OpenCASCADE的math包中实现了Gauss-Legendre积分算法。本文主要介绍其使用方法,进而对其应用进行理解。

2. The Gauss-Legendre Integration

Gauss型求积公式是数值稳定的,且对有限闭区间上的连续函数,Gauss求积的数值随节点数目的增加而收敛到准确积分值。

常用的Gauss型求积公式有Gauss-Legendre求积公式,Gauss-Chebyshev求积公式,Gauss-Laguerre求积公式和Gauss-Hermite求积公式等。

对于一般区间[a, b]上的Gauss型求积公式,可通过变量变换,由Gauss-Legendre求积公式得到:

其中:

OpenCASCADE中对应的类有math_GaussSingleIntegration,主要实现的函数为Perform(),计算过程如下:

v 查表求得Gauss点及求积系数;

//Recuperation des points de Gauss dans le fichier GaussPoints.
math::GaussPoints(Order,GaussP);
math::GaussWeights(Order,GaussW);

v 根据Gauss-Legendre求积公式计算;

// Changement de variable pour la mise a l'echelle [Lower, Upper] :
xm = 0.5*(Upper + Lower);
xr = 0.5*(Upper - Lower);
Val = .; Standard_Integer ind = Order/, ind1 = (Order+)/;
if(ind1 > ind) { // odder case
Ok1 = F.Value(xm, Val);
if (!Ok1) return;
Val *= GaussW(ind1);
}
// Sommation sur tous les points de Gauss: avec utilisation de la symetrie.
for (j = ; j <= ind; j++) {
dx = xr*GaussP(j);
Ok1 = F.Value(xm-dx, F1);
if(!Ok1) return;
Ok1 = F.Value(xm+dx, F2);
if(!Ok1) return;
// Multiplication par les poids de Gauss.
Standard_Real FT = F1+F2;
Val += GaussW(j)*FT;
}
// Mise a l'echelle de l'intervalle [Lower, Upper]
Val *= xr;

对比Gauss-Legendre求积公式来理解上述代码还是比较清晰的。下面给出使用此类的一个具体实例:

/*
* Copyright (c) 2014 eryar All Rights Reserved.
*
* File : Main.cpp
* Author : eryar@163.com
* Date : 2014-09-11 20:46
* Version : 1.0v
*
* Description : Demo for Gauss-Legendre Integration usage.
*
* Key words : OpenCascade, Gauss-Legendre Integration
*/ #define WNT
#include <math_Function.hxx>
#include <math_GaussSingleIntegration.hxx> #pragma comment(lib, "TKernel.lib")
#pragma comment(lib, "TKMath.lib") class Test_GaussFunction : public math_Function
{
public:
virtual Standard_Boolean Value(const Standard_Real x, Standard_Real &y)
{
y = x; return Standard_True;
} private:
}; void TestGaussIntegration(void)
{
Test_GaussFunction aFunction;
math_GaussSingleIntegration aSolver(aFunction, , , ); std::cout << aSolver << std::endl;
} int main(int argc, char* argv[])
{
TestGaussIntegration(); return ;
}

主要是从math_Function派生一个类来在虚函数Value()中重定义求积函数即可。上述实例中计算的是如下积分:

计算结果如下图所示:

Figure 2.1 Gauss-Legendre Integtation Result

3. Application

由高等数学知识可知,积分的应用主要用于计算图形面积,体积及曲线的弧长,功等。

积分在OpenCASCADE中的主要应用有计算曲线长度,曲面面积及实体的体积等。如下图所示:

Figure 3.1 Compute Area of a Surface

示例代码如下所示:

TopoDS_Shape S = BRepBuilderAPI_MakeFace(BSS, Precision::Confusion()).Face();

GProp_GProps System;
BRepGProp::SurfaceProperties(S,System);
gp_Pnt G = System.CentreOfMass ();
Standard_Real Area = System.Mass();
gp_Mat I = System.MatrixOfInertia();

4. Conclusion

OpenCASCADE中实现的Gauss-Legendre求积算法,由于是查表求得Gauss点及求积系数,所以计算速度快。唯一不足是对高斯点数有限制。

综上所述,可知数值计算在OpenCASCADE中重要作用。一个TKMath库相当于实现了一本《数值分析》课本中的大部分内容。所以有兴趣的朋友可结合《数值分析》或《计算方法》之类的书籍,来对OpenCASCADE的数学库TKMath进行理论联系实际的深入理解。

5. References

1. Wolfram MathWorld, Numerical Integration,

http://mathworld.wolfram.com/NumericalIntegration.html

2. 易大义,沈云宝,李有法编. 计算方法. 浙江大学出版社. 2002

3. 易大义,陈道琦编. 数值分析引论. 浙江大学出版社. 1998

4. 李庆杨,王能超,易大义.数值分析.华中理工大学出版社. 1986

5. 同济大学数学教研室. 高等数学(第四版). 高等教育出版社. 1996

OpenCASCADE Gauss Integration的更多相关文章

  1. OpenCASCADE 参数曲线曲面面积

    OpenCASCADE 参数曲线曲面面积 eryar@163.com Abstract. 本文介绍了参数曲面的第一基本公式,并应用曲面的第一基本公式,结合OpenCASCADE中计算多重积分的类,对任 ...

  2. OpenCASCADE 参数曲面面积

    OpenCASCADE 参数曲面面积 eryar@163.com Abstract. 本文介绍了参数曲面的第一基本公式,并应用曲面的第一基本公式,结合OpenCASCADE中计算多重积分的类,对任意参 ...

  3. OPEN CASCADE Multiple Variable Function

    OPEN CASCADE Multiple Variable Function eryar@163.com Abstract. Multiple variable function with grad ...

  4. OpenCASCADE Curve Length Calculation

    OpenCASCADE Curve Length Calculation eryar@163.com Abstract. The natural parametric equations of a c ...

  5. OpenCASCADE Interpolation - Lagrange

    OpenCASCADE Interpolation - Lagrange eryar@163.com Abstract. Power basis polynomial is the most simp ...

  6. Jenkins in OpenCASCADE

    Jenkins in OpenCASCADE eryar@163.com Abstract. Jenkins是一个开源软件项目,是基于Java开发的一个持续集成工具,用于监控持续复制的工作,旨在提供一 ...

  7. OpenCASCADE 3 Planes Intersection

    OpenCASCADE 3 Planes Intersection eryar@163.com Abstract. OpenCASCADE provides the algorithm to sear ...

  8. 在 Laravel 中使用图片处理库 Integration/Image

    系统需求 PHP >= 5.3 Fileinfo Extension GD Library (>=2.0) … or … Imagick PHP extension (>=6.5.7 ...

  9. OpenCASCADE AIS Manipulator

    OpenCASCADE AIS Manipulator eryar@163.com Abstract. OpenCASCADE7.1.0 introduces new built-in interac ...

随机推荐

  1. js基础2

    什么是DOM? DOM:文档对象模型(Document Object Model),又称为文档树模型.是一套操作HTML文档的API. 什么是文档对象模型? DOM将html文档看成了一个对象,htm ...

  2. web存储

    1. cookie: 如果想将cookie取到,可以通过document.cookie;取到的是所有的cookie数据 他是一直保存在网页中的:他有一个时间的限制,如果时间过期,则删除 写入:docu ...

  3. ThinkPHP框架的一些基础应用

    这是俺滴师傅给俺传授了的知识,特在此分享. TP框架,做PHP开发的都应该有所耳闻.下面,我们就来说说入口文件的生成: 创建新项目时,首先,在目录文件下创建一个新的文件夹.然后将Thinkphp框架文 ...

  4. redis的面试题,没答出来,直接被pass

    redis 内存数据集大小上升到一定大小的时候,就会施行数据淘汰策略.redis 提供 6种数据淘汰策略:  volatile-lru:从已设置过期时间的数据集(server.db[i].expire ...

  5. java分享第十四天(TestNG Assert详解)

     TestNG Assert 详解org.testng.Assert 用来校验接口测试的结果,那么它提供哪些方法呢? 中心为Assert测试类,一级节点为方法例如assertEquals,二级结点为参 ...

  6. UML类图关系--继承(泛化)、实现、关联、聚合、组合、依赖

    在UML类图中,常见的有以下几种关系:  泛化(Generalization),  实现(Realization),关联(Association),聚合(Aggregation),组合(Composi ...

  7. JVM垃圾收集器

    JVM中垃圾的回收由垃圾收集器进行,随着JDK的不断升级,垃圾收集器也开发出了各种版本,垃圾收集器不断优化的动力,就是为了实现更短的停顿. 下面是7种不同的分代收集器,如果两个收集器之间有连线,则表示 ...

  8. java并发之volatile

    volatile是轻量级的synchronized,它在多处理器应用开发中保证了共享变量的“可见性”(可见性指当一个线程修改共享变量后,其它线程可以看到这个修改). volatile如果使用合理会比s ...

  9. C#中的委托与事件并存的理由

    更多资源:http://denghejun.github.io 问题 有了委托为什么还要有事件? 理论上,事件能完成的事情委托完全可以胜任,但是我们思考的这一方面是功能性,我们必须从他们各自的特点分析 ...

  10. adb devices出现no permissions

    在做Android开发的时候,有些设备插入USB调试线,无法调试,出现以下错误: joey 09:58 $ adb devices List of devices attached 20080411 ...