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: for many years , the gold standard in minimization packages used by high energy physicists has been the fortran program minuit  @xcite . this code contains algorithms which are well - suited for many hep uses , especially when there is a large data set and minimization is needed to find a multi - parameter fit . the principal algorithm , called migrad , is a variable metric method . this exploits the idea that althought a costly second derivative computation is needed to move to the function minimum , that matrix need not be computed accurately for early steps ( since it will be changing significantly anyway ) . instead , it can be iteratively improved , so that by the time the method gets near the actual minimum , the matrix is accurate and rapid convergence occurs . another prime capability is a set of routines for analyzing the solution : the hesse method of supplying the parameter error correlation matrix , the minos method of nonparabolic chi - suqare error intervals , and contour , which depicts the shape of a fixed - sigma error curve around the solution point for functions which may not be well - described as bi - linear forms . minuit , however , is getting a bit stale . it is written in fortran , which is no longer the _ lingua franca _ of computing in hep ; it can be awkward to link c++ code with code from a different language . and modern users would like certain advantages relating to treating a minimization problem and its function , algorithms and domains , as objects . one such advantage is the option of performing multiple independant minimizations in tandem . the minuit code is difficult to maintain , as evidenced by the fact that there have been no enhancements since a major overhaul done by f. james about a decade ago , and no major algorithm improvements for some time before that . the hep community knows of promising alternative methods , such as fumili  @xcite and the related leamax  @xcite , and certainly the field has advanced in these years , with simulated annealing and genetic algorithms maturing . the fact that nobody has seen fit to incorporate these as options in minuit can be taken as evidence that such changes would be hard to make .    even assuming that minuit needs no maintenance because all bugs have been eliminated and that its capabilities are forever adequate , what happens when all the cernlib fortran modules start becoming awkward to bring forward ? there is reason to fear that minuit might fall into functional obsolecence . we asked a group of hep physicists who work heavily with c++ computation whether it would be a good idea to produce a c++ object - oriented minimization package with minuit s capabilities . about 40% reacted positively : `` it s about time we had this . '' some of these physicists had faced a tougher time accomplishing an actual project , for lack of a good stand - alone c++ minimizer . others reacted with some variant of `` are you out of your mind ? sombody must already have done this ! '' we will argue below that nobody actually has delivered an adequate object - oriented minimizer with the capabilities required . most commercial minimization packages fall short in capability aspects . in some cases , their algorthms are not as appropriate as migrad for the sort of problems which occur in hep ; more often they lack solution analysis methods equivalent to contour and minos . also , non - trivial licensing issues can hinder the adoption of commercial codes as _ de facto _ standards in the hep community . a more promising approach is to translate minuit directly into c or c++ . several analysis packages , including root  @xcite , provide fitting in this manner . another approach is to place a wrapper around the minuit code , to allow the interface to be expressed in c++ , and package the compiled fortran code along with the wrapper as an integrated library . gemini  @xcite does this , providing a shell over a choice of minuit or the commercial nag c minimizer . either approach solves the problem of inter - language awkwardness , but they do not address the issue of maintainability and ease of enhancement . and fundamentally , a minimization problem still can not be treated as a first - class object .    the root team , in particular , has done significant work both in cleaning up the _ f2c _ translation of minuit and in removing deficiencies involving multiple minimizations . but not everybody welcomes the requirement of linking to the large body of root code , when all that is needed is a minimizer . one other c++ object - oriented minimizer is definitely worth mentioning . matthias winkler , mentored by fred james at cern , is working on an `` object - oriented re - implementation of minuit in c++ '' under the auspices of the seal programme . clearly two nearly - identical good c++ minimizers are not needed . we have been in contact with dr . james , and intend to coordinate the two parallel efforts , hoping to emerge with a single package of superior design and implementation quality . it is worth saying a few things about the concept of a minimizer as opposed to a fitter . since most applications of minimization in hep are motivated by multi - parameter fits to data , the two can be confused . most of the packages which wrap minuit tend to be fitters , that is , they present a convenient way to provide a function form and a data set , and go about using minimization to optimize the parameters in a fit to that data . when the function @xmath0 is a chi - squared sum , the resulting minimization procedure is specialized , in the sense that the nature of the fitting problem assigns a meaning and a natural scale for differences in values of @xmath1 . algorithms ( such as migrad ) which yield a meaningful estimated distance above the true minimum can exploit a natural termination condition , since the value of @xmath1 assumed to be accurate only to some small fraction of a unit in chi - squared . and analysis of the estimated shape of @xmath1 can meaningfully be translated to statistical statements about correlations among the parameters . however , the kernel operation of minimizing an arbitrary function is a valid capability to ask for . by design , the c++ minimization package under development is a minimizer , not a fitter . one can automate the steps to turn a fitting task into a minimization problem . several others have done this in c++ , and though it may become convenient to package such a fitter adaptor with the minimizer , this adaptor need not be considered part of the development of the minimizer . the requirements guiding the developement of this package are driven by the goal of acceptance and widespread use in the hep community . this will require providing a significant set of capabilities , discussed below . the package must not to be deficient , with respect to minuit , in any capability . also , attention must be paid to certain performance considerations . the package must be clean to use and self - contained , and should not depend on other packages . and stringent testing will be required , to ensure that the package is robust .    at the same time , we wish to provide interfaces which will encourage code which is clear and readable and which does not deviate from the classical c++ coding idioms accepted by top practitioners . our goal is to provide a package that can be used for many years . it is impractical to provide something which is `` so good it will never need to change . '' instead , we strive to organize and code the package such that enhancements and maintenance changes will be as easy to incorporate as possible .    in order to stand the test of time , we feel the minimization package must include complete and easily understood user documentation , as well as full mathematical documentation of all algorithms used . the key requirement is that every capability present in minuit must be matched in our minimization package . it could be that minuit is the ideal minimizer for hep _ per se_. or perhaps the way physicists have grown accustomed to doing data analysis has evolved in tandem with minuit to take advantage of its capabilities . the point is moot  any package which can not match the capabilities present in minuit will be unacceptable to the hep community . these capabilities lie in several areas . clearly , the minuit algorithms ( migrad , simplex , minimize ) must be present . and user code must be able to control steps taken and selection of ( and switching among ) algorithms . code must be the able to establish limits on some or all of the function parameters , and fix and release their values . and the methods for error analysis , nonparabolic error intervals , and two - parameter contour lines ( analogous to hesse , minos , and contour ) must be provided . it is tempting , along the way , to make any tweaks that are noticed to be pure improvements . some of these are pure extensions , such as the ability to define parameters which have only an upper or lower limit . others might be minor algorithm improvements , such as a theoretically superior finite difference step size when numerically computing derivatives . in such cases , it is required that default behavior precisely mimic the procedures used by minuit . the minuit - matching requirement is a matter of capability , and not of interface . minuit assumes one of two usage paradigms :    1 . the minimizer is a static facility , which the user code explicitly initializes to embark on a problem . the user program then controls this facility by calling various subroutines . a commanding minimizer program is supplied with `` cards , '' each directing an aspect of initialization or an element of control .    clearly , both of these paradigms are obsolete in the context of a c++ , object - oriented package . the experienced minuit user may miss the mantra of calling mnexcm and the nine - argument calls to mnparm ; these will be replaced by methods which are more readable and which assigns one easily understood purpose to each method called . still , when the user constructs a minimization problem and invokes problem methods , she is logically employing capabilities which encompass the same set as was available in minuit .      to make best use of the object oriented abilities provided by c++ , it is not enough merely to ensure that a minimization problem can be used as an independant entity . we must use the tools at our disposal to make the setup and control of a minimization problem as clear and flexible as possible . the focus is on user code readability , and on avoiding traps where the user can do something which looks sensible but may have disasterous ( or worse , subtly incorrect ) consequences . this means isolating conceptually distinct areas of control , and providing ways of excercising control of one aspect of execution at a time . thus the user will be able to interact with termination conditions associated with a problem , and separately with the domain restricting ranges of parameters , and with the selection of algorithms . in each of these areas , the user controls behavior by interacting with specialized versions of objects . the minimization problem then interacts with these objects via by their base class interfaces . an example is in order :    to control how parameters are limited in range , the problem will own an instance of a domain . domain methods takes care of such necessities as translating points and gradients between internal unrestricted cartesian coordinates used by the algorithm , and external parameters known to the function . but domain by itself is merely a base interface class ; the actual domain attached to the problem would ( by default ) be a rectilineardomain , which provides the capabilities present in minuit : restricting the upper and/or lower limits on each parameter individually . rectilineardomain inherits from domain . the user who desires to alter these limits would obtain a pointer to the domain , and invoke methods of rectilineardomain such as setupperlimit(n , x ) . now suppose another user has a problem which should be restricted to a region delimited by arbitrary planes in parameter space . this user would provide a different sort of object inheriting from domain  say linearconstraintsregion . and then she might call methods of that class , such as forcelessthan(coefficients , maximum ) . the algorithm does not care which type of domain is in place .    in providing the user interface , it is tempting to parallel the control made available to minuit . this is good only insofar as there are nt opportunities to improve on the semantics provided by the minuit subroutines . we are finding that there is a tendency to rely too heavily on doing what minuit does , and we need to be careful to look for opportunities to provide interfaces which would be clearer to people who are not steeped in minuit usage . to paraphrase stroustrup , this package should be `` as close to minuit as possible , _ but no closer_. ''      the design of the package must be such that as new algorithms and anaylsis methods are developed , or as new types of termination conditions or domains are required , they can be added either to the package itself , or to one s local version of the package . we can not insist that the work associated with such an enhancement be small , because in principle coding some new algorithm might inherently involve a lot of effort . what we can , and do , insist on  and this drives the entire package design  is that the additional work imposed by being part of this minimizaztion package is a very small fraction of the typical work needed to create the new algorithm or domain itself . this `` additional work '' includes    * understanding all the steps one needs to do . * following the logic of existing package code if necessary . * discovering precisely where modifications will be needed , and what the nature of any new classes must be . * coding any c++ boilerplate needed to support the enhancement . these are the tasks which represent a barrier to putting an enhancement into any of the various c++ minuit translations , and we must ensure that this barrier is small .    moreover , we must see to it that the enhancer , who may well be an expert in algorithm coding , need not also be an expert in c++ .    the way to accomplish these goals is to carefully design the system such that any two concepts that need not be logically intertwined , are expressed as distinct classes which interact as little as possible . the overall package must be split into distinct subsystems , and dependencies carefully controlled . if this is done properly , then the enhancer need only be concerned about the small corner of the package this enhancement will involve . thus all the code to be examined in creating a new algorithm class would be in the algorithm base class , and all the boilerplate involved would be found in classes in the algorithm subsystem . the same philosophy that leads to ease of enhancement without undue tracing through an entanglement of code , will also greatly reduce the burden of maintenance . over the long term for a non - commercial package , one has to assume that any maintainer might well lose familiarity with the internals of the system . so it is crucial that proper coding idioms be employed , to avoid setting subtle traps . and good separation of functionality and control of dependencies are critical , so that when changes need to be introduced , they will be localized and not have tendrils of effect throughout the package . also , each class should have self - contained ( and easy - to - automate ) unit tests ; each subsystem should have integration tests only involving it and any subsystems it depends on ; and a good suite of regression integration tests should be provided to help catch any errors introduced . having set a goal of easy enhancement , we must recognize that there are different sorts of enhancements which will normally be performed by different levels of developers . a user might well want to create some new form of termination condition ; this must be very straightforward . more rarely , a more experienced coder might want to derive a new type of domain ( for example , a rectilineardomain but replacing the arcsin mapping function with some sort of sigmoid ) . still more rarely , adding a fundamentally new algorithm is the sort of ehancement which might be assisted by the principal package maintainers .      in the `` low - ceremony '' tight design iteration process so often sucessfully used in high energy physics code development , certain fundamental design considerations are often overlooked . two such considerations are the required levels of agility and efficiency under various circumstances . in the course of designing this package , there will be opportunities for tradeoffs among areas such as coding style , levels of agility when faced with minor changes in needs , and performance efficiency . an up - front grasp of the goals in these area will allow the developer to make intelligent choices . we have discussed above the agility considerations , in terms of extensibility and maintainability goals . as to efficiency , the typical first reaction is to say `` of course this should be as efficient as possible , '' but that is an inappropriate goal for this package . taken to the extreme , a total focus on efficiency might dictate unmaintainable coding tricks , avoidance of sound accepted c++ techniques , and ultimately abandoning the object - oriented paradigm . some aspects of performance efficiency are critical , but the design goals should permit other forms of efficiency to be sacrificed in favor of better modularity , semantic clarity , or other desirable features . the hep minimization problem has a natural set of performance goals , and these induce an efficiency philosophy which is implied in minuit , and which we state implicitly for our package : it is assumed that the work involved in an invocation of the function being minimized is very large compared to the overheads of bookkeeping , argument passing , and transformations of parameters . thus tradeoffs which involve attaining beneficial properties at the cost of incuring extra work in these bookkeeping areas are considered positive . furthermore , it is also assumed that the number @xmath2 of parameters for tractable minimization problems of the nature we are considering will not be too large . by this , we mean that even in an algorithm which requires operations ( such as inversion ) involving @xmath3 matrices , the time spent for such operations does not dominate the time needed to evaluate the fuction in forming those matrices .      this minimization package will perform well on problems in which the dominant time cost is that of evaluating the function . because of the tradeoffs between some forms of efficiency and other desirable features , this package may not perform efficiently for certain other classes of problems . it is tempting to say that such problems  which involve inexpensive function calls  are so low - cost that nobody should care about performance , so that this minimization package is suitable for all uses . we should realize that the above rationalization is a bit facile ; there are actually at least four cases in which efficiency considerations may become important :    1 . each call to evaluate the function is expensive . this is the case in most hep fitting applications , and is the case which our package will deal with well . some problems are potentially expensive not because each function call is costly , but because there are a large numbers of parameters , such that matrix operations or even variable transformations begin to become dominant . when the number of parameters is huge , the algorithm may require a huge number of function calls . if function calls are inexpensive relative to the number of parameters , then bookkeeping overhead can become significant . the minimization problems appearing in lattice qcd to find fermion propagators are of this nature . the overall scientific problem being solved may involve solving large numbers of easy minimization problems . thus , while no single minimization is costly , performance on each problem might remain a critical consideration . in that case , problem initialization costs can become important . the efficiency requirement laid on this minimization package , or for that matter on minuit , is that it is appropriate for efficient solution of problems in the first of the above categories . this has been found to meet the fitting and other minimization needs of most hep applications other than lattice qcd . the immediate preferences of potential hep users can come into tension with the desired use of classical c++ coding idioms accepted by top practitioners . concepts which are legacies from the fortran days lead physicsts to become comfortable with coding techniques which are today recognized as poor use of c++ . in such cases , it is better to use the accepted modern idiom , rather than to deliver inferior  but temporarily more familiar  interfaces . thus for example , it would be unwise to shun the use of auto_ptr when that is the appropriate way to convey the passing of responsibility for an object . the assumption is made that physicists who may at first be uncomfortable due to unfamiliarity with this idiom will quickly pick up on it . and the benefits in using concepts which have been battle - tested by the world of developers , and which match the understandings of skilled programmers , are considerable . few of the decisions discussed in this section were black and white , and some may change in the course of coordinating with the cern effort , but after some consideration , the following principles were adopted :      how does the user supply the function to be minimized ? the naive mechanism is that the user passes to the problem a global - scope function of some specified signature ( for example , taking a vector of coordinates and returning a value ) . in fact , this notion is so obvious to the physicist coming from fortran , that we must support it in some manner . but there are considerable advantages to being able to supply a functor ( a instance of a class that has an operator ( ) so that it looks like it can be called as a function ) instead . one advantage to suppying a functor is that unlike the global function , an instance of a functor may have internal state , and the class may provide methods to influence the function behavior . ( for example , one may wish to sum errors over fewer data points early in a minimization problem , and more later . ) this natural concept in c++ replaces the awkward nature of fcn in minuit , which contains six arguments . only two arguments are fundamental to the notion of a function to be minimized : the vector xval of coordinates , and the result fval . by allowing for a functor class , we obviate the need for iflag and futil ; and by separating the gradient method we eliminate the grad argument . the package has several situations in which behavior is defined by hooking into user code . two examples are the functor representing the function to be minimized , and a domain object implementing limits on cooordinates . the pattern of allowing a user to control behavior , by supplying a class which overrides one or more `` hook '' methods , resonates well with physicist coders . were we to say that the user constructs a problem supplying a reference or pointer to an instance of some class derived from function , the hep community would not complain . and if we further warned the user not to allow this function to go out of scope while the problem still needs it , most users would find that an obvious caution . however , the consequence of violating that stricture are dire : when the problem calls ( by pointer ) the out - of - scope method , confusion will reign , with no sensible diagnostic to point out why . good library developers protect against this by copying the vulnerable object . the problem , then , will _ own _ the function and domain , and not risk having those objects `` pulled out from under it . '' this notion is expressed by the user passing an auto_ptr to a new object ( on the heap ) ; the problem assumes responsibility for deleting that object when it no longer needs it . having made that decision , we must provide a way for the user to invoke controlling methods on the instance of the functor owned by the problem . the original object that the user supplied to the problem wo nt do , since the problem will be using a copy of it ; in fact , if the correct mantra was used to pass ownership to the problem , the user code should no longer be able to access the original object . instead , the correct recipe is :    * the user obtains a pointer to the base class ( in this example , function ) by calling a method of problem ( in this case , getfunction ( ) ) . * this pointer is dynamic cast into the class the user needs to control . the user can then invoke any control methods applicable for that class . * the pointer is now allowed to go out of scope , without the user calling delete on it .    the suggested syntax , which automatically obeys this recipe , involves using the returned pointer as a temprary variable which never is valid outside a single statement . the syntax can look like :    ....    dynamic_cast < rectilineardomain >      myproblem.getdomain()- >        setupperlimit(paramnumber , 35.1 ) ; ....    because this mantra may not be immediately familiar to our users , it is important that the package stick to the same rules ( _ never _ use delete on a pointer which problem has supplied ) for every case of object access . the powerful template facility in c++ would allow clever ways to improve the user interface and/or implementation of the minimization package . to give one example , methods giving a choice of supplying a user function or functor to act as the minimization target could be unified by expressing problem as a template , taking the function class as a template parameter . the use of templates has several adverse consequences which we perfer to avoid where practical : it makes it a bit tougher for a user to comprehend a header file , it can make maintenance more difficult , it makes compilation error messages much less readable , it introduces quirky linking issues , and it risks harming portability on certain weak compilers . there are valid answers to each of these objections ; still , absent compelling benefits , the package design avoids the use of templates . on the other hand , existing templated classes in the standard library are used freely . our philosophy concerning clever and subtle patterns to which convey minor benefits is that it is sometimes possible to be _ too _ clever . for example , the design of the function base class ( discussed in section [ function ] ) makes it possible to supply a gradient ( ) method but forget to override gradientavailable ( ) . that is a bit of a trap  the algorithm would ignore the user - supplied method when it needed a gradient . ( the same trap is present in minuit , if the user forgets to call setgradient ) . we could finesse this by using a variant on the visitor pattern : when the algorithm needs a gradient , it requests it from the function object but the signature of the gradient ( ) method also supplies the address of the algorithm s finite difference gradient - computing method . the base class gradient ( ) invokes this method ; if the user s functor class overrides gradient ( ) that routine replaces the finite difference method behavior . this pattern incurs some efficiency loss , but since that is merely a matter of an extra argument passed around , and an extra subroutine call , that cost would ( by the definition in section [ efficiency ] ) be acceptable . however , the pattern looks complicated and mysterious ; it might make maintenance more difficult and it might worry potential users who havent put in the time to understand it . since the gain is relatively slight , we have come down on the side of the simpler gradientavaliable ( ) solution .      in many cases , there may be multiple concepts for which the natural expression in bits and bytes will look like identical data structures . for example , both the set of coordinates identifying a point , and the set of values representing the components of the gradient of a function , are underneath indexed collections of doubles . the question is , should the design introduce distinct classes to represent the distinct concepts , or should it use a vector of doubles for both ? good c++ design takes the view that if two concepts are distinct , they should be represented by two distinct classes . the idea is that you ca nt be sure in advance that you will never alter the way you work with one of the concepts , and if they are represented by distinct classes , making such a change will be much less of a maintenance headache . introducing a distinct class to represent each distinct concept also provides some measure of type safety , in that you ca nt accidently use one concept when you mean the other . and naming each concept leads to superior code readability . in a package that emphasizes extensibility , these three advantages are significant . the disadvantage in introducing such classes is limited to the need to write some boilerplate code writing to provide their headers and near - trivial implementations . therefore , especially concerning minimization package internals , we tend to introduce distinct classes for distinct concepts , particularly in the area of geometric concepts ( see section [ geometric ] ) such as point and gradient . the minimization package must be a good c++ citizen . by that , we mean it must not step on potential user symbols ; we place everthing into a namespace to guard against this . other items we avoid are the use , in headers , of macro definitions ( which may clash with user defines ) and of using statements ( which defeat the purpose of the namespace protection if they appear in an included header ) . namespace use has been approved for the clhep library , which has severe portability requirements , so we are not concerned on that issue for the minimization package . the package is divided into several subsystems . a subsystem is defined as some group of classes and responsibilities , such that each subsystem depends on a minimal number of other subsystems , and such that these dependencies can be expressed as cleanly as possible . the first step in logical design was to identify the concepts present in minuit . it is no coincidence that each subsystem coresponds to one major concept in the minimization realm . the benefit of a good separation into subsystems comes during development of various classes , and later during enhancement implementation , because the coder can have in mind a smaller `` overall picture '' when deciding how to do things . good isolation of subsystems makes it possible to develop the more detailed package design one part at a time , without constantly having to consider the ramifications of each decision on every other part of the design . finally , the up - front thought about subsystem components has yielded immediate benefits in the form of solidifying the collection of concepts involved in the minimization package .    in principle a `` subsystem '' need not be tied to a group of classes . however , each identified subsystem in the minimization package is characterised by a key class embodying the interface of that subsystem to the user and/or to the other subsystems . this section will outline the responsibilities and dependancies of each subsystem , and discuss associated design issues . the problem subsystem has the responsibility for the user interface to a minimization process . it contains the problem class , which takes ownership of algorithm , domain , and termination objects , and contains the problemstate . minimizing involves the problem repeatedly requesting steps from the currently associated algorithm . after each step , it determines whether the user - supplied termination condition is satisfied , or whether the algorithm itself has declared that it ca nt profitably proceed further . if termination is reached , control returns to the user code , which may change the algorithm , the domain ( for example , releasing some parameters and fixing others ) , and/or the termination conditions , and again invoke the problem s minimize ( ) method . problemstate is a structure containing information about the state of solution which is not peculiar to any one algorithm or domain . for example , the concept of the current best guess is inherent in the nature of minimization ; it is applicable whehter the algorithm producing it was migrad or simplex . the user indirectly interacts with problemstate , via methods of problem , when she needs to inspect this information . subsystems of the minimization package interact with problemstate more directly . for example , termination conditions inspect its data to decide whether to terminate , and the algorithm subsystem is dependant upon problemstate , reading and altering its data . there is danger that problemstate can act as a catch - all for data that violates reasonable encapsulation goals , effectively replicating  with all its flaws  the notion of a big global common area . there is some art in restricting the contents to those items genuinely inherent in the minimization concept , and keeping those items which are associated with an algorithm in the specific algorithm classes . the algorithm base class provides interface between the problem and subclasses implementing specific minimization algorithms . this base class also standardizes ways for its derived classes to access the function and domain objects . for example , the mantra for obtaining a function gradient involves finding out whether the function object can provide it analytically , and if not , utilizing the algorithm s own finite difference methodology ; this mantra is encapsulated in algorithm . any given algorithm class has an iterate ( ) method , which `` takes a step , '' improving the state of solution . it also decides whether to declare that this algorithm is finished in the sense that no further meaningful improvement is forseen . between steps , the problem applies user - supplied termination conditions which may declare victory well before the algorithm declares exhaustion . this subsystem depends on the function and domain subsystems : algorithms , of course , repeatedly requests values of the function ( and if available its gradient and second derivatives ) , and utilize the domain to map its internal coordinates into external parameters to supply to those function calls . this subsystem has the responsibility for the various solution analysis tools such as minos , hesse , and contour . such tools have roughly the same needs and dependencies as algorithms , but a class derived from analysis can in addition also depend on algorithms which it needs to call . the function subsystem has the responsibility for providing function values and , if available , analytic gradient and second derivatives . the primary class is the abstract interface function . commonly , users will have an ordinary function ( taking a constant refernce to a point or to a vector of doubles , and returning a double ) to minimize . such users can make use of the free method userfunction ( ) , to supply this function when constructing the problem , as in :    ....    double f(const vector < double > v ) {      assert ( v.size()=nargs ;       // ... do work to return a value }    problem myprob ( userfunction ( nargs , f ) ) ; ....    an alternative signature would take both the function and a method to compute the gradient , if direct gradient computation is available . the user can also optionally supply global methods to compute second derivatives . userfunction ( ) constructs a function object on the heap , and returns an auto_ptr to it ; in the example code , this is immediately passed to the constructor of the problem . alternatively , a user can inherit from function to form his own `` functor '' class , and supply an instance of that when constructing a problem . in that case the functor class must override the method which returns the function value at a point . it may also override the gradient ( ) method , and if it does , it should override gradientavailable ( ) to return true  this replaces minuit s setgradient command , and more importantly , obviates the need for the user function to be passed an iflag argument to tell it what do do . similarly , the functor class may override the hessian ( ) and hessianavailable ( ) methods . the functor is to be created on the heap and an auto_ptr supplied to the problem so that there is no danger of the user getting mysterious errors due to the functor going out of scope while the problem is depending on it . the user can later recover a pointer to this function object and thus invoke its methods :    ....    class myfunct : public function {      myfunct ( ) : function(nargs ) { }      double operator()(const point & p ) ;      changestuff(int c ) ;      ... } ;    problem myprob ( auto_ptr < function >                    ( new myfunct )        ) ;    ( dynamic_cast < myfunct * >       ( myprob.getfunction()))->changestuff(3 ) ;       ....    the function subsystem depends only on the geometric concepts subsystem . the domain subsystem has the responsibility for providing mappings between the internal coordinates that an algorithm works with , and the external parameters that a function deals with . this solidifies the decoupling , already present in minuit , between the algorithm and any restrictions on the parameter values . that is , it would be too difficult to invent an algorithm that could deal with any arbitrary sort of variable restrictions ; instead , the algorithms themselves work with clean , unrestricted cartesian coordinates , and all the complexity involving restrictions is separated off into the domain classes . the domain base class provides the interface for algorithms and other parts of the package to transform variables . there is , of course , the complication that geometric objects of different kinds transform differently under the mappings provided by a domain . thus specific domain derived classes provide methods to map various geometric concepts : points , gradients , hessians , characteristicscales , ( error ) correlations , and whatever else algorithms and analyses may need . the base domain class supplies these methods in terms of virtual functions which supply the map s value and derivatives . domain encapsulates the minuit concepts of parameter limits and of fixing / releasing parameter values . half the interface of minuit is devoted to manipulating the way parameters can vary . most of this control belongs in the domain subsystem . the exemplar of a domain ( and the one associated with a problem by default ) is the rectilineardomain , which provides the mapping functions and capabilities present in minuit . a user could fairly easily create a class inheriting from rectilineardomain in order to modify the mapping functions used . rectilineardomain has the property that its individual components map in a separable way , that is , the map of each internal component depends on only one variable external parameter . algorithms can query a domain about this property . for example , when migrad utilizes diagonal second derivatives to avoid the need for computing the full hessian , it relies on this separablity property ; and in more complex domains , additional work may be needed . unlike functors derived from function , which need only supply the fundamental value - suppying method ( and may optionally supply gradient etc . ) , each domain subclass is required to supply all the geometric concepts mapping methods . while typical users are expected to create functors deriving from function , we can assume that those who embark on defining new types of domains are somehwhat more expert . the domain subsystem depends only on the geometric concepts subsystem . the termination subsystem has the responsibility for conditions conditions which the problem can apply , to determine whether to terminate the minimization iteration and return control to the user . the termination base class also defines an interface for potential user - created termination conditions . the subsystem contains termination classes corresponding to each of the criteria applicable in minuit , plus a time - based criterion . time - based termination is often preferable to a count of function calls , but has a potential disadvantage with respect to portability . however , the standard library does provide for timing at the one - second granulaty level , and almost all systems are compliant in this respect . termination also comes with methods to build `` compound '' termination objects as by logically combining ( and , or ) more basic termination objects . such a compound criterion could , without too much difficulty , be coded by each user as a class deriving from termination . but the desire to apply a group of criteria , terminating ( say ) if either the accuracy criterion is met , _ or _ the number of function calls is excessive , will be a very common case . so it was decided to make life easy for such users by supporting these combinations . we supply the logical combiners and and or , as this will satisfy the vast majority of applications . the termination subsystem depends on problemstate . the minimization process deals with a number of geometric concepts , which are logically distinct and not interchangeable . most of these behave like a vector of doubles . for example , a point in exterior coordinate space is not interchangeable with a gradient of the function at some point , even though they both have the same number of components . and neither is the same as an `` interiorpoint '' which is a set of coordinates in the cartesian space in which the algorithm is working . although the package code utilizes these geometric concept classes , there is good reason to provide user interfaces for common operations which deal with the familiar vector of doubles . for example , the userfunction ( ) method can accept a function of a std::vector where a point is logically called for :    ....    double f ( const std::vector < double > params ) ;    problem myprob ( userfunction ( nargs , f ) ) ; ....    the exterior concepts identified are point , gradient , secondderivatives , hessian , characteristicscale , and correlations . there are also classes corresponding to each of these , representing internal - coordinate concepts . the geometric concepts subsystem contains the collection of such classes . it depends on no other subsystem . the initial release of the minimization package is , by definition , to include all minuit functionality . of course , the object nature of a minimization problem will allow modes of usage which would be awkward with the fortran minuit , but basically the capabilities will at that point be those of minuit . however , certain extensions will be in the initial release . these appear in situations where the nature of a potential extension is obvious , and the implementation is straightforward . if the extension is clear in definition , and we have a firm idea of how it should by implemented , then there is no reason to delay that capability unitl a later release .    in the domain subsystem , the rectilineardomain class will immediately support half - open intervals for parameter limits ( that is , having only an upper or lower limit for a paramenter ) . it will also be easy to substitute for the arcsin mapping functions .    in the function subsystem , the package will allow supplying second derivatives ( either diagonal or full matrix ) along with the gradient .    in the termination subsystem , several new basic choices will be provided , including a time - based termination , and arbitrary logical combinations of termination conditions will be supported . a preliminary partial implementation of this package has been done , to obtain familiarity with the minuit algorithms and concepts , and to get a feel for the design issues . ( a spin - off of this work is a set of heavily - commented minuit fortran `` source '' code , and some mathematical documentation of the minuit algorithms . ) the cern effort ( as discussed with f. james ) is at the point where the minimization algorithms ( but not all of the analysis algorithms ) have first - pass implementions , again in a familiarization phase . since the two efforts are at comparable stages , the prospects for coordination seem reasonable .
Abstract:

Answer: a portion of the hep community has perceived the need for a minimization package written in c++ and taking advantage of the object - oriented nature of that langauge . to be acceptable for hep , such a package must at least encompass all the capabilities of minuit . aside from the slight plus of not relying on outside fortran compilation , the advantages that a c++ package based on o - o design would confer over the multitude of available c++ minuit - wrappers include : easier extensibility to different algorithms and forms opf constraints ; and usage modes which would not be available in the global - common - based minuit design . an example of the latter is a job persuing two ongoing minimization problems simultaneously . we discuss the design and implementation of such a package , which extends minuit only in minor ways but which greatly diminishes the programming effort ( if not the algorithm thought ) needed to make more significant extensions .


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: compact radio sources in nearby starburst galaxies often appear to be supernova remnants ( snr ) . snr in our galaxy and the large magellanic cloud ( lmc ) are sources of bright near - infrared [ feii ] line emission with luminosities ranging from @xmath5w to @xmath6w ( for the [ feii]@xmath0 m line , e.g. , oliva , moorwood , & danziger 1989 ; burton & spyromilio 1993 ) . this behavior suggests a connection between supernova activity and the [ feii ] emission in starbursts . this is supported by near - infrared spectroscopy ( moorwood & oliva 1988 ; greenhouse et al . 1991 ; vanzi & rieke 1997 ) , by high - resolution [ feii ] imaging studies ( forbes et al . 1993 ; greenhouse et al . 1997 ; alonso - herrero et al . 2000 ) , and by the tight correlation between the radio 6 cm and the [ feii ] emissions ( forbes & ward 1993 ) . further determination of the relation of the infrared [ feii ] emission with the supernova activity in nearby starforming galaxies would be helpful both for the understanding of these galaxies and for using the near - infrared [ feii ] luminosity as estimator of the supernova rate in more distant galaxies where individual radio supernovae can not be identified . there are indications of some complexities in this relation . for example , the bright [ feii ] emission may not be spatially coincident with the radio snr . lumsden & puxley ( 1995 ) conducted a spectroscopic survey for [ feii ] emission from radio snrs in m33 and found that in some cases the [ feii ] emission does not peak at the position of the radio snr , as well as that some radio snr show no [ feii ] emission . recently , morel , doyon , & st - louis ( 2002 ) have carried out an [ feii ] imaging survey of optically selected snr in m33 and found a similar result . both works have concluded that the [ feii ] emission from snr may be produced over a time scale of @xmath7 10@xmath8 years . greenhouse et al . ( 1997 ) found that none of the brightest sources of [ feii ] emission in m82 were coincident with the position of the radio snrs . m82 and ngc  253 are the nearest starburst galaxies ( distances @xmath9mpc and @xmath10mpc respectively ) . both galaxies show a large number of radio compact sources ( e.g. , @xmath11 in m82 kronberg et al . 1985 and huang et al . 1994 ; @xmath12 in ngc  253 turner & ho 1985 and ulvestad & antonucci 1991 , 1997 and references therein ) . their proximity makes these two galaxies ideal to study the relation between supernova activity and the [ feii ] emission in starburst galaxies . greenhouse et al . ( 1997 ) presented a fabry - perot [ feii]@xmath0 m map for m82 with a spatial resolution of @xmath13 . their main conclusions were that the compact [ feii ] sources contribute only @xmath14 of the total [ feii ] emission in m82 and that they trace a population of snr older than the radio snr . their work is somewhat inconclusive because of the limited angular resolution and also because strong extinction could hide [ feii ] emission from some of the radio snrs . this paper presents _ hst_/nicmos images of m82 and ngc  253 . the high spatial resolution of the nicmos [ feii]@xmath0 m images ( @xmath15 ; this corresponds to 5pc and 4pc for m82 and ngc  253 , respectively ) allows a detailed study of the location of the bright sources of [ feii ] and of the rle of supernova remnants in this line emission . _ hst_/nicmos observations were obtained on april 7 , 8 , 12 and 22 1998 ( for m82 ) and on august 13 and 20 1998 ( for ngc  253 ) using the nic2 and nic3 cameras . table  1 lists details of the observations . the observational strategy consisted of taking images in a spiral dither with a 5.5 pixel spacing with four positions for a given field . the plate scales for the nic2 and nic3 cameras are 0.076pixel@xmath16 and 0.20pixel@xmath16 , respectively . the fwhm of the point sources in the nic3 f166n continuum images is @xmath17 . lcclccc nic2 & f160w & @xmath18 & m82 & 192 & 4 & @xmath19 + & & & ngc  253 & 192 & 2 & @xmath20 + nic2 & f222 m & @xmath21 & m82 & 448 & 4 & @xmath19 + & & & ngc  253 & 512 & 2 & @xmath20 + nic2 & f187n & pa@xmath1 & m82 & 384 & 4 & @xmath19 + & & & ngc  253 & 384 & 2 & @xmath20 + nic2 & f190n & continuum & m82 & 384 & 4 & @xmath19 + & & & ngc  253 & 384 & 2 & @xmath20 + nic3 & f164n & [ feii]@xmath0 m & m82 & 4096 & 2 & @xmath22 + & & & ngc  253 & 2560 & 1 & @xmath23 + nic3 & f166n & continuum & m82 & 3840 & 2 & @xmath22 + & & & ngc  253 & 2560 & 1 & @xmath23 + nic3 & f215n & continuum & m82 & 2304 & 2 & @xmath22 + & & & ngc  253 & 2304 & 1 & @xmath23 +    the reduction of the nicmos images used routines from the package nicred ( mcleod 1997 ) . a master dark image was produced by combining between 10 and 20 darks for a given sample sequence after the subtraction of the first readout . the darks used were part of other proposals observed close in time , and the flatfield images were on - orbit data . the data reduction was performed using the following steps : subtraction of the first readout , dark current subtraction on a readout basis , correction for linearity and cosmic ray rejection ( using fullfit ) , and flatfielding . the individual dithered galaxy images were registered to a common position using fractional pixel offsets and cubic spline interpolation , and combined to produce the final image of each field . prior to the flux calibration of the images , the background needs to be subtracted ( only for wavelengths longwards of @xmath24 m ) . due to the large projected size of both galaxies , the field of view of the nic2 images is not large enough to enable measurements of the background on empty corners of the images , so we used background measurements taken during the servicing mission observatory verification ( smov ) program for the filter nic2 f222 m . the f215n images cover a larger field of view , and an estimate of the thermal background was obtained from the corners of the images . the flux calibration of the broad - band , on - line and off - line images was performed using the conversion factors ( from adu / s to mjy ) based on measurements of the standard star p330-e during smov ( marcia rieke 1999 private communication ) .    ideally one would like to measure the continuum on both sides of the emission line to estimate the continuum at the emission line . however , nicmos only provides continuum bands to the red of the emission line . in cases of objects such as m82 and ngc  253 , where the extinction to the nuclear regions is very high ( even at infrared wavelengths ) , a straight subtraction of the longer wavelength off - line image may result in an over correction of the continuum at the wavelength of the emission line . to estimate the continuum for the nic3 f164n images , we fitted the continuum on a pixel - by - pixel basis between @xmath25 m and @xmath26 m by linear regression of the flux as a function of the wavelength for the nic3 f166n and nic3 f215n line - free images . the extrapolated continuum was then subtracted from the nic3 f164n image to produce the final continuum - subtracted [ feii]@xmath0 m line image . this procedure does not correct for extinction , but it rather ensures the subtraction of the correct _ obscured _ continuum at each position of the galaxy . for the nic2 f187n image ( pa@xmath1 ) , since we did not have other narrow - band images taken with camera 2 , we used the nic2 f190n continuum image . fortunately , the continuum need not be subtracted as accurately for pa@xmath1 as for the [ feii ] emission line because the pa@xmath1 emission is much brighter . the fully - reduced images of each individual field were combined to produce the mosaics . the images presented in this paper were taken during different orbits . this resulted in slightly different orientations ( from @xmath27 to @xmath28 for m82 and from @xmath29 to @xmath30 for ngc  253 ) . before combining the images into mosaics , they were rotated by linear interpolation to north up , east to the left . this provided the same orientation for all the images , and more importantly allows comparisons with previously published data .    to construct extinction maps to correct the [ feii ] emission , and [ feii]/pa@xmath1 maps , the nic2 f187n , f190n , f160w and f222 m images were demagnified to match the plate scale of the nic3 images and filtered with a gaussian with @xmath31 to match the point spread function ( psf ) of the nic3 images . figure  1a and 1b are mosaics ( after rotation ) of the continuum at @xmath25 m and the continuum - subtracted [ feii]@xmath0 m line emission , respectively . both images are displayed on a logarithmic scale to stress the contrast between the diffuse extended continuum and the compact sources . the [ feii]@xmath0 m line map has been smoothed with a gaussian filter with @xmath32 . the line emission image shows some residuals which are due to both the under - sampling of the nic3 camera psf and the variation in the shape of the psf with wavelength . in addition , the right hand side of the line image shows a slightly lower s / n because the f164n and f166n images at that position had different orientations . the [ feii]@xmath0 m line emission exhibits two distinct components : a number of relatively compact regions distributed along the disk of the galaxy and an extended diffuse component in the form of filaments extending above and below the plane of the galaxy . the filaments seen in the [ feii]@xmath0 m emission line show a remarkable similarity with the chimneys at the base of the superwind detected in radio ( wills et al . 1999 ) , which are interpreted by these authors as expelled material from the central region of m82 .    in figures  2a and 2b we present contour maps of the @xmath25 m continuum emission and the [ feii]@xmath0 m line emission , respectively . comparison of the nicmos [ feii]@xmath0 m line map with that of greenhouse et al . ( 1997 ) reveals that some of the sources seen in the lower resolution map of greenhouse et al . ( 1997 ) break up into individual sources at higher spatial resolution . in particular , their source fe1 contains some 4 or 5 individual sources . source fe6 ( the most south - west source in greenhouse et al . 1997 map ) appears as a point source at our resolution in both the nic3 f164n and nic3 f166n images , and the signal to noise on it is compromised by the differing image orientations ( see above ) . therefore it is not surprising we do not detect line emission from it . the nicmos images detect a number of compact sources in addition to those reported in greenhouse et al . ( 1997 ) . in figures  2a and 2b we display the positions of 44 radio snrs with peak flux densities above 1mjy ( list provided by dr . karen wills , see more details in section  5.1 ) . the distribution of the radio snrs and their relation with the [ feii]@xmath0 m emission will be discussed in section  5 .    in figures  3a and 3b we display the mosaics of the @xmath33 m continuum ( nic2 f190n images ) and the pa@xmath1 line emission ( continuum subtracted nic2 f187n images ) respectively . the pa@xmath1 image is dominated by a ring of ionized gas with an approximate diameter of 400pc . the resemblance with the mid - infrared [ neii ] line emission map at @xmath34 m ( achtermann & lacy 1995 ) is remarkable . the south side of this ring of ionized material appears to be incomplete ( see achtermann & lacy 1995 ) , and it also shows a `` hole '' of emission at the center of the galaxy , which approximately coincides with the peak of the @xmath35 m surface brightness . ccccccc     + & ( ) & ( w m@xmath36 ) & ( w m@xmath36 ) & ( mag ) + pa@xmath1/h@xmath1 & 5.8 & @xmath37 & @xmath38 & @xmath39 & this work , oconnell & mangano ( 1978 ) + pa@xmath1/pa@xmath40 & 6 & @xmath37 & @xmath41 & @xmath42 & this work , mcleod et al . ( 1993 ) + pa@xmath1/pa@xmath43 & 3 & @xmath44 & @xmath45 & @xmath46 & this work , mcleod et al . ( 1993 ) + pa@xmath1/pa@xmath40 & 3 & @xmath44 & @xmath47 & @xmath42 & this work , mcleod et al . ( 1993 ) + br@xmath43/pa@xmath1 & 3.8 & @xmath48 & @xmath49 & @xmath50 & lester et al . ( 1991 ) , this work + br@xmath1/pa@xmath1 ( e1 ) & 6 ? & @xmath51 & @xmath52 & @xmath53 & achtermann & lacy ( 1995 ) , this work + br@xmath1/pa@xmath1 ( w1 ) & 6 ? & @xmath54 & @xmath55 & @xmath56 & achtermann & lacy ( 1995 ) , this work + br@xmath1/pa@xmath1 ( w2 ) & 6 ? & @xmath57 & @xmath58 & @xmath59 & achtermann & lacy ( 1995 ) , this work + @xmath60feii]1.64/[feii]1.25 & 3 & @xmath61 & @xmath62 & @xmath46 & this work , mcleod et al . ( 1993 ) +   + & ( ) & ( mjy ) & ( mjy ) & ( mag ) + f222m / f160w & 3 & 237 & 166 & @xmath56 & this work +    ccccccc + & ( ) & ( w m@xmath36 ) & ( w m@xmath36 ) & ( mag ) + pa@xmath1/br@xmath43 ( center ) & @xmath63 & @xmath64 & @xmath65 & @xmath66 & this work , harrison et al . ( 1998 ) + pa@xmath1/br@xmath43 ( center ) & @xmath67 & @xmath68 & @xmath69 & @xmath70 & this work , puxley & brand ( 1995 ) + pa@xmath1/br@xmath43 ( 3  ne ) & @xmath63 & @xmath71 & @xmath72 & @xmath73 & this work , harrison et al . ( 1998 ) + pa@xmath1/br@xmath43 ( 3.9  ne ) & @xmath67 & @xmath74 & @xmath75 & @xmath76 & this work , puxley & brand ( 1995 ) +   + & ( ) & ( mjy ) & ( mjy ) & ( mag ) + f222m / f160w ( center ) & @xmath67 & 213 & 142 & @xmath56 & this work +    it has been suggested that higher extinction in the south half of m82 may be causing the uneven distribution of ionized gas in the ring . however , the ring appears to be complete in [ feii]@xmath0 m line emission . if extinction were the cause , we should expect to see a similar morphology in the [ feii ] emission , which is more affected by extinction than the pa@xmath1 emission . it is therefore tempting to suggest that the south part of the ring of ionized gas is actually _ broken _ , perhaps due to sn explosions . this hypothesis would also explain the excess of [ feii ] emission in this region . the presence of the hole of gas emission at the center of m82 along with the anti - correlation between the morphology of the hydrogen recombination lines and the co index has been invoked as evidence for an outwards propagating starburst ( e.g. , satyapal et al . the diffuse component of pa@xmath1 emission is also in the form of `` fingers '' or `` filaments '' above and below the plane of the galaxy . it is now known from ground - based h@xmath1 imaging that this emission along the minor axis of the galaxy extends over 11kpc ( devine & bally 1999 ) . figure  1b and figure  3b show how there is a very close spatial correspondence between the morphologies of the extended [ feii]@xmath0 m and pa@xmath1 emissions . figures  4a and 4b are the @xmath25 m continuum and the [ feii]@xmath0 m line emission maps of ngc  253 . as seen in m82 , the [ feii]@xmath0 m emission is in the form of a number of compact regions along the disk of the galaxy superimposed on a diffuse emission component . figures  5a and 5b are contour plots of these two emissions . we also show the positions of the 6 cm radio sources with @xmath77mjy from ulvestad & antonucci ( 1991 ) , excluding those sources that are classified as hii regions in ulvestad & antonucci ( 1997 ) . ground - based [ feii]@xmath0 m and br@xmath43 maps at lower spatial resolution ( @xmath78 ) were previously observed by forbes et al . they found that the central [ feii ] emission of ngc  253 is dominated by three peaks of emission . our higher resolution nicmos image shows a larger number of [ feii ] compact sources . for instance , source b ( in the forbes et al . 1993 notation ) breaks up into three individual sources . in addition to the nucleus and sources a and b detected by forbes et al . ( 1993 ) , there are other sources along the disk of the galaxy , as well as diffuse emission perpendicular to the disk , more prominent in the south half of the galaxy . we believe that forbes et al . ( 1993 ) did not detect the diffuse component above and below the disk of the galaxy because of the lower s / n of their [ feii ] map . we also mark in these figures the position of the peak ( the filled square ) of the near - infrared emission ( also coincident with the peak of the pa@xmath1 emission ) which was identified as the nucleus of the galaxy by forbes et al . other authors ( e.g. , keto et al . 1999 ) claim that a bright non - thermal radio source located ne of the peak of the infrared continuum emission ( source 2 in turner & ho 1985 ) may be the true nucleus of ngc  253 . it is shown in these figures as a filled star symbol . the nicmos pa@xmath1 line map ( figure  6b ) shows very bright emission stemming from the peak of infrared continuum emission , together with more diffuse emission in a plume - like morphology along the plane of the galaxy , primarily to the ne of the nucleus . as for m82 , there is an extended component of the pa@xmath1 emission perpendicular to the disk , better seen in our nic3 images ( not shown here ) which are more sensitive to diffuse emission . this component is however not as strong as that in m82 . there is a remarkable resemblance of the central pa@xmath1 emission with both the [ neii]@xmath34 m emission line and the nearby continuum images presented in bker , krabbe , & storey ( 1998 ) and keto et al . both m82 and ngc  253 are known to be highly obscured ( see for instance the comparison between the optical _ hst_/wfpc2 f555w and the _ hst_/nicmos f222 m images of m82 in alonso - herrero et al . 2001 ) . to calibrate the supernova rate in terms of the [ feii]@xmath0 m emission and to obtain the total number of ionizing photons , it is necessary to correct the observed quantities for extinction . we can estimate the extinction to the gas in m82 using observations of the ratios of hydrogen recombination and near - infrared [ feii ] lines . we used published values of the line fluxes for the smallest apertures and matched those apertures on the nicmos line images when measuring the pa@xmath1 and [ feii]@xmath0 m fluxes . the intrinsic ratios of hydrogen recombination lines are from hummer & storey ( 1987 ) for an electron temperature of @xmath79k and electron density of @xmath80@xmath81 . the ratio between [ feii]@xmath82 m and [ feii]@xmath0 m does not depend upon the physical conditions of the gas . we used a polynomial fit to the near - infrared extinction law given in he et al . ( 1995 ) .    since we do not have nicmos images of two hydrogen recombination lines , to measure extinction at the full nicmos angular resolution we made use of the broad - band filters nic2 f160w and nic2 f222 m to construct a map of extinction to the stars in m82 . we assumed a color for the unreddened stellar population of @xmath83 ( which applies almost independently of the age of the population ) , and , for computational convenience , a foreground dust screen model following the extinction law of he et al . the contribution of dust emission at the @xmath21-band was assumed to be negligible . figure  7a is the extinction map ( @xmath84 ) from the @xmath85 color map of m82 , with the positions of the radio sources superimposed . the values of the extinction in the @xmath21-band derived from this map are between @xmath86mag and @xmath87mag . the m82 extinction map qualitatively looks very similar to those of satyapal et al . ( 1995 ) , but quantitatively the values of the extinction are different . the estimates of the extinction to the gas for m82 are presented in the first part of table  2 along with the line fluxes . the errors in table  2 are the uncertainties associated with the aperture size , aperture centering and flux calibration . for comparison , we also estimated the extinction through a small aperture ( 3 ) using the f222m / f160w continuum flux ratio , which gives an estimate of the visual extinction to the stars . as can be seen from this table the agreement between the extinctions to the gas and to the stars is excellent at infrared wavelengths . satyapal et al . ( 1995 ) reached a similar conclusion . a much lower level of extinction is indicated by comparing the h@xmath1 and pa@xmath1 line fluxes . this behavior is known to arise from the different optical depths probed by optical and near - infrared lines . this behavior is confirmed by the much higher extinction of @xmath88mag inferred by puxley et al . ( 1989 ) from the radio h53@xmath1 to mid - infrared br@xmath1 line ratio . we also computed the extinctions to the e1 , w1 and w2 regions in achtermann & lacy ( 1995 ) notation ( see table  2 for the positions ) . these were found to be bright in [ neii]@xmath34 m and br@xmath1 . we used their br@xmath1 measurements along with the mid - infrared extinction curve in rieke & lebofsky ( 1985 ) . since achtermann & lacy ( 1995 ) did not quote aperture sizes , we assumed that their fluxes were measured over the extent of these regions ( approximately @xmath89  in diameter ) . the values of the extinction for these three regions are slightly larger than the values we would measure from the @xmath85 extinction map , perhaps indicating minor differences in the optical depth . a further test of the extinction to the gas was based on the total number of ionizing photons from the extinction corrected pa@xmath1 line emission map . the number of ionizing photons obtained from the extinction - free h53@xmath1 radio recombination line is @xmath90s@xmath16 ( puxley et al . 1989 ) , and from the @xmath91 mm continuum is @xmath92s@xmath16 ( carlstrom & kronberg 1991 ) . we measured through a large aperture the pa@xmath1 flux on the observed and the extinction - corrected images . we obtained @xmath93s@xmath16 ( observed ) and @xmath94s@xmath16 ( corrected for extinction ) , for electron temperature @xmath95k . the total number of ionizing photons corrected for extinction is in excellent agreement with a number of estimates including those of puxley et al . ( 1989 ) , carlstrom & kronberg ( 1991 ) and frster schreiber et al . ( 2001 ) .      as for m82 , the extinction to the stars was computed from the @xmath85 color map of ngc  253 ( figure  7b ) . it shows values between @xmath96mag and @xmath97mag and is very similar to that of sams et al . ( 1994 ) derived from a @xmath98 color map . we find a similar value of the extinction ( @xmath99mag ) for the point with the highest value of the extinction which was identified by sams et al . ( 1994 ) as the true nucleus of ngc  253 .    to calculate the extinction to the gas in ngc  253 we used pa@xmath1 and br@xmath43 fluxes ( the latter from puxley & brand 1995 and harrison et al . 1998 ) . both works centered their apertures on the approximate coordinates of the secondary near - infrared peak and therefore we will assume that the largest br@xmath43 fluxes in these works correspond to the primary near - infrared peak . if their ( 0,0 ) position were to be identified with the peak of the continuum emission , their fluxes would not be consistent with forbes et al . ( 1993 ) br@xmath43 line emission map . this map and our pa@xmath1 line emission map ( figure  6b ) show that the brightest hydrogen line emission arises from the peak of the continuum emission . in table 3 we give estimates for the extinction for the primary near - infrared peak ( center ) and secondary near - infrared peak ( positions 3  and 3.9  to the ne of the primary peak for puxley & brand 1995 and harrison et al . 1998 works respectively ) . the errors take into account 10% and 20% combined flux calibration / centering uncertainties for the nucleus and ne source , respectively . the values of the gas extinction to these two regions as derived from the pa@xmath1/br@xmath43 line ratio are consistent ( to within the errors ) with the extinction to the stars from the @xmath85 color map . despite this apparent good agreement we caution the reader of the aforementioned centering uncertainties . we also measured the observed and corrected for extinction pa@xmath1 fluxes which yielded a number of ionizing photons of @xmath100 ( observed ) and @xmath101 ( corrected for extinction ) , for the field of view shown in figure  6b . the latter value is consistent with engelbracht et al . ( 1998 ) and keto et al . ( 1999 ) estimates of the number of ionzing photons from their large aperture br@xmath43 and [ neii ] measurements , but smaller than puxley et al . ( 1997 ) estimate ( @xmath102 ) from millimeter - wavelength hydrogen recombination lines . to determine whether there is a spatial correspondence between the compact radio sources and the [ feii ] sources in both m82 and ngc  253 , very precise coordinates are needed . unfortunately , due to problems with the _ hst _ telemetry during the observations of m82 , and the uncertainties when producing and rotating the mosaics , the coordinates of the nic3 images may not be completely reliable . in addition , it is expected that there are minor discrepancies between the radio and the optical / infrared coordinate reference systems . our approach was to determine the coordinates of the [ feii ] image by trying to find the highest possible number of coincidences between the positions of the radio snrs and the [ feii ] sources . we use daofind within the daophot package in iraf to find compact sources . the term `` compact '' is not used in the sense of point sources . in fact , most of these sources seem to be slightly resolved , although due to the under - sampling of the nic3 psf , some uncertainties remain . the sources were selected to be easily distinguishable from the extended diffuse emission detected above and below the plane of the galaxy . in m82 we identified 65 `` compact '' sources in the [ feii]@xmath0 m map ( not corrected for extinction ) . for the radio sources in m82 we used the positions of snrs identified by searching for compact features within radio maps at 408mhz , 1.4ghz , 5ghz and 8.4ghz , with peak flux densities above 1mjy ( dr . karen wills , private communication 2000 ) . we excluded from this list radio sources that are classified as hii regions based on their radio spectral index ( see mcdonald et al . 2002 for details ) , resulting in 44 radio snrs . we began the comparison of positions using the coordinates in the header of the nicmos nic3 f164n images and made small offsets of the [ feii ] image , cross - correlating the positions of the [ feii ] sources with the radio snrs . based on the resolution of the nic3 images and the average angular size of the radio snrs ( @xmath103 , huang et al . 1994 ) , we counted a coincidence when the angular separation between a radio source and an [ feii ] source was @xmath104 . with a small shift , we found 21 coincidences between radio snrs and [ feii ] sources , whereas if we used the original coordinates we could only find 7 coincidences . monte carlo simulations showed that the probability of finding 21 coincidences by pure chance is @xmath105 ( see whitmore & zhang 2002 for details on the method ) , so we believe that now the infrared and the radio images are properly aligned . for ngc  253 we found discrepancies between the coordinates derived for the near - infrared peak of emission from the nic2 and nic3 images ( which were taken during two different orbits ) . although no problems were reported for these two sets of observations , some error must be affecting the coordinates given in the headers of the images . again as with m82 , we cross - correlated the positions of the [ feii]@xmath0 m ( 36 sources ) and 35 radio sources taken from ulvestad & antonucci 1991 excluding those sources that are classified as hii regions in ulvestad & antonucci ( 1997) using the coordinates of the nic3 images . we find a total of 9 coincidences for a separation between the radio sources and [ feii ] sources @xmath106 . the probability of finding 9 coincidences in by pure chance in ngc  253 is approximately @xmath107 . one important conclusion , despite the uncertainties of the image alignment , is that we can not assign an [ feii ] source to each radio snr or vice versa . we find approximately @xmath3% of radio snrs in m82 and ngc  253 have an [ feii ] . this point is illustrated in the contour images of the [ feii]@xmath0 m line emission and continuum emission at @xmath25 m in figures  2a and 2b for m82 and figures  5a and 5b for ngc  253 . in all figures we show the positions of the radio snrs as asterisks . we also tried to determine if there is a relation between the presence of a coincidence and the age of the radio snr . such a correlation would be expected if the [ feii ] emission were brighter in older supernova remnants . for m82 we used the angular sizes of the radio snr as an indication of their age ( see huang et al . in general , we did not find that older radio snrs tend to have a coincidence with an [ feii ] source . we performed aperture photometry on the [ feii ] sources identified in m82 image registered to match the field of view and position of the pa@xmath1 image , the number of detected sources ( 71 ) is slightly different from that in section  5.1 ] and ngc  253 using a fixed @xmath108 diameter aperture . this aperture corresponds to a physical size of 12pc and 10pc for m82 and ngc  253 , respectively . figures  9a and 9b show the distributions of the measured [ feii]@xmath0 m luminosities ( not corrected for extinction ) . some of the sources may be larger than the aperture used for the photometry ; on the other hand , patchy extinction will make some of the sources appear smaller . from inspection of typical source sizes , however , the luminosities would be accurate on average . recent studies of galactic snr show that the [ feii ] emission is filamentary and have sizes of the order of 11pc ( e.g , ic443 rho et al . 2001 ; see also morel et al . 2002 for snr in m33 ) .    for the assumed distances to m82 and ngc  253 the observed [ feii]@xmath0 m luminosities of these sources range from @xmath109w to @xmath110w . as a consistency check for the photometry , we measured the [ feii]@xmath0 m fluxes for the fe1 , fe2 and fe3 sources in m82 identified by greenhouse et al . ( 1997 ) using the aperture diameters listed in their table  1 . the agreement with greenhouse et al . ( 1997 ) fluxes for these sources is excellent , to within @xmath111 . however , when we compared our photometry for the sources reported in forbes et al . ( 1993 ) for ngc  253 we found some discrepancies .    the corrected [ feii]@xmath0 m luminosities for compact sources in m82 and ngc  253 range from @xmath112w to @xmath113w . as a comparison , galactic snrs yield [ feii]@xmath0 m luminosities of @xmath114w to @xmath115w ( oliva et al . 1989 ; keller et al . 1995 ) , snrs in the lmc yield luminosities from @xmath116w to @xmath112w ( oliva et al . 1989 ) and snrs in m33 ( assuming a distance of 720kpc ) @xmath117 to @xmath112w ( lumsden & puxley 1995 and morel et al . 2002 ) . the individual sources account for some 22% and 27% of the total [ feii]@xmath0 m emission ( extinction corrected ) measured in m82 and ngc  253 , respectively . given the estimates that the [ feii ] emitting phase of a snr has a duration of @xmath118yr ( lumsden & puxley 1995 ; morel et al . 2002 ) and the supernova rates in m82 and ngc  253 are @xmath119yr@xmath16 ( huang et al . 1994 and ulvestad & antonucci 1997 ) , the detected [ feii ] sources account for @xmath120yr worth of supernovae . presumably much of the remaining unresolved [ feii ] arises from additional supernovae that are packed so densely that they are confused in our images . also , snrs in these two galaxies may expand and merge into a general hot ism within a few thousand years losing their identity as individual sources . photometry on the [ feii ] emission arising from the disk of m82 shows that as much as 70% of the total emission ( disk plus component perpendicular to the disk ) is associated with snrs . this together with the good agreement between the radio and [ feii ] morphologies of the components perpendicular to the disk gives further indication for a snr origin of the [ feii ] emission in starburst galaxies . lumsden & puxley ( 1995 ) pointed out that one possibility is that the [ feii ] lines reach maximum luminosity when the remnant as a whole undergoes the transition from adiabatic expansion to radiative expansion . according to shull & draine ( 1986 ) , this happens when radiative cooling occurs faster than dynamical times , around @xmath121yr after the supernova explosion . therefore if the brightest [ feii ] emission occurs at approximately @xmath121yr , it is not surprising that we do not find a correspondence for all the [ feii ] sources and all the radio snrs . muxlow et al . ( 1994 ) found that half of the radio sources in m82 are unresolved indicative of a young age , but the other resolved half are probably intermediate age . this together with the fact that we found that @xmath122 of the radio snr have an [ feii ] counterpart implies that the [ feii ] sources may be @xmath123yr . one possibility for the lack of one - to - one coincidences between the radio snrs and the [ feii ] sources would be that most of the [ feii ] emission in these two galaxies is produced in hii regions ( photoionization ) . to check this , we measured the pa@xmath1 fluxes for the sources detected in the [ feii ] extinction corrected map of m82 , and assumed an [ feii]@xmath0m / pa@xmath1 ratio for hii regions as that observed in orion . we found that only @xmath124 of the total [ feii]@xmath0 m arises from hii regions . this further supports the hypothesis that most of the [ feii]@xmath0 m emission along the plane of m82 is produced by snrs . the [ feii]/br@xmath43 ratios in hii regions and snr differ by a factor of a hundred or more , and therefore can be used to indicate the excitation mechanism of the [ feii ] line within our images . in alonso - herrero et al . ( 1997 ) we showed that the behavior of the [ feii]@xmath0m / br@xmath43 line ratio in starbursts and active galaxies can be understood as a progression from pure photoionization to pure shock excitation going from pure hii regions through starbursts , transition objects , seyferts to supernova remnants .    figures  8a and 8b are the observed [ feii]@xmath0m / pa@xmath1 maps for m82 and ngc  253 , respectively . the values of the [ feii]@xmath0m / pa@xmath1 ratio range from 0.02 in black ( mostly along the disk of the galaxy ) to 0.8 in white ( some of the bright [ feii ] sources ) . these ratios translate into [ feii]@xmath0m / br@xmath43 line ratios of between 0.24 and 10 ( see also the histograms in figure  10 ) . since these maps are not corrected for extinction , these numbers are actually lower limits . when the positions of radio sources are plotted in the [ feii]@xmath0m / pa@xmath1 maps we find again that there is a number of coincidences between the radio snrs and the regions of enhanced [ feii]@xmath0 m emission . this indicates that the relative contribution of the snrs is nearly as large in the diffuse background as it is in the radio snr . interestingly , we find that most of the coincidences with the radio snrs occur in regions of enhanced [ feii]@xmath0m / pa@xmath1 line ratio , or conversely we do not find a large number of coincidences in regions with high pa@xmath1 flux ( that is , hii regions ) . the supernova rate in m82 has been estimated from number counts of the radio supernova remnants and their ages to be @xmath125yr@xmath16 ( huang et al . the total [ feii]@xmath0 m luminosity in m82 is @xmath126([feii]@xmath127w . our value for the [ feii ] luminosity of m82 is approximately 1.5 and 2 times higher than the previous estimates by greenhouse et al . ( 1997 ) and vanzi & rieke ( 1997 ) , respectively . this difference is most likely due to better extinction estimates made possible by our nicmos images . if the supernova rate in m82 is @xmath128yr@xmath16 , then we can derive the following calibration for the supernova rate in terms of the [ feii]@xmath0 m luminosity :    @xmath129}}{10^{34 } { \rm w}}\,\,\ , ( { \rm yr}^{-1})\ ] ]    the total [ feii]@xmath0 m luminosity ( corrected for extinction ) measured in ngc  253 is @xmath130 . an upper limit on the supernova rate from radio supernova counts is @xmath131 ( ulvestad & antonucci 1997 ) . a rate of @xmath132 has been calculated assuming that the non - thermal radio flux is entirely due to synchrotron emission ( van buren & greenhouse 1994 ) ; this value is consistent with the similar number of compact radio sources in this galaxy compared with m82 . the calibration of supernova rate in terms of [ feii]@xmath0 m luminosity for this galaxy is then    @xmath133}}{10^{34 } { \rm w}}\,\,\ , ( { \rm yr}^{-1})\ ] ]    since m82 is the more thoroughly studied galaxy in the radio , we give higher weight to the calibration for it in averaging these results to estimate    @xmath134}}{10^{34 } { \rm w}}\,\,\ , ( { \rm yr}^{-1})\ ] ]    in the above calculations , we have included the entire [ feii ] luminosity , since that is the parameter usually measured in a distant galaxy . however , in these nearby examples , we can compare the [ feii ] emission produced in hii regions with that from supernovae . we assume that the ratio for a typical hii region is @xmath135([feii]@xmath0m@xmath136 ( e.g. , orion ) . we obtained that the [ feii ] luminosity produced in hii regions accounts for less than 8% the total [ feii ] luminosity in m82 . we have presented _ hst_/nicmos observations of the starburst galaxies m82 and ngc  253 . the morphology of the [ feii]@xmath0 m line emission in both galaxies displays two components : ( 1 ) a disk component with a number of superimposed compact sources and ( 2 ) a diffuse component in the form of filaments above and below the disk of the galaxy . in m82 we find a remarkable similarity between the extended component , the fingers or chimneys above and below the disk , of the radio and [ feii ] emissions .    in m82 , we have detected in both pa@xmath1 and [ feii ] a nuclear ring of star formation . this ring , with an approximate diameter of 400pc , has also been observed in other near- and mid - infrared hydrogen recombination lines and in the [ neii]@xmath34 m line . the peak of the @xmath35 m continuum brightness of m82 lies very close to ( within 1 ) the dynamical center , and corresponds to the center of the ring of star formation . in ngc  253 , bright pa@xmath1 emission stems from the near - infrared continuum peak and extends in a plume - like morphology along the disk of the galaxy . there is also a diffuse pa@xmath1 component above and below the plane of galaxy , however , much fainter than that in m82 .    prior to calibrating the surpernova rate in terms of the [ feii ] emission , we have derived the extinction to the central regions of ngc  253 and m82 . to assess the extinction correction , we have compared the number of ionizing photons infered from the extinction - corrected pa@xmath1 maps with other estimates based on longer wavelength indicators , and found a good agreement .    since the [ feii ] emission in starburst galaxies is thought to be associated with supernova activity , we have cross - correlated the positions ( with a spatial resolution of a few tenths of arcseconds ) of the radio snrs and the identified compact [ feii ] sources in m82 and ngc  253 . we have found that approximately @xmath3% of the radio snrs have an [ feii ] counterpart . however , we can not associate an [ feii ] source to each radio snr or vice versa . this is well understood in terms of the different ages traced by the radio snrs ( a few hundred years ) and the [ feii ] stage of snrs ( @xmath4yr ) . the compact [ feii ] emitting sources in m82 and ngc  253 only contribute to some @xmath137 of the total [ feii]@xmath0 m emission . however , much of the remaining emission in the disk of the two galaxies is most likely produced by snrs that expanded and merged into a general ism a few @xmath138yr ago . photometry on the [ feii]@xmath0 m emission arising from the disk of m82 shows that as much as @xmath139 of the total [ feii]@xmath0 m emission ( corrected for extinction ) is associated with snrs . the remaining diffuse [ feii ] emission perpendicular to the disk of m82 seems to be related to the superwind ( via expulsion of material from the central starburst ) , as argued from the radio emission morphology ( wills et al . 1999 ) . all these arguments support the supernova origin for the observed near - infrared [ feii ] line emission in starburst galaxies . thus we have recalibrated the supernova rate in terms of the [ feii ] luminosity relation for starburst galaxies using our newly measured ( corrected for extinction ) [ feii]@xmath0 m luminosities of m82 and ngc  253 . we would like to thank drs . karen wills , t. muxlow and a. pedlar for kindly providing us with the radio images for comparison with the nicmos images , and the list of radio snr . we are also grateful to dr . ernesto oliva for enlightening discussions , and an anonymous referee for constructive comments that helped improve the paper . during the course of this work aa - h was supported by the national aeronautics and space administration on grant nag 5 - 3042 through the university of arizona . the work was also partially supported by the national science foundation under grant ast-95 - 29190 . satyapal , s. , watson , d. m. , pipher , j. l. , forrest , w. j. , coppenbarger , d. , raines , s. n. , libonate , s. , pich , f. , greenhouse , m. a. , smith , h. a. , thompson , k. l. , fisher , j. , woodward , c. e. , & hodge , t. 1995 , apj , 448 , 611
Abstract:

Answer: we present _ hst_/nicmos [ feii]@xmath0 m , pa@xmath1 ( @xmath2 m ) and continuum images of the starburst galaxies m82 and ngc  253 at an unprecedented spatial resolution . in both galaxies we detect [ feii ] compact sources superimposed on a diffuse background in the disk of the galaxies together with a component above and below the plane of the galaxy . the radio and [ feii ] emissions perpendicular to the disk of m82 show a remarkable similarity to each other , suggesting that both emissions originate in shocks from supernova explosions . we find a spatial correspondence between bright _ compact _ [ feii ] emitting regions and the location of radio supernova remnants ( snr ) for approximately @xmath3% of radio snrs in m82 and ngc  253 . this lack of a one - to - one correspondence , more than being indicative of a different origin for the radio and [ feii ] emission in starbursts , suggests two populations of snrs : an older population ( @xmath4yr ) traced by the [ feii ] emission and a younger population ( a few hundred years old ) traced by the radio snrs . we therefore conclude that the [ feii ] emission in starburst galaxies provides a good estimate of the supernova activity . using our newly determined [ feii ] luminosities ( corrected for extinction ) of m82 and ngc  253 , we reevaluate the calibration of the supernova rate in terms of the [ feii ] luminosity for starburst galaxies .