beidou111的博客

原创 【持续更新】总结国内外图形学物理模拟相关学者和网站

总结图形学物理模拟相关学者和个人网站仅做个人笔记所用。欢迎在评论区补充！国外学者Matthias MüllerNV物理组的leader，1999年博士毕业于苏黎世理工博士论文是研究聚合物分子动力学模拟的，后研究图形学。XPBD作者之一。早期曾经研究SPH模拟变形体研究领域：PBD、SPHhttps://matthias-research.github.io/pages/Theodore Kim耶鲁大学副教授，曾任皮克斯高级研究员.曾参与寻梦环游记等动画制作PS: 2020年有一篇S

原创 Perfect Sampling in SPH （完美采样）

介绍下背景：SPH中粒子半径和支撑半径是有关系的，一般取粒子半径的4倍为支撑半径。其依据就是下面所谓的完美采样。由于SPH是靠粒子携带数据，再通过粒子平均得到场数据的。但粒子有时多有时少，就会出现欠采样和过采样的问题。而过采样或欠采样，总要有个衡量的标杆吧？完美采样说的就是认为采样充足的时候，粒子该是什么样的一个密度。PS:背景中画出的网格（实际上SPH本质上不需要网格，但用网格加速邻域搜索的时候也可以用网格），可以看成是MPM等方法借助的背景网格的一个间距，取粒子直径大小。核半径是粒子半径的四倍。

原创 【Houdini】Houdini输出attribute到csv的方法

3.csv如图所示。输出的是将所有属性flatten之后的。现在我还没搞清楚怎么给他不同属性分列，暂时可以先多连几个节点，一个属性一个属性地输出到不同的csv里面。attribute to write改为By attribute pattern。只不过是普通的top节点里默认连了一个geometry import。

原创 渲染出表面的网格(Houdini VEX)

思路：用vex增加polyline到prim当中。注意stage中sopimport 不要勾选导入N和P。用name分组出line和ball。

原创 Shortest path 代码

核心就rtcPointQuery，调用的是intel embree的API。用碰撞检测后的结果查找最近点。(Narrow Phase)在这个for loop当中查找每个Vert的最近点。用BVH做碰撞检测。（Broad Phase)40行之前是初始化（mesh和param)Test.cpp中是入口main函数。Fork:(在Win10上）

原创 SIGGRAPH投稿相关官方指导

author instruction使用latex模板。

原创 MINNE编译

【代码】MINNE编译。

原创 Houdini python code:参数指定文件路径

【代码】Houdini python code:参数指定文件路径。

原创 阅读li2019-DOT源码--逐步调试

solve函数-> fullyImplicit()函数。到这里就是主要的算法了。

原创 WSL git文件异常 所有文件均显示已修改

解决方案：删掉之前在windows下git clone的文件夹， 然后在WSL中重新git clone。原因：是因为在Windows系统中git clone的文件夹，在WSL中会显示冲突。如图，文件中没有任何修改，但是都显示多了一个^M。

原创 在windows WSL 安装DOT li2019 SIGGRAPH

然后在65行的include_directories补充"SuiteSparse/CHOLMOD/Include"，也就是改为。WSL版本：ubuntu22.04。首先安装 intel mkl。

原创 记录CUDA 编译器的很奇怪的一个BUG nvcc fatal : Could not set up the environment for Microsoft Visual Studio

这是因为在配置其他代码的时候，弄坏了环境变量。只要把当前多余的path变量情况就好了。

原创 【读文献】xian2019: A Scalable Galerkin multigird method

假如用多个粗网格点可以通过权重控制一个点，那么这个weight就是0到1的，是smoothed。假如没有任何的多重网格，也就是只有一层的话，那么U就是包含的material space(rest pos，X)信息。假如有多重网格，那么就用粗网格的点去控制某个细网格的点。其中i是被控制的细网格点，j是粗网格上的控制点。而这个weight就是用j控制i的权重，要么是0（不控制），要么是1（控制）。xi是3x1的位置坐标，Aj是3x4的矩阵，Xi是4x1的齐次坐标的rest pos。那么从上面的关系，改写成。

原创 houdini可视化distance约束的VEX代码

其中0号input连接sim geo 1号连cons geo。

原创 connect near points houdini

【代码】connect near points houdini。

原创 houdini肌肉刷pin点的方法

在视窗enter进入编辑模式，就出现笔刷，在两端刷上gluetoanimation属性（属性可以在节点的attribute选项卡更改）。笔刷用法见出现的说明（左键刷出来，control+中键抹除，滚轮放大笔刷大小等）。使用blast节点，在group输入@muscle_id=_mus_deltoidclavicular_L就能单独过滤出这块肌肉。按s，鼠标悬浮到某一个四面体，就能在右上角高亮其对应的group。使用attribute paint节点（或者muscle paint节点，两者基本是一样的）

原创 Houdini展示muscle id

原创 弹簧质点系统求Hessian

https://www.matrixcalculus.org/

原创 安装amgcl

你要把设置参数之类的放在main函数里面。每个case都是个单独的main函数。BOOST和MPI不需要安装也行。安装BOOST之后可以运行时指定参数。这是个header-only的库，不需要安装。编译poisson3Db example。下载数据（matrix maket格式）解压到amgcl根目录。

原创 用一维数组存储不等长的二维数组的方法

start存储每行起始的下标。start[0]必定是0。start[-1]必定是n(一维数组长度)。在numpy或者gpu中，通常会遇到不等长的二维数组，如果都用其最大长度存储，就会浪费许多内存。原本我希望取val2d[1,2] 即第1行第2个值（=e）的时候，等价于在val1d中求。我们只需要把2d数组flatten成一维数组，然后用辅助数组start即可存储。2维数组记为val2d，1维数组记为val1d。第1行：[c, d, e]即val1d[4] = e。

原创 安装pyamgx

将pyamgx/demo.py的前几部分改为下面的代码，并且运行。然后在pyamgx下。

原创 在当前终端加载“x64 Native Tools Command Prompt for VS 2022“

加入到系统PATH中。

原创 使用vscode debug cpp/python混合编程的程序（从python调用的C++编译的dll）

拷贝自 https://marketplace.visualstudio.com/items?采用了attach process的方式，大致是根据dll的process id将debugger挂载上去。

原创 bibtex引用网页的方式和wiki网页的方式

online也可以换成misc。

原创 Houdini 通过wedge来做模拟参数对比 (PDG TOP)

更改参数：attribute name处填入刚才那个变量的名字（如stretchstiffness）, 其他的就表示要从-10变动到10，然后共10个点。在打开的窗口中输入P@streachstiffnessexp （只要是@后面随便起个名字就行， P@表示PDG能够找到的变量）这是个简单的布料悬挂的例子。指定输出的文件路径（这里P@wedgeindex表示是第几个小绿点，例如0到9）我们现在的目的是想要对比不同的streach stiffness对模拟的影响。这个节点是用来输出结果到文件的。

原创 进入python时自动运行代码

方法2是设置PYTHONSTARTUP环境变量为一个启动脚本的路径。

原创 进入ipython自动引入numpy

首先找到ipython_config.py文件。更多请看ipython_config的用法请看。

原创 Houdini pbd_constraints.h的文件位置

Houdini安装目录下的houdini\vex\include文件夹。

原创 推导Hessian of XPBD

那么如何求C的二阶导数呢。

原创 推导 Primal Residual XPBD

原本的XPBD(Macklin 2016)Chen2023主张左边的红色二阶项可以舍掉，但右边的g（也就是Primal residual) 应该恢复。因为左端那个是Quasi Newton的常见做法。而舍弃右边需要保证lambda很小才成立。舍掉会导致XPBD残差下降停滞，仿真不稳定。

原创 pyamg Smoothed Aggregation算法流程

设置如下算例。分两步，第一步是setup，第二步是solve。第一步的目的就是得到P等。

原创 Newton法求解布料

【GAMES103-基于物理的计算机动画入门】 【精准空降到 1:25:38】 https://www.bilibili.com/video/BV12Q4y1S73g/?

原创 python局域网共享文件 http server

【代码】python局域网共享文件 http server。

原创 Debug Eigen

（下载地址： https://gitlab.com/libeigen/eigen/-/blob/master/debug/msvc/eigen.natvis ）（如果你git clone过Eigen文件夹，就发现它就在Eigen源码的Debug目录下面）假如我们直接想要检查Eigen中的VectorXf的值，我们会发现无法看到内部的值，如图。（其中%VSINSTALLDIR%是Visual studio的安装目录）只需要下载eigen.natvis这个文件，然后放到下面的目录即可。

原创 zeno编译

/step 0//step 1//step 2。

原创 WSL无法保存文件（权限不足）

其中username是你的用户名。

原创 git proxy

【代码】git proxy。

原创 Houdini vex清除过远的离散点（outliers）

filter num控制最小邻居数目。filter dist控制最小距离。

原创 Houdini中获取matrix的属性需要getcomp函数

matrix3x3的type是3@xxx。直接写mat[0][0]是错误的。

原创 houdini vex中的属性类型（attribute）

https://www.sidefx.com/docs/houdini/vex/snippets.html#parameters