Differences

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

Link to this comparison view

Next revision		Previous revision
documentation:language_reference:functions:chop [2016/10/09 21:27] – created Maurits W. Haverkortdocumentation:language_reference:functions:chop [2023/08/30 15:17] (current) – Maurits W. Haverkort
Line 1: Line 1:
 +====== Chop ======
  
 +###
 +Chop(//a//) replaces approximate real numbers in //a// that are close to zero by the exact integer 0. Chop(//a//, $\epsilon$) replaces numbers smaller in absolute magnitude than $\epsilon$ by 0. Chop uses a default tolerance of $\epsilon=10^{-10}$. Chop works on Real and Complex numbers as well as on operators, wavefunctions and Tight binding objects. 
 +###
 +
 +===== Input =====
 +
 +  * //a// : object that will be chopped, all types allowed
 +  * $\epsilon$ : real number setting the boundary between zero and finite
 +
 +===== Output =====
 +
 +  * //a'// : object of same type as //a//
 +
 +===== Example =====
 +
 +###
 +A small test script of chop is:
 +###
 +
 +==== Input ====
 +<code Quanty Chop.Quanty>
 +dofile("../definitions.Quanty")
 +
 +Opp = Opp1 + 0.000001*Opp2
 +print(Opp)
 +Opp3 = Chop(Opp,0.001)
 +print(Opp3)
 +
 +psi = psi1 + 0.00001 * psi2
 +print(psi)
 +psi3 = Chop(psi,0.001)
 +print(psi3)
 +
 +a = 0.000001 + 0.3 * I
 +b = Chop(a,0.001)
 +print("a = ",a,"\nChop(a) = ",b)
 +</code>
 +
 +==== Result ====
 +<file Quanty_Output Chop.out>
 +
 +Operator: Lx
 +QComplex         =          1 (Real==0 or Complex==1 or Mixed==2)
 +MaxLength        =          2 (largest number of product of lader operators)
 +NFermionic modes =          6 (Number of fermionic modes (site, spin, orbital, ...) in the one particle basis)
 +NBosonic modes   =          0 (Number of bosonic modes (phonon modes, ...) in the one particle basis)
 +
 +Operator of Length   2
 +QComplex      =          1 (Real==0 or Complex==1)
 +N             =          8 (number of operators of length   2)
 +C  3 A  1 |  7.071067811865476E-01 -7.071067811865476E-07
 +C  1 A  3 |  7.071067811865476E-01  7.071067811865476E-07
 +C  2 A  0 |  7.071067811865476E-01 -7.071067811865476E-07
 +C  0 A  2 |  7.071067811865476E-01  7.071067811865476E-07
 +C  5 A  3 |  7.071067811865476E-01 -7.071067811865476E-07
 +C  3 A  5 |  7.071067811865476E-01  7.071067811865476E-07
 +C  4 A  2 |  7.071067811865476E-01 -7.071067811865476E-07
 +C  2 A  4 |  7.071067811865476E-01  7.071067811865476E-07
 +
 +
 +
 +Operator: Lx
 +QComplex         =          0 (Real==0 or Complex==1 or Mixed==2)
 +MaxLength        =          2 (largest number of product of lader operators)
 +NFermionic modes =          6 (Number of fermionic modes (site, spin, orbital, ...) in the one particle basis)
 +NBosonic modes   =          0 (Number of bosonic modes (phonon modes, ...) in the one particle basis)
 +
 +Operator of Length   2
 +QComplex      =          0 (Real==0 or Complex==1)
 +N             =          8 (number of operators of length   2)
 +C  3 A  1 |  7.071067811865476E-01
 +C  1 A  3 |  7.071067811865476E-01
 +C  2 A  0 |  7.071067811865476E-01
 +C  0 A  2 |  7.071067811865476E-01
 +C  5 A  3 |  7.071067811865476E-01
 +C  3 A  5 |  7.071067811865476E-01
 +C  4 A  2 |  7.071067811865476E-01
 +C  2 A  4 |  7.071067811865476E-01
 +
 +
 +
 +WaveFunction: Wave Function
 +QComplex         =          0 (Real==0 or Complex==1)
 +N                =          3 (Number of basis functions used to discribe psi)
 +NFermionic modes =          6 (Number of fermions in the one particle basis)
 +NBosonic modes   =          0 (Number of bosons in the one particle basis)
 +
 +#      pre-factor         Determinant
 +   1   3.162277660168E-01       110000
 +   2   9.486832980505E-01       100010
 +   3   1.000000000000E-05       010010
 +
 +
 +
 +WaveFunction: Wave Function
 +QComplex         =          0 (Real==0 or Complex==1)
 +N                =          2 (Number of basis functions used to discribe psi)
 +NFermionic modes =          6 (Number of fermions in the one particle basis)
 +NBosonic modes   =          0 (Number of bosons in the one particle basis)
 +
 +#      pre-factor         Determinant
 +   1   9.486832980505E-01       100010
 +   2   3.162277660168E-01       110000
 +
 +
 +a = (1e-06 + 0.3 I)
 +Chop(a) = (0 + 0.3 I)
 +
 +</file>
 +
 +==== Options ====
 +
 +The last element of Chop can be a table of options. Possible options are:
 +  * "Distance" For tight binding objects one can chop according to the size of the hopping matrix element (keep all that are larger) or chop according to the distance between the atoms (keep all interactions that connect sites with a distance that is smaller than the given length).
 +
 +===== Table of contents =====
 +{{indexmenu>.#1}}

You are here: Quanty - a quantum many body script language » Documentation » Language Reference » Functions » Chop

Print/export