Differences

This shows you the differences between two versions of the page.

Link to this comparison view

forum:data:2020:problem_of_chop [2020/04/10 08:44] – Created from the form at forum:start Chang-Yang Kuoforum:data:2020:problem_of_chop [2020/04/17 10:17] (current) Maurits W. Haverkort
Line 19: Line 19:
  
 ++++++++   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)
Line 489: Line 489:
 )*(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~~
  
Print/export