Cauphenuny's Blog

绕轴 n\mathbf{n}n 旋转角度 α\alphaα 的旋转矩阵是： Rn,α=cos⁡(α)I+(1−cos⁡(α))nnT+sin⁡α(0−nznynz0−nx−nynx0)⏟N=I+sin⁡(θ)N+(1−cos⁡θ)N2\mathbf{R}_{\mathbf{n},\alpha}=\cos(\alpha)\mathbf{I}+(1-\cos(\alpha))\mathbf{n}\mathbf{n}^T+\sin\alpha \underbrace{ \begin{pmatrix}0 & -n_{z} & n_{y} \\ n_{z} & 0 & -n_{x} \\ -n_{y} & n_{x} & 0\end{pmatrix} }_{ \mathbf{N}...

又刷到个奇技淫巧，看完两眼一黑了 https://www.chiark.greenend.org.uk/~sgtatham/coroutines.html switch 是能跳转到任意控制流内部的，利用这个特性，可以在 C 中实现一个 generator 朴素的 generator 需要实现一个显式的状态机 123456789101112131415161718int func() { static int state = 0, i = 0; switch (state) { case 0: goto start; case 1: goto label1; case 2: goto label2; }start: state = 1; return 0;label1: for (i = 0; i < 10; i++) { state = 2; return i;label2:; } return...

https://danielchasehooper.com/posts/typechecked-generic-c-data-structures/ 比较可惜的是，在看到这篇博客之前，我已经写完了这学期数据结构的 oj 代码了 众所周知，用 C 写非侵入式泛型容器是很坐牢的事，要么全用 void* 一顿擦除，性能不高的同时还容易写错类型，要么写一个全是 define 的头文件，根据需要什么类型就 #define T int 然后 include 一下，产生很多重复代码，而且全 define 很难 debug 或代码提示。 这篇博客给出了一个令人耳目一新的方法 效果是这样的： 12345678910111213141516171819typedef struct { int value;} Foo;List(int) int_list = {0};list_prepend(&int_list, 3);List(Foo) foo_list =...

Fundamentals Theorem For any distribution over (X,Y)(X,Y)(X,Y), we have argminf E∣∣f(X)−Y∣∣2=E[Y∣X]\text{argmin}_{f} \space \mathbb{E} ||f(X)-Y||^2=\mathbb{E}[Y|X] argminf​ E∣∣f(X)−Y∣∣2=E[Y∣X] Proof argmin⁡ fE∣∣f(x)−y∣∣2=argmin⁡ fE∣∣f(x)−y∣∣2=argmin⁡ fEXEY∣X[f(x)−Y∣X=x]2=argmin⁡ f(x)EX[f(x)2−2f(x)EY∣X(Y∣X=x)+EY∣X(Y∣X=x)2]=EX(argmin⁡ f(x)f(x)2−2f(x)EY∣X(Y∣X=x)+EY∣X(Y∣X=x)2)\begin{aligned} &\mathop{\operatorname{argmin}\,}_{f} \mathbb{E} \lvert \lvert f(x)-y \rvert \rvert...

非原创 解决方案来源 在 Docker run 中加入 --cap-add=SYS_PTRACE --security-opt seccomp=unconfined 创建一个/usr/local/bin/gdb，因为PATH中/usr/local/bin比较靠前，会比/usr/bin/gdb先被找到 1234567891011121314151617181920#!/bin/bashSCRIPT_DIR=$( cd -- "$( dirname -- "${BASH_SOURCE[0]}" )" &> /dev/null && pwd )if [ "$#" != "1" ]; then echo "Usage: $0 <path/to/program to debug>" >&2 exit 1fiprog="$1"# Start program in...

最近在做一个编译原理lab ☝️🤓诶，我已经有结构化的AST和IR了，岂不是能不太费劲地地写出一些有意思的东西？ 就这样一拍脑袋给编译器加了个brainfuck后端 由于还没想好怎么实现栈，现在所有变量都是static的，也不支持递归函数调用 demo 12345678910/// file: test/int_size.cint main() { int bit = 0; int v = 1; while (v) { v = v * 2; bit = bit + 1; } return bit;} 1234567891011121314$ build/compiler -brainz...

没用的小玩具

计科导大作业 试玩链接: 原琴模拟器 6.14 音色库来源：@Colalala_冰阔落落落 试了半天，发现没法正常播放这个风物之诗琴的音源，决定先用 smplr 默认的钢琴音源把别的写了。 写好了键盘演奏函数、谱子 parser 以及自动播放功能。 version 0.1.0 6.15 加入 [] {} 表示时值变化。 添加临时升降记号 -/+ 以及临时高/低八度记号。 重写了教程 version 0.3.0 6.16 添加了固定调转调方式，更适合对着五线谱演奏。 version 0.4.0 6.17 优化了调号显示。 version 0.4.4 加了个键盘，可以看到哪些音被按下了。 version 0.5.4 把自动演奏和按键动画连接了起来。 version 0.6.0 加了个加载时的提示悬浮窗 version 0.6.1 6.18 重构了代码 version 0.7.0 给按钮加了点阴影，感谢Box-Shadow CSS Generator以及用filter: drop-shadow()给透明图片添加阴影 version 0.7.1 遇到了个...

定义 concept 12template<typename T>concept c = (bool_expression); requires(args) {...} 可以作为一个 bool 表达式，检查大括号中内容是否编译通过 e.g. 12template <typename T>concept printable = requires(std::ostream& os, T a) { os << a; }; requires 字句也能嵌套，这样就不用把两个无关的条件参数写到同一个括号里面。 12template <typename T>concept Field = <concepts> 头文件内有一些预定义的...

2024 年 3 月 屁事没干，破防了 对比一下