Re: Another AppendTo replacement problem
- To: mathgroup at smc.vnet.net
- Subject: [mg118208] Re: Another AppendTo replacement problem
- From: Iván Lazaro <gaminster at gmail.com>
- Date: Sat, 16 Apr 2011 07:36:07 -0400 (EDT)
I made a mistake in the code. Now it's fine. Sorry.
NumBasis = 10000;
q = matrA = ma = Table[0, {i, 2}];
M = RandomComplex[{-1 - I, 1 + I}, {NumBasis, 2, 2}];
M = Map[Orthogonalize, M];
matr = RandomComplex[{-1 - I, 1 + I}, {2, 2}]
Results = {};
Do[{ma[[k]] =
KroneckerProduct[M[[Nbase, k]], Conjugate[M[[Nbase, k]]]];
matrA[[k]] = Chop[matr.ma[[k]]];
matrA[[k]] = matrA[[k]]/Tr[matrA[[k]].matrA[[k]]] // Chop;
If[k == 2,
AppendTo[
Results, {M[[Nbase]], Total[Eigenvalues[matrA[[k]]]]}]];
}, {Nbase, 1, NumBasis}, {k, 1, 2}];
M = Sort[Results, #1[[2]] < #2[[2]] &][[1, 1]];
Thanks in advance!