Ada不同定点数类型数学运算报错的高效解决方法

该编译错误是Ada强类型规则的直接体现：不同命名的定点类型即使delta、取值范围完全一致，也属于完全独立的类型，彼此之间没有预定义的算术运算符。代码中Gain_Type和Offset_Type是两个独立定义的定点类型，Adc_Value_Type * Gain的运算结果属于Gain_Type类的定点运算结果，无法直接和Offset_Type类型的操作数做减法，因此触发报错。
盲目将所有运算数强制转换为最大位宽类型的方案确实会增加不必要的运算开销，在资源紧张的小型嵌入式目标上完全不可取。

不需要做全类型宽转，核心是利用Ada定点数的类型转换规则精准控制运算中间值类型，所有运算最终会被GNAT编译为纯整数运算，没有浮点开销，生成的代码和手写整数定点运算效率完全一致：

1. 首先确认量纲匹配：补偿公式(ADC原始值 * 增益 - 偏移) / 满量程值中，ADC原始值 * 增益和偏移的实际物理量纲完全一致，仅因为类型定义被编译器判定为不兼容。
2. 仅需在减法前，把偏移值做同delta的显式类型转换即可，不需要扩展位宽。
   修正后的Compensate函数代码如下：

function Compensate
 (Adc : in Adc_Encoded_Type; Gain : in Gain_Type; Offset : in Offset_Type)
  return Adc_Delta_Type
is
   -- 乘法结果与Gain_Type同delta，直接对偏移做同类型转换即可
   Scaled   : constant Gain_Type := Adc_Value_Type (Adc) * Gain;
   Adjusted : constant Gain_Type := Scaled - Gain_Type(Offset);
begin
   return Adc_Delta_Type(Adjusted / Adc_Mod);
end Compensate;

说明：此处Gain_Type(Offset)的转换是零开销的——两个类型delta完全一致、取值范围完全覆盖，GNAT编译时不会插入任何额外指令，仅做编译期的类型合法性校验。所有定点运算最终都会被映射为整数移位、加减乘操作，在Cortex-M0/M3这类无FPU的小型嵌入式目标上，运行效率和手写整数定点代码完全一致，不会引入浮点运算开销。

• 若不需要靠强类型做物理量区分，不要为同delta、同范围的量重复定义独立定点类型；如果确实需要靠类型约束取值范围，可以用子类型而非新类型，例如subtype Offset_Subtype is Gain_Type range Offset_Min .. Offset_Max;，这样子类型和父类型运算符兼容，不会触发类型不匹配错误，同时保留范围校验能力。
• 定点运算的中间值范围、delta不需要手动指定，GNAT会自动按照Ada标准选择满足精度、范围要求的最小整数类型作为底层实现，只要不用Float做中间中转，就不会生成浮点指令。
• 定点数乘除的缩放因子尽量使用2的整数次幂（如代码中Adc_Mod为2^10），编译器会自动将对应的乘除操作优化为移位操作，运算速度快于常量整数乘除。
• 不需要刻意把所有值转换为最大位宽类型：只要运算的范围、精度在编译期可推导，编译器会自动处理中间值位宽，既不会溢出，也不会浪费寄存器资源。

test2.adb:23:61: error: invalid operand types for operator "-"
test2.adb:23:61: error: left operand has type "Gain_Type" defined at line 11
test2.adb:23:61: error: right operand has type "Offset_Type" defined at line 12

procedure Test2 is
   Adc_Width   : constant Positive := 10;
   Adc_Delta   : constant Float    := 2.0**(-Adc_Width);
   Adc_Mod   : constant    := 2**Adc_Width;
   Error_Delta : constant          := 2.0**(-1);
   Gain_Min    : constant Float    := 1.0 - 2.0 * Adc_Delta;
   Gain_Max    : constant Float    := 1.0 + 2.0 * Adc_Delta;
   Offset_Min  : constant Float    := -0.5 * Adc_Delta;
   Offset_Max  : constant Float    := 2.0 * Adc_Delta;
   type Gain_Type is delta Adc_Delta * Error_Delta range Gain_Min .. Gain_Max;
   type Offset_Type is
      delta Adc_Delta * Error_Delta range Offset_Min .. Offset_Max;
   type Adc_Encoded_Type is mod Adc_Mod with
      Size => 16;
   subtype Adc_Value_Type is natural range 0 .. Adc_Encoded_Type'Modulus - 1;
   type Adc_Delta_Type is delta Adc_Delta range 0.0 .. 1.0 - Adc_Delta;
   function Compensate
    (Adc : in Adc_Encoded_Type; Gain : in Gain_Type; Offset : in Offset_Type)
     return Adc_Delta_Type
   is
   begin
      return Adc_Delta_Type (((Adc_Value_Type (Adc) * Gain) - Offset) / Adc_Mod);
   end Compensate;
begin
end Test2;

内容的提问来源于stack exchange，提问作者Torsten Knodt