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: the roadmap for fundamental physics in space issued by the european space agency ( esa ) in 2010  @xcite put an emphasis on gravitation tests in the solar system with missions to the outer planets . in this framework , it recommends the development of accelerometers compatible with spacecraft tracking at the 10 pm.s@xmath0 level at low frequencies . indeed , interplanetary probes can be used as test masses whose trajectories are to be compared to theoretical predictions . such a test has been performed by nasa with the pioneer 10 and 11 probes : the outcome was a discrepancy with respect to the predictions of general relativity @xcite . but these two probes lacked an instrument which could disentangle gravitational effects from non - gravitational ones . in order to improve the experiment made by the pioneer probes , several space missions have been proposed @xcite , many of them embarking an accelerometer . the oss mission @xcite relies , for its fundamental physics objectives , on the gravity advances package @xcite . it is an instrument composed of an electrostatic accelerometer called microstar , based on onera expertise in the field of accelerometry and gravimetry with champ , grace , and goce missions  @xcite , and a rotating platform , called bias rejection system  ( brs ) . in orbit technology is used with technological upgrades aiming at reducing power consumption , size and mass . this instrument aims at measuring the non - gravitational acceleration of the spacecraft with a precision compatible with esa requirement and with no bias . thus , it provides an additional observable and allows removing , during the orbit restitution process , the effect of the non - gravitational forces on the trajectory , enhancing deep space gravitation tests as well as gravity field recovery  @xcite . the measurements made by microstar or any other electrostatic accelerometer always include an intrinsic bias which may vary with time or temperature . it calls for an additional subsystem able to remove the bias through the modulation of the signal of interest . this is the role of the bias rejection system . in order to remove properly the bias of microstar , the bias rejection system rotates the accelerometer following a carefully designed periodical pattern  @xcite . in terms of measurement noise , this operation selects the noise of microstar around the modulation frequency . after post - processing and for a modulation period of @xmath1  min , it allows making absolute measurements with a white noise whose power spectrum density ( psd ) level is @xmath2  m.s@xmath0.hz@xmath3 with a cut - off frequency equal to @xmath4  hz . this corresponds , for an integration time of @xmath5  h , to a precision of @xmath6  pm.s@xmath0 and an exact measurement accuracy  @xcite . when taking into account the integration of the instrument in the spacecraft , and in particular the alignment accuracy ( @xmath7 @xmath6  mrad ) , the positioning accuracy and the spacecraft self - gravity , a global precision of 10  pm.s@xmath0 is expected .    the spacecraft of the oss mission has been specifically designed to reach this level of precision . the three main drivers taken into account are to ( i ) provide the lowest and most axisymmetrical gravitational field as viewed from gap , ( ii ) make coincide the dry mass center of gravity , the propellant center of gravity , the radiation pressure force line and the gap , and ( iii ) ensure a stable and reliable alignment between the gap and the high gain antenna to ensure consistency between radio science data and non - gravitational acceleration measurements . this led to the architecture presented in  @xcite : the platform is built as a flat ring and the gap is on a settable plate at the center , held by thermally stable struts . in addition , four hydrazine tanks are distributed symmetrically around the gap , with fluidic interconnections . it allows balancing the quantity of propellant and thus controlling the position of the center of mass of the spacecraft . the goal of the research project presented in this article is to validate experimentally the data processing method developed in @xcite . first , it will be shown that it is experimentally possible to completely remove the bias from the measurements . then , the precision on the unbiased acceleration obtained after post - processing will be compared to the predicted precision given the noise of the devices used in this experiment . these results will demonstrate the possibility to make unbiased acceleration measurements in real conditions using in space technology . therefore , they are a key milestone toward embarking the gravity advanced package on an interplanetary mission .    since the gravity advanced package is not available yet , commercial devices were used to mimic it : the `` measurement device '' which replaces the gravity advanced package and the overall experimental setup is described in section  [ section : materials ] . in section [ section : methods ] , the experimental methodology as well as the data processing is presented . finally , the results are discussed in section  [ section : results ] . the overall experimental setup is shown in fig .  [ fig : experimentalsetup ] shows . the `` measurement device '' , used to mimic the gravity advanced package , is described in section  [ subsection : accelerometer ] . a pendulum controlled by an electrostatic accelerometer , called dm1 , is used to impose a monitored external acceleration to the measurement device ( cf . section  [ subsection : pendulum ] ) . the measurement device is composed of :    * a rotating platform rgv100bl from newport @xcite which acts as the bias rejection system . the axis of rotation is called @xmath8 and is vertical . the rotation is parameterized by an angle called @xmath9 . the angle repeatability is equal to @xmath10. given the noise of the q - flex ( cf . [ fig : psd_qflex ] ) , the rotating platform does not introduce uncertainty in the measurement process . * two q - flex qa-700 accelerometers from honeywell @xcite . they are mono - axial accelerometers and are positioned orthogonally so as to have bi - axial measurements . the q - flex a and b measure local gravity along , respectively , the @xmath11 and @xmath12 axis of the frame attached to the rotating part of the rotating platform . .measured scale factor of the q - flex a [ cols= " < , < " , ]     the effect of the sampling frequency has also been investigated . because the harmonics of the calibration signal in the frequency domain decreases rapidly , it is expected that the precision does not depend on the sampling frequency as far as it is two order of magnitude larger than the modulation frequency . that is why the sampling frequency does not appear in the approximated formula . in these experiments , the larger modulation frequency is @xmath13  hz and the explored sampling frequencies are 10 , 100 and 200  hz , which are two orders of magnitude larger than 0.0167  hz . the results are summarized in table  [ table : frequency ] and show that the precision does not depend significantly on the sampling frequency for the range explored in these experiments . the set of experiments presented in this article was designed to validate the data processing scheme developed in order to remove the bias from the measurements made with an electrostatic accelerometer . to do so , two q - flex accelerometers mounted on a rotating platform were placed on a pendulum whose inclination was precisely controlled via a high - precision electrostatic accelerometer . the data gathered provided an experimental validation of the method . first , the ability to recover , after post - processing , the mean and the slope of the bias and of the unbiased acceleration has been demonstrated with experimental data . then , concerning the precision of the unbiased acceleration measurement , the experimental precisions have been compared to the theoretical predictions and a solid agreement has been obtained .    the experimental validation of this method devised to make unbiased acceleration measurement with an electrostatic accelerometer opens new opportunities . it will allow improved orbit restitution of interplanetary probes using an instrument with a strong heritage . we express our deep appreciation of the support from our colleagues at onera : damien boulanger , dominique chauvin , patrick flinoise , bernard foulon , yoann lelay and franoise liorzou . we are also grateful to cnes ( centre national dtudes spatiales ) for its financial support . b.  christophe , p.  h. andersen , j.  d. anderson , s.  asmar , p.  brio , o.  bertolami , r.  bingham , f.  bondu , p.  bouyer , s.  bremer , j.  m. courty , h.  dittus , b.  foulon , p.  gil , u.  johann , j.  f. jordan , b.  kent , c.  lmmerzahl , a.  lvy , g.  mtris , o.  olsen , j.  pramos , j.  d. prestage , s.  v. progrebenko , e.  rasel , a.  rathke , s.  reynaud , b.  rievers , e.  samain , t.  j. sumner , s.  theil , p.  touboul , s.  turyshev , p.  vrancken , p.  wolf , and n.  yu . . , 23(2):529547 , 2009 . http://dx.doi.org/10.1007/s10686-008-9084-y [ ] . b.  christophe , l.  j. spilker , j.  d. anderson , n.  andr , s.  w. asmar , j.  aurnou , d.  banfield , a.  barucci , o.  bertolami , r.  bingham , p.  brown , b.  cecconi , j.  m. courty , h.  dittus , l.  n. fletcher , b.  foulon , f.  francisco , p.  j.  s. gil , k.  h. glassmeier , w.  grundy , c.  hansen , j.  helbert , r.  helled , h.  hussmann , b.  lamine , c.  lmmerzahl , l.  lamy , b.  lenoir , a.  lvy , g.  orton , j.  pramos , j.  poncy , f.  postberg , s.  v. progrebenko , k.  r. reh , s.  reynaud , c.  robert , e.  samain , j.  saur , k.  m. sayanagi , n.  schmitz , h.  selig , f.  sohl , t.  r. spilker , r.  srama , k.  stephan , p.  touboul , and p.  wolf . . , 34(2):203242 , 2012 . http://dx.doi.org/10.1007/s10686-012-9309-y [ ] . u.  johann et  al . . in h.  dittus , c.  lmmerzahl , and s.  g. turyshev , editors , _ lasers , clocks and drag - free control _ , volume 349 of _ astrophysics and space science library _ , pages 577604 . springer , 2008 . http://dx.doi.org/10.1007/978-3-540-34377-6_26 [ ] . b.  lenoir , b.  christophe , a.  lvy , b.  foulon , s.  reynaud , j.  m. courty , b.  lamine , h.  dittus , t.  van zoest , c.  lmmerzahl , h.  selig , s.  lon - hirtz , r.  biancale , g.  mtris , f.  sohl , and p.  wolf . . in _ 61st international astronautical congress _ , prague , czech republic , 10 2010 . http://publications.onera.fr/exl-doc/doc397477_s1.pdf[iac-10.a3.6.5 ] . http://arxiv.org/abs/1107.2316 [ ] . b.  lenoir , b.  christophe , and s.  reynaud . . in _ journes 2011 de la socit franaise dastronomie & dastrophysique _ , paris , france , 2011 . http://lesia.obspm.fr / semaine - sf2a/2011/proceedings/2011/2011sf2a.conf.% .0663l.pdf[http://lesia.obspm.fr / semaine - sf2a/2011/proceedings/2011/2011sf2a.conf.% .0663l.pdf ] . http://arxiv.org/abs/1110.0342 [ ] . b.  lenoir , a.  lvy , b.  foulon , b.  lamine , b.  christophe , and s.  reynaud . electrostatic accelerometer with bias rejection for gravitation and solar system physics . , 48(7):12481257 , 2011 . http://dx.doi.org/10.1016/j.asr.2011.06.005 [ ] . newport . . online [ 2011/06/28 ] : http://www.newport.com / rgv100bl - high - speed - precision - rotation - stage/542% 416/1033/catalog.aspx[http://www.newport.com / rgv100bl - high - speed - precision - rotation - stage/542% 416/1033/catalog.aspx ] .
Abstract:

Answer: in order to test gravitation in the solar system , it is necessary to improve the orbit restitution of interplanetary spacecrafts . the addition of an accelerometer on board is a major step toward this goal because this instrument measures the non - gravitational acceleration of the spacecraft . it must be able to perform measurements at low frequencies with no bias to provide an additional observable of interest . since electrostatic accelerometers suffer a bias , a technological upgrade has been proposed by onera . it consists in adding to an electrostatic accelerometer a rotating platform which allows modulating the signal of interest and retrieving it without bias after post - processing . using this principle , a measurement method and a post - processing method have been developed . the objective of this article is to validate these methods experimentally . to do so , a horizontally controlled pendulum was used to apply a known signal to an accelerometer mounted on a rotating platform . the processing of the experimental data demonstrates the ability to make acceleration measurements with no bias . in addition , the experimental precision on the unbiased acceleration obtained after post - processing corresponds to the precision predicted theoretically . [ [ keywords ] ] keywords + + + + + + + +    electrostatic accelerometer ; bias rejection ; modulation ; data processing ; precision ; colored noise . = 1


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: the generic matrix - matrix multiplication ( gemm ) is given by the equation : @xmath0 where @xmath1 , @xmath2 and @xmath3 are matrices , and @xmath4 and @xmath5 are scalars . gemm is arguably the most popular computational kernel of the 20th century . the apparent simplicity of gemm has haunted generation after generation of researchers who have evaluated its performance on generation after generation of computer systems , while uncovering layer after layer of its hidden complexity . for example , discovering the beneficial effects of cache blocking on gemm performance  @xcite has fuelled research on locality optimizations in compilers for many years .    yet , surprisingly , no common methodology for evaluating gemm performance has been established over the many decades of using this kernel for comparing architectures , compilers and ninja - class programmers . consequently , the reader of a report presenting gemm results is often left wondering :    * was the kernel specialized , for example , to @xmath6 ? ( in other words , @xmath7 and @xmath8 . ) * which of the data types were used : single precision ( sgemm ) , double precision ( dgemm ) , complex single precision ( cgemm ) , or complex double precision ( zgemm ) ? * which data layouts were used : normal ( n ) or transposed ( t ) ? typically results in better locality , because @xmath9 is read row - wise . ] if transposed , did the execution time include the overhead for transposition ? * which data shapes were used : square or rectangular ? if rectangular , did the execution time depend on the ratio between the dimensions ? * which data sizes were used : small or large ? * on a system with caches , did ` large ' result in cache thrashing ; did ` small ' result in good locality ( no thrashing ) ? * on a heterogeneous system equipped with a discrete accelerator , did the execution time include the overhead for copying the data to the accelerator and back , or only the kernel execution time ? * did the evaluation include power or energy measurements ? * if a diesel generator was used to get the system running , how many megaflops per gallon were they getting ? * more seriously , have we achieved significant improvements in energy efficiency of floating - point operations over the last decade ? * how much human effort and ingenuity was involved in writing the kernel or in implementing the compiler that generated the kernel ? * can we compare the generators , for example , based on polyhedral compilation  @xcite and functional expression rewriting  @xcite in a fair way ( including code quality , code generation time and robustness ) ? * can we evaluate the generators against ninja - class programmers  @xcite or vendor libraries ? * have we used all the tricks up our sleeves to get the fastest gemm implementation for our hardware and problem at hand ? * can we _ adapt _ our gemm implementations to work well across a range of architectures , data types , data sizes , etc . ? given that we are discussing gemm , a simple kernel intended to give us insights for solving more complex ` real - world ' problems , it is essential to start getting some of the answers right to facilitate our learning and knowledge sharing . we introduce gemmbench , a framework and methodology for evaluating performance of gemm implementations . gemmbench is implemented on top of collective knowledge ( ck ) , a lightweight framework for reproducible and collaborative r&d in computer systems . our initial implementation supports hand - written opencl kernels operating on matrices consisting of single- and double - precision floating - point values , and producing single or multiple output elements per work - item ( via thread coarsening and vectorization ) .    over time , we plan to involve the community to add further hand - written and generated kernels ( e.g. from @xcite ) , and , importantly , to collectively study the gemm performance across multiple platforms , data sizes and data types . the gemmbench framework reads from a json file the metadata describing a kernel . the json file specifies the data type ( s or d ) , the layout of the matrices ( n or t ) , the thread - coarsening configuration ( di for the number of rows and dj for the number of columns in a block computed by a single work - item ) , and so on . for example , the sgemm kernel that assumes that @xmath1 is non - transposed and @xmath2 is transposed and outputs a single element per work - item :    .... kernel void gemm (      global float const * restrict a ,      global float const * restrict b ,      global float        * restrict c ,      float alpha , float beta , uint n ) {      const uint j = get_global_id(0 ) ;      const uint i = get_global_id(1 ) ;        float abij = 0.0f ;      for ( uint k = 0 ; k < n ; k + = 1 )      {          abij + = a[i*n + k ] * b[j*n + k ] ;      }      c[i*n + j ] = alpha * abij + beta * c[i*n + j ] ; } ....    is described by the following metadata :    .... {      " name "    : " sgemm_nt_1x1 " ,      " file "    : " sgemm_nt_1x1.cl " ,      " type "    : " s " ,      " transa " : " n " ,      " transb " : " t " ,      " dj "      : 1 ,      " di "      : 1 } ....    see further examples in the dataset entries of the gemmbench repository . install collective knowledge ( ck ) _ e.g. _ : .... $ export ck_repos=~/ck $ mkdir $ ck_repos $ export ck_root=$ck_repos / ck $ git clone https://github.com/ctuning/ck.git $ ck_root $ export path=$ck_root / bin:$path $ ck status your version is up - to - date : v1.6.12 ....    for using ck in python scripts ( _ e.g. _  with ipython notebook ) :    .... $ cd $ ck_root & & sudo python setup.py install $ python -c " import ck.kernel as ck ;               print ck.version({})['version_str ' ] " 1.6.12 ....      install gemmbench along with other ck repositories it depends upon : ck - autotuning , ck - env , _ etc . _ .... $ ck pull repo : gemmbench \ --url = https://github.com / dividiti / gemmbench ....      register an ( already installed ) opencl driver _ e.g. _ :    .... $ ck find soft : lib.opencl * ~/ck / ck - env / soft / lib.opencl.linux ... $ ck setup soft : lib.opencl.linux ....    enter the opencl driver version ( a string to identify it later ) . enter the path to ` libopencl.so ` without ` lib ` _ e.g. _ ` /usr ` for ` libopencl.so ` in ` /usr / lib ` . take a look at one of the available scripts for disabling frequency and voltage scaling ( dvfs ) and fixing the frequencies :    .... $ ck find platform.init : * ~/ck / ck - autotuning / platform.init / generic - android ~/ck / ck - autotuning / platform.init / chromebook - ubuntu ~/ck / ck - autotuning / platform.init / generic - odroid ~/ck / ck - autotuning / platform.init / generic - linux ....    customize the scripts called from ck - set - performance ( you may leave them blank initially ) and add them to the system path _ :    .... $ cp ~/ck / ck - autotuning / platform.init / generic - linux \ ~/ck / ck - autotuning / platform.init / my - linux - platform ... $ export path=$path:\    ~/ck / ck - autotuning / platform.init / my - linux - platform ....    ( ensure the scripts are executable . )      to compile gemmbench :    .... $ ck compile program : gemmbench - cl - launcher-1.0 ....    ( the cjson and xopenme libraries should be installed automatically the first time you compile gemmbench . )      to run gemmbench with the default parameters :    .... $ ck run program : gemmbench - cl - launcher-1.0 ....    select the default command ( press `` enter '' ) , one of the four currently supported `` flavours '' ( sgemm nn , sgemm nt , dgemm nn , dgemm nt ) , and finally one of the kernel variants .    to override the default parameters , use _ e.g. _    .... $ ck run program : gemmbench - cl - launcher-1.0 \ --extra_run_cmd ="- p 1 -d 1 -n 512 -lws 4,16 " ....    to run on platform 1 , device 1 , with the matrix order of 512 and the local work size of ( 4,16 ) . .... $ ck find script : sgemm * ~/ck / gemmbench / script / sgemm_nt $ cd ~/ck / gemmbench / script / sgemm_nt $ ./_clean_program_pipeline.sh $ ./_setup_program_pipeline.sh ... pipeline is ready ! $ ./explore - f - n.sh $ ./explore - n - lws.sh ....      to replay our experiments , obtain our experimental data for this paper :    .... $ ck pull repo : gemmbench - adapt16 \ --url = https://github.com / dividiti / gemmbench - adapt16 $ ck find experiment : sgemm_nt * ~/ck / gemmbench - adapt16/experiment / sgemm_nt - explore - f - n ~/ck / gemmbench - adapt16/experiment / sgemm_nt - explore - n - lws ....    start the ck web server :    .... $ ck start web .... open http://localhost:3344 in a web browser . select gemmbench - adapt16 in the `` repository '' dropdown menu . open sgemm_nt - explore - f - n or sgemm_nt - explore - n - lws . select gemmbench - view in the `` select experiment view '' menu under a qr code . the table shows one experiment ( with a number of statistical repetitions ) per row . click on a `` copy to clipboard '' button in the rightmost column to obtain a command to replay the corresponding experiment . explains the meaning of most other columns . ] for example , to replay the intermittently failing experiment mentioned in section  [ sec : validating ] , run :    .... $ ck replay \ experiment:8bcbe025bd8803c2 --point = ca4a5dbe25613c7d .... we demonstrate using the gemmbench framework for evaluating 3 sgemm nt opencl kernels on a hardkernel odroid  xu3 board ( table  [ odroid ] ) . the odroid xu3 board has 4 integrated power consumption sensors :    * for the lpddr3 ram ; * for the cpu cluster 0 comprised of 4 arm cortex - a7 ( `` little '' ) cores ; * for the cpu cluster 1 comprised of 4 arm cortex - a15 ( `` big '' ) cores ; * for both the gpu cluster 0 and the gpu cluster 1 comprised respectively of 4 and 2 arm mali - t628 cores . we reused the _ pipeline _ functionality of the underlying collective knowledge framework to conduct experiments under controlled conditions . we ran the sgemm kernels on the gpu cluster 0 ( opencl device 0 ) ; we disabled dynamic voltage and frequency scaling ( dvfs ) and set the frequency to the maximum of 600 mhz . likewise , we set the cpu governors to the performance mode and the cpu frequencies to the maximum .      the sgemm nt kernels are contained in 3 separate files :    * ` sgemm_nt_1x1.cl ` : a nave version shown in section  [ sec : implementation ] which computes a single element of matrix @xmath3 per work - item ; * ` sgemm_nt_4x1.cl ` : a vectorised version which computes a vector of four adjacent elements of matrix @xmath3 per work - item ; * ` sgemm_nt_4x1_barrier.cl ` : a similarly vectorised version which synchronises work - items in a work - group with a barrier to improve cache utilisation  @xcite .    [ cols="<,<",options="header " , ]     for the ` sgemm_nt_1x1.cl ` program and the matrix order of 64 , no results matched . in fact , the maximum absolute differences suggest a possible bug in either the opencl or the reference implementation . ( neither has been around for long or code - reviewed . )    for the ` sgemm_nt_4x1.cl ` program and the matrix order of 64 , the results are mixed : they twice matched and twice did not . the failures may be difficult to debug , since they are intermittent . luckily , replaying an experimental point under the collective knowledge framework is a matter of running a single command for the command to replay this very experiment . ] which would help investigate the failures and quickly test potential solutions .    for the ` sgemm_nt_4x1_barrier.cl ` program , the results show repeatable failures for the orders of 64 , 96 , 192 . this may give us a clue to what goes wrong here . returning to choosing the value of @xmath10 , we note that these results would pass under @xmath11 . it is likely , however , that even this @xmath10 would need to be changed had the input values been drawn from a different distribution _ e.g. _  the uniform distribution over the range @xmath12 .    rather than making an arbitrary choice of @xmath10 for an arbitrary choice of the random distribution , we could look into defining representative datasets . ideally , they would come from real - world problems along with precision requirements . we have presented gemmbench , a framework and methodology for systematically evaluating performance of matrix multiplication implementations . our initial implementation supports hand - written opencl kernels , producing single or multiple output elements per work - item ( via thread coarsening and vectorization ) . our goal is to involve the community to extend gemmbench to evaluate performance of compiler - generated opencl kernels , non - opencl implementations , library implementations and so on , across many target platforms . to this end , the underlying collective knowledge framework provides unique opportunities for the community to gradually gather and share valuable knowledge for optimizing performance of matrix multiplication and other programs , as well as of compilers and processors . we will build upon other strengths of the collective knowledge framework including support for multiple operating systems ( windows , linux , android and macos ) , compilers ( llvm , gcc , icc , msvc , _ etc . _ ) and interfaces to packages for data mining and predictive analytics .    where do we start ? first , we encourage the interested reader help us investigate the failures reported in section  [ sec : validating ] . eric s.  raymond s proposition that `` given enough eyeballs , all bugs are shallow '' should apply well in this case . indeed , some of the failures may have nothing to do with numerical ( in)stability . programming errors ( such as tacit assumptions ) may be detected with static and dynamic analysis tools such as gpuverify  @xcite and oclgrind  @xcite . in the final version of this article , we will acknowledge those who help us explain the existing failures and perhaps find new ones , and of course fix them . second , we welcome contributions in the form of experimental data in the collective knowledge format . the contributors should acknowledge their compliance with the terms of use of their systems , specifically that they do not breach confidentiality .. as far as we are aware , the standard hardkernel bsp is distributed without such limitations . ] third , we welcome contributions in the form of improvements for the core gemmbench code _ e.g. _  support for rectangular matrices and complex floating - point numbers . fourth , we welcome contributions in the form of opencl kernels or other implementations . the initial set of kernels optimised for the arm mali - t600 architecture is intentionally small . we would like to see contributed kernels optimised for different architectures . the contributors should acknowledge their copyright in the source code and specify the licensing terms . fifth , we welcome contributions in the form of representative datasets or their descriptions . we envision gemmbench to inspire community - driven development of other representative workloads for use in performance evaluation and optimisation of computer systems . we thank grigori fursin , cto of dividiti and chief scientist of the ctuning foundation , for designing and implementing the collective knowledge framework , on top of which we implemented gemmbench . we look forward to the public discussion of this draft and hope to get useful feedback and contributions from the community , which will be acknowledged in the final version of this article .    1    u.  beaugnon , a.  kravets , s.  van haastregt , r.  baghdadi , d.  tweed , j.  absar , and a.  lokhmotov . : a vehicle for optimized basic linear algebra . in _ proceedings of the 2014 sigplan / sigbed conference on languages , compilers and tools for embedded systems _ , lctes 14 , pages 115124 , new york , ny , usa , 2014 . a.  betts , n.  chong , a.  donaldson , s.  qadeer , and p.  thomson . gpuverify : a verifier for gpu kernels . in _ proceedings of the acm international conference on object oriented programming systems languages and applications _ , oopsla 12 , pages 113132 , new york , ny , usa , 2012 . g.  fursin , a.  memon , c.  guillon , and a.  lokhmotov . : towards performance- and cost - aware software engineering as a natural science . g.  fursin , r.  miceli , a.  lokhmotov , m.  gerndt , m.  baboulin , d.  malony , allen , z.  chamski , d.  novillo , and d.  d. vento . : towards practical and collaborative auto - tuning . , 22(4):309329 , july 2014 . k.  goto and r.  a. v.  d. geijn . anatomy of high - performance matrix multiplication . , 34(3):12:112:25 , may 2008 . j.  gronqvist and a.  lokhmotov . optimizing opencl kernels for the arm mali - t600 gpus . in _ gpu pro 5 : advanced rendering techniques_. 2014 . m.  d. lam , e.  e. rothberg , and m.  e. wolf . the cache performance and optimizations of blocked algorithms . in _ proceedings of the fourth international conference on architectural support for programming languages and operating systems _ , asplos iv , pages 6374 , new york , ny , usa , 1991 . j.  price and s.  mcintosh - smith . oclgrind : an extensible opencl device simulator . in _ proceedings of the 3rd international workshop on opencl _ , 2015 . m.  steuwer , c.  fensch , s.  lindley , and c.  dubach . generating performance portable code using rewrite rules : from high - level functional expressions to high - performance opencl code . in _ proceedings of the 20th acm sigplan international conference on functional programming _ , icfp 2015 , pages 205217 , new york , ny , usa , 2015 . sharing gemmbench experimental results is easy ! optionally , provide a license file and update further details in .cm / info.json and .cm / meta.json in your experimental entries . ( if not provided , the standard ck license and copyright details will apply . )        please do not forget to send us a link to your new repository ! we will maintain a list of such repositories for the community to access , similarly to this paper s repository ( section  [ sec : reproduce ] ) .
Abstract:

Answer: the generic matrix - matrix multiplication ( gemm ) is arguably the most popular computational kernel of the 20th century . yet , surprisingly , no common methodology for evaluating gemm performance has been established over the many decades of using gemm for comparing architectures , compilers and ninja - class programmers . we introduce gemmbench , a framework and methodology for evaluating performance of gemm implementations . gemmbench is implemented on top of collective knowledge ( ck ) , a lightweight framework for reproducible and collaborative r&d in computer systems . using ck allows the r&d community to crowdsource hand - written and compiler - generated gemm implementations and to study their performance across multiple platforms , data sizes and data types . our initial implementation supports hand - written opencl kernels operating on matrices consisting of single- and double - precision floating - point values , and producing single or multiple output elements per work - item ( via thread coarsening and vectorization ) .


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: long - range beam - beam interactions are known to cause emittance growth or beam loss in the tevatron @xcite and are expected to deteriorate beam quality in the lhc . a possible remedy to reduce their effects is to increase the crossing angle . however , it has a side effect of reducing the luminosity due to the reduction in geometric overlap . compensation of long - range beam - beam interactions by applying external electromagnetic forces has been proposed for the lhc @xcite . at large beam - beam separation , the electromagnetic force which a beam exerts on individual particles of the other beam is proportional to @xmath0 , which can be generated and canceled out by the magnetic field of a current - carrying wire . direct - current wires were installed in sps @xcite , da@xmath1ne @xcite , and later in rhic @xcite for tests . results of the sps wire excitation experiments have been reported earlier @xcite . during the kloe experiment in da@xmath1ne , an improvement of beam lifetime without luminosity loss due to the compensation of parasitic collisions has been observed @xcite . experiments in rhic supplement the sps and da@xmath1ne tests . the beam lifetime in rhic is typical for a collider storage ring and better than in the other two machines . unlike the sps , two beams circulate and interact in rhic and compensation of a single long - range interaction is possible . the energies and other beam conditions are also different from the sps and da@xmath1ne . so the results reported here will result in simulations being benchmarked under a wider range of operating conditions . two current carrying wires , one for each beam , have been installed between the magnets @xmath2 and @xmath3 of ip6 in the rhic tunnel . their impact on a beam was measured during the physics run with deuteron and gold beams . no attempt was made to compensate the beam - beam interaction since so far only ion beams were available which have a beam - beam parameter about a factor of 3 smaller than proton beams . however , the experimental results help to understand the beam - beam effects because the wire force is similar to the long - range beam - beam force at large separations . in this report we discuss the results of numerical simulations of a wire acting on a beam in rhic using a multi - particle tracking code and compare with measurements in rhic run 7 and 8 @xcite . the organization of the paper is as follows : the physical models used in the simulation code are described in section [ sec : model ] . the beam and wire parameters are given in section [ sec : parameters ] . section [ sec : results ] presents the simulation results of the effects of the wire on the beam dynamics and the comparison between measurement and simulation results . we summarize our results in section [ sec : summary ] . in the collider simulation , the two beams moving in opposite direction are represented by macroparticles of which the charge to mass ratio is that of each beam . the number of macroparticles is much less than the bunch intensity of the beam because even with modern supercomputers it is too time consuming to track 10@xmath4 particles for the number of revolutions of interest . the macroparticles are generated and loaded with an initial distribution for a specific simulation purpose , for example , six - dimensional gaussian distribution for long - term beam evolution . the transverse and longitudinal motion of particles is calculated by transfer maps which consist of linear and nonlinear maps . in the simulations , the following nonlinearities are included : head - on and long - range beam - beam interactions , external electromagnetic force by current carrying wire , multipole errors in the interaction region ( ir ) quadrupole triplets , and sextupoles for chromaticity correction . the finite bunch length effect of the beam - beam interactions is considered by slicing the beam into several sections in the longitudinal direction and by a synchro - beam map @xcite . each slice in a beam interacts with particles in the other beam in turn at the collision points . in addition , the magnetic field of a finite length of wire is calculated using the biot - savart law at the particles betatron amplitude . the six - dimensional accelerator coordinates @xmath5 are applied , where @xmath6 and @xmath7 are horizontal and vertical coordinates , @xmath8 and @xmath9 the trajectory slopes of each coordinates , @xmath10 the longitudinal distance from the synchronous particle , and @xmath11 the momentum deviation from the synchronous energy . coupling between the transverse planes is included in the transfer map between elements . we adopt the weak - strong model to treat the beam - beam interactions where one beam is strong and is not affected by the other beam while the other beam is weak and experiences a beam - beam force due to the strong beam during the interactions . the density distribution of the strong beam is assumed to be gaussian . it is well known that for a large separation distance @xmath12 at parasitic crossings , the strength of long - range interactions is inversely proportional to the distance . its effect on a test beam can be compensated by current carrying wires which create just the @xmath0 field . the advantage of such an approach consists of the simplicity of the method and the possibility to deal with all multipole orders at once . for a finite length of a wire embedded in the middle of a drift length @xmath13 and tilted in pitch and yaw angles , the transfer map of a wire can be written as @xcite @xmath14 where @xmath15 represents the tilt of the coordinate system by horizontal and vertical angles @xmath16 to orient the coordinate system parallel to the wire , @xmath17 is the drift map with a length @xmath18 , @xmath19 is the wire kick integrated over a drift length , and @xmath20 represents a shift of the coordinate axes to make the coordinate systems after and before the wire agree .    at the start of the simulation , the particles in the weak beam are distributed over the phase space @xmath5 . the number of particles @xmath21 is limited by the computational power . in order to make the best use of a relatively small number of simulation particles compared to the bunch intensity , the initial distribution should be optimized . indeed the initial distribution is very important because a proper choice can reduce the statistical noise in the physical quantities . the simulation particles are generated in two steps :    1 . the particle coordinates @xmath22 of particles can be directly generated from the spatial gaussian distribution , @xmath23 , where @xmath24 . since the particle coordinates are not correlated , one can generate them by inverse mapping of each cumulative distribution function of horizontal , vertical , and longitudinal gaussian distributions , and together with the bit - reversed sequence to minimize nonphysical correlations @xcite . 2 .   an equilibrium distribution in transverse phase space e.g. in the @xmath25 plane is @xmath26 . since the spatial coordinate @xmath6 is determined at the first step , the slope @xmath8 can be obtained from the random variate @xmath27 of a univariable gaussian , i.e. , @xmath28 .    following the above physical model , a beam - beam simulation code ` bbsimc ` has been developed at fnal over the past few years to study the effects of the machine nonlinearities and the beam - beam interactions @xcite . if required , time dependent effects such as tune modulation and fluctuation , beam offset modulation and fluctuation , dipole strength fluctuations to mimic rest - gas scattering etc can be included in the model . the code is under continuous development with the emphasis being on including the important details of an accelerator and the ability to reproduce observations in diagnostic devices . at present , the code can be used to calculate tune footprints , dynamic apertures , beam transfer functions , frequency diffusion maps , action diffusion coefficients , emittance growth and beam lifetime . calculation of the last two quantities over the long time scales of interest is time consuming even with modern computer technology . in order to run efficiently on a multiprocessor system , the resulting model was implemented by using parallel libraries which are mpi ( interprocessor message passing interface standard ) @xcite , state - of - the - art parallel solver libraries ( portable , extensible toolkit for scientific calculation , petsc ) @xcite , and hdf5 ( hierarchical data format ) @xcite . the transverse electric field of a round gaussian bunch with transverse rms size @xmath29 is @xmath30 where @xmath31 and @xmath32 are the number of particles and the electric charge in the opposite bunch respectively , and @xmath27 is the radial distance from the center of the beam . at a small radius @xmath33 , the field is proportional to @xmath27 and shifts the betatron tune . this tune shift is characterized by the beam - beam parameter @xmath34 , where @xmath35 is the classical radius of the particle , @xmath36 is the normalized emittance . while particles at small @xmath27 undergo a linear tune shift , the particles with @xmath37 experience a @xmath0 force . the long - range effect is nonlinear and may vary from bunch to bunch if the the bunch pattern is asymmetric . a current carrying wire generates a magnetic force which is @xmath38 , the same as the long - range beam - beam force at large separations . the wire current required to compensate a long - range interaction is @xmath39 , where @xmath40 is the wire current , and @xmath41 its length . in addition to the wire strength , the phase advance from the location of the wire to the location of the long - range interaction should be small for an effective compensation of the long - range interaction @xcite . in rhic , a possible parasitic collision point is near the separation dipole dx magnets near ip6 and ip8 . there is sufficient space in the straight sections between quadrupoles @xmath42 and @xmath43 for the wires . the phase advance between the dx magnet and the wire location is 5.7@xmath44 for @xmath45 1 m and is expected to provide effective compensation @xcite . figure [ fig : layout ] ( a ) shows the overall layout of rhic which has two head - on beam - beam collisions at ip6 and ip8 . the yellow beam revolves in a counter clockwise direction , and the blue beam revolves clockwise . the elements and the location of the wire at ip6 are shown in fig . [ fig : layout ] ( b ) . two wires are currently installed in rhic : one for each beam . in these simulations the effects of a current carrying wire on a beam have been studied for different rhic runs : gold beam at both injection and store energies , and deuteron beam at store energy . the parameters of rhic used in our simulations are shown in table [ tab : rhic ] .    [ cols="^,^,^,^,^ " , ]     the simulations have studied smaller separations than would be practical in rhic - the large loss rates seen at the smallest separation would have quenched the machine . however , the simulations show the steep climb in the loss rate beyond a threshold separation and consequently will have a higher power law behavior . the power law at injection is somewhat higher than at collision but the power laws for gold and deuteron beams are very close despite the large difference in loss rates . we expect that the power laws depends on the details of the machine and not to be universal . in fact different power laws have been reported for the sps and the tevatron @xcite . in order to study the effects of the machine nonlinearities and the beam - beam interactions , including strong localized long - range beam - beam interactions , we have developed a six - dimensional weak - strong code ` bbsimc ` . machine nonlinearities , beam - beam collisions and the field of a current carrying wire can be included in the model . we have studied the effects of the wire on the beam dynamics in three different cases : gold beam at injection , gold beam at store , and deuteron beam at store . results show that the betatron tune change due to the wire is well tracked by the simulation , and that the stability of particle motion is strongly influenced by the wire which causes a significant increase in tune spread and diffusion for both gold and deuteron beam . tune diffusion with the gold beam and the wire is appreciable only at amplitudes larger than 2 @xmath29 while the tune diffusion with the deuteron beam and the wire is significantly larger extending to the beam core at amplitudes down to 0.5 @xmath29 .    at injection energy , the dynamic aperture is largely determined by the wire , while the wire and the ir multipoles have a major impact on the stability boundary at collision . with the beam - wire separation only along the vertical axis , it is found that the stability boundary near the vertical axis is linearly proportional to the beam - wire separation . we observe , from the tune scan of dynamic aperture , that at injection the largest and the smallest dynamic apertures are distributed along a band parallel to the zone along @xmath46 and @xmath47 respectively at all wire separations while at collision the largest dynamic apertures are distributed nearly along the diagonal between @xmath48 and @xmath49 . by modeling the stripline kicker we obtain the amplitude and phase angle of the transverse btf . results show that the betatron tune is well identified by the simulation , and that the wire changes characteristics of the beam response significantly in the horizontal plane , similar to the result for the tune shift .    in simulations the wire enhances the diffusion . this effects is particularly pronounced in the plane of the wire , and at amplitudes of about 3 @xmath29 , where the wire can increase the diffusion rate by more than an order of magnitude . the action diffusion for the deuteron beam is larger than for the gold beam , similar to the result for frequency diffusion . simulations of the beam loss rate when the wire is present are in good agreement with the experimental observations . the threshold separation at which there is a sharp rise in the loss rates agree to better than 1 @xmath29 . the loss rates with the deuteron beam are found in measurements and simulations to be nearly an order of magnitude higher than with the gold beam . comparisons of tune footprints and dynamic apertures of these two cases show no indication that the loss rates may be higher with the deuteron beam and in fact the tune footprint with the deuteron beam is free of low order resonances . however frequency diffusion and action diffusion with the deuteron beam are substantially higher than with the gold beam implying that these diffusion measures may be better indicators of loss rates .    in this paper simulation results of the beam - wire interactions were compared with measurements taken when the long - range beam - beam interactions were not present . it is expected that compensation of parasitic interactions using the wires will be tested during the rhic run in 2009 . we plan to compare simulations with these forthcoming experimental results and also to determine the effectiveness of the compensation of long - range interactions in the lhc . we thank y. luo and r. calaga for help with the rhic lattice . v. boocha , b. erdelyi and v. ranjbar made significant contributions to the development of the code . some of the parallel computations were performed at the nersc facility at lbl . this work was supported by the us - larp collaboration which is funded by the us department of energy .
Abstract:

Answer: the beam - beam interaction is one of the dominant sources of emittance growth and luminosity lifetime deterioration . a current carrying wire has been proposed to compensate long - range beam - beam effects in the lhc and strong localized long - range beam - beam effects are experimentally investigated in the rhic collider . tune shift , beam transfer function , and beam loss rate are measured in dedicated experiments . in this paper , we report on simulations to study the effect of beam - wire interactions based on diffusive apertures , beam loss rates , and beam transfer function using a parallelized weak - strong beam simulation code ( ` bbsimc ` ) . the simulation results are compared with measurements performed in rhic during 2007 and 2008 .


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: in both spoken and written language , word occurrences are not random but vary greatly from document to document . indeed , the field of information retrieval ( ir ) relies on the degree of departure from randomness as a discriminative indicator . ir systems are typically based on unigram statistics ( often referred to as a `` bag - of - words '' model ) , coupled with sophisticated term weighting schemes and similarity measures  @xcite . in an attempt to mathematically realise the intuition that an occurrence of a certain word may increase the chance that the same word is observed later , several probabilistic models of word occurrence have been proposed . much of this work has evolved around the use of ( a mixture of ) the poisson distribution  @xcite . recently , church and gale have demonstrated that a continuous mixture of poisson distributions can produce accurate estimates of variable word rate  @xcite . lowe has introduced a beta - binomial mixture model which was applied to topic tracking and detection  @xcite . although a constant word rate is an unlikely premise , it is nevertheless adopted in many areas including @xmath0-gram language modelling . in order to address the problem of variable word rate , several adaptive language modelling approaches have been proposed with a moderate degree of success . typically , some notion of `` topic '' is inferred from the text according to the `` bag - of - words '' model . information from different language model statistics ( , a general model and/or models specific to each topic ) are then combined using methods such as mixture modelling  @xcite or maximum entropy  @xcite . the _ dynamic cache model _ @xcite is a related approach , based on an observation that recently appearing words are more likely to re - appear than those predicted by a static @xmath0-gram model . it blends cached unigram statistics for recent words with the baseline @xmath0-grams using an interpolation scheme . theoretically , it should not be necessary to rely on an _ ad hoc _ device such as a cache in order to model variable word occurrences . all the parameters of a language model may be completely determined according to probabilistic model of word rate , such as a poisson mixture .    in this paper , we outline the theoretical background for modelling the variable word rate , and illustrate a key observation that word rates are not static using spoken data transcripts . the constant word rate assumption is then eliminated , and we introduce a variable word rate @xmath0-gram language model . an approach to estimating relative frequencies using prior information of word occurrences is presented . it is integrated with standard @xmath0-gram modelling that naturally involves discounting and smoothing schemes for practical use . using the darpa / nist 4e north american broadcast news task , the approach demonstrates the reduction of perplexity up to 10% . in this section , we illustrate how the assumption of a constant word rate fails to capture the statistics of word occurrence in spoken ( or written ) documents . we show that the word rate is variable and may be modelled using a poisson distribution or a continuous mixture of poissons . the poisson distribution is one of the most commonly observed distributions in both natural and social environments . it is fundamental to the queueing theory : under certain conditions , the number of occurrences of a certain event during a given period , or in a specified region of space , follows a poisson distribution ( a _ poisson process _  @xcite ) .    by assuming randomness in a poisson process , word rate is no longer uniform . firstly , we provide a loose definition of a document as a unit of spoken ( or written ) data of a certain length that contains some topic(s ) , or content(s ) . we consider a model in which a word occurs at random in a fixed length document . for a set of documents we assume that each document produces this word independently and that the underlying process is the poisson with a single parameter @xmath1 . formally , a poisson distribution is a discrete distribution ( of a random variable @xmath2 ) which is defined for @xmath3 such that @xmath4}(x )      = { \cal p}(x = x;\ \lambda )      = \frac{e^{- \lambda } \lambda^x}{x ! } \label{eq : poisson}\end{aligned}\ ] ] whose expectation and variance are given by @xmath5 = \lambda$ ] and @xmath6 = \lambda$ ] , respectively  @xcite .      a less constrained model of variable word rate is offered by a multiple of poissons , rather than a single poisson . suppose the parameter @xmath7 of the pdf  ( [ eq : poisson ] ) is distributed according to some function @xmath8 , then we define a continuous mixture of poisson distributions by @xmath9}(x ) \phi(\lambda ) d\lambda \ . \label{eq : mixture}\end{aligned}\ ] ] in particular , if @xmath8 is a gamma distribution , , @xmath10 for @xmath11 and @xmath12 , then the integral  ( [ eq : mixture ] ) is reduced to a discrete distribution for @xmath3 such that @xmath13}(x ) \hspace{-0.1 in } & = & \hspace{-0.1 in }      { \cal nb}(x = x;\ \alpha , \beta ) \nonumber \\    & = & \hspace{-0.1 in }      \left ( \begin{array}{c } \alpha + x - 1 \\ x \end{array } \right )      \frac{\beta^x}{(1 + \beta)^{\alpha + x } } \ . \label{eq : nb}\end{aligned}\ ] ] this @xmath14}(x)$ ] is a negative binomial distribution be @xmath15 in  ( [ eq : mixture ] ) . this integration is straightforward using the definition of the gamma function , @xmath16 , and the recursion , @xmath17 . the resultant pdf  ( [ eq : nb ] ) has a slightly unconventional form in comparison to that in most of standard textbooks ( ,  @xcite ) , but is identical by setting a new parameter @xmath18 with @xmath19 . ] and its expectation and variance are respectively given by @xmath5 = \alpha\beta$ ] and @xmath6 = \alpha\beta ( \beta + 1)$ ] . the histograms in figure  [ fig : uhist ] show the number of word ( unigram ) occurrences in spoken news broadcast , taken from transcripts of the 4e broadcast news acoustic training data ( 199697 ) . these transcripts were separated into documents according to section markers and those with less than 100 words were removed , resulting 2583 documents containing slightly less than 1.3 million words in total . in the following , the number of word occurrences were normalised to 1000-word length documents .    and  appeared approximately the same number of times across all the transcripts . using a constant word rate assumption , they would have been assigned a probability of around 0.0086 . however their occurrence rates varied from document to document ; about 11% and 33% of all documents did not contain  and  ( respectively ) , while 1% and 3% contained these words more than 30 times . this seems to indicate that occurrences of  is less dependent on the content of the document . a negative binomial distribution was used to model the variable word rate in each case ( the solid line in figure  [ fig : uhist ] ) . the negative binomial seems to model word occurrence rate relatively well for most vocabulary items , regardless of frequency . figure  [ fig : uhist ] illustrates this for one of the most frequent words  ( probability of 0.023 according to the constant word rate assumption ) and the less frequently occurring  ( less than 0.00029 ) . in particular ,  appeared only in 93 out of 2583 documents , but 28 of them contained more than 10 instances , suggesting strong correlation with document content .    . they are fitted by negative binomial distributions ( solid lines).,width=307 ]    we also collected statistics of bigrams appearing in the broadcast news transcripts . figure  [ fig : bhist ] show histograms and their negative binomial fits for bigrams  and . although very sparse ( , they appeared in 127 and 6 documents , respectively ) , this suggests that variable bigram rate can also be modelled using a continuous mixture of poissons . taking word occurrence rate into account changes a probabilistic language model from a situation akin to playing a lottery , to something closer to betting on a horse race : the odds for a certain word improve if it has come up in the past . in this section , we eliminate the constant word rate assumption and present a variable word rate @xmath0-gram language model . let @xmath20 denote a relative frequency after we observe @xmath21 occurrences of word @xmath22 . it is calculated by @xmath23 the function is defined for @xmath24 , where @xmath25 is a fixed document length ( , @xmath25 is normalised to 1000 in figures  [ fig : uhist ] and  [ fig : bhist ] ) . @xmath26 is the occurrence rate for word @xmath22 in an @xmath25-length document ( , poisson , negative binomial ) , satisfying @xmath27 in particular , @xmath28 which corresponds to the case with no prior information of word occurrence . for the conventional approach with the constant word rate assumption , this @xmath29 is not modified regardless of any word occurrences . further , function  ( [ eq : rf_var ] ) satisfies our intuition ; the value of @xmath20 increases monotonically as the number of observation @xmath21 accumulates ( easy to verify ) , and it reaches a unity ( ` 1 ' ) when @xmath30 .    ) . the right figure demonstrates relative frequencies after a certain number of word occurrences . circles ( ` ' ) correspond to relative frequencies under the constant word rate assumption ( 0.023 for  and 0.00029 for ) .,width=307 ]    the characteristics of function  ( [ eq : rf_var ] ) are illustrated in figure  [ fig : rfrequency ] . the right hand figure shows relative frequencies for  and  after a certain number of previous observations of the word . it indicates that the first few instances of the frequent word ( ) do not modify its relative frequency very much , but have a substantial effect on the relative frequency of the less common word ( ) . as the number of observations increases , the former is caught up by the latter .    finally , in order to convert this relative frequency model to any type of probabilistic model for language , normalisation is required . this is achieved by dividing @xmath20 by @xmath31 , where @xmath32 implies a set of vocabulary . variable relative frequencies for bigrams can also be calculated in a similar fashion .      for any practical application , smoothing of the probability estimates is essential to avoid zero probabilities for events that were not observed in the training data . let @xmath33 denote a bigram entry ( a word @xmath34 followed by @xmath22 ) in the model . further , @xmath35 implies a relative frequency after we observe @xmath36 occurrences of the bigram . a bigram probability @xmath37 may be smoothed with a unigram probability @xmath38 . using the interpolation method  @xcite : @xmath39 where @xmath40 implies a `` discounted '' relative frequency ( described later ) and @xmath41 is a non - zero probability estimate ( , the probability that a bigram entry @xmath33 exists in the model ) . alternatively , the back - off smoothing  @xcite may be applied : @xmath42 in  ( [ eq : backoff ] ) , @xmath43 is a back - off factor and is calculated by @xmath44 a unigram probability @xmath38 can be obtained similarly by smoothing with some constant value .    finally , a number of standard discounting methods exist for constant word rate models ( see , ,  @xcite ) . analogous discounting functions for variable word rate models may be @xmath45 for the absolute discounting , and @xmath46 for the good - turing discounting . discounting factors ( @xmath47 and @xmath48 ) may be obtained using zero prior information case  , @xmath49 s of all bigrams in the model  and the rest should be referred to , ,  @xcite or  @xcite .      as noted in section  [ sc : word ] , we extracted 2583 documents from the transcripts of the broadcast news acoustic training data , each with a minimum of 100 words . a vocabulary of 19885 words was selected and 390000 bigrams were counted . in these experiments , the absolute discounting scheme  ( [ eq : absolute ] ) was applied , followed by interpolation smoothing  ( [ eq : interpolation ] ) . figure  [ fig : perplex ] shows perplexities for the reference ( key ) transcription of the 1997 4e evaluation data , containing three hours of speech and approximately 32000 words . using conventional modelling with a constant word rate assumption , unigram and bigram perplexities were 936.5 and 237.9 , respectively . for the variable word rate models , the poisson distribution was adopted because of simplicity in calculation . the number of word occurrences were normalised to @xmath25-word length document with @xmath25 being between 200 and 50000 , and the model parameters were modified ` on - line ' during the perplexity calculation . for each occurrence of a word ( bigram ) in the evaluation data , a histogram of the past @xmath25 words ( bigrams ) was collected and their relative frequencies were modified according to the poisson estimates ( appropriate normalisation applied ) , then discounted and smoothed . as figure  [ fig : perplex ] indicates , the variable word rate models were able to reduced perplexities from the constant word rate models . a unigram perplexity of 843.4 ( 10% reduction ) was achieved when @xmath50 , and a bigram perplexity of 219.0 ( 8% reduction ) when @xmath51 . the difference was predictable because bigrams were orders of magnitude more sparse than unigrams . in this paper , we have presented a variable word/@xmath0-gram rate language model , based upon an approach to estimating relative frequencies using prior information of word occurrences . poisson and negative binomial models were used to approximate word occurrences in documents of fixed length . using the broadcast news task , the approach demonstrated a reduction of perplexity up to 10% , indicating potential although the technique is still premature . because of the data sparsity problem , it is not clear if the approach can be applied to language model components of current state - of - the - art speech recognition systems that typically use 3/4-grams . however , we believe this technique does have application to problems in the area of information extraction . in particular , we are planning to apply these methods to the named entity annotation task , along with further theoretical development . k.  sprck jones , s.  walker , and s.  e. robertson . a probabilistic model of information retrieval : development and status . technical report tr446 , university of cambridge , computer laboratory , 1998 . available from . stephen  a. lowe . the beta - binomial mixture model for word frequencies in documents with applications to information retrieval . in _ proceedings of eurospeech-99 _ , volume  6 , pages 24432446 , budapest , september 1999 . f.  jelinek and r.  l. mercer . interpolated estimation of markov source parameters from sparse data . in _ proceedings of the workshop : pattern recognition in practice _ , pages 381397 , amsterdam , may 1980 .
Abstract:

Answer: the rate of occurrence of words is not uniform but varies from document to document . despite this observation , parameters for conventional @xmath0-gram language models are usually derived using the assumption of a constant word rate . in this paper we investigate the use of variable word rate assumption , modelled by a poisson distribution or a continuous mixture of poissons . we present an approach to estimating the relative frequencies of words or @xmath0-grams taking prior information of their occurrences into account . discounting and smoothing schemes are also considered . using the broadcast news task , the approach demonstrates a reduction of perplexity up to 10% . ( 0,0 ) ( 0,0 )


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: the study of exactly solved lattice models in statistical mechanics has been a richly rewarding area of activity over many decades . among the pantheon of exactly solved models is the two - dimensional ising model , which is fundamental to the theory of phase transitions and critical phenomena @xcite . two - dimensional lattice models with local interactions like the ising model possess remarkable properties at criticality . the fluctuations in such systems are conformally invariant at the critical point . one of the triumphs of statistical physics is that this property entails the classification of universality classes of second order phase transitions in terms of conformal field theories @xcite . roughly speaking , there is an exactly solved model for each universality class . the beauty of exactly solved models is that their critical boltzmann weights satisfy the star - triangle or yang - baxter equation @xcite . a remarkable link has been uncovered between yang - baxter integrability and the condition of discrete holomorphicity @xcite . this involves the establishment of a local identity which is satisfied by an operator defined on the lattice . this operator is an example of a discrete parafermionic observable appearing in conformal field theory . the observable obeys a lattice version of the cauchy - riemann equations and is thus known to be discretely holomorphic . the condition of discrete holomorphicity leads to a system of linear equations which can be solved to yield the boltzmann weights of the lattice model . these are the boltzmann weights which are known to satisfy the yang - baxter equations at criticality . this connection has been observed for a number of exactly solved models , @xcite including the z@xmath0 model,@xcite the potts model @xcite and the o@xmath1 loop model @xcite .    here we will discuss these developments in terms of the o@xmath1 loop model . the o@xmath1 loop model on the square lattice was originally constructed by nienhuis @xcite . it is exactly solved by means of the bethe ansatz , with periodic @xcite and open @xcite boundaries . as an integrable model , its underlying @xmath2-matrix is that of the izergin - korepin , @xcite or @xmath3 vertex model @xcite . the partition function is given by @xmath4 which is a sum over all possible loop configurations @xmath5 . here @xmath6 is the number of occurrences of boltzmann weight @xmath7 and @xmath8 is the total number of closed loops of fugacity @xmath9 in any given configuration . the nine non - zero boltzmann weights @xmath7 are shown in fig . [ fig : weights ] . these weights are given by @xcite @xmath10 with @xmath11 , @xmath12 and @xmath13 . the parafermion operator is defined on the lattice as @xcite @xmath14 here @xmath15 is the probability of configuration @xmath5 , @xmath16 is a set of loop configurations for which points @xmath17 and @xmath18 belong to the same loop , @xmath19 is the winding angle of this loop from points @xmath17 and @xmath18 , and @xmath20 is the spin of the parafermion @xmath21 . for the o@xmath1 loop model @xmath22 for some path @xmath23 across a domain @xmath24 . what ikhlef and cardy @xcite have shown is that the parafermion @xmath21 is discretely holomorphic , i.e. , it satisfies @xmath25 on each vertex as defined in fig . [ fig : vertex ] . this condition is equivalent to requiring that @xmath26 is divergence and curl free .    to establish this , one chooses a particular mid edge @xmath27 . there are four possible different external connectivities , as shown in fig . [ fig : loops ] . these possibilities result in four equations @xmath28 where @xmath29 and @xmath30 . the remarkable point is that ikhlef and cardy @xcite were able to solve this set of linear equations to yield the boltzmann weights @xmath31 and the conformal spin @xmath32 provided @xmath33 . in particular , at @xmath34 , @xmath35 , the isotropic point . substitution of the boltzmann weights ( [ weights ] ) into equations ( [ relations ] ) shows that these equations are indeed satisfied . the surprising point is that the yang - baxter equations are _ cubic _ in the boltzmann weights , whereas imposing discrete holomorphicity led to a set of equations which are _ linear _ in the boltzmann weights . there thus appears to be a simpler route to uncovering integrable lattice models at criticality . why is this ? and what is the precise link between discrete holomorphicity and integrability ? these questions will be taken up further in section 4 in the context of the z@xmath0 model , which can be discussed from an algebraic point of view . we now turn to the related developments in the theory of self - avoiding random walks . the random walk , also known as the drunkard s walk , was introduced over 100 years ago @xcite . the study of random walks has since developed into the central pillar of the theory of probability and random processes @xcite . a key feature of the random walk is that , like the drunkard , it has no memory of where it has been . this absence of memory is the defining property of a markovian process . a more challenging class of mathematical problems arise when a random walk is not allowed to revisit a site it has already visited . such self - avoiding random walks were originally proposed as models of long polymer chains in a good solvent . despite the simplicity of their definition , the non - markovian nature of self - avoiding random walks makes them notoriously difficult to study . nevertheless much is known about the theory of self - avoiding walks , particularly in higher dimensions @xcite . much is also known in two and three dimensions from extensive computer enumeration . the self - avoiding constraint is more restrictive in two dimensions and thus self - avoiding walks on planar lattices are arguably more interesting to study .    in 1982 nienhuis @xcite was able to obtain some key information for self - avoiding walks by considering the o@xmath1 loop model on the honeycomb ( hexagonal ) lattice depicted in fig . [ fig : lattices](a ) . in general , the o@xmath1 model has a connection to self - avoiding walks in the limit @xmath36 @xcite . the critical boltzmann weights of the model on the honeycomb lattice were obtained from an exact mapping to the potts model @xcite . the boltzmann weights defining the critical line of the o@xmath1 loop model on the honeycomb lattice also follow as a condition for the model to be solved exactly in terms of the bethe ansatz @xcite . more generally the o@xmath1 loop model on the honeycomb lattice follows as a limiting case of the exactly solved o@xmath1 loop model on the square lattice . or @xmath37 in ( [ weights ] ) , corresponding to the two orientations in fig .  [ fig : orient ] . ]    the critical boltzmann weights of the model give the value for the connective constant @xmath38 , which is a measure of the growth in the number of allowed self - avoiding walks with increasing length . the nienhuis result , @xmath39 was generally believed to be exact , having been confirmed to high precision from extensive enumeration of self - avoiding walks on the honeycomb lattice . no such exact results are known for the other regular planar lattices . similar exact results are known for self - avoiding walks on related semi - regular lattices . these are the 3 - 12 lattice ( fig . [ fig : lattices](b ) ) and the martini lattice ( fig . [ fig : lattices](c ) ) . for the 3 - 12 lattice , the exact value @xcite @xmath40 follows from the solution of @xmath41 and for the martini lattice , the exact value @xcite @xmath42 follows from the solution of @xmath43 more generally , the o@xmath1 model on the honeycomb lattice can be mapped to the o@xmath1 model on the 3 - 12 @xcite and martini lattices @xcite .    up until recently , a rigorous mathematical proof of the result ( [ mu ] ) had seemed out of reach . in a remarkable paper , duminil - copin and smirnov @xcite have given an ingenious proof of the nienhuis result . the first part of the proof involves defining a parafermionic observable on the honeycomb lattice which is shown to be discretely holomorphic . the boltzmann weights of the model follow from the condition of discretely holomorphicity . in this sense the value ( [ mu ] ) also appears as a consequence of discrete holomorphicity . the second part of the proof builds on previously hard won results in the general theory of self - avoiding walks to rigorously establish that this value is indeed the connective constant . these arguments have been extended @xcite to provide a rigorous proof of the critical surface fugacity @xmath44 of self - avoiding walks on the honeycomb lattice in the presence of a boundary . this is the critical adsorption temperature at which the polymer becomes attached to the surface . for @xmath45 the polymer is adsorbed and for @xmath46 the polymer is desorbed . the critical value depends on the orientation of the lattice ( see fig .  [ fig : orient ] ) .    for the orientation shown in fig . 2(a ) , the exact value @xcite @xmath47 had been obtained earlier using the underlying integrability of the o@xmath1 loop model with a boundary @xcite . for the orientation shown in fig . 2(b ) , the exact value @xcite @xmath48 also followed from the integrable boundary weights of the o@xmath1 loop model @xcite . building on the approach of duminil - copin and smirnov , these results have been shown by beaton _ _ @xcite to also follow from the discrete holomorphicity of a parafermionic observable . their mathematical proof that these values are indeed the critical adsorption temperatures is an impressive _ tour de force _ , involving the adaption of previous rigorous work on self - avoiding walks and the rigorous treatment of self - avoiding bridges in a finite strip . we turn now to further examine the connection between discrete holomorphicity and integrability . the z@xmath0 model is an @xmath49-state model in which the spins @xmath50 at vertex @xmath51 take values from the @xmath49th roots of unity , @xmath52 , with @xmath53 and @xmath54 . the model contains a number of well known models as special cases , namely the ising model ( @xmath55 ) , the three - state potts model ( @xmath56 ) and the ashkin - teller model ( @xmath57 ) . the z@xmath0 model has a well known set of critical points obtained by fateev and zamolodchikov ( fz ) as solutions to the star - triangle relations @xcite . by defining lattice parafermions as products of neighbouring order and disorder variables with a suitable phase factor , it was shown by rajapour and cardy that the fz boltzmann weights ensure that the lattice parafermions obey the discrete version of the cauchy - riemann equations , i.e. , they enjoy the property of discrete holomorphicity @xcite . these results were extended to lattices with anisotropic interactions . the underlying motivation for this work on the z@xmath0 model is the construction of discretely holomorphic parafermions , which are expected to become the holomorphic parafermions of the conformal theory . here our interest is in the connection with integrability . from that perspective , the key finding is that the contour sum around each elementary face or rhombi of the lattice vanishes , the key ingredient for which is the set of fz critical boltzmann weights . this approach was extended by explicitly considering the condition of discrete holomorphicity on two and three adjacent rhombi @xcite . for two rhombi this leads to a quadratic equation for the boltzmann weights which implies the known inversion relations for the z@xmath0 model . for three rhombi it leads to a cubic equation from which the star - triangle relation follows . the simplicity of the discrete holomorphic approach is that the two- and the three - rhombus equations are equivalent to the one - rhombus equation . as alluded to in the previous section , the original motivation behind the recent progress on various aspects of discrete holomorphicity covered in this article was to set up the necessary discrete building blocks for the passage to the continuous scaling limit . there is ample motivation for doing this . for example , the geometry of the curves traced out by self - avoiding walks is conjectured to be described by the mathematical theory of schram - loewner evolution ( sle@xmath58 ) @xcite . basically , if the scaling limit of two - dimensional self - avoiding walks exists and has a certain conformal invariance property , then the scaling limit must be described by the particular stochastic process sle@xmath58 . however , establishing the existence of the scaling limit ( and therefore also its conformal invariance ) is a difficult open problem . nevertheless , significant progress in this direction has been made for other two - dimensional models , most notably percolation and the ising model . @xcite . although the programme to establish the scaling limit of self - avoiding walks in a mathematically rigorous way is not yet complete , the first step , involving establishing analytic properties on the lattice via the discrete cauchy - riemann equations , has led to some remarkable developments . these include the surprising mathematical proofs of the results ( [ mu ] ) , ( [ kappa1 ] ) and ( [ kappa2 ] ) for the connective constant and critical surface fugacity for self - avoiding walks on the honeycomb lattice @xcite . there are also ongoing developments with regard to the connection with integrability , particularly at a boundary . a boundary condition can also be imposed on a discretely parafermionic observable @xcite . this approach has led to a set of new integrable boundary weights for the z@xmath0 model @xcite . several new sets of boundary weights have also been obtained in this way for the o@xmath1 and @xmath59 loop models @xcite . more developments are eagerly awaited . this work has been partly supported by the australian research council .
Abstract:

Answer: over the past few years it has been discovered that an `` observable '' can be set up on the lattice which obeys the discrete cauchy - riemann equations . the ensuing condition of discrete holomorphicity leads to a system of linear equations which can be solved to yield the boltzmann weights of the underlying lattice model . surprisingly , these are the well known boltzmann weights which satisfy the star - triangle or yang - baxter equations at criticality . this connection has been observed for a number of exactly solved models . i briefly review these developments and discuss how this connection can be made explicit in the context of the z@xmath0 model . i also discuss how discrete holomorphicity has been used in recent breakthroughs in the rigorous proof of some key results in the theory of planar self - avoiding walks .


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: once fluidized a vibrated granular medium looks like a gas of particles that can be described using kinetic theory of usual gases . however , granular gases basically differ from ordinary ones mostly due to the inelasticity of collisions , i.e. , nonconservation of energy . while over the years many attempts based on kinetic theory @xcite have been made to describe such dissipative granular gases , no agreement has been found so far both with experiments @xcite and numerical simulations @xcite , for the dependence of the granular temperature , i.e. , the mean kinetic energy per particle , on the parameters of vibration @xcite . the aim of this study is to guess possible gas - like state equations for such a dissipative granular gas and to observe new kinetic behaviors which trace back to the inelasticity of collisions , e.g. , the tendency of such media to form clusters . although this feature has probably been known since the early observation of planetary rings @xcite and although various cluster types in granular flows have been observed numerically @xcite , there exist only a few recent laboratory experiments . one experiment , with a horizontally shaken 2d layer of particles , displayed a cluster formation , but the coherent friction force acting on all the particles was far from being negligible @xcite . we performed a similar experiment by exciting a 3d granular medium with a vertically vibrated piston . we did observe clustering , but we could not rule out a lock - in mechanism involving the time scale connected with gravity and the period of vibration @xcite . we thus repeated this experiment in a low - gravity environment , where inelastic collisions were the only interaction mechanism , and observed a motionless dense cluster that confirms that the inelasticity of collisions alone can generate clustering @xcite . this paper is devoted to the study of the low - density situation , where clustering does not occur . it is organized as follows . section 2 is devoted to the presentation of our experimental setup . the experimental results are presented in secs . we report in sec . 3 ( resp . 4 ) the measurements of pressure ( resp . volume ) of a gas of spherical particles excited by a vibrating piston and undergoing inelastic collisions . at constant external driving , we show that the pressure passes through a maximum for a critical number of particles before decreasing for large @xmath0 . the density of particles as a function of the altitude is studied in sec . 5 , where we observe an exponential atmosphere far enough from the piston . from measurements of secs . 35 , we show in sec . 6 that the dependence of the granular temperature , @xmath1 , on the piston velocity , @xmath2 , is of the form @xmath3 , where @xmath4 is a decreasing function of @xmath0 . finally , we discuss our results in the light of previous works @xcite and give our conclusions in sec . 7 . the experiment consists of a transparent cylindrical tube , with an inner diameter of @xmath5  mm , filled from @xmath6 up to @xmath7 stainless steel spheres , @xmath8  mm in diameter , roughly corresponding to @xmath9 up to about @xmath10 particle layers at rest . an electrical motor , with eccentric transformer from rotational to translational motion , drives the particles sinusoidally with a @xmath11 or @xmath12 mm amplitude , @xmath13 , in the frequency range from @xmath14 to @xmath6 hz . a lid in the upper part of the cylinder , is either fixed at a given height , @xmath15 ( constant  volume experiment ) or is mobile and stabilized at a given height @xmath16 due to the bead collisions ( constant  pressure experiment ) . the heights @xmath15 and @xmath16 are defined from the lower piston at full stroke . when the vibration is strong enough and the number of particles is low enough , the particles display ballistic motion between successive collisions like molecules in a gas ( see fig .  [ fig01]a ) . when the density of the medium is increased , the gaslike state is no longer stable but displays the formation of a dense cluster bouncing like a nearly solid body ( see fig .  [ fig01]b ) . this paper is devoted to the study of the gaslike state . its aim is to determine experimentally a state equation  between pressure , volume , particle number and the vibration parameters ( amplitude and frequency ) .    ; ( b ) @xmath17 , respectively corresponding to roughly @xmath18 and @xmath19 particle layers at rest . the parameters of vibration are @xmath20  hz and @xmath21  mm . the driving piston is at the bottom ( not visible ) , the inner diameter of the tube being @xmath22 mm ] as a function of @xmath0 . from the lower ( @xmath23marks ) to the upper ( hexagrams ) curve , vibration frequency @xmath24 varies from @xmath25 to @xmath6  hz with a @xmath18  hz step . for all these experiments , @xmath26 mm and @xmath27  mm , @xmath28 being the bed height at rest . one single layer of particles at rest corresponds to @xmath29 . lines join the data points ]    time averaged pressure measurements have been performed as follows . initially , a counterwheight of mass @xmath30  g balances the lid mass . the piston drives stainless steel spheres in erratic motions in all directions ( see fig .  [ fig01]a ) . particles are hitting the lid all the time , so that to keep it at a given height @xmath15 we have to hold the lid down by a given force , @xmath31 , where @xmath32 is the mass of a weight we place on the lid and @xmath33 the acceleration of gravity . at a fixed @xmath15 , _ i.e. _ at a constant  volume , fig .  [ fig02 ] shows the time averaged pressure @xmath34 exerted on the lid as a function of the number @xmath0 of beads in the container , for different frequencies of vibration , @xmath35 being the area of the tube cross - section . at constant external driving , _ i.e. _ at fixed @xmath24 and @xmath13 , the pressure passes through a maximum for a critical value of @xmath0 roughly corresponding to @xmath36 particle layers at rest . this critical number is independent of the vibration frequency . a further increase of the number of particles leads to a decrease in the mean pressure since more and more energy is dissipated by inelastic collisions . note that gravity has a small effect in these measurements that are performed for @xmath37 , where @xmath38 is the maximum velocity of the piston . for @xmath0 such that one has less than one particle layer at rest , most particles perform vertical ballistic motion between the piston and the lid . thus , the mean pressure increases roughly proportionally to @xmath0 . when @xmath0 is increased such that one has more than one particle layer at rest , interparticle collisions become more frequent . the energy dissipation is increased and thus the pressure decreases . , as a function of @xmath0 , for various frequencies @xmath24 of vibration . from the lower ( @xmath39marks ) to the upper ( @xmath40marks ) curve , @xmath24 varies from @xmath41 to @xmath6  hz with a @xmath18  hz step and @xmath27  mm . one single layer of particles at rest corresponds to @xmath29 . lines join the data points ] we now consider the bed expansion under the influence of collisions on a circular wire mesh lid placed on top of the beads leaving a clearance of about 0.5 mm between the edge of the lid and the tube . due to the bead collisions , the lid is stabilized at a given height @xmath16 from the piston at full stroke . although the lid mass is roughly @xmath42 times smaller than the total mass of the beads , the lid proves to be quite stable and remains horizontal . the expansion , @xmath43 , of the bed is displayed in fig . [ fig03 ] as a function of @xmath0 for different vibration frequencies . @xmath28 is the bed height at rest . at fixed @xmath24 , the expansion passes through a maximum for a critical value of @xmath0 roughly corresponding to @xmath44 particle layers at rest . this critical number is independent of the vibration frequency . when @xmath0 is further increased , the expansion decreases showing , as for pressure measurements , an increase in dissipated energy by inelastic collisions . note that the height @xmath16 of the granular gas is much larger than for the pressure measurements of fig . [ fig02 ] . consequently , gravity is important here . as already found experimentally in 1d @xcite and 3d @xcite and numerically @xcite , we observe that , at @xmath0 fixed , the granular medium exhibits ( see fig .  [ fig04]a ) a sudden expansion at a critical frequency corresponding to a bifurcation similar to that exhibited by a single ball bouncing on a vibrating plate @xcite . moreover , fig . [ fig04]a shows that this critical frequency depends on the number of layers , @xmath45 . when @xmath45 increases above 0.4 , a transition from the 1d  like behavior to a 3d one is observed : the expansion at the critical frequency becomes less abrupt and tends to increase regularly with @xmath24 ( see figs . [ fig04]a  b ) .    [ cols="^,^ " , ] the aim of this study was an experimental determination of the state equation of dissipative granular gases . it is known that for a fixed number of granular layers at rest , one has in the low - density limit @xmath46 @xcite , where @xmath47 is the mean pressure , @xmath48 the volume and @xmath1 the granular temperature . however , the dependence of @xmath1 on the vibration amplitude , @xmath13 , and frequency , @xmath24 , of the piston and on the number of particles , @xmath45 , is still a matter of debate @xcite . kinetic theory @xcite or hydrodynamic models @xcite show @xmath49 , whereas numerical simulations @xcite or experiments @xcite give @xmath50 , with @xmath51 and @xmath52 , where @xmath38 is the maximum velocity of the piston . these previous conflicting results may be explained by the @xmath2dependence of @xmath1 that we have found from our measurements : @xmath53 , with @xmath4 continuously varying from @xmath54 when @xmath55 , as expected from kinetic theory , to @xmath56 for large @xmath45 . we emphasize that we have only considered the dependence of the granular temperature on the vibration velocity @xmath2 . for the dependence on the number of particles , another term of the form @xmath57 should exist , that we can not determine from these measurements . j.  t.  jenkins , s.  b.  savage , a theory for the rapid flow of identical , smooth , nearly elastic , spherical particles . fluid mech . * 130 * , 187202 ( 1983 ) ; c.  s.  campbell , rapid granular flows . fluid mech . * 22 * , 5792 ( 1990 ) s.  warr , j.  m.  huntley , g.  t.  h.  jacques , fluidization of a 2d granular system : experimental study and scaling behavior . e * 52 * , 55835595 ( 1995 ) s.  luding , e.  clment , a.  blumen , j.  rajchenbach , j.  duran , studies of columns of beads under external vibrations . e * 49 * , 16341646 ( 1994 ) s.  luding , h.  j.  herrmann , a.  blumen , simulations of 2d arrays of beads under external vibrations : scaling behavior . e * 50 * , 31003108 ( 1994 ) s.  luding , granular materials under vibration : simulations of rotating spheres . e * 52 * , 44424457 ( 1995 ) h.  j.  herrmann , s.  luding , modeling granular media on the computer . . thermodyn . * 10 * , 189231 ( 1998 ) j.  lee , scaling behavior of granular particles in a vibrating box . physica a * 219 * , 305326 ( 1995 ) s.  mcnamara , s.  luding , energy flows in vibrated granular media . e * 58 * , 813822 ( 1998 ) v.  kumaran , temperature of a granular material fluidizedby external vibrations . e * 57 * , 56605664 ( 1998 ) j.  m.  huntley , scaling laws for a 2d vibro - fluidized granular material . e * 58 * , 51685170 ( 1998 ) p.  goldreich , s.  tremaine , the dynamics of planetary rings . . astrophys . * 20 * , 249283 ( 1982 ) m.  a.  hopkins , m.  y.  louge , inelastic microstructure in rapid granular flows of smooth disks fluids a * 3 * , 4757 ( 1991 ) ; s.  mcnamara , w.  r.  young , inelastic collapse and clumping in a 1d granular medium . fluids a * 4 * , 496504 ( 1992 ) ; i.  goldhirsch , g.  zanetti , clustering instability in dissipative gases . phys * 70 * , 16191622 ( 1993 ) a.  kudrolli , m.  wolpert , j.  p.  gollub , cluster formation due to collisions in granular material . * 78 * , 13831386 ( 1997 ) . falcon , s.  fauve , c.  laroche , cluster formation , pressure and density measurements in a granular medium fluidized by vibrations . j. b * 9 * , 183186 ( 1999 ) . falcon , r.  wunenburger , p.  vesque , s.  fauve , c.  chabot , y.  garrabos , d.  beysens , cluster formation in a granular medium fluidized by vibrations in low gravity . lett . * 83 * , 440443 ( 1999 ) c.  e.  brennen , s.  ghosh , c.  r.  wassgren , vertical oscillation of a bed of granular material . j. appl . mech . * 63 * , 156161 ( 1996 ) m.  l.  hunt , s.  hsiau , k.  t.  hong , j. fluid eng . * 116 * , 785791 ( 1994 ) y.  lan , a.  d.  rosato , macroscopic behavior of vibrating beds of smooth inelastic spheres . phys . fluids * 7 * , 18181831 ( 1995 ) d.  a.  kurtze , d.  c.  hong , effect of dissipation on density profile of 1d gas . physica a * 256 * , 5764 ( 1998 ) s.  mcnamara ( 1998 ) private communication
Abstract:

Answer: we report experimental results on the behavior of an ensemble of inelastically colliding particles , excited by a vibrated piston in a vertical cylinder . when the particle number is increased , we observe a transition from a regime where the particles have erratic motions ( granular gas ) to a collective behavior where all the particles bounce like a nearly solid body . in the gaslike regime , we measure the pressure at constant volume , and the bed expansion at constant external pressure , as a function of the number @xmath0 of particles . we also measure the density of particles as a function of the altitude , and find that the atmosphere  is exponential far enough from the piston . from these three independent measurements , we determine a state equation  between pressure , volume , particle number and the vibration amplitude and frequency .