- 博客(109)
- 收藏
- 关注
RSS订阅原创 高斯过程处理大型数据集时返回值为 0 的问题
但是当我对一个 100 维的问题，用 1000 个训练样本训练高斯过程模型的时候，得到的预测值几乎接近于 0。我不知道这是我代码的问题还是模型本身的问题，于是我去用了 sklearn 的库。这是我在网上搜的高斯过程的代码，挺容易理解的。
2025-02-04 18:33:51
211
原创 Using HPC for Computational Fluid Dynamics 笔记
但是按列分配有他固有的问题，就是对于某一个列处理的话，那么就只能按元素计算，或者是按向量计算，但是没有矩阵，没法利用 BLAS 的第二级和第三级，也就是矩阵和向量、矩阵和矩阵。他就是说的取某个对角线元素，要更新他右下角的块矩阵，然后接下来要取下一个对角线元素，就是刚刚更新的块矩阵的左上角，所以会有这样一种嵌套的感觉。指的是专门拿一个处理器去通信，同时其他处理器在计算。这个 A(i,k) * A(k,j) 还真的是某种意义上的多余的值，只要每次迭代的时候都减去它，那么行就是 a，列就可以直接除。
2025-02-04 18:33:20
838
原创 VDB：具有动态拓扑的高分辨率稀疏体积 VDB: High-Resolution Sparse Volumes with Dynamic Topology 论文阅读笔记
我们的“VDB”之所以如此命名，是因为它是一个体积 Volumetric 动态 Dynamic 网格，与 B+ 树共享多个特征，利用时变数据的空间一致性来单独、紧凑地编码数据值和网格拓扑。VDB 模拟了一个几乎无限的 3D 索引空间，允许对高分辨率的稀疏卷进行高速缓存一致和快速的数据访问。它对体数据的稀疏性没有拓扑限制，并且在插入、检索或删除数据时支持快速（平均 O(1)）随机访问模式。这与大多数现有的稀疏体积数据结构形成鲜明对比，后者采用静态或流形拓扑，并需要特定的数据访问模式来补偿缓慢的随机访问。
2025-02-04 18:32:58
835
1
原创 分层 RLE 水平集：紧凑且多功能的可变形表面表示 Hierarchical RLE Level Set: A Compact and Versatile ... 论文阅读笔记
结合了 DT-Grid 和 RLE 水平集H-RLE 水平集以维度递归方式采用 RLE。RLE 方案允许连续非窄带区域的紧凑存储，同时沿每个轴的维度递归编码有效地压缩非窄带平面和体积。本文除了介绍 H-RLE 水平集数据结构及其高效的核心算法之外，还描述了从我们使用该结构中受益的众多应用：我们的统一隐式对象表示、高效且鲁棒的网格到水平集转换、快速射线跟踪、水平集变形、碰撞检测和完全稀疏流体模拟（包括 RLE 向量和矩阵表示）。
2025-02-04 18:32:35
771
1
原创 共轭传热和浸没边界耦合相关的论文的阅读笔记 2
在这篇文献中，浸没边界法（IBM）与共轭传热的耦合通过多材料问题的处理实现。文中使用的IBM基于浸没不连续网格的方法，界面由不同材料区域的相邻笛卡尔面组成。通过将这些面投影到真实材料界面，并计算速度重构系数，使用二阶插值来满足材料界面的物理条件。对于共轭传热问题，在流体和固体界面处，IBM方法要求温度场的连续性（Dirichlet边界条件）和热流的平衡（Neumann边界条件）。
2025-02-04 18:31:56
580
原创 UE UMG 多级弹出菜单踩坑
这个 Get widget > As Canvas Slot - set Size / Position 的方法确实可以。所以这里也是涉及到了 slot，只是有一个函数使得 canvas 作为 slot？但是一般情况下我怎么获取到一个 widget 然后把他视为 slot 变量？对比不同方法，BP 中的设置保持相同，比如我的 anchor 始终没有变过。只有让新添加的 widget 是上一级菜单的 child 才能对齐。随便找的 slot 变量也不是这个类型，怪。其中我的子菜单的 UI 是这样的。
2024-12-18 10:06:08
537
原创 UE5.4 Release 源码构建的坑
重新构建之前记得 Clean，因为这是构建配置方面的错误，构建设置写在 makefile 里面，如果不清除构建的话，那么就会沿用上一次的构建设置。于是安装 MSVC v143 - VS 2022 C++ x64/x86 build tools (v14.38-17.8)于是发现我这里都有提示。
2024-10-18 11:36:31
661
原创 浸没边界 直接强迫法 圆球绕流验证 阅读笔记
我们还引入了一种改进的体积力分布函数，它将IB点离散体积中的体积力完全传递到周围的笛卡尔网格。所得到的求解器的总体精度是二阶的，这在涉及流经球体的层流的两个测试案例中得到了证明。在每个时刻，通过沿着物体表面的局部法线重建解，在位于物体外部但紧邻物体的笛卡尔网格节点处应用边界条件，来考虑物体对流动的影响。注意到振荡是由直接强迫点处数值描述的瞬时变化引起的，我们提出了一种允许这些点处数值描述平滑过渡的公式。使用 VPM 方法。研究圆球绕流对平面边界的影响的，研究的是速度分布、边界层厚度，湍流强度等。
2024-10-06 23:16:24
1361
1
原创 Git unable to resolve reference ‘refs/remotes/origin/HEAD‘: reference broken 的解决方法
清除本地仓库中不存在的引用，也没有效果。设置远程仓库的默认分支到本地仓库。于是尝试一些解决方法。
2024-09-22 23:16:55
440
原创 制造解法 Manufactured Solutions 相关的论文的阅读笔记
对于本文提出的纳维-斯托克斯模拟，绝对粘度被选择为一个大的常数值（μ = 10 N · s/m^2），以获得对流和扩散之间的平衡。导热系数 k 由粘度通过指定常数确定。他这里是直接把制造解代入偏微分方程式，并且把所有项移到一边，得到的就是源项。解结构也必须足够复杂，以便测试所测试的控制方程中的所有项都得到执行。他的写法似乎也是，直接写了 NS 方程全部移到一边作为源项。然后 sci hub 下载的都是公式不全的。但是他没有讲制造解的源项具体是怎么构造的。用制造解法还可以找出编程错误。
2024-09-22 22:20:10
598
原创 共轭传热和浸没边界耦合相关的论文的阅读笔记
我有点不知道这是哪里体现了曲线网格了可能是体现在，他在开篇说，它使用的 IBM 方案是曲线坐标系的但是我是不知道他之后提出的那些算法有什么是特定需要曲线坐标系才能完成的。
2024-09-21 19:21:47
1068
原创 Efficient monolithic immersed boundary projection method for incompressible... 阅读笔记
lambda* 通过 A1 lambda* = R1 解出来，A1 是 A 的左上角，R1 是 R 的前四项，那 lambda* 似乎也应该具有类似的物理意义。他这个的构造的第一行是能量方程，第二行是动量方程，第三行是能量外力的插值方程，第四行是动能外力的插值方程，都是 IBM 的经典核函数插值。我大概理解 bc 是 boundary condition 的缩写，m 是动量的缩写，e 是能量的缩写。第三行是能量外力的插值方程，第四行是动能外力的插值方程，都是 IBM 的经典核函数插值。
2024-09-19 07:27:49
790
原创 使用 MPI 做 3D 带状矩阵的转置
一开始我大概有一个想法，假设两个进程的话，就把整个大矩阵分成 4 * 4 的块，进行分配，但是我不知道怎么分，就算是把转置之前的数据分布和转置之后的数据分布写出来了，也似乎找不到规律。后面通过 chatgpt 做出来了。每个进程的输入：大矩阵的。每个进程的输出：大矩阵的。
2024-05-25 14:55:34
297
原创 SolidWorks 创建零件时，阵列出来的实体是新实体，需要和已有实体合并
解决方法是插入一个组合特征，把所有实体都选中，组合起来，就不会有这个问题了。通过阵列实体创建出来零件之后，发现零件里面有些实体没有合并到一起。视角进入到内部来看也可以看到它们不是一个东西。
2024-04-02 10:37:51
1574
原创 vulkan 渲染模型加载找错，stride 不对
回去看一下顶点数据，就是 6 个 6 个 float 的读进顶点着色器的。可以看到，顶点之间都是紧密挨着的，就是正常的三角形。加载下来的顶点数组，我跟教程是一样的。但是输入到顶点着色器的数据就不一样了。一看就是 stride 不对。
2024-03-17 01:09:00
477
原创 matlab 使用 contourf 或 scatter 绘制动图输出成视频，避免帧与帧之间重合
一般的读取场的每个时刻的状态，然后输出成视频的 matlab 代码类似这样。后来我想用 scatter 画的时候，发现帧与帧之间的图像重合了。最后发现是我在循环里面创建 axes 的缘故。
2024-03-13 00:43:50
713
原创 NUMERICAL SIMULATION OF THE PROPELLER/WING INTERACTIONS FOR TRANSPORT AIRCRAFT 笔记
运输机螺旋桨/机翼相互作用巡航条件（跨音速飞行速度，中等推力）起飞情况（低速和高推力）对这两种情况进行数值实验致动盘模型预测螺旋桨滑流与不同运输机机翼配置的相互作用的能力已经针对必须考虑这种相互作用的两个主要重要飞行条件进行了评估：起飞和巡航配置。在第一种情况下，欧拉流计算显示了与高升力系统的复杂流动相互作用，这导致升力能力显着增加。为了更现实地优化螺旋桨位置，还应计算机舱并考虑粘性效应。对于第二种情况，模型已准确预测了机身设计所用的主要气动参数（整体升力系数、翼展载荷变化和压力分布）。
2024-02-29 21:23:16
1170
原创 用于模拟颗粒流的直接强迫浸没边界法 An immersed boundary method with direct forcing 笔记
使用固定且均匀的计算网格来计算悬浮刚性颗粒周围的不可压缩粘性流主要思想是将 Peskin 的正则化 delta 函数方法 [Acta Numerica 11 (2002) 1] 纳入流固相互作用力的直接公式中，以便允许欧拉和拉格朗日表示之间的平滑转换，同时避免时间步长的严格限制。该技术是在有限差分和分数步环境中实现的。我们提出了一种改进的浸入边界方法，直接公式化了流固相互作用力。Peskin [23] 的正则化 delta 函数用于任意拉格朗日位置和离散欧拉位置之间的关联。
2024-02-26 11:26:10
1718
原创 动态管状网格：高分辨率水平集的有效数据结构和算法 Dynamic tubular grid: An efficient data structure... 论文阅读笔记
XN−1X_{N-1}XN−1表示对于 N 维情况，把网格点投影到前 N - 1 维张成的空间里得到的点集称为 p-column (projection column 的简称) number所以对于 2 维情况，p-column 就是一维的N 维的连通分量 connected component 被定义为 p-column 内相邻网格点的最大集合例如，在 2D 中，p-column x 被定义为管状网格中通过正交投影到 X 轴而投影到 (x, 0) 的网格点集。
2024-02-20 01:29:12
1178
原创 使用动态网格的流体动画 Fluid Animation with Dynamic Meshes 论文阅读笔记
原文：Klingner, Bryan M., et al. “Fluid animation with dynamic meshes.” ACM SIGGRAPH 2006 Papers. 2006. 820-825.使用 [Alliez et al., 2005] 的方法动态生成不规则的四面体网格根据边界的位置、边界的形状、基于流体和速度场的视觉重点部分的标准来构建一个尺寸场。这个尺寸场表明要生成的四面体网格在某点处的尺寸使用不规则的网格，而不是轴对称的，因为不规则的网格更加适应弯曲边界和不规则边界并且不
2024-02-19 00:50:04
1314
原创 使用八叉树模拟水和烟雾 Simulating Water and Smoke with an Octree Data Structure 论文阅读笔记
后面的都无关紧要了总之他这个文章的核心就是，现在八叉树的两个网格就只考虑彼此，暴力忽略其他网格，就这么得到了对称的模板忽然发现称为模板隐含着一个很优雅的事情，就是它可以包括某个点是流体或者是固体的情况这样，我们只需要最后讨论一下边界条件怎么放进来就好了，一开始推公式的时候就不用想边界好酷啊总之好羡慕啊，我现在的心情就跟我第一次看到半拉格朗日平流的 stable fluid 一样，感觉这么简单但是有效的东西，他们就能够研究到，然后发出来，就很优雅，很有应用上的美感的这么一件事。
2024-02-18 01:17:25
1249
原创 Windows 编译 yangfengzzz/fluid-engine-OpenVDB
我想将 OpenVDB 接入 doyubkim 的流体引擎然后搜到已经有人做过这件事了。
2024-02-17 10:04:50
1193
原创 OpenVDB Cookbook 笔记
如果不是对应同一个世界坐标，那就是两个网格的变换不同，一般的方法是先将一个网格重新采样到另一个网格的索引空间。可以获得元数据的 string 形式的类型名称，进而得知某个元数据的数据类型，进而进行类型转换。两层循环，外层迭代树结构的所有叶子节点，内层对于每一个叶子节点，根据属性的名字构建。遍历到的元素的 SDF 值如果在这个区间内，也就是在窄带内，那么就赋值，否则跳过。为了单线程调试 TBB，可以把输入的范围做成没有拆分的，对于。，这个迭代器的范围是由 tbb 划分的，原始迭代范围的子区间。
2024-02-15 15:51:27
1442
1
原创 并行计算导论 笔记 1
比如对于一个 if else，所有处理器都用相同的这个代码处理不同的数据，假设一半的数据能满足 if，另外一半只能满足 else，那么在执行这段代码的时候，在所有处理器处理 if 的时候，有一半的处理器忙碌而另一半空闲，else 也是一样。书中说的”虽然没有对共享变量进行更新，但是系统并不进行检测“应该说的是，假设讨论的处理器有 A、B、C，那么 A 对自己的共享变量进行更新的时候，A 更新的并不是 B、C 的共享变量，所以无法设计一个程序来检测假共享？总之这话听起来挺怪的。
2024-02-11 15:07:21
1345
原创 GDC2023 Photon Water System Reading Note
【代码】GDC2023 Photon Water System Reading Note。
2024-02-01 10:16:36
1219
原创 Fluid Engine Development 创建 PIC, FLIP, Level Set 演示 Demo
然后他这个 Marching Cube 的算法里面写了对于边界的处理，原仓库的代码里面特意设置了不创建底部的面的标志，我不知道为什么要这么做。一开始我还不知道有这个标记，于是等那个溃坝的模拟，一帧五分钟，等了五十帧，之后才开始做，这个时候才开始发现模型没有底面……用 Blender 渲染的话，为了做光线追踪，液体网格应该封闭才能达成正确的光线反射，现在这个不创建底面直接就没有那个正确的效果了。提供了自己创建求解器的模板，但是现在暂时应该用不到，现在想输出他的求解器的图像。项目的 CMakeList。
2024-01-30 14:58:16
454
原创 Blender Object 的 name 字段不是 string
之前在获取 object 的 name 之后直接当作 string，对他调用 string 的函数，结果导致了 blender 崩溃，完全不知道哪错了。后面发现 object 的 name 字段似乎不是 string，需要先转到 string。感觉 blender 这点不好，调用了错误的函数之后就直接崩溃了……Blender 版本 4.0.2。
2024-01-29 22:41:16
585
原创 Fluid Engine Development PIC/FLIP 代码分析
把 Fluid Engine Development 看完了，但是仍然感觉不懂感觉还是应该了解整体代码怎么写的，所以做个总结看着看着，感觉还是从底层开始看起从底层重新开始看的时候，感觉就来了而且作者也有很多注释，感觉能够体会到别人的思路他这里也有很多内容，我选择从 PIC/FLIP 开始看起然后我选择它的 hybrid_liquid_sim 工程它的 main 文件就很容易找到了，在 src\examples\hybrid_liquid_sim\main.cpp。
2024-01-24 10:09:17
1011
原创 Fluid Engine Development Windows 编译出现 xmemory ‘_invalid_parameter‘: is not a member of ‘`global 解决方法
Windows 编译《Fluid Engine Development》配套的代码时出现了 ‘_invalid_parameter’: is not a member of ‘`global namespace’’ 类似的错误。在代码中找到 pybind11 的 common.h 去更改如下部分的代码。
2024-01-23 14:28:10
460
原创 计算机图形学流体模拟 blender 渲染脚本
做流体模拟的时候，想要复现别人的成果，但是别人的代码都是每帧输出 ply 格式的文件，渲染部分需要自己完成看了一下，似乎用 blender 是最简单的，于是记录一下过程中用到的代码Blender 版本 4.0。
2024-01-09 17:00:37
1224
原创 最优化理论期末复习笔记 Part 1
数学基础线性代数从行的角度从列的角度行列式的几何解释向量范数和矩阵范数向量范数矩阵范数的更强的性质的意义几种向量范数诱导的矩阵范数1 范数诱导的矩阵范数无穷范数诱导的矩阵范数2 范数诱导的矩阵范数各种范数之间的等价性向量与矩阵序列的收敛性函数的可微性与展开一维优化问题牛顿莱布尼茨公式对多维的拓展Lipschitz 连续中值定理凸优化问题凸函数的判断f 在 D 一阶可微正定矩阵f 在 D 二阶可微无约束问题的最优性条件
2024-01-03 10:21:31
764
2
原创 VS Code ssh-remote 连接时由于 clash 导致 XHR failed 报错的一个解决方法
所以出现 XHR failed 的错误之后，上网一搜，发现是网络相关的问题，于是我就把 clash for windows 的 proxy 关掉，还是有报错，于是我就把 Profiles 里面的 yaml 切换回默认的 config.yaml 就好了。我个人不太懂 clash for windows 怎么运行的。这只是我个人情况的一个解决方法，可能别人的错误原因是不一样的。但是我知道我所有跟网络代理相关的操作只有用它。
2023-12-28 11:02:39
957
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人