Conjugate transpose: Difference between revisions

 

Complex!T[][] conjugateTranspose(T)(in Complex!T[][] m)

return iota(m[0].length)

.map!(i => m.transversal(i).map!conj.array)

}

 

const Bt = B[0].length.iota.map!(i=> B.transversal(i).array).array;

return A.map!(a => Bt.map!(b => a.dotProduct(b)).array).array;

 

bool isNormal(T)(in Complex!T[][] m, in Complex!T[][] ct)

return [matMul(m, ct), matMul(ct, m)].areEqual;

}

 

.map!(r => side.iota.map!(c => complex(r == c)).array)

.array;

 

bool isUnitary(T)(in Complex!T[][] m, in Complex!T[][] ct)

const mct = matMul(m, ct);

const ident = mct.length.complexIdentitymatrix;