public void posdef2()
{
var A = BlockMatrix.Define(new float[] { 4, 12, -16, 12, 37, -43, -16, -43, 98 }, n3);
var B = A.Replicate();
Claim.yea(mkl.potrf(B));
}
public void posdef1()
{
var notPD = BlockMatrix.Define(new float[] { 1, 2, 2, 1 }, n2);
var eval1 = mkl.posdef(notPD);
Claim.nea(eval1);
var PD = BlockMatrix.Define(new float[] { 2, -1, 0, -1, 2, -1, 0, -1, 2 }, n3);
var eval2 = mkl.posdef(PD);
Claim.yea(eval2);
}