next | previous | forward | backward | up | top | index | toc | directory | Macaulay 2 web site

clean -- Set to zero elements which are approximately zero

Synopsis

Description

If the input is a matrix or a ring element, then the result has the same type, where each real or complex number coefficient which is less than epsilon in absolute value is replaced with zero.

i1 : e = 1e-11;
i2 : M = random(RR^4,RR^5)

o2 = | .346961  .121086  .558462 .965621 .437911 |
     | .462035  .234583  .601796 .989548 .14228  |
     | .0809594 .959007  .20566  .209385 .643739 |
     | .396497  .0974429 .787142 .867366 .941691 |

                4          5
o2 : Matrix RR    <--- RR
              53         53
i3 : M * M^-1

o3 = | 1           4.44089e-16 -1.11022e-16 -4.44089e-16 |
     | 8.88178e-16 1           -1.11022e-16 -6.66134e-16 |
     | 1.11022e-16 2.22045e-16 1            -3.88578e-16 |
     | 0           0           0            1            |

                4          4
o3 : Matrix RR    <--- RR
              53         53
i4 : clean_e oo

o4 = | 1 0 0 0 |
     | 0 1 0 0 |
     | 0 0 1 0 |
     | 0 0 0 1 |

                4          4
o4 : Matrix RR    <--- RR
              53         53
Cleaning a polynomial is a way to get rid of small terms.
i5 : CC[x];
i6 : f = product(5,j -> x - exp(2*pi*j*ii/5))

      5                                  4                   3     
o6 = x  + (2.22045e-16 - 1.11022e-16*ii)x  - 1.11022e-16*ii*x  + (-
     ------------------------------------------------------------------------
                                   2
     1.11022e-16 - 1.11022e-16*ii)x  - 1.11022e-16*ii*x - 1 + 5.55112e-16*ii

o6 : CC  [x]
       53
i7 : clean_e f  

      5
o7 = x  - 1

o7 : CC  [x]
       53

See also

Ways to use clean :