夜雨聆风字数 3143，阅读大约需 16 分钟
第十一章：数学协处理器模拟：软件实现的浮点运算
80386 CPU本身不包含浮点运算单元，浮点运算由独立的80387数学协处理器完成。但并非所有系统都配备了协处理器，因此Linux 0.12实现了协处理器模拟，通过软件模拟浮点指令的执行，让没有协处理器的系统也能运行浮点程序。
11.1 协处理器架构
80387协处理器有8个80位浮点寄存器（ST0-ST7），组织为栈结构。浮点指令从内存或寄存器栈加载操作数，执行运算，将结果压入栈。协处理器还有控制字寄存器（设置舍入模式、精度、异常屏蔽）和状态字寄存器（记录异常标志、栈指针、条件码）。
80387寄存器与状态
80387的完整状态保存在i387_state结构中，包含控制字CW（舍入模式、精度、异常屏蔽位）、状态字SW（异常标志、栈顶指针TOP、条件码）、标签字TW（标记每个寄存器的类型），以及8个80位扩展精度浮点寄存器ST0-ST7。每个进程的task_struct中保存一份i387状态，支持上下文切换时保存和恢复。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
/* include/linux/math_emu.h - 协处理器状态结构 */
/* 80位扩展精度浮点数: bit79=符号, bit78-64=指数(偏置16383), bit63-0=尾数 */
typedefstruct {
unsigned short exp; /* 符号+指数 */
unsigned long hi; /* 尾数高32位 */
unsigned long lo; /* 尾数低32位 */
} FP_REG;
struct i387_state {
unsigned short control_word; /* CW: 舍入/精度/异常屏蔽 */
unsigned short status_word; /* SW: 异常标志/栈指针TOP/条件码 */
unsigned short tag_word; /* TW: 每个寄存器类型(有效/零/特殊/空) */
unsigned long fip; /* 指令指针 */
unsigned short fcs; /* 指令段选择符 */
unsigned short opcode; /* 操作码 */
unsigned long foo; /* 操作数偏移 */
unsigned short fos; /* 操作数段选择符 */
FP_REG st[8]; /* ST0-ST7浮点寄存器栈 */
};
协处理器指令以ESC转义前缀（0xD8-0xDF）开头，后跟操作码和操作数。例如FADD指令（浮点加法）、FMUL指令（浮点乘法）、FSQRT指令（平方根）。协处理器与主CPU协同工作：主CPU执行到ESC指令时，将指令和操作数传递给协处理器，协处理器执行浮点运算，主CPU继续执行下一条指令（异步并行）。
如果系统没有协处理器，主CPU执行ESC指令时会触发device_not_available异常（INT 7）。Linux 0.12的异常处理程序检测到这是浮点指令，就调用协处理器模拟器math_emulate()来软件模拟指令执行。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
┌─────────────┐
│ 执行FADD指令 │
└──────┬──────┘
│
▼
┌──────────────────┐
│ 有协处理器? │
└────────┬─────────┘
│
┌────────────┴──────────────┐
│ 是 │ 否
▼ ▼
┌────────┐ ┌─────────────┐
│ 硬件执行 │ │ 软件模拟执行 │
│ (快速) │ │ (慢速，透明) │
└────────┘ └─────────────┘
11.2 协处理器模拟器
math_emulate.c是协处理器模拟器的入口。它接收异常时保存的指令地址和操作数，解析浮点指令的操作码，根据指令类型调用相应的处理函数。
device_not_available异常处理
当CPU执行浮点指令时，如果CR0的TS位（任务切换时置1）或EM位（无协处理器时置1）被设置，就触发INT 7异常。处理程序先清除TS位，然后根据情况恢复或初始化当前进程的浮点状态。Linux采用"延迟保存"策略——只在进程实际使用浮点时才保存/恢复协处理器状态，避免无谓的开销。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
/* kernel/traps.c - 设备不可用异常处理 */
/* CR0: TS位(位3)=任务切换标志, EM位(位2)=协处理器模拟标志 */
void device_not_available(void)
{
__asm__("clts"); /* 清除CR0.TS位 */
if (current->flags & PF_USEDFPU) {
/* 进程曾使用浮点: 恢复其协处理器状态 */
math_state_restore();
return;
}
/* 第一次使用: 初始化浮点状态 */
current->flags |= PF_USEDFPU;
init_i387(¤t->i387);
}
math_state_restore与延迟保存
Linux采用"延迟保存"策略处理协处理器状态，即只有当新任务第一次使用浮点运算时才恢复状态，而不是在每次任务切换时保存/恢复。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
/* kernel/math/math_emulate.c - 协处理器状态恢复 */
struct task_struct *last_task_used_math = NULL;
void math_state_restore(void)
{
if (last_task_used_math == current) return; /* 同一进程,无需操作 */
if (last_task_used_math)
save_i387(&last_task_used_math->i387); /* 保存上一个任务的FPU状态 */
if (current->used_math)
restore_i387(¤t->i387); /* 恢复当前任务的FPU状态 */
else {
init_i387(¤t->i387); /* 首次使用,初始化 */
current->used_math = 1;
}
last_task_used_math = current; /* 更新FPU拥有者 */
}
math_emulate入口
当CR0.EM=1时，所有浮点指令都会产生INT 7。math_emulate()从触发异常的eip处读取ESC前缀字节(0xD8-0xDF)确定指令组，再读取ModR/M字节解析寻址模式，然后分派到do_esc0~do_esc7八个处理函数。每个处理函数再根据Reg字段确定具体操作（加法、乘法、比较等）。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
/* kernel/math/math_emulate.c - 模拟器入口 */
void math_emulate(long edi, long esi, long ebp, long esp, long ebx,
long edx, long ecx, long eax,
unsigned short ds, es, fs, gs,
unsigned long eip, unsigned short cs, unsigned long eflags)
{
unsigned char opcode = get_fs_byte((unsigned char *)eip);
unsigned char modrm = get_fs_byte((unsigned char *)(eip + 1));
unsigned long addr;
decode_ea(modrm, &addr); /* 解析ModR/M计算有效地址 */
switch (opcode) {
case 0xD8: do_esc0(modrm, addr); break; /* 单精度运算 */
case 0xD9: do_esc1(modrm, addr); break; /* 控制/状态/常数 */
case 0xDC: do_esc4(modrm, addr); break; /* 双精度运算 */
/* 0xDA, 0xDB, 0xDD, 0xDE, 0xDF 类似处理 */
default:
printk("math_emulate: unknown ESC %02x\n", opcode);
do_exit(SIGFPE);
}
/* 根据指令长度调整eip */
}
/* ModR/M字节: bit7-6=Mod(寻址模式), bit5-3=Reg(子操作码), bit2-0=R/M */
模拟器维护一个模拟的协处理器状态，保存在进程的task_struct->i387字段中。这个字段包含8个浮点寄存器、控制字、状态字、标志字等，与真实的80387状态完全对应。
浮点指令的模拟分为几个步骤：首先解析指令获取操作码和地址模式，然后从内存或寄存器栈加载操作数（可能需要格式转换，80387支持单精度、双精度、扩展精度等多种格式），调用相应的运算函数（加减乘除、开方、三角函数等），将结果写回寄存器栈或内存，更新状态字（异常标志、栈指针、条件码）。
由于浮点运算的复杂性，模拟器的实现非常庞大。Linux 0.12的math子目录包含多个C文件：add.c实现加法，mul.c实现乘法，div.c实现除法，compare.c实现比较，get_put.c实现加载和存储，convert.c实现格式转换。每个文件都有数百行代码，处理各种精度、舍入模式、特殊值（无穷大、NaN）。
11.3 浮点数表示
80387使用IEEE 754浮点标准。扩展精度格式（80位）包含1位符号、15位指数、64位尾数。双精度格式（64位）包含1位符号、11位指数、52位尾数。单精度格式（32位）包含1位符号、8位指数、23位尾数。
IEEE 754格式细节
IEEE 754定义了三种浮点格式。单精度32位（1位符号+8位指数+23位尾数，偏置127）、双精度64位（1位符号+11位指数+52位尾数，偏置1023）、扩展精度80位（1位符号+15位指数+64位尾数，偏置16383）。特殊值通过指数全0或全1编码：零、非规格化数、无穷大、NaN。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
/* 单精度浮点数（32位） */
typedefunion {
float f;
struct {
unsigned int mantissa:23; /* 尾数 */
unsigned int exponent:8; /* 指数(偏置127) */
unsigned int sign:1; /* 符号 */
} parts;
} float32;
/* 双精度浮点数（64位） */
typedefunion {
double d;
struct {
unsigned long mantissa:52;
unsigned long exponent:11; /* 偏置1023 */
unsigned long sign:1;
} parts;
} float64;
/* 特殊值编码:
* 指数全0+尾数全0 = ±0 指数全0+尾数非0 = 非规格化数
* 指数全1+尾数全0 = ±∞ 指数全1+尾数非0 = NaN
* 其他 = 规格化数(隐含最高位1) */
模拟器中的浮点表示
模拟器内部使用80位扩展精度运算以保证精度。加载单精度数到扩展精度寄存器时，需要处理零、非规格化数（需规范化）、无穷大/NaN等特殊情况，并调整指数偏置（127→16383）和尾数对齐。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23
/* 将单精度数加载到扩展精度寄存器中 */
void load_single(FP_REG *st, float32 src)
{
int e = src.parts.exponent;
unsigned int m = src.parts.mantissa;
if (e == 0 && m == 0) {
st->exp = 0; st->hi = st->lo = 0; /* 零 */
} else if (e == 0) {
/* 非规格化数: 规范化 */
while (!(m & 0x400000)) { m <<= 1; e--; }
e++;
st->exp = (e - 127 + 16383) | (src.parts.sign << 15);
st->hi = (unsigned long)m << 8; st->lo = 0;
} else if (e == 255) {
st->exp = 0x7FFF | (src.parts.sign << 15); /* ∞/NaN */
st->hi = st->lo = 0xFFFFFFFF;
} else {
m |= 0x800000; /* 隐含的1 */
st->exp = (e - 127 + 16383) | (src.parts.sign << 15);
st->hi = (unsigned long)m << 8; st->lo = 0;
}
}
浮点运算需要处理这些特殊情况。例如正数除以零得正无穷，零除以零得NaN，无穷大加无穷大得无穷大，无穷大减无穷大得NaN。模拟器必须正确实现这些语义，保证与硬件协处理器的行为一致。
11.4 舍入与精度
浮点运算结果可能无法精确表示（例如1/3），需要舍入。80387支持四种舍入模式：向最近舍入（默认，四舍五入）、向零舍入（截断）、向正无穷舍入（向上）、向负无穷舍入（向下）。控制字寄存器的RC字段设置舍入模式。
舍入模式与精度控制
控制字的RC字段（bit10-11）控制舍入方向，PC字段（bit8-9）控制运算精度（24/53/64位）。模拟器在运算后根据这两个设置调整结果尾数，处理舍入引起的精度异常。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
/* 控制字中的舍入模式(RC)和精度控制(PC) */
#define RC_NEAR 0x0000 /* 向最近舍入(默认) */
#define RC_DOWN 0x0400 /* 向负无穷 */
#define RC_UP 0x0800 /* 向正无穷 */
#define RC_CHOP 0x0C00 /* 向零截断 */
#define PC_24 0x0000 /* 单精度(24位) */
#define PC_53 0x0200 /* 双精度(53位) */
#define PC_64 0x0300 /* 扩展精度(64位) */
/* 舍入流程(伪代码):
* 1. 根据PC确定保留位数(24/53/64)
* 2. 分析尾数多余低位和舍入位/粘滞位
* 3. 根据RC选择舍入方向
* 4. 检查是否触发精度/溢出/下溢异常 */
11.5 性能考虑
软件模拟浮点运算比硬件协处理器慢得多，大约慢100倍甚至更多。一条简单的FADD指令，硬件协处理器几个时钟周期即可完成，而软件模拟需要成千上万条指令。因此对于浮点密集型程序（如科学计算、图形渲染），没有协处理器的系统性能会严重下降。
Linux 0.12在实现时尽量优化模拟器性能，例如使用查表法加速某些运算，使用整数运算模拟浮点运算（因为整数运算快得多）。但即使如此，软件模拟仍然是性能瓶颈。好在到90年代中期，浮点单元已经集成到CPU内部（如Pentium），协处理器模拟逐渐不再需要。
Linux 0.12的协处理器模拟展示了操作系统的适应性：通过软件补偿硬件的不足，让系统能在各种配置下运行。虽然性能不佳，但总比完全无法运行浮点程序要好。这种权衡在早期计算机时代非常常见。
通过本章的学习，我们理解了Linux 0.12的协处理器模拟实现：从80387架构到模拟器工作原理，从浮点数表示到舍入和精度，以及性能考虑。协处理器模拟虽然复杂且性能不佳，但它让Linux能够在没有浮点硬件的系统上运行，体现了操作系统的兼容性设计。
11.6 参考资料
本章内容基于以下资料编写：
Intel官方文档
• Intel 80387 Programmer's Reference Manual (1987) • 详细描述了80387协处理器架构、指令集、浮点格式
IEEE标准
• IEEE 754-1985 Standard for Binary Floating-Point Arithmetic • 定义了浮点数表示和运算语义
Linux 0.12源代码文件
• kernel/math/math_emulate.c - 协处理器模拟器入口 • kernel/math/add.c - 浮点加法 • kernel/math/mul.c - 浮点乘法 • kernel/math/div.c - 浮点除法 • kernel/math/compare.c - 浮点比较 • kernel/math/get_put.c - 加载和存储 • kernel/math/convert.c - 格式转换 • kernel/math/ea.c - 有效地址计算 • kernel/traps.c - device_not_available异常处理 • include/linux/math_emu.h - 协处理器模拟器数据结构
"精诚所至，金石为开。" —— 《后汉书》
