> restart:
> Digits:=14:
> st:= time():
> with(LinearAlgebra):
> with(PolynomialTools):
> with(DiscreteTransforms):
> with(ImageTools):
> read"../Src/One_GreyImage.mpl":
> infolevel[S_stls]:=1:
> interface(rtablesize=690):
> printlevel:=0:
> Ph:=Read("../Result/Figure.3/OGI_blur.jpg"):
> Pho1:=RGBtoGray(Ph):
> View([Ph,Pho1]);
Initializing Java runtime environment.
> Write("../Result/Figure.3/OGI_blur.jpg",Pho1):
> e:=3:ydeg:=25:
> Phoo:=SD_ImageDeconvolution(Pho1,e,ydeg):
SD_ImageDeconvolution:

           Before we compute gcd  using computing time is


                             st1 = 0.110

Ap1D_gcd:

                    the time for gcd from cofactor


                              t = 0.015

SD_Dividen:

   we get Polynomials' coefficient matrix using computing time is


                               st1 = 0.

SD_Dividen:

                        The first FFT time is


                             st1 = 0.094

SD_Dividen:

                       The second FFT time is


                             st1 = 0.125

SD_Dividen:

                         The divide time is


                             st1 = 0.578

SD_Dividen:

                       The inverse FFT time is


                             st1 = 0.156

SD_ImageDeconvolution:

            After we compute gcd  using computing time is


                             st1 = 1.203

> View([Ph,Pho1,Phoo]);
> st:=time()-st;#the time is

                             st := 3.672

> Write("../Result/Figure.3/OGI_deblur.jpg",Phoo):
> 
>