Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: compressed sensing ( cs ) aims at recovering high dimensional sparse vectors from considerably fewer linear measurements . the problem of sparse recovery through @xmath6 norm minimization is not tractable . e. candes @xcite , d. donoho @xcite , have made pioneering contributions reposing the problem as a simple linear programming problem ( lpp ) . they have then established the conditions that ensure the stated equivalence between the original @xmath6 problem and its reposed version . it is known that restricted isometry property ( rip ) is one sufficient condition to ensure the equivalence . random matrices with gaussian or bernoulli entries have been shown to satisfy rip with high probability @xcite .    in the recent literature on cs @xcite , the deterministic construction of cs matrices has gained momentum . a. amini et . @xcite , r. devore @xcite , p. indyk @xcite and s. li et . @xcite have constructed deterministic binary sensing matrices using ideas from algebra and graph theory . in the present work , however , we attempt to construct deterministic binary sensing matrices using euler squares . the advantages of the proposed methodology are as follows :    * matrices of general row size ( different from a prime and its square ) and small density are constructed . * simplicity in construction is achieved    here , density is defined as the ratio of number of nonzero entries to the total number of entries of the matrix . sparse sensing matrices may contribute to fast processing with low computational complexity in compressed sensing @xcite . in the recent past , gaussian and @xmath7 bernoulli random matrices ( which are shown to satisfy rip @xcite ) have been used to project data into lower dimension for the purpose of classification @xcite . there are , however , the following advantages of using deterministic binary matrices for the stated purpose :    * binary matrices being sparse and possessing @xmath8 as elements provide multiplier - less and faster dimensionality reduction operation , which is not possible with gaussian matrices * there is a nonzero probability of non - compliance with rip for gaussian matrices    the problem ( @xcite ) of searching for similar images in a large image repository based on content is called content based image retrieval ( cbir ) . cbir has several multimedia applications like text based retrieval ( such as google search ) and retrieval from medical databases . as images of large size typically involve more complexity , one needs to project them to lower dimensional spaces . it is demonstrated that the proposed binary sensing matrices project data into lower dimensional spaces in such a way that the reduced vectors are useful for the purpose of cbir . moreover , the proposed dimensionality reduction technique through binary sensing matrices allows for reconstruction , which is important for tele - medicine ( @xcite ) . the other dimensionality reduction techniques and even the sparsity seeking dictionary based methods @xcite in general do not provide this advantage . the paper is organized into several sections . in section [ sec:2 ] , basic cs theory and the equivalence between @xmath6  norm and @xmath9 norm problems are given . in section [ sec:3 ] , the deterministic construction procedure of binary sensing matrices using euler squares is presented . section [ sec:4 ] gives comparison with existing constructions . in section [ sec:5 ] , an application to content based image retrieval is demonstrated . concluding remarks are given in section [ sec:7 ] . as stated already , cs refers to the problem of reconstruction of an unknown vector @xmath10 from linear measurements @xmath11 with @xmath12 being the inner - product between @xmath13 and @xmath14 . the basic objective in cs is to design a recovery procedure based on the sparsity assumption on @xmath13 when the number of measurements @xmath0 is very small compared to @xmath15 . sparse representations have merit for various applications @xcite in areas such as image / signal processing and numerical computation .    a vector @xmath10 is @xmath16sparse if it has at most @xmath17 nonzero coordinates . let @xmath18 stand for @xmath19 . the problem of obtaining the sparse vector from its linear measurements may be posed as @xmath20 this @xmath21minimization problem is np - hard @xcite in general . there are greedy algorithms for solving @xmath22 problem and orthogonal matching pursuit ( omp ) is one of the popular methods @xcite . several researchers ( for example , d. donoho @xcite and e. candes @xcite ) have established the conditions that ensure the recovery of solution to @xmath22 from : @xmath23 here @xmath24 denotes the @xmath25-norm of the vector @xmath26 . denote the solution to @xmath27 by @xmath28 and solution to @xmath29 by @xmath30 . the coherence @xmath31 of a given matrix @xmath32 is the largest absolute inner - product between different normalized columns of @xmath32 . denoting the @xmath17-th column in @xmath32 by @xmath33 , one defines @xmath31 , the coherence , as @xmath34 for a matrix of size @xmath35 , @xmath4 satisfies : @xmath36 , called the welch bound . it is known @xcite that for @xmath4-coherent matrices @xmath32 , one has @xmath37 , provided @xmath13 is @xmath16sparse with @xmath38    the restricted isometry property ( rip ) plays an important role @xcite in cs as it establishes the equivalence between the @xmath22 and @xmath39 problems . an @xmath35 matrix @xmath32 ( which we refer to a cs matrix ) is said to satisfy the restricted isometry property(rip ) of order @xmath17 with constant @xmath40 if for all @xmath16sparse vectors @xmath41 , we have @xmath42 it is known @xcite@xcite that the rip along with some conditions on @xmath43 imply the equivalence between the @xmath22 and @xmath39 problems . one of the important problems in cs theory deals with constructing cs matrices that satisfy rip for the largest possible range of @xmath17 . it is known that the random constructions satisfy the rip for the largest possible range on sparsity , which is @xmath44 @xcite .      to the best of our knowledge , designing good deterministic constructions of rip matrices is still a very interesting problem . r. devore @xcite has constructed deterministic binary sensing matrices of size @xmath45 with coherence @xmath46 , where @xmath2 is a prime power and @xmath47 . later on , s. li , f. gao et . al . @xcite have generalized the work in @xcite , constructing binary sensing matrices of size @xmath48 , where @xmath49 is any prime power and @xmath50 is the set of all rational points on algebraic curve @xmath51 over finite field @xmath52 . p. indyk @xcite has constructed binary sensing matrices using hash functions and extractor graphs with sizes @xmath53 , where @xmath54 . a. amini et . @xcite have constructed binary sensing matrices using ooc codes . in all these constructions , cs matrices are given for specific sets of row sizes . j. bourgain et . @xcite have constructed rip matrices of order @xmath55 , for some @xmath56 and @xmath57 using additive combinatorics . it is remarkable that this construction overcomes the natural barrier @xmath58 for those based on coherence . j. l. nelson et . al @xcite have given lower and upper bounds on maximal possible column size in terms of fixed row size @xmath0 and coherence @xmath4 . they have given a class of deterministic cs matrices of size @xmath59 with coherence @xmath60 , where @xmath2 is a prime number . r. calderbank et.al . @xcite have constructed cs matrices of size @xmath61 with coherence @xmath62 , using the delsarte - goethals codes , where @xmath63 is an odd number and @xmath64 is a constant integer . m. a. herman et.al . @xcite have constructed gabor frame of size @xmath65 with coherence @xmath66 , using alltop sequence , where @xmath67 is any prime number . the constructions , possessing sparsity @xmath68 have been given by calderbank et . al . in @xcite based on tight frames such as @xmath69 chirp matrices and alltop gabor frames @xcite , here @xmath0 is prime or prime power . the authors in @xcite have constructed matrices which achieve the square root bottleneck asymptotically .    in the present work , binary sensing matrices are constructed using euler squares . in particular , given any positive integer @xmath70 ( for prime @xmath2 ) , the procedure allows to construct binary cs matrix of size @xmath35 , where @xmath15 is @xmath71 , with @xmath4 being the coherence parameter of the matrix ( which is dependent on @xmath0 as explained in sub - sequent sections ) and @xmath5 . as the proposed methodology does not involve any function evaluations and gives matrices of small density , this method involves less complexity as compared to some of the existing ideas stated above . the applicability of construction is demonstrated through an application to cbir . an euler square of order @xmath72 , degree @xmath17 and index @xmath73 is a square array of @xmath74 @xmath16ads of numbers , @xmath75 , where @xmath76 and @xmath77 for @xmath78 and @xmath79 for @xmath80 and @xmath81    harris f. macneish @xcite has constructed euler squares for the following cases : + 1 . index @xmath82 , where @xmath2 is a prime number . index @xmath83 , for @xmath2 prime . index @xmath73 , where @xmath84 for distinct odd primes @xmath85 . here , @xmath86 equals the least of the numbers @xmath87 + for example , the euler square of index @xmath88 is as follows :    @xmath89 @xmath90 @xmath91    @xcite[lem : e ] let @xmath92 . then the existence of the euler square of index @xmath73 implies that the euler square of index @xmath93 exists .      using the concept of euler square of index @xmath73 , deterministic binary cs matrices @xmath32 that possess small coherence can be constructed . in particular , for @xmath94 a cs matrix of size @xmath95 , which may be treated as a block matrix consisting of @xmath17 number of @xmath96 blocks can be obtained . each column in @xmath32 corresponds to a @xmath17-ad in the euler square . the matrix @xmath32 is defined as follows : for @xmath97 , @xmath98 where @xmath99 is the @xmath100 @xmath16ad , @xmath101 is @xmath102 element in @xmath100 @xmath16ad and @xmath103 is the largest integer not greater than @xmath104 . to summarize , for @xmath105 the @xmath100 column @xmath14 of @xmath32 is generated as an @xmath106 binary vector from the @xmath100 @xmath17-ad @xmath99 with 1 occurring at the positions @xmath107 for @xmath108 . there are exactly @xmath17 ones in each column of @xmath32 and size of the matrix is @xmath109 . for example , the @xmath110 matrix given by the euler square of index @xmath88 is as follows :    @xmath111    the following lemma finds the bound on the coherence of @xmath32 . +    [ thm : lem ] the coherence of @xmath32 is at most equal to @xmath112 . suppose there exist two columns @xmath113 of @xmath32 such that the overlap between them is at least 2 . that is , if @xmath114 are the support vectors of @xmath113 , then @xmath115 , which implies that there exist two @xmath16ads @xmath116 with @xmath117 and @xmath118 such that @xmath119 * case 1 : * suppose @xmath120 and @xmath121 , from ( [ eqn : ci ] ) , @xmath122 and @xmath123 for @xmath121 , which is a contradiction from the definition of euler square . + * case 2 : * suppose @xmath80 and @xmath124 . from ( [ eqn : ci ] ) , @xmath125 and @xmath126 for @xmath80 , which is again a contradiction . + * case 3 : * suppose @xmath80 and @xmath121 , from ( 3 ) , @xmath127 and @xmath128 . since , @xmath129 for all @xmath130 , we have @xmath131 and @xmath132 , which implies that @xmath133 for @xmath80 and @xmath121 , which is a contradiction from the definition of euler square . hence there are no two such @xmath16ads satisfying ( [ eqn : ci ] ) . therefore , the overlap between any two columns of @xmath32 is at most 1 . since each column in @xmath32 contains fixed number of @xmath17 ones , it follows that the coherence of @xmath32 , @xmath31 , is at most equal to @xmath112 .    * remark 1 : * the maximum possible column size of any binary matrix is @xmath134 @xcite , where @xmath0 is the row size , @xmath17 is the number of ones in each column and @xmath135 is the overlap between any two columns . euler square of index @xmath73 results in a binary sensing matrix of size @xmath95 . in this matrix each column contains @xmath17 number of ones and coherence is at most @xmath112 . the maximum possible column size is thus @xmath136 for the euler square based matrix @xmath137 and @xmath138 , and hence @xmath139 , which is in the order of maximum possible column size . hence the aspect ratio is also in the maximum possible order . the density for these matrices is @xmath140 . the following proposition @xcite relates the rip constant @xmath141 and @xmath142    [ thm : pro ] suppose that @xmath143 are the unit norm columns of the matrix @xmath32 possessing the coherence @xmath4 . then @xmath32 satisfies the rip of order @xmath144 with constant @xmath145 .    from lemma [ thm : lem ] and proposition [ thm : pro ] , it follows that the matrix @xmath32 so constructed satisfies rip . [ thm : main ] the matrix @xmath146 satisfies rip with @xmath147 for any @xmath148 . * remark 2 : * @xmath149 and @xmath148 gives @xmath150 which implies that @xmath151 . consequently @xmath152    the previous procedure can be generalized for any row size @xmath0 different from a prime and its square . the following theorem summarizes the main result . [ thm : main1 ] suppose @xmath0 is any positive integer different from @xmath1 for a prime @xmath2 . then there exists a binary sensing matrix of size @xmath35 with coherence @xmath153 . * case-1 : * if @xmath154 then @xmath0 can be written as @xmath155 . since @xmath156 , we have @xmath157 . from @xcite , it is known that the euler square of index @xmath158 exists . since @xmath157 , from lemma  [ lem : e ] euler square of index @xmath159 exists . if we construct the binary matrix using this euler square , then the row size of it becomes @xmath160 . + * case-2 : * let @xmath0 be any integer such that @xmath161 for @xmath162 . from the fundamental theorem of arithmetic @xmath0 is factorized as @xmath163 . let @xmath164 and @xmath165 . with out loss of generality assume that @xmath166 ( the following arguments hold true even if @xmath167 for some @xmath168 ) . from the construction of macneish @xcite , it is known that the euler square of index @xmath169 exists . let @xmath170 and without loss of generality assume that @xmath171 . since @xmath172 , euler square of index @xmath173 exists , which in view of lemma  [ lem : e ] implies that the euler square of index @xmath174 exists . the binary cs matrix constructed through this euler square has a row size of @xmath175 , which is @xmath0 .    * remark 3 : * in both the cases above @xmath176 , which is in the order of maximum possible column size as stated in remark-1 and they allow for successful signal recovery of sparsity @xmath177 . +    the afore - described construction methodology results in a binary cs matrix of size @xmath95 from an euler square of index @xmath73 . as explained in remark 1 , the maximum possible column size for our construction in terms of row size @xmath178 and coherence @xmath112 is @xmath179 . since , @xmath180 , we have @xmath181 and @xmath182 . in theorem 5 , we give the construction for @xmath183 and in the next section , when @xmath72 is of the form @xmath184 ( that is , @xmath72 is not a power of single prime ) we try to increase the column size of above construction to @xmath185 , where @xmath186 , leaving the row size and the coherence intact . the procedure discussed above uses smallest factor , while in increasing the column size , other factors are also used .      as before , let @xmath84 , @xmath144= min@xmath187 and @xmath165 . suppose @xmath188 is the matrix of size @xmath95 that we obtain from the euler square of index @xmath73 . without loss of generality , assume that @xmath189 , let @xmath190 max@xmath191 . clearly , @xmath192 . since the euler square of index @xmath193 exists and using it we construct @xmath194 matrix , say @xmath195 , with coherence at most @xmath112 . now using these two matrices @xmath188 and @xmath195 , a new matrix @xmath196 is obtained as detailed below :    the matrix @xmath188 may be viewed as a block matrix possessing @xmath17 blocks with each block being of size @xmath96 . since @xmath197 and @xmath198 , we have @xmath199 the suitably zero - padded @xmath195 may be added as an additional column block to @xmath188 to generate a matrix of size @xmath200 . one may obtain @xmath17 number of zero - padded column blocks by placing @xmath195 in the corresponding row - block - locations of @xmath188 . using all these column blocks , we obtain a new binary matrix of size @xmath201 , say @xmath196 . from lemma [ thm : lem ] , the inner - product between any two columns in @xmath188 and @xmath196 is at most 1 . further the inner - product between any two columns one from @xmath188 and other from @xmath196 is at most 1 , this is because the @xmath17 nonzero entries of a column from @xmath196 fall only in one of the blocks of @xmath188 which has only one nonzero entry . the elements of the matrix @xmath196 may be generated as follows : for @xmath202 , also for @xmath203 , @xmath204 and for @xmath205 , where @xmath99 is the @xmath100 @xmath16ad in the euler square of index @xmath193 .    similarly , using @xmath206 max@xmath207 , @xmath208 , a cs matrix of size @xmath209 ( say @xmath210 ) is obtained . the matrix @xmath195 may be viewed as a block matrix possessing @xmath17 blocks with each block being of size @xmath211 . as shown in fig . 1 , the matrix @xmath196 contains @xmath17 such @xmath195 matrices , therefore there are @xmath212 blocks each of size @xmath211 in @xmath196 . since @xmath213 and @xmath214 , we have @xmath215 the suitably zero - padded @xmath210 may be added in these @xmath212 blocks as additional column blocks to @xmath196 as shown in fig . 1 and generate a new matrix , say @xmath216 , of size @xmath217 . continuing this way till euler square of index @xmath218 , a matrix , say @xmath219 , of size @xmath220 is obtained at the @xmath102 stage . finally , concatenating @xmath221 , results in a matrix with coherence at most @xmath112 and column size @xmath222 $ ] . this methodology is explained diagrammatically in figure [ fig : block ] for @xmath223 , that is , @xmath224 . the column extended matrix has rip compliance as summarized by the following theorem . [ thm : ma ] the matrix @xmath225 $ ] satisfies the rip with @xmath147 for any @xmath148 .    * remark 4 : * combining the results in sections [ subsec : eu ] and [ subsec : col ] , it can be concluded that the constructed matrix attains column size @xmath226 for some @xmath227 in particular , when @xmath0 is power of a single prime ( that is , @xmath228 ) , @xmath229 , and in other cases @xmath230 can be made to lie in @xmath231 , as justified in section [ subsec : col ] . as discussed already in the section [ subsec:1 ] , existing methods provide binary cs matrices for specific set of integers . the present work , however , provides a procedure for constructing binary matrices for a large class of row sizes . from remark 1 , the maximum possible column size is between @xmath74 and @xmath232 . since @xmath233 , it follows that @xmath234 and hence @xmath235 . in the constructions provided by @xcite , @xmath235 holds , albeit for specific set of values of @xmath0 . while in the construction of @xcite , @xmath236 holds . it is known @xcite that the welch bound is not sharp in this case . consequently the gain obtained in terms of increased column size does affect the coherence , which inturn restricts the sparsity of the solution to be recovered . in addition to providing general row size , the present construction is simple in the sense that it does not involve function evaluations like in @xcite and gives matrices with small density , which support algorithms with low computational complexity . for example , to generate an euler square matrix of size @xmath82 , it is only required to store two cyclic permutations of length @xmath2 and @xmath237 respectively and for index @xmath238 , it is sufficient to store at most @xmath239 permutations @xcite .    to the best of our knowledge , the constructions possessing sparsity @xmath240 ( that is , coherence @xmath241 ) exist for non - binary matrices with row size @xmath0 being prime or prime power @xcite . while in our construction methodology , if we use the euler squares of index @xmath242 for @xmath243 , then the size of the matrix @xmath32 is @xmath244 and coherence @xmath4 is @xmath245 . since the row size @xmath0 of @xmath32 is @xmath246 , we have @xmath247 . hence , this matrix provides guarantees for signals of sparsity up to @xmath248 . this is true for any euler square of index @xmath249 , where @xmath2 is a prime number , @xmath250 and @xmath251 hence we construct the binary matrices that provide guarantees for signals of sparsity up to @xmath252 for different class of row sizes such as @xmath253 , where @xmath2 is prime @xmath250 and @xmath251    for an arbitrary binary matrix , if the inner - product between any two columns is at most 1 , every column contains fixed ( @xmath254 ) number of ones ( that is coherence is at most @xmath255 ) and row size is @xmath0 , then the maximum possible column size @xmath256 as mentioned earlier .    using euler squares , we have constructed matrices that provide guarantees for signals of sparsity up to @xmath248 with @xmath257 . through an extension of our construction methodology , it is possible to generate ternary matrices , for which @xmath258 holds ( that is , coherence is in the order of @xmath259 ) with column size @xmath260 , which is explained in the following :    suppose @xmath32 is a matrix constructed from the euler square of index @xmath261 for @xmath262 1 or 2 and @xmath263 a hadamard matrix of size @xmath264 or @xmath265 . to get a ternary cs matrix ( i.e. matrix containing @xmath266 as elements ) @xmath267 , for each column of @xmath32 , we replace each of its @xmath268-valued entries with a distinct row of @xmath263 and @xmath269-valued entries with a zero row of same size . the size of the matrix @xmath267 becomes @xmath270 . in this construction , it is very easy to check that the coherence of the matrix @xmath267 is @xmath271 , which is in the order of @xmath259 ( welch bound ) with column size @xmath260 .    with a view towards comparing the numerical performance against the standard gaussian ( with entries drawn from @xmath272 and bernoulli ( with entries @xmath273 , each with probability @xmath274 ) random matrices , binary matrices of size @xmath275 and @xmath276 are generated using the euler squares of indices @xmath277 and @xmath278 . the omp algorithm is used to solve the @xmath279 minimization problem for a signal @xmath104 . let @xmath280 denote the recovered solution . from the reconstruction , the signal - to - noise ratio ( snr ) of @xmath104 is computed using @xmath281 for recovery at each sparsity level , 1000 input @xmath16sparse signals @xmath104 ( the nonzero indices are chosen uniformly randomly and the nonzero entries are drawn with @xmath282 ) are considered . the recovery is considered to be good if @xmath283 . simulation results on these matrices indicate ( shown in figure [ fig : comp_matrix ] ) that the euler square matrix gives better performance than the gaussian and bernoulli random matrices for certain higher sparsity levels . while for other sparsity levels , these matrices give the same recovery performance . it is observed that the cs matrices constructed from the euler squares of indices @xmath284 for a prime @xmath2 and @xmath285 to @xmath237 , give relatively superior performance than the other choices of @xmath286 against their gaussian and bernoulli counterparts . ( top plot ) and ( b ) @xmath276 ( bottom plot ) . these plots indicate that the euler square based matrix shows superior performance for some sparsity levels , while for other levels all matrices result in the same performance . the @xmath104 and @xmath287 axes in both plots refer respectively to the sparsity level and the success rate ( in % terms ) . for matrices in ( a ) and ( b ) , the coherences 0.5226 , 0.34 ( for gaussian matrices ) 0.52 , 0.28 ( for bernoulli matrices ) and 0.2 , 0.1 ( for euler based matrices ) . the values of coherence of euler matrices are equal to their theoretical bound @xmath288.,title="fig:",width=377 ]   ( top plot ) and ( b ) @xmath276 ( bottom plot ) . these plots indicate that the euler square based matrix shows superior performance for some sparsity levels , while for other levels all matrices result in the same performance . the @xmath104 and @xmath287 axes in both plots refer respectively to the sparsity level and the success rate ( in % terms ) . for matrices in ( a ) and ( b ) , the coherences 0.5226 , 0.34 ( for gaussian matrices ) 0.52 , 0.28 ( for bernoulli matrices ) and 0.2 , 0.1 ( for euler based matrices ) . the values of coherence of euler matrices are equal to their theoretical bound @xmath288.,title="fig:",width=377 ]      the phase transition diagrams depict the largest @xmath17 ( with fixed @xmath0 and @xmath15 ) for faithfully recovering @xmath16sparse vectors via @xmath9-norm minimization . in figure [ fig : phase ] , the region above the curve is the one in which successful reconstruction is not possible , while the region below consists of points at which successful recovery is possible . again , recovery is considered to be successful if snr is greater than @xmath289 . given a set of points @xmath290 , we have generated phase transition by finding the largest sparsity @xmath17 such that the successful recovery of 90 percent is achieved by considering the average recovery over 1000 iterations at each point by using the matrix which gives that particular point . figure  [ fig : phase ] provides phase transition for euler , gaussian and bernoulli random matrices for different values of @xmath291 which were given by the matrices of size @xmath35 , where @xmath292 and m=121 . it may be inferred from figure  [ fig : phase ] that the euler square based recovery is very competitive .     and @xmath293 respectively . this plot is generated for @xmath294 and @xmath295,width=377 ]      the efficacy of euler square based matrix is demonstrated using image reconstruction from lower dimensional patches , where the patches are generated via the sensing matrices ( as explained in more detail in the next section and figure  [ fig : prop ] ) . since the cs theory relies on the concept of sparsity , the key property needed of the image is enough sparsity in its original or some transform domain .    as compressed sensing allows for the reconstruction of sparse vector @xmath104 from its linear measurements @xmath296 , we demonstrate reconstruction performance via a medical image . the reconstructions shown in figure  [ fig:1 ] correspond to different down - sampling factors , viz 2.6 and 1.6 . the associated reconstruction errors in term of snr are shown in table  [ fig : tab1 ] . here by down - sampling factor , we mean the ratio of original patch size to reduced patch size which is same as @xmath297 ( where , @xmath35 is the size of the matrix used for projecting data to the lower dimensional space ) . from figure  [ fig:1 ] and table  [ fig : tab1 ] , it may be concluded that the euler square based matrices provide better reconstruction performance . , the images in ( b ) and ( d ) are those reconstructed via the euler square based matrices with down - sampling factors 2.6 and 1.6 respectively . the images in ( c ) and ( e ) are those obtained via the corresponding gaussian matrices . this figure states that euler square based cs matrices provide competitive reconstruction performance . the corresponding reconstruction errors are reported in table [ fig : tab1].,width=264,height=132 ]    , the images in ( b ) and ( d ) are those reconstructed via the euler square based matrices with down - sampling factors 2.6 and 1.6 respectively . the images in ( c ) and ( e ) are those obtained via the corresponding gaussian matrices . this figure states that euler square based cs matrices provide competitive reconstruction performance . the corresponding reconstruction errors are reported in table [ fig : tab1].,title="fig:",height=132 ] , the images in ( b ) and ( d ) are those reconstructed via the euler square based matrices with down - sampling factors 2.6 and 1.6 respectively . the images in ( c ) and ( e ) are those obtained via the corresponding gaussian matrices . this figure states that euler square based cs matrices provide competitive reconstruction performance . the corresponding reconstruction errors are reported in table [ fig : tab1].,title="fig:",height=132 ]    .a comparative error analysis of reconstruction by euler based and gaussian matrices for different down - sampling factors @xmath298 4 , 2.6 , 2 , 1.6 . the average error over 1000 iterations is reported for gaussian matrices . [ cols="^,^,^ " , ]     the classification performances of different classes by binary and gaussian matrices are shown in figures  [ fig : cone ] and [ fig : cong ] in terms of a confusion matrix . the confusion matrix gives the accuracy of the classification results . the diagonal elements of the confusion matrix indicate if each class is classified correctly . the corresponding nonzero column entries indicate the presence of mis - classification by the cbir scheme . from table [ tab:2 ] , it may be concluded that the binary based retrieval gives the performance comparable to that of its gaussian counter part . the retrieval performances by both matrices may be improved by fine - tuning the associated parameters ( like patch size ) and by taking different sparsifying basis . since the main focus of the present work is to construct cs matrices of general size , we do not go into the details of improvement of performance of cbir and the reconstruction of database member from the reduced dimensional vectors any further . in this paper , euler squares of index @xmath299 are used to construct binary sensing matrices of size @xmath109 . using these and with the help of prime factorization , matrices of general row size with asymptotically optimal column size and good coherence have been constructed . further , using block wise extension , the column size of the constructed matrices has been enhanced . simulation results on two test cases show that the generated matrices capture the support of the unknown vector and give better performance than gaussian matrices . it has also been demonstrated that the proposed binary sensing matrices exploit the pattern of inherent sparsity present in natural images and project data into lower dimensional spaces in such a way that the reduced vectors are useful for the purpose of cbir . the first author gratefully acknowledges the support ( ref no . 20 - 6/2009(i)eu - iv ) that he receives from ugc , govt of india . the third author is grateful to csir , gov . of india ( no . 25(219)/13/emr - ii ) for the financial support that he received . authors are thankful to dr . s. jana , dr . m. srinivas , mr . sandip chandra and mr . pradip sasmal for helpful discussions . 99 c. b. akgul , d. l. rubin , s. napel , c. f. beaulieu , h. greenspan , and b. acar ,  content based image retrieval in radiology : current status and future directions , " j. digital imag.vol . 208 - 222 , 2010 . l. applebauma , s. d. howardb , s. searlec and r. calderbank ,  chirp sensing codes : deterministic compressed sensing measurements for fast recovery , " appl . , 283 - 290 , mar . 2009 . w. u. bajwa , r. calderbank , and s. jafarpour ,  why gabor frames ? two fundamental measures of coherence and their role in model selection , " j. commun . 289 - 307 , 2010 . r. calderbank , howard s , jafarpour s.  construction of a large class of deterministic sensing matrices that satisfy a statistical isometry property , " ieee journal on selected topics in signal processing , 358 - 374 , 2010 .        y. c. chen , c. s. sastry , v. m. patel , p. j. phillips , r. chellappa . ,  in - plane rotation and scale invariant clustering using dictionaries , " ieee tans . image processing , vol . 22 , no . 6 , pp.2166 - 2180 , 2013 . a. gilbert and p. indyk ,  sparse recovery using sparse matrices , " proc . ieee , vol . 98 , no . 6 , pp . 937 - 947 , 2010 . m. a. herman and t. strohmer ,  high - resolution radar via compressed sensing , " ieee transactions on signal processing , vol . 2275 - 2284 , 2009 . p. indyk ,  explicit constructions for compressed sensing of sparse signals , " in proc . 19th annu . acm - siam symp . algorithms , 30 - 33 , 2008 . t.m , guld . o ,  deselaers . thomas , keysers . daniel ,  schubert . henning ,  spitzer . klaus ,  ney . hermann ,  wein.b.b , `` automatic categorization of medical images for content - based retrieval and data mining , '' computerized medical imaging and graphics . , vol.9 , no . 2 , pp.143 - 155 , 2005 . s. li , f. gao , g. ge , and s. zhang ,  deterministic construction of compressed sensing matrices via algebraic curves , " ieee trans . theory , vol . 58 , 5035 - 5041 , 2012 . s. li and g. ge ,  deterministic sensing matrices arising from near orthogonal systems , " ieee trans . theory , vol . 2291 - 2302 , 2014 . li zhuo , bo cheng and jing zhang ,  a comparative study of dimensionality reduction methods for large - scale image retrieval , " neurocomputing , pp . 202 - 210 , 2014 . r. prasad , c. r. murthy , and b. d. rao ,  joint approximately sparse channel estimation and data detection in ofdm systems using sparse bayesian learning , " ieee transactions on sig . proc . 3591 - 3603 , 2014 . shai shalev - shwartz and shai ben - david ,  understanding machine learning : from theory to algorithms , " cambridge university press , 2014 . n. y. yu and y. li ,  deterministic construction of fourier - based compressed sensing matrices using an almost difference set , " eurasip journal on advances in signal processing , vol . 2013 , no . 1 , pp . 1 - 14 , 2013 .
Abstract:

Answer: in compressed sensing the matrices that satisfy the restricted isometry property ( rip ) play an important role . but to date , very few results for designing such matrices are available . for applications such as multiplier - less data compression , binary sensing matrices are of interest . the present work constructs deterministic and binary sensing matrices using euler squares . in particular , given a positive integer @xmath0 different from @xmath1 for a prime @xmath2 , we show that it is possible to construct a binary sensing matrix of size @xmath3 , where @xmath4 is the coherence parameter of the matrix and @xmath5 . the matrices that we construct have small density ( that is , percentage of nonzero entries in the matrix is small ) with no function evaluation in their construction , which support algorithms with low computational complexity . through experimental work , we show that our binary sensing matrices can be used for such applications as content based image retrieval . our simulation results demonstrate that the euler square based cs matrices give better performance than their gaussian counterparts . [ theorem]definition [ theorem]remark [ theorem]example [ theorem]lemma [ theorem]proposition [ theorem]corollary    compressed sensing , coherence , rip , binary sensing matrices , euler squares , cbir .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: wilson lines ( also known as gauge links or eikonal lines ) can be naturally introduced in any gauge field theory . these objects are generically defined via traces of path - ordered exponentials of a gauge field evaluated along a given trajectory @xmath0}\ dz^{\mu } \ { \cal a}_{\mu}(z ) \ \right ]    $ ] . the path @xmath1 is a curve along which the gauge field @xmath2 gets transported from the initial point to the final one . wilson lines defined on closed contours are called wilson loops . they are path - dependent non - local functionals of the gauge field , invariant under gauge group transformations . putting the matter of question more mathematical , one can construct a space with its elements being wilson loops defined on an infinite set of contours . reformulation of qcd in terms of the elements of a generic loop space would allow one to use gauge - invariant quantities as fundamental degrees of freedom instead of the quarks and gluons from the standard qcd lagrangian @xcite . observables can then be obtained via correlation functions of wilson loops : @xmath3 \ . \label{eq : wl_def}\ ] ] complete information on the quantum dynamical properties of the loop space is accumulated in the schwinger - dyson equations : @xmath4 where @xmath5 stands for an arbitrary functional of the gauge fields . let the functionals @xmath5 be the wilson exponentials @xmath6 ( [ eq : wl_def ] ) . then eqs . ( [ eq : sch_dy_ym ] ) turn into the makeenko - migdal ( mm ) equations @xcite : @xmath7 where the basic operations are the area- @xmath8 and the path- @xmath9 derivatives @xcite : @xmath10 and the contour @xmath11 is obtained from the initial one by means of the infinitesimal area deformation @xmath12 at the point @xmath13 , while the path variation without changing the area gives rise to the path derivative @xmath14 the area derivative can be written as well in the so - called polyakov form  see , e.g. , @xcite for a discussion of an alternative approach . note that the derivation of the mm equations from the schwinger - dyson equations is grounded on the mandelstam formula @xmath15\ ] ] and/or on the stokes theorem , so that the wilson functionals which do not satisfy the corresponding restrictions ( such as , e.g. , cusped light - like loops ) apparently can not be straightforwardly treated within the same scheme . there are several other issues limiting the predictive power of the mm equations . namely , there exists an interesting class of wilson loops which possess very specific singularities originating , in particular , from the cusps and/or self - intersections of the contours and , in addition , from the light - like segments of the integration paths . the simplest example is given by a wilson exponential evaluated along a cusped contour with two semi - infinite light - like sides , fig . [ fig:0 ] . @xmath16{cusp_1loop}\ ] ]    already the leading order contribution to this wilson exponential possesses all the peculiar singularities : the pure ultraviolet , the infrared ( due to the infinite lengths of the sides ) , and the light - like cups divergences . this simple contour will arise in what follows as a building unit of many important wilson loops and correlation functions . physically it corresponds to the soft part of the factorized quark form factor , which has been studied in detail in @xcite .    in the present work we propose and discuss a new approach to these issues , having in mind , as an instructive example , a very special type of wilson loops  planar rectangles with light - like sides . considerable interest to cusped light - like wilson polygons has arisen thanks to the recently conjectured duality between the @xmath17gluon planar scattering amplitudes in the @xmath18 super - yang - mills theory and the vacuum average of planar wilson loops formed , correspondingly , by @xmath19 light - like segments connecting space - time points @xmath20 , so that their `` lengths '' @xmath21 are chosen equal to the external momenta of the @xmath17gluon amplitude ( see , e.g. , @xcite and references therein ) . it has been demonstrated that the infrared singularities of the former corresponds to the ultraviolet singularities of the latter , and the cusp anomalous dimension is the crucial constituent of the evolution equations @xcite . wilson exponentials possessing light - like segments ( or that are fully light - like ) have been studied also in a different context @xcite . the main observation is that the renormalization properties of these wilson loops are more intricate than those of cusped wilson loops defined on off - light - cone integration contours . namely , the light - cone cusped wilson loops are not multiplicatively renormalizable because of the additional light - cone singularities ( besides the standard ultraviolet and infrared ones ) . it is possible , however , to construct a combined renormalization - group equation taking into account ultraviolet as well as light - cone divergences . the cusp anomalous dimension , which is the principal ingredient of this equation , is remarkably universal : it controls , e.g. , the infrared asymptotic behavior of such important quantities as the qcd and qed sudakov form factors , the gluon regge trajectory , the integrated ( collinear ) parton distribution functions at large-@xmath13 , the anomalous dimension of the heavy quark effective theory , etc . @xcite . another interesting field of application of cusped light - cone wilson lines could be transverse - momentum dependent parton densities ( tmds ) @xcite . the latter are introduced to describe the intrinsic transverse momentum of partons inside the nucleon , which is needed in the study of semi - inclusive processes within the ( generalization of ) the qcd factorization formalism @xcite . let us discuss the emergent singularities arising in tmds beyond tree - approximation . at one - loop level , the following three classes of divergences appear : @xmath22 standard ultraviolet poles , which are removable by a normal renormalization procedure ; @xmath23 pure rapidity divergences , which depend on an additional rapidity cutoff , but do not violate renormalizability of tmds ; they can be resummed by means of the collins - soper evolution equation ; @xmath24 very specific overlapping divergences : they contain the ultraviolet and rapidity poles simultaneously and thus break down the standard renormalizability of tmds . this situation resembles the problems with renormalizability of the light - like wilson loops discussed above . however , the structure of wilson lines is quite involved already in the tree - approximation . the most straightforward definition of `` a quark in a quark '' tmd , which meets the requirement of the parton number interpretation , reads @xmath25 with @xmath26 . here we define the semi - infinite wilson lines evaluated along a four - vector @xmath27 as @xmath28 \ , \ ] ] where the vector @xmath27 can be light - like @xmath29 , or transverse @xmath30 . formally , the integration of ( [ eq : general ] ) over @xmath31 is expected to give the collinear ( also called integrated ) pdf @xmath32 however , this is only justified in tree approximation . it is worth noting that the normalization of the above tmd @xmath33 can be most easily obtained by making use of the canonical quantization procedure in the light - cone gauge , where longitudinal wilson lines become equal to unity and where equal - time commutation relations for creation and annihilation operators @xmath34 immediately yield the parton number interpretation @xmath35 the usage of `` tilted '' gauge links in the operator definition of tmds does not meet this requirement . we visualize the geometrical layout of various wilson lines in the operator definition of tmds in fig . [ fig : geo1 ] , [ fig : geo2 ] , [ fig : geo3 ] and discuss relevant issues in their captions . geometry of the contours in unsubtracted tmds with light - like ( upper panel ) and off - light - cone ( lower panel ) longitudinal wilson lines and their reduction to integrated pdfs in tree approximation . in the former case , the transverse wilson lines vanish after @xmath31-integration , while the longitudinal wilson lines turn into an one - dimensional connector @xmath36 . in the off - light - cone schemes , the mutual compensation of transverse wilson lines at infinity is not visible . moreover , the integrated configuration contains two non - vanishing off - light - cone wilson lines , which apparently are not equivalent to the collinear connector @xmath36 . the interrogation marks next to the transverse wilson lines symbolize the lacking of any consistent treatment in tmd formulations with off - light - cone ( shifted ) wilson lines . in contrast , the transverse wilson lines appear naturally in `` light - cone '' schemes.,scaledwidth=45.0% ]    comparative layout of wilson lines in unsubtracted soft factors and visualization of the reduction to the collinear case . the upper panel shows the soft factor in momentum space , as proposed in refs . the lower panel presents the tilted off - light - cone integration paths in impact parameter space , as well as the result of the reduction to the collinear @xmath37 configuration.,scaledwidth=45.0% ]    comparative layout of wilson lines in subtracted soft factors . the upper panel corresponds to the soft factor of the tmd distribution function which enters the factorization with pure light - like wilson lines . the lower panel has the same setup , but with the longitudinal wilson lines shifted off the light - cone.,scaledwidth=45.0% ]    beyond tree - approximation , the virtual diagrams producing terms with overlapping singularities are shown in fig . [ fig:1 ] . the typical extra divergency stems from the one - loop vertex - type graph fig . [ fig:1](a ) in covariant gauges or from the self - energy graph fig . [ fig:1](b ) in the light - cone gauge ( in the large-@xmath38 limit ) and reads @xmath39^{\epsilon }   \ \delta ( 1-x ) \delta^{(2 ) } ( \bm{k}_\perp ) \ { \int_0 ^ 1\ ! \frac{x^{1-\epsilon}}{(1-x)^{1+\epsilon } } } \ .\ ] ] the standard ultraviolet pole in the gamma - function @xmath40 is accompanied by an additional singularity in the integral . the latter is due to the integration over infinite gluon rapidity and can not be treated by dimensional regularization , calling for an extra ( rapidity ) cutoff . the reason for renormalizability violation in the leading order contribution to tmds is that light - like wilson lines ( or the `` standard '' quark self - energy in light - cone gauge ) produce more singular terms than usual green functions do . @xmath16{wl_lc_1loop}\ ] ]    to solve the problems with extra singularities and renormalizability in tmds , a variety of ( possibly non - equivalent ) methods has been proposed . working in the covariant feynman gauge , ji , ma and yuan proposed a scheme which utilizes tilted ( off - light - cone ) longitudinal wilson likes directed along the vector @xmath41 @xcite . transverse wilson lines at the light - cone infinity cancel in covariant gauges , while the rapidity cutoff @xmath42 marks the deviation of longitudinal wilson lines from a pure light - like direction . a subtracted soft factor then contains non - light - like wilson lines as well . obviously , such off - light - cone unsubtracted tmds with the light - like vector @xmath43 replaced by the vector @xmath44 do not obey the equation ( [ eq : u_to_i ] ) , not even at tree level . however , it is possible to formulate a `` secondary factorization '' method which allows one to express off - light - cone tmds ( in impact parameter space @xmath45 ) as a convolution of integrated pdfs and perturbative coefficient functions in the perturbative region ( that is , at small @xmath46 ) , see @xcite .    in publications @xcite it was proposed to explore the renormalization - group properties of unsubtracted tmds ( [ eq : general ] ) and to make use of their anomalous dimension as a tool to discover the minimal layout of wilson lines in the soft factor that provides a cancelation of overlapping dependent terms . it has been demonstrated ( in the leading @xmath47-order ) that the extra contribution to the anomalous dimension is exactly the cusp anomalous dimension @xcite , which is a crucial element of the investigation of non - renormalizible cusped light - like wilson loops . making use of specially chosen soft factors , one can get rid of the extra divergences in the operator definition of the tmds , however paying a price in the form of significant complication of the structure of the wilson lines in the above definition . in the present work we discuss another approach to the problems of light - cone cusped wilson loops @xcite . to this end , it appears instructive to study those properties shared by such apparently different quantities as tmds , light - like wilson polygons , etc . , which originate in their light - cone structure and arise in the form of the `` too singular '' non - renormalizable terms . we made use of the observation that in the large-@xmath38 limit , in the transverse null - plane , for the light - like planar dimensionally regularized ( not renormalized ) wilson rectangles , the area derivatives introduced in the previous sections can be reduced to the normal ones . the area variational equations in the coordinate representation describe the evolution of light - like wilson polygons and represent , therefore , the `` equations of motion '' in loop space , valid for a specific class of its elements . as a result , the obtained differential equations give us a closed set of dynamical equations for the loop functionals , and can in principle be solved in several interesting cases . let us start with the quantum dynamical principle proposed by schwinger @xcite : the action operator @xmath48 defines variations of arbitrary quantum states , so that @xmath49 the area variations ( [ eq : area_derivative ] ) of field exponentials @xmath50 yield @xmath51 where @xmath52 is yet to be defined . the loop space consists of scalar objects with different topological and geometrical features , hence the equations of motion in this space must be the laws which state how those objects change their shape . it means that `` motion '' in loop space is equivalent to variation of the shapes of integration contours in the wilson loops @xcite . therefore , we have to find the proper operator @xmath52 , which governs the shape variations of the light - like cusped loops ( wilson planar polygons ) . following the standard method , one makes use of eq . ( [ eq : principle ] ) in the form ( [ eq : sch_dy_ym ] ) and gets the system of the mm eqs . ( [ eq : mm_general ] ) . alternatively , we will try to get rid of the operations which tacitly assume the property of smoothness of the wilson loops of interest . consider , e.g. , a generic wilson loop @xmath53 without mentioning if it is smooth or not . the two leading terms of its perturbative series are given by @xmath54 where @xmath55 is the dimensionally regularized ( @xmath56 ) free gluon propagator @xmath57^{1- \epsilon } } \ . \label{eq : gluon_prop_dr}\ ] ] for convenience s sake , we work in the feynman covariant gauge and separate out the scalar part of the propagator @xmath58 . the issues related to gauge- and regularization independence of the calculations will be considered elsewhere . then , if the l.h.s . of the eq . ( [ eq : principle_mod ] ) acts on the wilson exponential ( [ eq : w_generic_pert ] ) , one gets @xmath59 using the stokes theorem ( let us assume for a moment that we are allowed to do so ) , we obtain @xmath60 where @xmath1 is considered as the boundary of the surface @xmath61 . we get then the leading perturbative term of the mm equation ( [ eq : mm_general ] ) : @xmath62 we must treat this result with due caution : in the course of the derivation , we assumed that the stokes theorem is valid for all wilson loops which we consider . however , the last statement is not true in general , that is why we mark the `` good '' contours with a special index @xmath63 . it is interesting that in @xmath64 qcd , the area differentiation turns into the ordinary derivative , by virtue that the gluon propagator ( [ eq : gluon_prop_dr ] ) for @xmath65 gets logarithmic in @xmath66 : @xmath67 \ , \ \sigma = \ { \rm area \ inside\ } \gamma_\bigcirc \ ,    \label{eq:2d_area}\ ] ] so that @xmath68 evaluating , in the same way , the nlo terms , one obtains the full mm eq . ( [ eq : mm_general ] ) . however , we interrupt here and return a bit backward , since we are mostly interested in those wilson functionals that do not satisfy ( or , at least , do not satisfy straightforwardly ) the conditions of the applicability of the stokes theorem . we will , therefore , continue with the study of the shape variations of wilson loops without relying upon the stokes theorem , but keeping in mind an explicit form of the free gluon propagator ( which possesses a specific light - cone / rapidity divergence ) , see eq . ( [ eq : gluon_prop_dr ] ) . we are now in a position to extend the schwinger approach to a more complicated case and to try to derive the corresponding area evolution equations . the calculation of cusped light - cone wilson loops beyond tree approximation in different gauges and the justification of gauge independence calls for a careful treatment of a variety of divergences already in leading order . special attention must by paid to the separation of the rapidity divergences and the standard ultraviolet ones @xcite . in the t hooft ( large-@xmath38 ) limit one obtains @xcite @xmath69^\epsilon + \[{2 n^+ n^-\mu^2 + i0}\]^\epsilon   \ ) \\ & &   + \frac{\alpha_s n_c}{2\pi } \ ( \frac{1}{2 } \ln^2 \frac{n^+n^-+ i0}{-n^+n^-+i0 } + { \rm finite\ terms }   \ ) + o(\alpha_s^2 ) \ , \nonumber \label{eq : wl_lc_1loop}\end{aligned}\ ] ] with the mandelstam variables in the momentum space , @xmath70 and @xmath71 , map onto the area variables in the coordinate space , so that @xmath72 . we will show separately that the result ( [ eq : wl_lc_1loop ] ) is not only gauge invariant , but is independent of any regularization of light - cone and ultraviolet divergences and of the way they are separated . the latter issue is of considerable importance to the study of the operator structure of transverse - momentum dependent parton densities , the jet and soft functions in the soft - collinear effective theory , the infrared properties of the high - energy scattering amplitudes , etc . ( see , e.g. , @xcite and references therein ) . the transverse null - plane is defined by the condition @xmath73 ; therefore , the shape variations ( which give rise to the infinitesimal area changes ) are defined as @xmath74 we assume that these operations do make sense only at the corner points @xmath20 , and introduce the `` left '' and `` right '' variations , as shown in fig . [ fig:2 ] . @xmath75{wl_lc_area}\ ] ]    rectangular planar wilson loop @xmath76 is known to lack multiplicative renormalizability @xcite . in order to decrease the power of singularity that violates the renormalizability , one can follow the scheme proposed in @xcite . having in mind eq . ( [ eq : delta_area ] ) , we define the area logarithmic derivative @xmath77 and act with this operator on the r.h.s . of the eq . ( [ eq : wl_lc_1loop ] ) : @xmath78^\epsilon + \[{2n^+n^-\mu^2 + i0}\]^\epsilon \ ) \ . \label{eq : log_der}\ ] ] then , the result is finite ( after additional logarithmic differentiation in the ultraviolet scale @xmath79 ) and is given by the cusp anomalous dimension @xmath80 therefore , we have obtained the result ( [ eq : full_der ] ) given that the infinitesimal shape variations are defined as in ( [ eq : area_derivative ] ) . eq.([eq : full_der ] ) describes then the dynamics of the cusped planar light - like wilson loops @xcite . we established , therefore , the connection between the geometry of the loop space ( in terms of the area / shape differentials ) and the dynamical properties of the fundamental degrees of freedom  the gauge- and regularization - invariant planar light - like wilson loops . the very possibility to obtain a finite result by means of eqs . ( [ eq : log_der ] , [ eq : full_der ] ) is a direct consequence of the geometrical properties of loop space , when considering the non - renormalizable cusped light - like wilson loops . to show this explicitly , we restrict ourselves to shape variations ( [ eq : delta_area ] ) , and apply the area derivative to the planar wilson rectangle @xmath81^{1 -\epsilon } } \ ,   \label{eq : wl_rect_area}\end{aligned}\ ] ] where the sides are parameterized as @xmath82 with the vectors @xmath83 having the dimension @xmath84 $ ] @xcite . a peculiar feature of the planar light - like contours is that the area dependence separates out from the integrals and can be calculated explicitly ( taking into account that @xmath85 , see eq . ( [ eq : delta_area ] ) ) @xmath86^{1 -\epsilon } }    \ . \label{eq : wl_rect_area_nn}\ ] ] moreover , light - like wilson lines with @xmath87 develop an extra singularity , which shows up in the form of a second - order pole @xmath88 , while the cusps violate conformal invariance of the wilson loop because the `` skewed '' scalar products @xmath89 replace the conformal ones @xmath90 . hence , performing the area @xmath91 and the uv - scale logarithmic differentiation of eq . ( [ eq : wl_rect_area_nn ] ) and summing up all relevant terms , we obtain @xmath92 }   =   - \sum { \gamma_{\rm cusp } } \ ,   \label{eq : mod_schwinger}\ ] ] which has been foreseen in eq . ( [ eq : full_der ] ) and which is derived now directly from the schwinger quantum dynamical principle . it is natural that this result is akin , in some sense , to the situation in @xmath64 qcd mentioned above . the area derivative becomes the ordinary derivative for the same reason : one has effectively planar wilson loops , thus the mm eqs . becomes a closed system and , in principle , solvable @xcite . let us emphasize that the r.h.s . of eq . ( [ eq : mod_schwinger ] ) is regulated by the cusp anomalous dimension , which is a universal quantity that is independent of the shape of particular contour and which is known perturbatively up to the @xmath93 order . it is therefore worth analyzing if the above result is only a leading order approximation , or if it also will be valid for higher order approximations . let us take into consideration the property of linearity of the ( angular - dependent ) cusp anomalous dimension in the infinitely large angle asymptotics with respect to the logarithm of the cusp angle @xmath94 @xcite : @xmath95 with the `` maximally non - abelian '' coefficients given by @xmath96 and @xmath97 are cusp - independent factors . this asymptotic regime coincides with the light - cone case with the angular - dependent logarithms being transformed into extra pole terms in @xmath98 : @xmath99 , see @xcite . more specifically , the area variable @xmath100 enters the regularized cusp anomalous dimension in the light - like limit as @xmath101 and , after logarithmic differentiation , one obtains a perturbative expansion of the cusp anomalous dimension @xmath102 which justifies the previous result ( [ eq : mod_schwinger ] ) in the nlo by virtue that @xmath103 this implies that the result ( [ eq : mod_schwinger ] ) is , in fact , an all - order one , akin to the mm eq . ( [ eq : mm_general ] ) : they both are exact and non - perturbative , while the r.h.ss of each one can be calculated perturbatively . it is worth noting that eq . ( [ eq : full_der ] ) agrees with the non - abelian exponentiation theorem for the dimensionally regularized wilson loops @xmath104 \ , \ ] ] with the summation going over all two - particle irreducible diagrams , whose contribution is given by the so - called `` web '' functions @xmath105 @xcite . thus , eq . ( [ eq : full_der ] ) can be used in calculation of the higher - order terms in the cusp anomalous dimension , by virtue that one has a closed recursive system of the perturbative equations . now we can apply the methodology described above to the tmds with the light - like longitudinal wilson lines @xmath106 , eq . ( [ eq : general ] ) , what yields @xmath107    =    2 { \gamma_{\rm cusp } } \ ,    \label{eq : tmd_combined}\ ] ] with the `` area '' being encoded in the rapidity cutoff parameter @xmath108 @xcite . another interesting example is provided by the @xmath109-shape semi - loop with one finite light - like segment @xcite . in the one - loop order , we have in the large-@xmath38 limit @xmath110 \ , \\    & & l(nn^- )    = \frac{1}{2}\(\ln ( \mu n n^- + i0 ) + \ln ( \mu nn^- + i0 ) \)^2 \ , \nonumber    \label{eq : pi_1loop }    \end{aligned}\ ] ] where the area is given by the product of the light - like @xmath111 and non - light - like @xmath112 vectors , see fig . [ fig:5 ] . @xmath16{wl_lc_pi_shape}\ ] ]    the @xmath109-shaped wilson loop ( [ eq : pi_1loop ] ) also satisfies eq . ( [ eq : mod_schwinger ] ) : @xmath113    =    - 2 { \gamma_{\rm cusp } } \ , \ ] ] the latter controls for the renormalization - group evolution of the integrated pdfs in the large-@xmath13 limit , as well as the anomalous dimensions of conformal operators with large lorentz spin @xcite . the @xmath109-shape contour can be split and moved apart to separate two planes by the transverse distance @xmath114 . the wilson loop obtained in such a way is expected to be `` dual '' to the tmd . this duality implies that both the double - planar @xmath109-shaped wilson loop , fig . [ fig:7 ] , and the tmd ( [ eq : general ] ) obey the same combined evolution equation ( [ eq : tmd_combined ] ) . the detailed analysis of this configuration will be presented elsewhere . @xmath75{tmd_dual}\ ] ] the universal quantum dynamical principle by schwinger provides a proper approach to the description of the dynamics of the loop space . the gauge invariant wilson loops are considered then as the only degrees of freedom , and the makeenko - migdal equations ( [ eq : mm_general ] ) stem from the schwinger - dyson equations applied to the renormalizable loops . in general case , the system of the mm eqs . is not closed and can not be immediately used in practical calculations in qcd . @xmath16{wl_lc_general}\ ] ]    in the present work we showed that it is possible to design a relevant system of equations of motion valid for the cusped planar light - like wilson loops , taking into account that the latter possess a very specific singularity structure compared to their off - light - cone analogues . general solution of this problem has not been found yet , but we have managed to demonstrate that some simplifications make it possible to propose a new potentially fruitful method to deal with such wilson exponentials . namely , in the large-@xmath38 limit , the planar rectangular light - like contours at @xmath115 enable us to reduce the area functional derivative to the normal derivative for dimensionally regularized wilson loops . as the result , the equations which describe the infinitesimal shape variations in coordinate space appear to be dual to the energy evolution equations for cusped wilson loops in the space of the mandelstam momentum variables . within the framework we proposed , the dynamics of elements of loop space is introduced by means of obstructions of the initially smooth wilson loops , which play , therefore , the role of _ sources _ in schwinger s `` fields and sources '' picture . we have argued , therefore , that the schwinger quantum dynamical principle can be used as an effective tool to study ( at least ) one important class of the elements of loop space  the cusped planar wilson polygons on the light - cone . we implemented the program only in one of the simplest situations , the planar rectangle . in fig . [ fig:3 ] , a more involving configuration is visualized , the arbitrary quadrilateral path , the area evolution of which is far from being trivial and deserves a separate study . we appreciate stimulating discussions and useful critical remarks by i.v . anikin , y.m . makeenko and n.g . we thank the participants of the jlab qcd evolution workshop , the summer meetings in ect * , trento , and the joint seminar at the university of antwerp for useful discussions . s. mandelstam , phys . rev .  175 ( 1968 ) 1580 ; k.g . wilson , phys . rev . * d10 * ( 1974 ) 2445 ; y. nambu , phys . lett . * b80 * ( 1979 ) 372 ; a.m. polyakov , phys . lett . * b82 * ( 1979 ) 247 . a.m. polyakov , nucl . * b164 * ( 1979 ) 171 ; v.s . dotsenko , s.n . vergeles , nucl . b169 * ( 1980 ) 527 ; i.ya . arefeva , phys . lett . * b93 * ( 1980 ) 347 ; n.s . craigie , h. dorn , nucl .  phys . * b185 * ( 1981 ) 204 ; s. aoyama , nucl .  phys . * b194 * ( 1982 ) 513 ; n.g . stefanis , nuovo cim . * a83 * ( 1984 ) 205 ; i.i . balitsky , v.m . braun , nucl . * b311 * ( 1989 ) 541 . makeenko , a.a . migdal , phys .  lett . * b88 * ( 1979 ) 135 , [ erratum - ibid . * b89 * ( 1980 ) 437 ] ; nucl .  phys . * b188 * ( 1981 ) 269 ; a.a . migdal , annals phys . * 109 * ( 1977 ) 365 ; annals phys . * 126 * ( 1980 ) 279 ; phys .  rept . * 102 * ( 1983 ) 199 . brandt , f. neri , m .- a . sato , phys . * d24 * ( 1981 ) 879 ; r.a . brandt , a. gocksch , m .- a . sato , f. neri , phys . * d26 * ( 1982 ) 3611 . a. kernemann , n.g . stefanis , phys .  rev . * d40 * ( 1989 ) 2103 ; n.g . stefanis , hep - th/9607063 ; g.c . gellas , a.i . karanikas , c.n . ktorides , n.g . stefanis , phys .  lett . * b412 * ( 1997 ) 95 ; a.i . karanikas , c.n . ktorides , jhep * 11 * ( 1999 ) 033 ; a.i . karanikas , c.n . ktorides , n.g . stefanis , eur . j. * c26 * ( 2003 ) 445 . korchemskaya , g.p . korchemsky , phys .  lett . * b287 * ( 1992 ) 169 ; a. bassetto , i.a . korchemskaya , g.p . korchemsky , g. nardelli , nucl .  phys . * b408 * ( 1993 ) 62 ; g.p . korchemsky , mod .  phys .  lett . * a4 * ( 1989 ) 1257 ; g.p . korchemsky , g. marchesini , phys .  lett . * b313 * ( 1993 ) 433 ; g.p . korchemsky , j.m . drummond , e. sokatchev , nucl . b795 * ( 2008 ) 385 . ivanov , g.p . korchemsky , a.v . radyushkin , yad .  fiz . * 44 * ( 1986 ) 230 ; g.p . korchemsky , a.v . radyushkin , phys .  lett . * b171 * ( 1986 ) 459 ; g.p . korchemsky , phys .  lett . * b217 * ( 1989 ) 330 ; phys .  lett . * b220 * ( 1989 ) 629 . alday , j.m . maldacena , jhep * 06 * ( 2007 ) 064 ; l.f . alday , b. eden , g.p . korchemsky , j. maldacena , e. sokatchev , jhep * 99 * ( 2011 ) 123 ; n. beisert _ et al . _ , lett .  math * 99 * ( 2012 ) 3 ; m. giordano , r. peschanski , s. seki , acta phys .  polon . * b43 * ( 2012 ) 1289 . korchemsky , a.v . radyushkin , nucl .  phys . * b283 * ( 1987 ) 342 . karanikas , c.n . ktorides , n.g . stefanis , phys .  rev . * d52 * ( 1995 ) 5898 ; g.c . gellas , a.i . karanikas , c.n . ktorides , n.g . stefanis , phys .  lett . * b412 * ( 1997 ) 95 ; a.i . karanikas , c.n . ktorides , n.g . stefanis , s.m.n . wong , phys . b455 * ( 1999 ) 291 . j.c . collins , d.e . soper , nucl . b193 * ( 1981 ) 381 ; nucl .  phys . * b194 * ( 1982 ) 445 ; j.c . collins , f. hautmann , phys . b472 * ( 2000 ) 129 ; f. hautmann phys .  lett . * b655 * ( 2007 ) 26 ; j.c . collins , acta phys . pol.*b34 * ( 2003 ) 3103 . cherednikov , n.g . stefanis , phys .  rev . * d77 * ( 2008 ) 094001 ; nucl . b802 * ( 2008 ) 146 ; phys .  rev . * d80 * ( 2009 ) 054008 ; n.g . stefanis , i.o . cherednikov , mod . * a24 * ( 2009 ) 2913 ; i.o . cherednikov , a.i . karanikas , n.g . stefanis , nucl .  phys . * b840 * ( 2010 ) 379 . collins , d.e . soper , g. sterman , adv . high energy phys . * 5 * ( 1988 ) 1 ; a.v . belitsky , a.v . radyushkin , phys .  rept . * 418 * ( 2005 ) 1 ; a.v . belitsky , x. ji , f. yuan , nucl . b656 * ( 2003 ) 165 ; d. boer , p.j . mulders , f. pijlman , nucl . b667 * ( 2003 ) 201 . x. ji , j. ma , f. yuan , phys .  rev . * d71 * ( 2005 ) 034005 . cherednikov , t. mertens , f.f . van der veken , phys . * d86 * ( 2012 ) 085035 ; int . j. mod * 20 * ( 2012 ) 109 . schwinger , phys .  rev . * 82 * ( 1951 ) 914 . collins , int . j.  mod . ser . * 4 * ( 2011 ) 85 ; i.o . cherednikov , n.g . stefanis , int . j.  mod . .  ser . * 4 * ( 2011 ) 135 ; j .- y . chiu , a. jain , d. neill , i.z . rothstein , jhep * 05 * ( 2012 ) 084 ; phys . * 108 * ( 2012 ) 151601 ; e. avsar , arxiv:1203.1916 [ hep - ph ] ; i.o . cherednikov , pos * qnp * ( 2012 ) 061 ; n.g . stefanis , arxiv:1211.7218 [ hep - ph ] ; j. collins , arxiv:1212.5974 [ hep - ph ] . a. idilbi , i. scimemi , phys . b. 2011 . v. 695 . p. 463 ; m.g . garca - echevarra , a. idilbi , i. scimemi , phys .  rev . * d84 * ( 2011 ) 011502(r ) ; jhep * 1207 * ( 2012 ) 002 ; y. li , s. mantry , f. petriello , phys .  rev . * d84 * ( 2011 ) 094014 ; s. mantry , f. petriello , int . j.  mod . phys .  conf . * 04 * ( 2011 ) 106 . gatheral , phys . b133 * ( 1993 ) 90 ; j. frenkel , j.c . taylor , nucl .  phys . * b246 * ( 1984 ) 231 ; a. mitov , g. sterman , i. sung , phys .  rev . * d82 * ( 2010 ) 096010 ; o.  erdogan , g.  sterman , arxiv:1112.4564 [ hep - th ] . korchemsky , g. marchesini , nucl .  phys . * b406 * ( 1993 ) 225 ; y. makeenko , jhep * 01 * ( 2003 ) 007 . n. drukker , d.j . gross , h. ooguri , phys . * d60 * ( 1999 ) 125006 ; y. makeenko , p. olesen , g.w . semenoff , nucl .  phys . * b748 * ( 2006 ) 170 .
Abstract:

Answer: we propose and discuss a new approach to the analysis of the correlation functions which contain light - like wilson lines or loops , the latter being cusped in addition . the objects of interest are therefore the light - like wilson null - polygons , the soft factors of the parton distribution and fragmentation functions , high - energy scattering amplitudes in the eikonal approximation , gravitational wilson lines , etc . our method is based on a generalization of the universal quantum dynamical principle by j. schwinger and allows one to take care of extra singularities emerging due to light - like or semi - light - like cusps . we show that such wilson loops obey a differential equation which connects the area variations and renormalization group behavior of those objects and discuss the possible relation between geometrical structure of the loop space and area evolution of the light - like cusped wilson loops .