| |
forum:data:2020:problem_of_chop [2020/04/10 08:44] – Created from the form at forum:start Chang-Yang Kuo | forum:data:2020:problem_of_chop [2020/04/17 10:17] (current) – Maurits W. Haverkort |
---|
| |
++++++++ The operator I am using ++++++++++++++ | ++++++++ The operator I am using ++++++++++++++ |
| <code> |
Udd = 0.0166631955291370*(0.707106781186548*NewOperator("Cr",NF,0)+I*0*NewOperator("Cr",NF,0) + 0*NewOperator("Cr",NF,2)+I*0*NewOperator("Cr",NF,2) + 0*NewOperator("Cr",NF,4)+I*0*NewOperator("Cr",NF,4) + 0*NewOperator("Cr",NF,6)+I*0*NewOperator("Cr",NF,6) + 0*NewOperator("Cr",NF,8)+I*0.707106781186548*NewOperator("Cr",NF,8) | Udd = 0.0166631955291370*(0.707106781186548*NewOperator("Cr",NF,0)+I*0*NewOperator("Cr",NF,0) + 0*NewOperator("Cr",NF,2)+I*0*NewOperator("Cr",NF,2) + 0*NewOperator("Cr",NF,4)+I*0*NewOperator("Cr",NF,4) + 0*NewOperator("Cr",NF,6)+I*0*NewOperator("Cr",NF,6) + 0*NewOperator("Cr",NF,8)+I*0.707106781186548*NewOperator("Cr",NF,8) |
)*(0.707106781186548*NewOperator("Cr",NF,7)+I*0*NewOperator("Cr",NF,7) + 0*NewOperator("Cr",NF,1)+I*0*NewOperator("Cr",NF,1) + 0*NewOperator("Cr",NF,3)+I*0*NewOperator("Cr",NF,3) + 0*NewOperator("Cr",NF,5)+I*-0.707106781186548*NewOperator("Cr",NF,5) + 0*NewOperator("Cr",NF,9)+I*0*NewOperator("Cr",NF,9) | )*(0.707106781186548*NewOperator("Cr",NF,7)+I*0*NewOperator("Cr",NF,7) + 0*NewOperator("Cr",NF,1)+I*0*NewOperator("Cr",NF,1) + 0*NewOperator("Cr",NF,3)+I*0*NewOperator("Cr",NF,3) + 0*NewOperator("Cr",NF,5)+I*-0.707106781186548*NewOperator("Cr",NF,5) + 0*NewOperator("Cr",NF,9)+I*0*NewOperator("Cr",NF,9) |
)*(0*NewOperator("An",NF,0)+I*0*NewOperator("An",NF,0) + 0*NewOperator("An",NF,4)+I*0*NewOperator("An",NF,4) + 0*NewOperator("An",NF,6)+I*0*NewOperator("An",NF,6) + 0*NewOperator("An",NF,8)+I*0*NewOperator("An",NF,8) + 1.00000000000000*NewOperator("An",NF,2)+I*0*NewOperator("An",NF,2) | )*(0*NewOperator("An",NF,0)+I*0*NewOperator("An",NF,0) + 0*NewOperator("An",NF,4)+I*0*NewOperator("An",NF,4) + 0*NewOperator("An",NF,6)+I*0*NewOperator("An",NF,6) + 0*NewOperator("An",NF,8)+I*0*NewOperator("An",NF,8) + 1.00000000000000*NewOperator("An",NF,2)+I*0*NewOperator("An",NF,2) |
) | ) |
| </code> |
| |
| <code> |
StartRestrictions = {NF, NB, {"1111111111 000000",2,2}} | StartRestrictions = {NF, NB, {"1111111111 000000",2,2}} |
Eigensystem(Udd, StartRestrictions, 10, {{"Zero", 0.0000000001},{"Epsilon",0.0000000001}}) --it works | Eigensystem(Udd, StartRestrictions, 10, {{"Zero", 0.0000000001},{"Epsilon",0.0000000001}}) --it works |
Eigensystem(Chop(Udd), StartRestrictions, 10, {{"Zero", 0.0000000001},{"Epsilon",0.0000000001}}) --it doesn't work | Eigensystem(Chop(Udd), StartRestrictions, 10, {{"Zero", 0.0000000001},{"Epsilon",0.0000000001}}) --it doesn't work |
| </code> |
| |
</WRAP> | </WRAP> |
| |
~~DISCUSSION|Answers~~ | ~~DISCUSSION|Answers~~ |
| |