OpenCASCADE 平面与球面求交 eryar@163.com OpenCASCADE提供了类IntAna_QuadQuadGeo用来计算两个二次曲面quadric(球面.圆柱面.圆锥面及平面,平面是二次曲面的特例)之间的交线.他们之间可能的结果有: l 一个点 l 一条或两条直线 l 一个点和一条直线 l 圆 l 椭圆 l 抛物线 l 双曲线 将源码结合<高等数学>.<解析几何>等书,可以来学习如何将理论付诸实践.本文主要介绍这个类中平面与球面求交的源码实现.  将源码列出如…

OpenCASCADE 平面求交 eryar@163.com OpenCASCADE提供了类IntAna_QuadQuadGeo用来计算两个二次曲面quadric(球面.圆柱面.圆锥面及平面,平面是二次曲面的特例)之间的交线.他们之间可能的结果有: l 一个点 l 一条或两条直线 l 一个点和一条直线 l 圆 l 椭圆 l 抛物线 l 双曲线 将源码结合<高等数学>.<解析几何>等书,可以来学习如何将理论付诸实践.本文主要介绍这个类中两个平面求交的源码实现.从源码中也可以看出Ope…

OpenCASCADE直线与平面求交 在<解析几何>相关的书中都给出了直线和平面的一般方程和参数方程.其中直线的一般方程有点向式形式的. 由于过空间一点可作且只能作一条直线平行于已知直线,所以当直线上一点(x0, y0, z0)和它的一方向向量(m,n,p)为已知时,直线就完全确定了.所以在OpenCASCADE中直线类gp_Lin有一个构造函数: gp_Lin (const gp_Pnt &P, const gp_Dir &V) 即通过点和方向来构造直线.由直线的点向式方程容…

OpenCASCADE圆与平面求交 eryar@163.com 在 解析几何求交之圆与二次曲面中分析了OpenCASCADE提供的类IntAna_IntConicQuad可以用来计算圆与二次曲面之间的交点,这个算法是将平面Plane作为二次曲面的一个特例来处理,最后主要是对三角函数方程进行求解. 当直接使用圆和平面作为参数时,IntAna_IntConicQuad重载了函数Perform来对圆和平面进行求交计算,这时的算法与前面解三角函数不同,代码如下: void IntAna_IntConic…

光线求交 光线定义:position \(a(t)\) = \(o\) + \(t\vec{d}\); 球定义: center p, radius r; 平面定义:normal \(\vec{n}\) , offset t; 三角形定义:position \(a_1\), \(a_2\), \(a_3\), normal \(\vec{n}\); 光线与球相交 (Ray/Sphere Intersection) c++代码 : bool HitTest(const Ray& ray, HitTe…

不管是哪种全局光照算法,最根本的都要落实到光线与物体的求交.主要分为光线与参数曲面和非参数曲面的求交,典型的参数曲面有球.盒.圆柱等基本体及基本体的组合体,以及一些更为复杂的参数曲面.非参数曲面就是所谓的网格(Mesh),通过一族顶点集和三角面的集合来近似表示物体曲面.一个简单的渲染器至少应包含如下几种求交算法:光线与任意三角形的求交,光线与平面求交,光线与盒相交(更特殊的是光线与轴对齐盒的求交),光线与球体相交等.这里先分享一下常见的光线与三角形求交的算法: —————————————————…

引擎中,ray与quad求交,算法未细看,但有求解二次方程,不解.ray与triangle求交,使用的是97年经典算法,仔细看过论文,多谢小武同学指点,用到了克拉默法则求解线性方程组.想模仿该方法,做ray与quad的求交,发现方程里不仅有u和v,还有uv,没法变换成线性方程组的形式.本以为引擎中quad中四个点可以不共面,看过接口,不然,“不共面和退化的多边形不保证正确结果“.而后又有两个问题,一是,与一个quad求交比与两个三角形求交快吗?二是,如果前面的问题答案为否,即两个三角形更快,为何…

任意多边形求交: private void btnPolygon_Click(object sender, EventArgs e) { , , , , , , , , , , , , , }; , , , , , , , }; ; var coords = new Coordinate[m]; var coordscheck = new Coordinate[m]; ; i < Xs.Length / ; i++) { ]; + ]; coords[i] = new Coordinate(x,…

pro:给定一枚蛋糕,蛋糕上某个位置有个草莓,寿星在上面切了N刀,最后寿星会吃含有草莓的那一块蛋糕,问他的蛋糕占总蛋糕的面积比. sol:显然需要半平面交求含有蛋糕的那一块,然后有圆弧,不太方便求交. 所以我们可以直线构成的边界,求出平面交: 然后用这个多边形去和圆求交. (百度了一下很多人都没过,好像是这题很卡精度,反正我每个地方都改过,还是WA,大概wa了4个小时了,要不以后再回来改. 当然也不排除有其他问题. #include<bits/stdc++.h> #define rep(i,a…

hdu 5111 树上求交(树链剖分 + 主席树) 题意: 给出两棵树,大小分别为\(n1\),\(n2\), 树上的结点权值为\(weight_i\) 同一棵树上的结点权值各不相同,不同树上的结点权值可以出现重复 每次查询 \(u1\) \(v1\) \(u2\) \(v2\) 第一棵树上\(u1\) 到 \(v1\)的路径上所有结点权值组成的集合\(S1\) 第二棵树上\(u2\) 到 \(v2\)的路径上所有结点权值组成的集合\(S2\) 求\(S1\) 与 \(S2\) 的交集 \(1…

题意:https://ac.nowcoder.com/acm/contest/1109/C 问你有几个x满足A,B集合都能XOR出x. 思路: 就是线性基求交后,有几个基就是2^几次方. #define IOS ios_base::sync_with_stdio(0); cin.tie(0); #include <cstdio>//sprintf islower isupper #include <cstdlib>//malloc exit strcat itoa system(&…

题意: 给你n个基,q个询问,每个询问问你能不能 l~r 的所有基都能表示 x . 思路: 建一颗线性基的线段树,up就是求交的过程,按照线段树区间查询的方法进行check就可以了. #define IOS ios_base::sync_with_stdio(0); cin.tie(0); #include <cstdio>//sprintf islower isupper #include <cstdlib>//malloc exit strcat itoa system(&qu…

[练习3.4] 给定两个已排序的表L1和L2,只使用基本的表操作编写计算L1∩L2的过程. [练习3.5] 给定两个已排序的表L1和L2,只使用基本的表操作编写计算L1∪L2的过程. 思路比较简单,测试代码如下,两道题比较相似,测试代码就放一起了. #include <iostream> #include "linklist.h" using namespace std; using namespace linklist; template class List<in…

Github链接传送:JasonCeng/MultipartyPSI-Pro 大家好,我是阿创,这是我的第29篇原创文章. 今天是一篇纯技术性文章,希望对工程狮们有所帮助. 向大家推荐一个我最近改造的算法库(JasonCeng/MultipartyPSI-Pro) 主要涉及隐私计算.联邦学习领域,是一个支持多方的隐私求交算法库,主要是是基于osu-crypto/MultipartyPSI****的增强实现. 详细介绍就是下方的README啦~ 基于可编程不经意伪随机数的多方隐私求交算法库 Pro…

Intersection Time Limit: 4000/4000 MS (Java/Others) Memory Limit: 512000/512000 K (Java/Others) Total Submission(s): 602 Accepted Submission(s): 247 Problem Description Matt is a big fan of logo design. Recently he falls in love with logo made up by…

Description Given one triangle and one circle in the plane. Your task is to calculate the common area of these two figures. Input The input will contain several test cases. Each line of input describes a test case. Each test case consists of nine flo…

教学班级:周三上午三四节 项目地址:https://github.com/875571216/- PSP表格 psp2.1 Personal Software Process Stages 预估耗时(分钟) 实际耗时(分钟) Planning 计划 60 30 Estimate 估计这个任务需要多少时间 10 10 Development 开发 Analysis 需求分析 (包括学习新技术) 180 240 Design Spec 生成设计文档 30 10 Design Review 设计复审…

Problem Intersecting Lines (POJ 1269) 题目大意 给定两条直线,问两条直线是否重合,是否平行,或求出交点. 解题分析 主要用叉积做,可以避免斜率被0除的情况. 求交点P0: 已知P1 P2 P3 P4 运用 P0P1 X P0P2 = 0 和 P0P3 X P0P4 = 0 C++ 用%.2lf g++ 用 %.2f!!! C++ 用%.2lf g++ 用 %.2f!!! C++ 用%.2lf g++ 用 %.2f!!! 参考程序 #include <cstd…

http://acm.hdu.edu.cn/showproblem.php?pid=1086 You can Solve a Geometry Problem too Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 8861    Accepted Submission(s): 4317 Problem Description Many…

Atlantis 题目连接 http://poj.org/problem?id=1151 Description here are several ancient Greek texts that contain descriptions of the fabled island Atlantis. Some of these texts even include maps of parts of the island. But unfortunately, these maps describ…

题意: 传送门 给\(n\)个集合,每个集合有一些数.给出\(m\)个询问,再给出\(l\)和\(r\)和一个数\(v\),问你任意的\(i \in[l,r]\)的集合,能不能找出子集异或为\(v\).简单点说,\(v\)能用\([l,r]\)任意一个集合的子集异或和表示. 思路: 子集异或和显然是用线性基.我们用线段树维护任意区间的线性基交集即可. 代码: /** 求交集 O(logn * logn) **/ LBasis intersection(const LBasis &a, const…

B样条曲线生成 Bezier曲线缺点:改变任一控制点的位置,将影响整条曲线的形状. B样条曲线是对Bezier曲线的改进,可进行局部控制,生成的曲线与控制多边形的外形更接近,将Bezier曲线作为一特例. 三次B样条曲线的矩阵表达式 \[\begin{aligned}P(t) &= \frac 1 6 \begin{bmatrix} t^3 & t^2 &t & 1\end{bmatrix}\begin{bmatrix}-1 & 3 & -3 & 1…

Bezier曲线生成 法国工程师Pierre Bezier在雷诺公司使用该方法来设计汽车.一条Bezier曲线可以拟合任何数目的控制点. 公式 设\(n+1\)个控制点\(P_0,P_1--P_n\),其中$P_k=(X_k,Y_k,Z_k),0≤k≤n $ 则\(n\)次Bezier曲线为: \[P(t)=∑P_iB_{i,n}(t)\qquad 0≤t≤1 \] 其中,\(B_{i,n}(t)\)是Bernstein基函数,即 \[B_{i,n}(t)=c_n^it^i(1-t)^{n-i}…

Matlab geom3d函数解析 geom3d函数库 geom3d库的目的是处理和可视化三维几何原语,如点.线.平面.多面体等.它提供了操作三维几何原语的底层功能,使得开发更复杂的几何算法变得更加容易. | 序号 | 函数名 | 作用 | 输入 | 输出 | | --- | --- | --- | --- | --- | | 1 | area = triangleArea3d(pt1, pt2, pt3) | 获取三角形面积 | 三个点坐标 [10 10 10], [30 10 10], [1…

题面 题意:给出小于10个点形成的凸多边形 和一个半径为r 可以移动的圆 求圆心在何处的面积交最大,面积为多少 题解:三分套三分求出圆心位置,再用圆与多边形面积求交 #include<bits/stdc++.h> #define inf 1000000000000 #define M 100009 #define eps 1e-12 #define PI acos(-1.0) using namespace std; struct node { double x,y; node(){} nod…

OpenCASCADE BRep vs. OpenNURBS BRep eryar@163.com Abstract. BRep short for Boundary Representation. First give the definition of the BRep, then compare the BRep mode between OpenCASCADE and OpenNURBS. There are 3 main representation method: use face/…

OpenCASCADE 基础 转载▼ 一直在用OCC作项目,但这方面的中文资料很少,看来OCC在中国还不是十分普及: 后来,项目中使用OCC和DirectX结合使用,取得了很好的效果: 随着OCC6.3版本的推出,Open CASCADE在速度方面已有了很大的改变.以下为一些OCC的基础知识,愿与各位OCC爱好者共同学习: 一:OCC中的基础类: gp_Pnt 在OCC中,gp_Pnt表示一个顶点,gp_Vec表示一个向量,可以用两个顶点来生成一个向量. 比如: gp_Pnt P1(0,0,0)…

Two analytical 2d line intersection in OpenCASCADE eryar@163.com Abstract. OpenCASCADE geometric tools provide algorithms to calculate the intersection of two 2d curves, surfaces, or a 3d curve and a surface. Those are the basis of the Boolean Operat…

OpenCASCADE Make Face With Holes eryar@163.com OpenCASCADE提供了构造Face的类BRepBuilderAPI_MakeFace,使用这个类可以构造出带孔的面.如下图所示: 当然,要得到上图所示的结果,还可以使用Boolean操作,用一个面去Cut几个圆柱.当使用布尔操作就会涉及到一些复杂算法,如求交,重构Topo体等,比较耗时.既然可以直接在生成面的时候挖孔,这个不涉及复杂算法,速度.稳定性都比使用布尔操作要好.本文主要来介绍如何使用BR…

OpenCASCADE 麻花钻头造型实例分析 eryar@163.com Abstract. OpenCASCADE provides a simple twist drill bit modeling example in Tcl script. The blog will give a details of some key points, when you understand the key points, you can modeling the shape like that. Ke…