十年磨一剑

原创 Scale-Dependent 3D Geometric Features

这篇文章主要讲如何提取3D模型中的corner和edge. 但是是scale-dependent的我理解的scale就是文中平滑的时候影响的范围大小

原创 Non-Rigid Registration Under Isometric Deformations

2. OverviewCorrespondence Computation.这部分讲matching过程不仅考虑最近点, 而且还要考虑feature最相似的点3. Downsampled RepresentationFeatures算principle curvatures作为它的feature discriptor5.1 算cluster, 注意s

原创 non-rigid shape registration using similarity-invariant differential coordinates

这篇其实是as rigid as possible的变种, 只不过把R换成了sR(2)->(3) 同时除个标量说得过去, 因为这里取的是向量变后的方向需要一致, 同时除个标量不影响方向(5)->(6) 公式(6)原文写错了应该为∑j∈Ni4wij(e′ij−12(siRi+sjRj)eij)\sum_{j \in N_i}4w_{ij}(e'_{ij}-\frac{1}{2}(s_iR_i+s_

原创 Optimal Step Nonrigid ICP Algorithms for surface registration

这篇paper的约束有三个:距离, stiffness, landmark

原创 Deformation Transfer for Triangle Meshes

1.where S refers to the set of triangle indices for the source mesh.这里S是指source三角形的个数2.M集合里面的pair个数是指, source里的s1三角形的形变应与target里t1三角形的形变接近3.对于公式7, |T|代表target里面三角形的个数. 公式7的几何意义, source中的

原创 Projective Dynamics: Fusing Constraint Projections for Fast Simulation

1.(3)->(4)qn+1−qn−hvn−h2Mfext=h2Mfint(qn+1)q_{n+1}- q_n-hv_n-\frac{h^2}{M}f_{ext}= \frac{h^2}{M}f_{int}(q_{n+1}) 1h2M(qn+1−sn)=fint(qn+1)\frac{1}{h^2}M(q_{n+1}- s_n)= f_{int}(q_{n+1})右边项左移然后对qn+1q_{n+

原创 矩阵积分

假设XX为未知矩阵,那么对于矩阵XX的积分 ∫XdX=12∥X∥2F\int X dX = \frac{1}{2}\|X\|_F^2另外一种方式理解对于XX中的每个元素xijx_{ij}, 对其积分∫xijdxij=12x2ij\int x_{ij}dx_{ij} = \frac{1}{2}x_{ij}^2 所以 ∫XdX=∑i∑j∫xijdxij=12∑i∑jx2ij=12∑i∑jxij∗

原创 arrayfun用法

arrayfun(@(n) norm(A(n,:)), 1:size(A,1))注意当返回不是标量的时候记得将'UniformOutput'参数改为0arrayfun(@(n) norm(A(n,:)), 1:size(A,1), 'UniformOutput', 0);

转载 GPU CUDA and MEX Programming

https://uk.mathworks.com/help/distcomp/gpu-cuda-and-mex-programming.htmlmxGPUCopyFromMxArrayCopy mxArray to mxGPUArraymxGPUCopyGPUArrayDuplicate (deep copy) mxGPU

原创 matlab中cuda编程中分配grid和block dimension的时候的注意事项

这两天在MATLAB里面用CUDA, 结果老是不对, 找了半天才发现, 原来是分配的问题所以提醒小伙伴们, 分配的时候先看下你的GPU性能>> gpuDeviceans = CUDADevice (具有属性): Name: 'Quadro K2000' Index:

原创 CUDA里面GRID, BLOCK 边界检测

dim3 block(BLOCK_X, BLOCK_Y);dim3 grid((roi.width + block.x - 1) / block.x, (roi.height + block.y - 1) / block.y);int y = blockIdx.y * blockDim.y + threadIdx.y;int x = blockIdx.x * blockDim.x + th

原创 mex 里面调用matlab函数

mxArray *rhs[1]; rhs[0] = (mxArray *)prhs[4];//phrs[4]为mexfunction传过来的arguement mxArray *nneighbors_array; mexCallMATLAB(1, &nneighbors_array, 1, rhs, "numel"); int nneighbors = (int)

原创 mexcuda输入nvcc中的参数

mexcuda -v COMPFLAGS="-arch=sm_30" show_neighbor.cu

原创 mexcuda中矩阵数据的传输

数据的传输以0开始,然后以列序传输/* * Example of how to use the mxGPUArray API in a MEX file. This example shows * how to write a MEX function that takes a gpuArray input and returns a * gpuArray output, e.g. B

转载 MATLAB Primitive Types

原创 laplacian,degree,adjacency and oriented incidence matrix, differential and laplacian coordinates

注意购造incident matrix时, 边上点索引小的为-1,大的为1这里的L没有标准化处理,要标准化处理的话不需要对称的话,则可以表示为laplacian coordinates