Groups > comp.soft-sys.math.mathematica > #1722

Re: Another AppendTo replacement problem

From	Iván Lazaro <gaminster@gmail.com>
Newsgroups	comp.soft-sys.math.mathematica
Subject	Re: Another AppendTo replacement problem
Date	2011-04-16 11:37 +0000
Organization	Steven M. Christensen and Associates, Inc and MathTensor, Inc.
Message-ID	<iobv13$bfp$1@smc.vnet.net> (permalink)

Show all headers | View raw

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!

Back to comp.soft-sys.math.mathematica | Previous | Next — Next in thread | Find similar | Unroll thread

Thread

Re: Another AppendTo replacement problem Iván Lazaro <gaminster@gmail.com> - 2011-04-16 11:37 +0000
  Re: Another AppendTo replacement problem ADL <alberto.dilullo@tiscali.it> - 2011-04-17 11:51 +0000
  Re: Another AppendTo replacement problem "Sjoerd C. de Vries" <sjoerd.c.devries@gmail.com> - 2011-04-17 11:51 +0000
    Re: Another AppendTo replacement problem Ray Koopman <koopman@sfu.ca> - 2011-04-17 23:15 +0000

csiph-web