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: astronomers are limited to simply observing the physical universe , and one of the key tools is recording the brightness of some object of interest as a function of time . this technique of _ time - series photometry _ has yielded many important results over past decades , particularly in the area of optical astronomy . these successes include the study of eclipsing binary stars , pulsating stars and other types of stellar variability .    astronomy research at victoria university of wellington ( vuw ) over the past few decades has involved the study of pulsating white dwarfs using time - series photometry techniques . a succession of photomultiplier - based instruments@xcite for obtaining light curve data has evolved into an efficient and fast ccd ( charge coupled device ) camera system that we call _ _ puoko - nui__@xcite . there are in fact two flavours of the instrument : _ puoko - nui south _ , based at mt john university observatory in new zealand and _ puoko - nui north _ based at the university of texas in austin , usa . as part of the development of this ccd photometer , a light box with controllable light emitting diodes ( leds ) as artificial stars was constructed for bench testing the instrument . a very useful addition was the provision of microprocessor control of the led intensities , which enabled the creation via software of a variety of simulation light curves . it is basically a modified version of this testing apparatus that we describe here . for the laboratory instrumentation , we have replaced the research grade frame transfer ccd system with a much cheaper commercial ccd camera from the sbig@xcite corporation ( sbig st-402me ; the same variety of ccd that we use for autoguiding the telescope with _ puoko - nui south _ attached ) , and the custom led controller was redesigned using an arduino@xcite board and through - hole electronic components in order to simplify reproduction possibilities . the particular class of pulsating star of interest here , and the primary motivation for our laboratory simulation experiment , is that of the pulsating white dwarf stars@xcite . these stars exhibit luminosity variations with timescales in the range of @xmath0s to @xmath1s . monitoring these variations and extracting frequency information from the resulting light curve enables the observer to identify pulsation ( normal ) modes of the star , and thereby provide information about its otherwise hidden internal structure . the technique is called _ _ asteroseismology__@xcite due to its similarity to seismic studies of the earth s internal structure . there are a number reasons why simulating the photometric study of the pulsating white dwarf ( wd ) stars yields a quality laboratory experiment . first , isolated pulsating wds exhibit extremely coherent frequencies , and the analysis of the resulting light curves provides an excellent introduction to the usefulness of fourier techniques in such scenarios . second , the hundreds - of - seconds periodicities observed in real stars means that the laboratory experiment can be operated with timescales matching or nearly matching an actual observing run at a telescope . third , simulating the light output variations from a real astrophysical object , as distinct from simply creating a variable signal to analyze , yields a less artificial experiment and permits natural connections to the underlying physics of the object . and , fourth , use of a modern ccd camera provides a direct introduction to a key instrument that dominates the current world of optical observational astronomy .        in the following sections we describe in detail the experimental apparatus and the analytical tools required to extract the periodicities that are present in the detected time - series data . we also provide a few insights into the physics of pulsating wds . the experimental apparatus ( pictured in fig .  [ figure : experimentalhardware ] ) consists of a light - tight box with the sbig ccd camera and a small lens at one end imaging a pair of leds attached to a plate on the inside of the box at the other end . this configuration allows the ccd to be removed for general inspection , and to adjust the relative position of the lens in order to optimize the focus . a key part of the apparatus is the ` star simulator ' unit ( pictured in fig . [ figure : starsimphoto ] ) that drives the leds with simulated intensity profiles corresponding to the behaviour of chosen astronomical objects . it is attached to a panel of the box at the opposite end to the ccd and connects to the leds which are on the inside . the core of the simulator electronics is an arduino@xcite microcontroller , with a custom daughter - board ( made from an off - the - shelf prototyping ` shield ' and through - hole electronic components ) that holds the led drivers . a dc - dc converter outputs @xmath212v to power the drivers , and an output relay is included such that one of the output channels can be redirected to an external led , which can be viewed directly by eye outside the box . [ figure : starsimblockdiagram ] gives a schematic overview of the star simulator hardware . each of the two led drivers are built around an ad711 operation amplifier ( op - amp ) configured as a constant current source . in this configuration , the output current is determined by the resistance between the anode of the led and ground . the output driver provides four software - selectable currents ( 5@xmath3a , 50@xmath3a , 500@xmath3a , 5ma ) by connecting four parallel resistors ( with different values ) in series with individual bs107 field - effect transistors . these transistors are connected to digital output pins on the arduino , and act as software - controlled switches that can be switched set the effective resistance to ground . all four paths run through a fifth transistor , which is connected to one of the pulse width modulation ( pwm ) outputs on the arduino . this transistor is used to rapidly switch the led on and off ( using a 10-bit timer at 250hz ) to provide 1024 fine - grained intensity levels between zero and the maximum intensity set by the current source . the combination of current and pwm control provide the large dynamic range that is required for our simulations . the software programming ( firmware ) inside the arduino is relatively straightforward . when the unit is powered on it reads the desired simulation parameters from read - only memory , and configures the current selector transistors and the pwm timers . the unit then enters a mode where it will evaluate and update the led intensities every 16ms and checks for new commands from the usb connection to the control pc . we have developed four types of simulation profiles for our experiment : ( a ) a constant intensity ; ( b ) a sum of multiple sinusoids ( with arbitrary amplitude , frequency and phase ) ; ( c ) a sum of multiple gaussians ( with fixed repeating period , but arbitrary amplitude , width , and offset within a cycle ) ; ( d ) a repeating linear ramp ( with configurable period ) . these profiles provide the building blocks to define a range of simulations ( such as the examples described in section [ section : example ] ) and apparatus tests . in addition to these profile types , a cloud simulation can be enabled on a per - profile basis to attenuate both output channels in a semi - random manner . a simple random number generator ( based on the relative drift of two on - board oscillators ) seeds control point values for a spline interpolation algorithm , which results in a smoothly varying attenuation as a factor of time . additional per - point noise could be included to provide a more realistic simulation , but we omit this for the sake of experimental simplicity . the profile parameters for each simulation are defined in the firmware , and uploaded to the unit via its usb connection . the active profile can be changed at any time by running a simple configuration tool on the acquisition computer . the main steps required in the experiment are ( a ) acquisition of a sequence of ccd frames imaging the illuminated leds , ( b ) extraction of the led intensities from each of the frames . ( c ) conversion of the led intensities to form an appropriate time series for fourier analysis , and ( d ) performing the fourier analysis on various time series light curves and interpreting the results . we discuss each of these steps in turn . after selection of one of the simulation light curve options , the first step is to acquire a sequence of @xmath4 ccd frames using a chosen sampling cycle time of @xmath5s . each frame will have some exposure length @xmath6s , which will typically be shorter than @xmath5s by some amount due to the ccd readout time and possibly other hardware constraints . with the sbig ccd ( st-402me ) we use in our laboratory experiment , the minimum value for @xmath5 is about 1s , but the exposure times can be much shorter as the device uses a mechanical shutter ( which is standard ) to block the light input during the readout phase .    clearly some thought needs to be put into deciding on the these values for a chosen simulation run such that the nyquist criterion ( a minimum of two samples for the shortest period ) is met for the highest frequency in the light curve , and the run length @xmath7s is adequate to provided sufficient frequency resolution in the fourier analysis .    in order to provide uniformity , a conventional ccd chip will measure the stored charge ( formally electron - hole pairs in the si semiconductor substrate ) in each pixel serially , using a single amplifier and adc ( analogue to digital converter ) . this conversion can take quite some time ( seconds or more ) , especially when operated at slow speeds in order to reduce measurement noise ; consequently , exposure intervals even in the many seconds domain become rather inefficient and one can end up closing the telescope to scarce photons for a significant fraction of an observing run .    the _ _ puoko - nui__@xcite photometers employ frame transfer ccds , which have adjacent masked pixel arrays of the same size that act as storage regions . the integrated charges in the exposure frames can be rapidly shifted to the storage frames ( in milliseconds ) and the frames can then be digitized and read out in parallel with the next exposure . thus the ccd photometers can be operated without mechanical shutters , and the deadtime losses become negligible .    the laboratory simulation experiment of course does nt suffer from a deficit of photons so the standard ccd is perfectly adequate for the task in hand . the sbig ccd has a fast ( @xmath8ms ) readout time , and so the readout dead time is not of practical concern for this experiment . thinking in terms of led brightness variations with periods of @xmath0s to match the real white dwarf variations , exposure times of @xmath9s or less are appropriate and provide sample rates significantly higher than the minimum nyquist rate corresponding to the highest frequency variation in the light curves .    for their use in astronomical observing , ccds need to be cooled below ambient temperatures in order to minimize the thermal creation of electron - hole pairs that contribute to the background . this can now be mostly achieved thermoelectrically ; for example when observing , we cool our _ puoko - nui south _ primary ccd to @xmath10c . for the laboratory experiment , with control over the brightness of the leds , this can be avoided , but in the interests of alerting students to the issue it is sensible to implement cooling . our sbig ccd can be thermoelectrically cooled and we typically operate it at @xmath11c .    another factor which should nt remain hidden from students is the issue of _ dark _ frames . these are exposures with the shutter closed and employing the same time interval as the observation frames : they are subtracted from the latter to produce the _ science _ frames . with our sbig ccd , the manufacturer - provided acquisition software includes options to automate this procedure . in real astronomical observing , after dark frame subtraction , the remaining background is dominated by the sky brightness , which will vary with the conditions ; for the experiment in a light - tight box , this should in theory be zero . in practice , however , a finite charge is added to each pixel before digitization in order to ensure better linearity . zero time exposures to measure these values are known as _ bias _ frames in the ccd literature . the dark subtraction also removes this bias charge . this fact should at least be mentioned in any notes in order to satisfy questions from any inquisitive students . it is perhaps relevant to point out that when using our frame transfer ccd in shutterless mode , we need to accumulate a separate library of suitable ccd dark frames outside an observing run . although it is not a required step in the laboratory experiment , it would be appropriate to at least mention the part played by `` flat fields '' in ccd astronomical photometry . invariably because of the telescope optics , an external uniformly bright field does not produce a uniform pixel count across the ccd . this can be accounted for by observing a physical flat field in high precision and scaling the science frame pixel counts accordingly . in practice , a flat field is obtained in astronomy by observing the bright clear sky at dusk or dawn ( and removing any stellar images ) or a uniformly illuminated board in the telescope dome . ( and we will simply mention the debate concerning which is better ) .    finally , the brightest intensity of the leds should be monitored , and the exposure time should be chosen such that the image intensity will not saturate . the adcs used to digitize the pixel charges have a finite number of digital bits . in the case of our sbig ccd , a 16 bit adc is used yielding a maximum count of @xmath12 adu ( analog to digital units ) . in practice it is best to maintain the signal of interest at a level below this maximum value . the number of exposures collected for the time series should ensure that a reasonable number of complete cycles of the longest period are covered and understanding the significance of this can also be part of the experiment .      for the real observing situation with _ puoko - nui _ we have written code that employs circular synthetic apertures to extract stellar intensities from each exposure frame . the software package , ` tsreduce`@xcite , uses two concentric apertures centered on each stellar image of interest to obtain the stellar intensity above the sky background , and the developing light curves are presented in real time as the observing proceeds . the top section of fig . [ figure : reduction ] illustrates the technique using a real observational science frame of the pulsating white dwarf qu tel@xcite ( also known as ec20058@xmath135234 ) . it was obtained as part of an observing campaign in 2011 with the _ puoko - nui south _ photometer attached to the mt john one metre telescope in new zealand . the top encircled image corresponds to the pulsating white dwarf , while the other object acts as a comparison .    for the simulation experiment the optics are much simpler and one need not be unduly concerned with the circular symmetry of telescope optics and the resulting point spread function ( including the effects of astronomical `` seeing '' at the observing site ) . instead it is both simpler and adequate to employ synthetic square apertures to determine the integrated led flux and the background value ( refer to the bottom section of fig . [ figure : reduction ] ) . a close to optimum estimate of the led intensity per frame is obtained by subtracting the summed pixels in the background aperture from the summed pixels in the led enclosed aperture . this approach is straightforward to implement in a computing environment such as matlab or python , and if deemed appropriate the students themselves can develop their own extraction code as part of the experimental method . it is relevant to point out that the synthetic aperture technique ( square or circular ) for extracting intensities is only really suitable for sparse stellar fields ; with the aim of making connections where appropriate with real observing , we simply mention here the techniques of `` point spread function '' fitting and `` difference imaging '' for conducting photometry in more crowded fields . we refer the interested reader to the extensive literature for explanations of these techniques . the information of interest in the light curves is encoded in the variations about the mean level . the actual detected intensity is rarely of interest in the world of astronomy and one is almost always required to account for the effects of flux variations created by cloud interference and changing airmass . one approach to making these corrections is to compute the normalized ratio of the variable star intensity to a nearby ( constant ) comparison star . this situation can be simulated if required in the experiment . the appropriate form for the time - series data prior to performing a fourier analysis is to compute modulation values : percentage variation in intensity units are an obvious scale to use , but millimodulation intensity ( mmi ) units are often used in the white dwarf pulsation literature as they are suitable for the typical @xmath14% intensity variations of these stars : 10mmi corresponds to a 1% modulation in the time domain . units consistent with these in fourier amplitude space ( see below ) are millimodulation amplitude ( mma ) values . furthermore , for periodicities in the @xmath15s region , @xmath16hz or mhz are suitable for frequency values . we will use these units in this paper .      starting with the time series data featuring modulation intensities versus time , this step involves `` simply '' computing a discrete fourier transform ( dft ) of the data . however , given the simulation data to be analyzed can include multiple frequencies ( with close spacings ) and has high frequency stability and coherence ( that match the observed white dwarf pulsations ) , there are a number of features and techniques that can be exploited that will give students good insights into the use of fourier analysis . we will simply summarize the main features here and allude to their use in the world of pulsating white dwarf asteroseismic research . first , it is important to note that the analysis will be using a discrete version of the fourier transform , and to be aware of its properties and limitations . the phase information is of no interest in wd asteroseismic research , so computing the sinusoidal power ( the power spectrum ) will yield the amplitude of the periodicities ; and , as it is best to think in terms of amplitudes rather than power ( see below ) we require the square root of the power spectrum , which we will call the fourier amplitude spectrum or simply dft henceforth . no doubt in most circumstances the dft will be obtained by a `` black box '' function ( by using matlab routines for example ) , but it is important to make clear in student minds the difference between a dft and the fft ( fast fourier transform ) . the fft is nothing more than an efficient way of computing a dft  with definite limitations . for an fft algorithm the samples must be evenly spaced , and in order to achieve maximum efficiency the number of samples needs to be a power of 2 ( which is often achieved in practice by zero padding ) . furthermore , the number of points in the frequency domain is dictated by the number of time domain samples ; this can be a limitation when inspecting a fft of real data , particularly with noise included .    given the speed of modern computers we have found in practice that it is much more useful to use the most simple ( but inefficient ) dft algorithm , with direct control over the number and spacing of the frequency domain estimates . this is especially so and essential when combining different data sets that have different sampling intervals and times , as sometimes occurs in multi - site observing campaigns run by groups such as the whole earth telescope ( wet)@xcite . second , we need to introduce the concept of the dft window and , for real life observing campaigns , generalize the concept somewhat . in the general fourier literature the dft standard window function is the sinc ( @xmath17 ) function ( actually sinc@xmath18 for power ) , which is the fourier transform of a rectangular `` top - hat '' function that defines the finite length of the time domain sample . the dft of the actual data is then a convolution of this sinc function and the fourier transform of the theoretically infinite sinusoids in the data . in some applications , tapering the `` top - hat '' ends in order to reduce the amplitude of the wiggles in the sinc function ( at the expense of broadening the central peak ) is considered effective . this is called _ apodisation _ in the literature . however , in the world of pulsating white dwarfs in particular , with the occurrence of unavoidable data gaps due to such things as cloud interference and the rising of the sun , what is needed is a direct transformation of the signal turning on and off ( often multiple times ) into the fourier domain . this is achieved by computing the dft of a single noise free sinusoid with a suitable period calculated for each sample time for the real data . this function can be compared with various peaks in the dft of the real data in order to see if they are consistent with only one frequency . even if the experiment is limited to investigating the noise - free simulation light curves , it is instructive to compute the window function when examining light curves with multiple gaps and closely spaced frequencies .    finally , in the interests of completeness and making connections with real observational data , it is relevant to discuss the impact of noise on the frequency identification process . as part of generating the light curves , observational uncertainties or noise can easily be added using for example an algorithmic ( pseudo-)random number generator . the issue of identifying as real low level peaks in the dft computed from an observational light curve is an important part of the analysis in observational astronomy . a technique that has proved to be useful , especially in the case of low amplitude periodicities that are close to large amplitude ones , is referred to as _ _ prewhitening__@xcite . one identifies the periods of the larger amplitude peaks in the dft and then uses these to compute a synthetic light curve . this synthetic light curve is then subtracted from the real time - series data and another dft calculated for the now `` prewhitened '' data . this new dft should provide a clearer indication of the reality of any remaining peaks above any nearby noise peaks . if the periodicities @xmath19 of the peaks being subtracted are known accurately enough , they can be fixed in a straight forward _ linear _ least squares procedure by fitting values for @xmath20 and @xmath21 in functions of the following form    @xmath22         + b_{i}\cos\left[\frac{2\pi t}{t_{i}}\right ]        = c_{i}\sin\left[\frac{2\pi t}{t_{i } } + \phi_{i}\right]\ ] ]    alternatively , the @xmath19 values can also be optimized using _ nonlinear _ least squares techniques , which require iterative procedures . in this case , it would more convenient to also optimize the parameters @xmath23 and @xmath24 in the above functions . it is instructive to compare the amplitudes @xmath23 of the least squares derived sinusoids with those in the dft amplitude spectrum . of course something would be radically wrong if they were not very similar , but making the comparison reinforces the basic fact that the least squares estimates and the dft amplitudes _ should _ be very similar even though they are derived from apparently different procedures . with this comparison in mind , it is useful to view the dft as a set of cross correlations of the time - series data with a sequence of sinusoids with different frequencies . when deciding whether a peak in a dft corresponds to a real periodicity in the light curve , it is important to have a statistical model for the time domain noise and its effect on the uncertainties in the dft . with the speed of modern computers , simulation procedures can make this task largely model - free . we are attracted to a simulation implementation@xcite of the false alarm probability ( fap ) concept@xcite for characterizing uncorrelated noise levels in dfts . this method is intuitively straight forward and essentially free of complicated statistical theory . it is certainly preferable to simply quoting some formula for the dft noise level . one first _ prewhitens _ ( removes ) the main periodicities from the light curve , leaving residuals that characterize the time - domain noise level . the time stamps for each of the sample values are then randomly shuffled , thus destroying any remaining coherent signals in the data , but leaving the random noise level intact . a dft of the time - shuffled data is calculated and the height of the highest peak in the region of interest is noted . this procedure is repeated n times and a record kept of each highest peak , resulting in an ensemble of highest peaks . one then asserts that there is a one in @xmath4 chance of the highest ensemble peak occurring through a noise conspiracy in the data . for @xmath25 say the chance or fap is 1 in 1000 . thus , if a peak is visible in the dft at this level , there is 0.1% chance of it simply being noise . it is instructive to construct a histogram of the ensemble of highest peaks obtained in the simulation  some examples can be found in the references@xcite . the implementation and effectiveness ( along with some explanations ) of the above techniques can be found in a papers describing observational programs on pulsating white dwarfs@xcite . in particular , an implementation of a fap simulation procedure should really make use of the much more efficient fft algorithm given the requirement to compute a large sample of dfts . in this section we provide three example cases with different light curves . the first one involves sampling a light curve that features beating between two sinusoidal signals . this provides a good starting example as it provides a relatively simple signal while the students are familiarizing themselves with the general operation of the camera and the data reduction procedures . it also illustrates the main features associated with sampling a signal with two adjacent sinusoidal signals over a finite interval . [ figure : beating ] depicts both the set of modulated time series samples covering an interval of 600s and the resulting dft . the dft window function for this sample is also shown in the insert . the second example employs a light curve generated using multiple frequencies corresponding to those observed in the pulsating white dwarf qu tel@xcite ( also known as ec20058@xmath135234 ) . cloud simulation ( affecting both led intensities ) is also added to the mix and the correct light curve variations recovered by computing intensity ratios after background estimates have been subtracted ( refer to fig.[figure : ec20058cloud ] ) . [ figure : ec20058dft ] compares a 10h run using the qu tel simulated light curve with a real @xmath810h 2004 observing run on this pulsating white dwarf carried out at mt john observatory in new zealand @xcite . included in the figure is a comparison of the two dfts .     to accelerate the data collection time to only 1 hour . ] the final example illustrates the principles of signal folding , and provides an impressive example of the ability to reconstruct periodic signals from sparsely sampled data when there is a high level of signal coherency . the simulation profile is based on optical observations of the pulsar ( psrb1919 + 21 ) in the crab nebula that were acquired in 2013 using the _ puoko - nui north _ photometer@xcite . a pulsar is a rotating neutron star that emits two beams of radiation such that if one of the beams intersects the earth it is detected ( primarily ) in the radio frequency domain as an extremely stable periodic signal . the crab pulsar is one of a handful that have also been detected optically . the first crab pulsar optical detections were obtained via photometry synchronized to the known radio ephemeris @xcite , but one can also achieve the detection using signal folding as shown in this simulation example . because the sbig camera has a maximum frame rate of 1 frame per second , the 33ms period of the simulation has been increased by a factor of 100@xmath26 to a more manageable @xmath27 ; this yields roughly three frames sampled per pulsation cycle . figure [ figure : crabdft ] shows the raw time series acquired in a 10minute collection run and the dft of this data . the pulsation profile is non - sinusoidal ( it is modelled using two gaussians ) , and so the main features in the frequency domain are the fundamental pulsation period and its harmonics . the shape of the pulsation profile in the time domain are encoded in the relative amplitude and phase of the harmonics in the frequency domain . the dft plot shows another important aspect of fourier analysis : frequencies that are greater than the nyquist frequency ( 500mhz for these 1s samples ) are ` aliased ' back inside the dft range by reflection around either the nyquist frequency ( at the high end ) or zero ( at the low ) . it can be difficult to distinguish aliases from true signals , but in this case we can uniquely identify the fundamental frequency based on the known time - domain period ( which can be timed manually using the external led on the lightbox ) , and it is then a straightforward exercise to identify the remaining peaks by taking integer multiples of the fundamental frequency and applying the appropriate reflections . the precise period of the pulsar signal can be identified from the dft plot , and this period is then used to ` fold ' the signal over itself in order to create a detailed form for a single pulsation cycle . this can be represented mathematically using the equation    @xmath28    where @xmath29 denotes the integer division operator ( divide and discard the remainder ) . in other words , we are converting the time - based signal to a phase ( which may be itself measured as a time , angle , or fractional number ) . this is shown in figure [ figure : crabfold ] , which compares the simulated pulsar photometry against our real observational data . we have outlined a range of topics connected with the observation of pulsating white dwarf stars that can be explored using our simulation apparatus . depending on the course level that the experiment might be offered in , one would want to make a selection of the topics discussed . during the early development of this experiment we used it very successfully as the basis of a one semester project in 2013 in our 400-level honours course in physics , and in 2014 it was run as a 200-level experiment . ideally , a laboratory experiment using the simulator should endeavor to draw as many links as possible to the real world of astronomical observing and the physics of the target objects . the pulsating white dwarfs are certainly a suitable subject for this task . when white dwarfs pulsate ( at several well - defined surface temperature regimes ) they are usually multiperiodic and exhibit very stable frequencies . this is largely a consequence of their compact very slowly evolving structure  about the mass of the sun compressed into a sphere about the size of the earth . they are prevented from gravitational collapse by the temperature - independent degeneracy pressure resulting from the essentially separate electron gas in the plasma being gravitationally confined at very high number densities . the observed luminosity variations result from nonradial _ gravity mode _ pulsations  the more common and intuitively simpler radial pulsations due to pressure variations ( seen in many other stellar types ) have never been observed in white dwarfs . the underlying driving mechanism for stellar pulsation is the interaction of the escaping core radiation with a partially ionized chemical element ( usually h or he ) in a layer near the star s surface . this mechanism was first postulated by eddington@xcite early last century .    for white dwarfs , the very high surface gravities clearly inhibit radial material motion and hence limit any consequent luminosity variations . in addition , the large bulk modulus ( high `` stiffness '' ) of the material leads to high sound speeds and hence to expected periodicities of seconds or less for the pressure modes . these oscillations , if they exist , are likely to be very difficult to detect@xcite .    the gravity modes , however , are the result of material density gradients in the local gravitational field and feature much longer periods : the so - called brunt - visl frequency@xcite determines these periods and its calculated values are consistent with the hundreds of seconds luminosity variations seen in the pulsating white dwarfs . a simplified and instructive physical picture for the phenomenon is to envisage the oscillations of an object floating at the interface between two fluids of different density ( eg water and air ) . however , in white dwarfs the gravity modes involve largely nonradial material motion and a rich array of possible normal modes . it is interesting to note that the gravity mode pulsations in white dwarfs were discovered@xcite before a theoretical model was proposed . a final note on the asteroseismic study of white dwarfs is worth making and should be at least mentioned in any experimental notes . the detected frequencies in a pulsating white dwarf are used to infer normal modes of the star and these are used to guide and constrain theoretical models of the structure . including some of these steps in an experiment is possible , but it would take us beyond the scope of this article .    obviously , the apparatus could also be used to simulate other observational programs by generating appropriate synthetic light curves : eclipsing binary stars , planetary transit events , or even gravitational microlensing events involving extrasolar planets orbiting stars@xcite are interesting possibilities .    from an observer s perspective , simulating the observation of the short period binary object ny vir ( also known as pg1336@xmath13018)@xcite provides an interesting experience : it has a @xmath82.5h period , a deep primary eclipse , a smaller secondary eclipse combined with a ` reflection effect ' and one component is also a pulsating ( subdwarf b ) star . a number of classic binary star light curve characteristics are therefore in evidence , but more detailed analysis would require some advanced modelling . also , the detection of extrasolar planets via transits is certainly of current topical interest@xcite and monitoring and analyzing a simulation light curve corresponding to a relatively rapid planetary transit@xcite could make a useful laboratory experiment . however , we think that the optimum and least artificial use of the equipment is to produce and analyse light curves corresponding to pulsating white dwarfs , develop expertise in fourier techniques and also explore at least some of the physics of these exotic yet ubiquitous astrophysical objects . d.  j.  sullivan , t.  s.  metcalfe , d.  odonoghue et al . , `` whole earth telescope observations of the hot helium atmosphere pulsating white dwarf ec20058@xmath135234 '' , mon.not.r.astron.soc . , * 387 * , 137152 ( 2008 ) . r.  silvotti , g.  fontaine , m.  pavlov , t.  r.  marsh , v.  s.  dhillon , s.  p.  littlefair and f.  getman , `` search for p - mode oscillations in da white dwarfs with vlt - ultracam . i. upper limits to the p - modes '' , astron . & astrophys . , * 525 * , a64 ( 2011 ) . d.  kilkenny , b.  y.  welsh , c.  koen , a.  a.  s.  gubis and m.  m.  kotze , `` a search for p - mode pulsations in white dwarf stars using the berkeley visible imaging tube detector '' , mon.not.r.astron.soc . , * 437 * , 18361839 ( 2014 ) . an interesting gravitational microlensing planetary event : b.  s.  gaudi , d.  p.  bennett , a.  udalski et al . , `` discovery of a jupiter / saturn analogue with gravitational microlensing '' , science , * 319 * , 927930 ( 2008 ) .
Abstract:

Answer: we have designed a realistic simulation of astronomical observing using a relatively low - cost commercial ccd camera and a microcontroller - based circuit that drives leds inside a light - tight box with time - varying intensities . as part of a laboratory experiment , students can acquire sequences of images using the camera , and then perform data analysis using a language such as matlab or python to : ( a ) extract the intensity of the imaged leds , ( b ) perform basic calibrations on the time - series data , and ( c ) convert their data into the frequency domain where they can then identify the frequency structure . the primary focus is on studying light curves produced by the pulsating white dwarf stars . the exercise provides an introduction to ccd observing , a framework for teaching concepts in numerical data analysis and fourier techniques , and connections with the physics of white dwarf stars . ( submitted to the american journal of physics . )


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 this paper we consider the problem of online update of a pareto archive composed of mutually non - dominated solutions when a new candidate solution shows up . this problem is also referred to as dynamic non - dominance problem  @xcite . online update of a pareto archive is typically used in multiobjective metaheuristics ( momhs ) , e.g. multiobjective evolutionary algorithms , whose goal is to generate a good approximation of the pareto front . most of the nowadays momhs use an external archive of potentially efficient solutions , i.e. pareto archive containing solutions not dominated by any other solutions generated so far . we consider here momhs that generate iteratively new candidate solutions and use them to immediately update pareto archive . updating pareto archive with a new solution @xmath0 means that :    * all solutions dominated by @xmath0 are removed from pareto archive , * @xmath0 is added to pareto archive if it was non - dominated w.r.t . to any solution in pareto archive . online update of a pareto archive may also be used as an intermediate step of other methods . for example drozdk et al . @xcite propose to use the online updated pareto archive to speed - up non - dominated sorting procedure used in many multiobjective evolutionary algorithms . the time needed to update pareto archive generally grows with growing number of objectives and growing number of solutions . in some cases it may become a crucial part of the total running time of a momh . the simplest data structure for storing pareto archive is a plain list . when a new solution @xmath0 is added , @xmath0 is compared to all solutions in pareto archive until either all solutions are checked or a solution dominating @xmath0 was found .    in order to speed - up the process of updating a pareto archive some authors proposed the use of specialized data structures and algorithms , e.g. quad - trees . however , the results of computational experiments reported in the literature are not conclusive and in some cases such data structures may in fact increase the update time compared to simple list .    the direct motivation for us to revisit the issue of updating pareto archive comes from special properties of pareto local search ( pls ) method  @xcite . pls works directly with the pareto archive . for each solution in the archive its neighborhood is searched to find new potentially efficient solutions and update the archive . standalone pls starting from random solutions is very inefficient since it spends a lot of time generating large numbers of solutions being still very far from the pareto front . pls , however , is used as a crucial component in some of the best methods for multiobjective knapsack  @xcite , biobjective traveling salesperson problem ( btsp )  @xcite and set covering problem  @xcite . the general idea of such methods is to start pls from a set of high quality solutions generated by some other methods , e.g. the powerful lin - kernighan heuristic for tsp  @xcite . pls has also been used in other contexts , e.g. for solving scheduling problems  @xcite , multi - agent problems  @xcite or multiobjective markov decision processes  @xcite . the importance of pls may be explained by the fact that unlike most other momhs it is especially good in the search for new solutions along pareto front . in order to generate a good approximation of the pareto front a momh has to achieve two different goals : the solutions generated by such method should both approach the pareto front and cover all areas along this front . in other words , the momh has to search both towards and along pareto front  @xcite . pls seems to be especially good in search along pareto front , i.e. when started from a sample of solutions located very close to the true pareto front , pls may generate a large number of new potentially efficient solutions in a very short time .    from the point of view of updating pareto archive , pls has some special characteristics :    * pls generates large number of new candidate solutions in a very short time ( as the neighborhood solutions for many problems may be quickly generated ) . in particular , before a neighbor solution is added to the archive only the new values of objectives are needed . these values may often be efficiently calculated using changes of the objectives . for example , in the case of tsp , the typical two edges exchange move requires only removing two edges and adding two other edges to the solution . so , the calculation of change of each objective requires just few arithmetic operations . this means that the time needed to generate a new solution may be much shorter than the time needed to update the archive . * since the new candidate solutions are neighbors of potentially efficient solutions they are usually quite good , i.e. even if the new solution is dominated there are only relatively few dominating solutions . in this case the plain list is especially inefficient since the first dominating solution will on average be found after checking large fraction of solutions . * pareto archive is a crucial element of pls . new solutions are accepted or rejected based on the comparison to the archive . so pareto archive has to be updated online , while in some other momhs it is at least potentially possible to store all candidate solutions and only at the end to filter out the dominated ones . based on these motivations , we propose a new method , called nd - tree , for updating a pareto archive . a thorough experimental study will show that we can obtain substantial computational time reductions comparing to state - of - the - art methods . the remainder of the paper is organized as follows . basic definitions related to multiobjective optimization are given in section 2 . in section 3 , we present a state of the art of the methods used for online updating a pareto archive . the new method nd - tree is described in section 4 . finally , computational experiments are reported and discussed in section 5 . we consider in this section a general multiobjective combinatorial optimization problem , defined as follows : @xmath1 @xmath2    a combinatorial structure is associated to this problem , which can be path , tree , flow , tour , etc . we denote by @xmath3 the feasible set in decision space , defined by @xmath4 . the image of the feasible set in objective space is called @xmath5 and is defined by @xmath6 . due to the contradictory features of the objectives , there does not exist a feasible solution simultaneously minimizing each objective ( that is why the word is placed between quotation marks ) but a set of feasible solutions called _ efficient solutions_. we present below some definitions that characterize these efficient solutions . we first define dominance relations :    pareto dominance relation : we say that a vector @xmath7 _ dominates _ a vector @xmath8 if , and only if , @xmath9 . we denote this relation by @xmath10 . we can now define an efficient solution , a non - dominated point , the efficient set and the pareto front . efficient solution : a feasible solution @xmath11 is called _ efficient _ if there does not exist any other feasible solution @xmath12 such that @xmath13 . for the sake of simplicity we will use the dominance relation w.r.t . solutions as well , i.e. @xmath14 .    non - dominated point : the image @xmath15 in objective space of an efficient solution @xmath16 is called a non - dominated point . efficient set : the efficient set denoted by @xmath17 contains all efficient solutions . pareto front : the image of the efficient set in @xmath5 is called the pareto front ( or non - dominated frontier / set ) , and is denoted by @xmath18 .    since in many places we will need dominance or equality relation we define it as coverage relation . coverage relation : we say that a point @xmath19 _ covers _ a point @xmath20 if @xmath10 or @xmath21 . we denote this relation by @xmath22 . we will also use the coverage relation w.r.t . solutions as well , i.e. @xmath23 .    please note that coverage relation is sometimes referred to as weak dominance  @xcite . the terminology , however , is not consistent and other authors use weak dominance as the synonym of pareto dominance  @xcite ( in contrast to strong dominance where a solution is better on all objectives ) . thus to avoid confusion we use the name `` coverage '' .    to define a pareto archive , we need to introduce the mutually non - dominated relation :    mutually non - dominated relation : we say that two solutions are _ mutually non - dominated _ or _ non - dominated w.r.t . each other _ if none of the two solutions covers the other one . we can now define the pareto archive :    pareto archive : set of solutions such that any pair of solutions of the set are mutually non - dominated .    in the context of momhs , the pareto archive contains thus the mutually non - dominated solutions generated so far ( i.e. at a given iteration of a momh ) . we will denote @xmath24 the pareto archive generated by a momh . in other words @xmath24 contains solutions that are potentially efficient at a given iteration . the new method nd - tree is based on _ local _ ideal and nadir points that we define below . local ideal point of a subset @xmath25 denoted as @xmath26 is the point in the objective space composed of the best coordinates of all solutions belonging to @xmath27 , i.e. @xmath28 . naturally , local ideal point covers all solutions in @xmath27 .    local nadir point of a subset @xmath25 denoted as @xmath29 is the point in the objective space composed of the worst coordinates of all solutions belonging to @xmath27 , i.e. @xmath30 . naturally , local nadir point is covered by all solutions in @xmath27 . approximations of these two points will be also used in nd - tree :    approximate local ideal point of a subset @xmath25 denoted as @xmath31 is a point in the objective space such that @xmath32 . naturally , approximate local ideal point also covers all solutions in @xmath27 . approximate local nadir point of a subset @xmath25 denoted as @xmath33 is a point in the objective space such that @xmath34 . naturally , approximate nadir ideal point is also covered by all solutions in @xmath27 . we present here a number of methods for the online update of pareto archive presented in the literature and used in the comparative experiment . this review is not supposed to be exhaustive . other methods can be found in  @xcite and reviews in  @xcite . we describe two popular methods : linear list and quad - tree , the composite points method  @xcite and one recent method , m - front  @xcite , when used as a part of the non - dominated sorting procedure gave excellent performances compared to jensen - fortin s algorithm  @xcite , one of the fastest non - dominated sorting algorithm .        in this structure , a new solution is compared to all solutions of the list until a covering solution is found or all solutions are checked . the solution is only added if it is non - dominated w.r.t . all the solutions in the list , that is we need to browse whole list before adding a solution . the complexity in terms of number of pareto dominance check is thus in @xmath35 with @xmath36 the size of the list . when only two objectives are considered , we can use the following special property : if we sort the list according to one objective ( let s say the first ) , the non - dominated list is also sorted according to the second objective . therefore , roughly speaking , updating of the list can be efficiently done in the following way . let us consider two minimized objectives . the list is kept sorted according to the first objective . we first determine the potential position of the new candidate solution in the sorted list according to its value for the first objective , with a binary search . please note that in this case , it is more efficient to code the list with a dynamic array to get a constant access to the solutions of the list . once the place of the new solution is found and it is not equal to some existing solution , it is enough to use these two procedures :    * we compare the new solution to the preceding one ( if there is one ) in the sorted list : this solution has a better value according to the first objective compared to the new solution . therefore , if the preceding solution has also a better value according to the second objective , the new solution is dominated and can not be added . * if the solution is not dominated by the preceding solution , the solution can be added since the next solutions have a higher value for the first objective and can not thus dominated the solution . we then need to check if there are some dominated solutions : we browse the next solutions of the list , until finding a solution that has a better evaluation according to the second objective compared to the new solution . all the solutions found that have a worse evaluation according to the second objective have to be removed since there are dominated by the new solution . the worst - case complexity is still in @xmath35 since it can happen that a new solution has to be compared to all the other solutions ( in the special case where we add a new solution in the first position and all the solutions of the sorted list are dominated by this new solution ) . but on average , experiments show that the behavior of this structure for handling two objectives updating problems is much better than a simple list since the worst - case scenario rarely occurs . the use of quad - tree for storing potentially efficient solutions was proposed by habenicht  @xcite and further developed by sun and steuer  @xcite , and mostaghim and teich  @xcite . in quad - tree , solutions are located in both internal nodes and leafs . each node may have @xmath37 sons corresponding to each possible combination of results of comparisons on each objective where a solution can either be better or not worse . in the case of mutually non - dominated solutions in fact @xmath38 sons are possible since the combinations corresponding to dominating or covered solutions are not possible . quad - tree allows for a fast checking if a new solution is dominated or covered . a weak point of this data structure is that when an existing solution is removed its whole sub - tree has to be re - inserted to the structure . thus , removal of a dominated solution is in general costly . in this paper we use quad - tree2 version as described by mostaghim and teich  @xcite . m - front has been proposed relatively recently by drozdk et al . @xcite as a part of method for non - dominated sorting problem . the idea of of m - front is as follows . assume that in addition to the new solution @xmath0 a reference solution @xmath39 relatively close to @xmath0 and belonging to the pareto archive @xmath24 is known . the authors define two sets :    @xmath40    @xmath41    and prove that if a solution @xmath42 is dominated by @xmath0 then it belongs to @xmath43 and if @xmath44 dominates @xmath0 then it belongs to @xmath45 . thus , it is sufficient to compare the new solutions to sets @xmath43 and @xmath45 only . to find all solutions with objective values in a certain interval m - front uses additional indexes one for each objective . each index sorts the pareto archive according to one objective .    to find a reference solution being close to @xmath0 m - front uses the k - d tree data structure . the k - d tree is a binary tree in which each intermediate node divides the space into two parts based on a value of one objective . while going down the tree the algorithm cycle over particular objectives , selecting one objective for each level . drozdk et al . @xcite suggest to store references to solutions in leaf nodes only , while intermediate nodes keep only split values . to insert a solution k - d tree selects either lower or upper - or - equal half - space . please note , that the paper of drozdk et al . @xcite does not give the precise algorithm of k - d tree . in our implementation when a leaf is reached a new division is made using the average value of the current level objective . the split value is average between the value for new solution and the solution in the leaf . also alike drozdk et al . @xcite approximate nearest neighbor is found exactly as in the standard exact nearest neighbor search , but only four evaluations of the distance are allowed . we do not use the rebalancing of the k - d tree described in  @xcite since it would never be activated in the case of all data sets used in this paper and is very unlikely to be used in the case of solutions generated by pls .      while implementing m - front for the purpose of computational experiment we noticed that a number of elements can be improved to reduce further the running time . the improvement is in our opinion significant enough to call the new method m - front - ii ( see algorithm  [ mfrontupdate ] ) . the modifications we introduced are as follows :    * in original m - front the sets @xmath43 and @xmath45 are built explicitly and only then the solutions contained in these sets are compared to @xmath0 . in m - front - ii we build them only implicitly , i.e. we immediately compare the solutions that would be added to the sets to @xmath0 . * we analyze the objectives such that we start from objectives for which @xmath46 . in other words we start with solutions from set @xmath45 . since many new solutions are dominated it allows to stop the search immediately when a solution dominating @xmath0 was found . note that a similar mechanism is in fact used in original m - front but only after the sets @xmath43 and @xmath45 are explicitly built . * the last modification is more technical . m - front uses linked lists ( ` std::list ` in c++ ) to store the indexes and a hash - table ( ` std::unordered_map ` in c++ ) to link solutions with their positions in these lists . we observed , however , that even basic operations like iterating over the list are much slower with linked lists than with dynamic arrays ( like ` std::vector ` in c++ ) . thus we use dynamic arrays for the indexes . in this case , however , there is no constant iterator that could be used to link the solutions with their positions in these indexes . so , we use the fast binary search to locate a position of a solution in the sorted index whenever it is necessary . the overhead of the binary search is anyway smaller than the savings due to use of faster hash - table . parameter @xmath47 : a pareto archive @xmath24 with @xmath48 indexes @xmath49 , i.e. lists containing references to all solutions sorted incrementally according to objective @xmath50 and @xmath51 a k - d tree containing references to all solutions parameter @xmath52 : new candidate solution @xmath0    add @xmath0 to @xmath24 and its reference to all indexes @xmath53 and to @xmath51 find reference solution @xmath39 using approximate nearest neighbor search in @xmath51 stop : @xmath0 is rejected -@xmath54-@xmath55 comparison to set @xmath45 find @xmath56 , i.e. position of @xmath39 in index @xmath53 using binary search decrement @xmath56 to the last position where @xmath57_k = ref_k$ ] stop : @xmath0 is dominated it++ insert @xmath0 to @xmath53 right before @xmath56    -@xmath54-@xmath55 comparison to set @xmath43 find @xmath56 , i.e. position of @xmath39 in index @xmath53 using binary search increment @xmath56 to the last position where @xmath57_k = ref_k$ ] remove @xmath57 $ ] from @xmath24 , indexes @xmath53 and @xmath51 it insert @xmath0 to @xmath53 right after @xmath56    @xmath58 add @xmath0 to @xmath51      the composite points method has been developed by fieldsend  _ et al _ @xcite in 2003 . they introduce two new data structures called dominated tree and non - dominated tree . the dominated tree allows to detect if a new solution should be included into the pareto archive ( that is the new solution is not covered ) and the non - dominated tree allows to determine which solutions of the pareto archive are dominated by the new solution . some particular points , called _ composite points _ , are used to build the trees . these points are based on the maximal and minimal coordinates of the points of the pareto archive , and are built such that a total order between the composite points is obtained , according to the coverage relation . a binary search among the composite points of each tree allows then to update the pareto archive when a new candidate point is added . this method has been compared to the simple list for standard multiobjective continuous problems ( zdt test functions  @xcite ) but no significant difference in speed has been achieved . the new method proposed in this paper is based on the following idea that allows to reduce the number of comparisons to be made . consider a subset @xmath59 composed of mutually non - dominated solutions and a new candidate solution @xmath0 . assume that some approximate local ideal @xmath31 and approximate local nadir points @xmath33 are known . we can define the following simple properties :    1 . if @xmath0 is covered by @xmath33 , then @xmath0 is covered by each solution in @xmath27 and thus can be rejected . + this property is a straightforward consequence of the transitivity of the coverage relation . if @xmath0 covers @xmath31 , then each solution in @xmath27 is covered by @xmath0 . + this property is also a straightforward consequence of the transitivity of the coverage relation . 3 .   if @xmath0 is non - dominated w.r.t . both @xmath33 and @xmath31 , then @xmath0 is non - dominated w.r.t . each solution in @xmath27 . + if @xmath0 is non - dominated w.r.t . @xmath33 then there is at least one objective on which @xmath0 is worse than @xmath33 and thus worse than each solution in @xmath27 . if @xmath0 is non - dominated w.r.t . @xmath31 then there is at least one objective on which @xmath0 is better than @xmath31 and thus better than each solution in @xmath27 . so , there is at least one objective on which @xmath0 is better and at least one objective on which @xmath0 is worse than each solution in @xmath27 . 4 .   if none of the above holds , i.e. @xmath0 is neither covered by @xmath33 , does not cover @xmath31 , nor is non - dominated w.r.t . both @xmath33 and @xmath31 , then all situations are possible , i.e. @xmath0 may either be non - dominated w.r.t . all solutions in @xmath27 , covered by some solutions in @xmath27 or dominate some solutions in @xmath27 . + this property can be proven by showing examples of each of the situation . consider for example a set @xmath60 with @xmath61 and @xmath62 . a new solution @xmath63 dominates a solution in @xmath27 , a new solution ( 1 , 1 , 2 ) is dominated ( thus covered ) by a solution in @xmath27 , and solutions @xmath64 and @xmath65 are non - dominated w.r.t . all solutions in @xmath27 . the properties are graphically illustrated for the biobjective case in figure  [ fig : properties ] . as can be seen in this figure , in the biobjective case , if @xmath0 is covered by @xmath31 and @xmath0 is non - dominated w.r.t . @xmath33 then @xmath0 is dominated by at least one solution in @xmath27 . note , however , that this does not hold in the case of three and more objectives as shown in the example used in the proof - the point @xmath66 is covered by @xmath67 , non - dominated w.r.t . @xmath68 and @xmath66 is not dominated by any points in @xmath27 . in fact it is possible to distinguish more specific situations if property 4 holds , e.g. situation when a new solution may be covered but can not dominate any solution , but since we do not distinguish them in the proposed algorithm we do not define them formally .     based on comparisons to @xmath31 and @xmath33 only . ] the above properties allow in some cases to quickly compare a new candidate solution @xmath0 to all solutions in a set @xmath27 without the need for further comparisons to individual solutions belonging to @xmath27 . such further comparisons are necessary only if the situations described by property 4 hold . intuitively , the closer the approximate local ideal and nadir points the stronger are these properties , i.e. it is more likely that the further comparisons can be avoided . to obtain close approximate local ideal and nadir points we should :    * split the whole set of potentially efficient solutions into subsets of solutions located close in the objective space . * have good approximations of the exact local ideal and nadir points . on the other hand calculation of the exact points may be computationally demanding and a reasonable approximation may assure the best overall efficiency . based on these properties , we can now define the nd - tree structure . nd - tree data structure is a tree with the following properties :    1 .   with each node @xmath69 is associated a set of solutions @xmath70 . 2 .   each leaf node contains a list @xmath71 of solutions and @xmath72 . 3 .   for each internal node @xmath69 , @xmath70 is the union of sets associated with all sons of @xmath69 . each node @xmath69 stores an approximate ideal point @xmath73 and approximate nadir point @xmath74 . 5 .   if @xmath75 is a son of @xmath69 , then @xmath76 and @xmath77 . the algorithm to update a pareto archive with nd - tree is given in algorithm  [ algoupdate ] ( ` update ` ) . the idea of the algorithm is as follows . we start by checking if the new solution @xmath0 is dominated , covered or non - dominated w.r.t . all solutions in @xmath24 by going through the nodes of nd - tree and skipping sons ( and thus their sub - trees ) for which property 4 does not hold . this procedure is presented in algorithm  [ algoupdatenode ] ( ` updatenode ` ) . the new solution is first compared to the approximate ideal point @xmath78 and nadir point ( @xmath74 of the current node . if the new solution is dominated by @xmath79 it is immediately rejected ( * property 1 * ) . if @xmath80 is covered , the node is deleted and its whole sub - tree as well ( * property 2 * ) . otherwise if @xmath81 or @xmath82 ( * property 4 * ) and if @xmath69 is a leaf node , @xmath0 may be dominated by or dominate some solutions of @xmath69 and it is necessary to browse the whole list @xmath71 of the node @xmath69 . if a solution dominating @xmath0 is found , @xmath0 is rejected , and if a solution dominated by @xmath0 is found , the solution is deleted from @xmath71 . if after checking nd - tree the new solution was found to be non - dominated it is inserted by adding it to a close leaf ( algorithm  [ algoinsert ] ` insert ` ) . to find a proper leaf we start from the root and always select a son with closest distance to @xmath0 . as a distance measure we use the euclidean distance to the middle point between approximate ideal and approximate nadir points . the middle point is a point with individual coordinates equal to the average of corresponding coordinates in approximate ideal and approximate nadir points . such distance measure is very fast to compute based on the two points only . since we do not need the value of the euclidean distance but just the order induced by this distance , we used squared euclidean distance to avoid calculation of the root . once we have reached a leaf node , we add the solution @xmath0 to the list @xmath71 of the node and possibly update the ideal and nadir points of the node @xmath69 ( algorithm  [ algoupdateidealnadir ] ` updateidealnadir ` ) . however , if the size of @xmath71 became larger than the maximum allowed size of a leaf set , we need to split the node into a predefined number of sons ( algorithm  [ algosplit ] ` split ` ) . we first select one solution of @xmath71 for each new subnode . the first solution selected is the one with the highest euclidean distance to all other solutions in @xmath71 . then we select the solution with the highest average euclidean distance to all the solutions contained in all subnodes already created . we continue this procedure until the required number of subnodes is created . finally , we add each remaining solution of @xmath71 to the closest subnode . we use as distance measure the euclidean distance to the middle point between approximate ideal and approximate nadir points , as done in the ` insert ` algorithm . the approximate local ideal and nadir points are updated only when a solution is added . we do not update them when solution(s ) are removed since it is a more complex operation . this is why we deal with approximate ( not exact ) local ideal and nadir points . please note that our goal is to propose a practically fast method . we do not present any analysis of worst case complexity but we will show that this method performs well in computational experiments . parameter @xmath47 : a pareto archive @xmath24 organized as nd - tree parameter @xmath52 : new candidate solution @xmath0    create a leaf node @xmath69 with an empty list set @xmath71 and use it as a root ` \mathcal{l}(n ) \leftarrow \mathcal{l}(n ) + x ` @xmath83 root node ` updatenode`(@xmath84,@xmath85 ) ` insert`(@xmath84,@xmath85 )    parameter @xmath47 : a node @xmath69 parameter @xmath52 : new candidate solution @xmath0 compare @xmath0 to @xmath73 and @xmath74 -@xmath54-@xmath55 property 1 stop : @xmath0 is rejected -@xmath54-@xmath55 property 2 remove @xmath69 and its whole sub - tree -@xmath54-@xmath55 property 4 stop : @xmath0 is rejected remove @xmath44 ` updatenode ` ( @xmath86,@xmath85 ) remove @xmath75    parameter @xmath47 : a node @xmath69 parameter @xmath52 : new candidate solution @xmath0    ` \mathcal{l}(n ) \leftarrow \mathcal{l}(n ) + x ` ` updateidealnadir ` ( @xmath87,@xmath85 ) ` split ` ( @xmath87 ) find subnode @xmath75 of @xmath69 being closest to @xmath0 ` insert`(@xmath86,@xmath85 )    parameter @xmath47 : a node @xmath69    find the solution @xmath88 with the highest average euclidean distance to all other solutions in @xmath71 create a new subnode @xmath75 with an empty list set @xmath89 @xmath90 ` updateidealnadir ` ( @xmath91,@xmath92 ) @xmath93 create a new subnode @xmath75 with an empty list set @xmath89 @xmath90 ` updateidealnadir ` ( @xmath91,@xmath92 ) @xmath93 @xmath94 first solution in @xmath71 find subnode @xmath75 of @xmath69 being closest to @xmath44 @xmath90 ` updateidealnadir ` ( @xmath91,@xmath92 ) @xmath93    parameter @xmath47 : a node @xmath69 parameter @xmath52 : new candidate solution @xmath0    check in any component of @xmath0 is lower than corresponding component in @xmath73 or greater than corresponding component in @xmath74 and update the points if necessary @xmath95 parent of @xmath69 ` updateidealnadir ` ( @xmath96,@xmath85 ) we will show results obtained with nd - tree and other methods in three different cases :    1 . results for sets artificially generated , which allow us to easily control the number of solutions of the sets and the quality of the solutions . results when nd - tree is integrated into pls for solving the multiobjective traveling salesperson problem . results for sets coming from solutions generated by pls for solving the multiobjective traveling salesperson problem .    to avoid the influence of implementation details all methods were implemented from the scratch in c++ and c in as much homogeneous way as possible , i.e. when possible the same code was used to perform the same operations like pareto dominance checks . the code , as well as test instances and data sets , are available from the authors upon request . most of the results have been obtained on an intel core i7 - 5500u cpu at 2.4 ghz . the results involving pls ( section  [ secpls ] ) were obtained on an i5 - 450 m cpu at 2.4 ghz . the artificial sets are composed of @xmath69 points with @xmath48 objectives to minimize . the sets are created as follows . we generate randomly @xmath69 points @xmath97 in @xmath98 with the following constraint ( otherwise there can be few non - dominated points , all located near to the point @xmath99 ) : @xmath100 . with these constraints , all the non - dominated points will be close to the hypersphere of center @xmath101 and with a radius of length equal to @xmath102 . in order to control the quality of the points generated , we also add a quality constraint : @xmath103 . with a small @xmath104 , only high - quality solutions will be generated and this simulates well the behavior of a momh . we have generated data sets composed of 100 000 and 200 000 points , with @xmath105 , and for @xmath106 to 6 . in the main experiment we use data sets with 100 000 points because for the larger sets running times of some methods became very long . for each value of @xmath48 , five different quality sets are considered :    * quality = 1 , @xmath107 * quality = 2 , @xmath108 * quality = 3 , @xmath109 * quality = 4 , @xmath110 * quality = 5 , @xmath111    the fraction of non - dominated solutions grows both with increasing quality and number of objectives and in extreme cases all solutions may be non - dominated ( see table  [ tab : numbersnd ] ) .    [ tab : numbersnd ]    .numbers of non - dominated solutions in artificial sets [ cols="^,^,^",options="header " , ] according to the results of the computational experiments nd - tree should be a method of choice for storing and updating a pareto archive in the case of three and more objectives problems . in biobjective case the best choice remains the sorted list . we believe that with the proposed method for updating a pareto archive , new state - of - the art results could be obtained for many moco problems ( with @xmath112 ) and some momhs could be adapted to this special data structure in order to improve their results . in fact in this paper we were able to apply pls to motsp instances with up to 6 objectives while according to our knowledge so far pls was applied to biobjective and three - objective instances only . an interesting direction for further research is to adapt nd - tree to be able to deal with archives of a relatively large but limited size . finally , it would be interesting to study how nd - tree could be modified to update archives with stronger dominance relations than the pareto dominance ( e.g. sorted - pareto dominance  @xcite or lorenz dominance  @xcite ) . in this way , preferences of a decision maker could be integrated and methods such as pls will not be limited by the high number of pareto efficient solutions when solving many - objective combinatorial optimization problems . the research of andrzej jaszkiewicz was funded by the the polish national science center , grant no .  umo-2013/11/b / st6/01075 . a.lara , sanchez , g. , coello , c.  c. , schtze , o. , feb 2010 . hcs : a new local search strategy for memetic multiobjective evolutionary algorithms . ieee transactions on evolutionary computation 14  ( 1 ) , 112132 .      angel , e. , bampis , e. , gourvs , l. , 2004 . a dynasearch neighborhood for the bicriteria traveling salesman problem . in : gandibleux , x. , sevaux , m. , srensen , k. , tkindt , v. ( eds . ) , metaheuristics for multiobjective optimisation . lecture notes in economics and mathematical systems vol . 535 , berlin , pp . 153176 . brockhoff , w. , 2010 . theory of randomized search heuristics : foundations and recent developments . world scientific publishing company , ch . theoretical aspects of evolutionary multiobjective optimization , pp . 101139 . coelho , v.  n. , souza , m. j.  f. , coelho , i.  m. , guimares , f.  g. , lust , t. , cruz , r.  c. , 2012 . multi - objective approaches for the open - pit mining operational planning problem . electronic notes in discrete mathematics 39 , 233240 . fortin , f .- a . , grenier , s. , parizeau , m. , 2013 . generalizing the improved run - time complexity algorithm for non - dominated sorting . in : proceedings of the 15th annual conference on genetic and evolutionary computation . gecco 13 . acm , new york , ny , usa , pp . 615622 . habenicht , w. , 1983 . essays and surveys on multiple criteria decision making : proceedings of the fifth international conference on multiple criteria decision making , mons , belgium , august 913 , 1982 . springer berlin heidelberg , ch . quad trees , a datastructure for discrete vector optimization problems , pp . 136145 . inja , m. , kooijman , c. , de  waard , m. , roijers , d. , whiteson , s. , september 2014 . queued pareto local search for multi - objective optimization . in : ppsn 2014 : proceedings of the thirteenth international conference on parallel problem solving from nature . . 589599 . kooijman , c. , de  waard , m. , inja , m. , roijers , d. , whiteson , s. , april 2015 . pareto local policy search for momdp planning . in : esann 2015 : proceedings of the 23rd european symposium on artificial neural networks , special session on emerging techniques and applications in multi - objective reinforcement learning . . 5358 . mostaghim , s ;  teich , j. , 2004 . quad - trees : a data structure for storing pareto sets in multiobjective evolutionary algorithms with elitism . in : abraham , a ;  jain , l. g.  r. ( ed . ) , evolutionary multiobjective optimization , book series : advanced information and knowledge . springer - verlag london ltd , pp . 81104 . mostaghim , s. , teich , j. , tyagi , a. , 2002 . comparison of data structures for storing pareto - sets in moeas . in : evolutionary computation , 2002 . proceedings of the 2002 congress on . vol .  1 . 843848 . omahony , c. , wilson , n. , 2013 . sorted - pareto dominance and qualitative notions of optimality . in : gaag , l.  c. ( ed . ) , symbolic and quantitative approaches to reasoning with uncertainty : 12th european conference , ecsqaru 2013 , utrecht , the netherlands , july 8 - 10 , 2013 . springer berlin heidelberg , pp . 449460 . paquete , l. , chiarandini , m. , sttzle , t. , 2004 . pareto local optimum sets in the biobjective traveling salesman problem : an experimental study . in : gandibleux , x. , sevaux , m. , srensen , k. , tkindt , v. ( eds . ) , metaheuristics for multiobjective optimisation . lecture notes in economics and mathematical systems vol . 535 , berlin , pp . 177199 . schtze , o. , 2003 . a new data structure for the nondominance problem in multi - objective optimization . in : proceedings of the 2nd international conference on evolutionary multi - criterion optimization . springer - verlag , berlin , heidelberg , pp . 509518 . schtze , o. , mostaghim , s. , dellnitz , m. , teich , j. , 2003 . covering pareto sets by multi - level evolutionary subdivision techniques . in : proceedings of the second int . conf . on evolutionary multi - criterion optimization . springer , berlin , pp . 118132 . teixeira , c. , covas , j. , sttzle , t. , gaspar - cunha , a. , 2009 . application of pareto local search and multi - objective ant colony algorithms to the optimization of co - rotating twin screw extruders . in : viana , a. , et  al . ( eds . ) , proceedings of the eu / meeting 2009 : debating the future : new areas of application and innovative approaches . . 115120 .
Abstract:

Answer: in this paper we propose a new method called nd - tree for fast online update of a pareto archive composed of mutually non - dominated solutions . nd - tree uses a tree structure in which each node represents a subset of solutions contained in a hypercube defined by its local approximate ideal and nadir points . a leaf is a subset of solutions organized as a simple list , and an internal node is subset of solutions composed of the union of all its sub - nodes . using heuristic rules we build subsets , either leafs or internal nodes , containing solutions located close in the objective space . using basic properties of local ideal and nadir points we can efficiently avoid searching many branches in the tree . nd - tree may be used in any multiobjective metaheuristics e.g. in an multiobjective evolutionary algorithm to update the external archive of potentially efficient solutions . we experimentally compare nd - tree to simple list , quad - tree , and m - front methods using artificial and realistic benchmarks . finally we apply nd - tree within two - phase pareto local search for traveling salesperson problems instances with up to 6 objectives . we show that with this new method substantial reduction of the computational time can be obtained . multiobjective optimization , pareto archive , quad - tree , pareto local search , many - objective optimization