AVR如何执行浮点运算

我正在尝试在没有FPU的CPU上实现对double和的支持float以及相应的基本算法。

我知道这一切都有可能AVR ATmega controllers。ATmega也没有FPU。因此出现了一个问题：它是如何工作的？是否有任何文献建议或带有解释和实例的链接？

在最好的情况下，我将为以下代码提供支持：

double twice ( double x )
{
  return x*x;
}

提前谢谢，亚历克斯