用c语言写矩阵开平方的函数

怎么用c语言实现矩阵开平方的函数呢？类似matlab的sqrtm函数

c语言里 有 sqrt函数,在math.h头文件里
玩矩阵 用boost库吧。