From owner-reliable_computing [at] interval [dot] louisiana.edu Wed May 1 10:49:46 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g41Fnja09912 for reliable_computing-outgoing; Wed, 1 May 2002 10:49:45 -0500 (CDT) Received: from marnier.ucs.louisiana.edu (root [at] marnier [dot] ucs.louisiana.edu [130.70.132.233]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g41Fnch09908 for ; Wed, 1 May 2002 10:49:39 -0500 (CDT) Received: from liberty (liberty.louisiana.edu [130.70.46.171]) by marnier.ucs.louisiana.edu (8.11.3/8.11.3/ull-ucs-mx-host_1.6) with SMTP id g41FnbZ18273 for ; Wed, 1 May 2002 10:49:37 -0500 (CDT) Message-Id: <2.2.32.20020501155857.015688ac [at] 130 [dot] 70.132.231> X-Sender: rbk5287 [at] 130 [dot] 70.132.231 X-Mailer: Windows Eudora Pro Version 2.2 (32) Mime-Version: 1.0 Content-Type: text/plain; charset="us-ascii" Date: Wed, 01 May 2002 10:58:57 -0500 To: reliable_computing [at] interval [dot] louisiana.edu From: "R. Baker Kearfott" Subject: New preprint available: verification of minimax Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Colleagues, I would like to call your attention to a new short preprint, "On Verification of Solutions to Nonlinear Minimax Problems" that you can download from http://interval.louisiana.edu/preprints/2002_nonlinear_minimax.ps or http://interval.louisiana.edu/preprints/2002_nonlinear_minimax.dvi or http://interval.louisiana.edu/preprints/2002_nonlinear_minimax.pdf In it, I show that, if Lemarechal's suggestions are used to convert the problem into a smooth problem with constraints, the corresponding linear problem must have Fritz-John matrix that is multiply singular at the solution. Thus, that avenue cannot be used for solution verification. I propose an alternate technique, based on Jansson's work. If anyone knows of similar work in the literature, I would appreciate knowing of it. We're presently investigating the "alternate technique" :-) Best regards, Baker --------------------------------------------------------------- R. Baker Kearfott, rbk [at] louisiana [dot] edu (337) 482-5346 (fax) (337) 482-5270 (work) (337) 981-9744 (home) URL: http://interval.louisiana.edu/kearfott.html Department of Mathematics, University of Louisiana at Lafayette Box 4-1010, Lafayette, LA 70504-1010, USA --------------------------------------------------------------- From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 2 06:50:14 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g42BoDW12209 for reliable_computing-outgoing; Thu, 2 May 2002 06:50:13 -0500 (CDT) Received: from mailbox.univie.ac.at (mailbox.univie.ac.at [131.130.1.27]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g42Bo6h12205 for ; Thu, 2 May 2002 06:50:07 -0500 (CDT) Received: from univie.ac.at (hektor.mat.univie.ac.at [131.130.16.21]) by mailbox.univie.ac.at (8.12.2/8.12.2) with ESMTP id g42BnwKH1338958; Thu, 2 May 2002 13:50:00 +0200 Message-ID: <3CD127E7.AA54A4A5 [at] univie [dot] ac.at> Date: Thu, 02 May 2002 13:49:59 +0200 From: Arnold Neumaier Organization: University of Vienna X-Mailer: Mozilla 4.78 [en] (X11; U; Linux 2.4.9-31 i686) X-Accept-Language: en, de MIME-Version: 1.0 To: "R. Baker Kearfott" CC: reliable_computing [at] interval [dot] louisiana.edu Subject: Re: verification of minimax References: <2.2.32.20020501155857.015688ac [at] 130 [dot] 70.132.231> Content-Type: text/plain; charset=iso-8859-15 Content-Transfer-Encoding: 7bit Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk "R. Baker Kearfott" wrote: > > I would like to call your attention to a new short preprint, > "On Verification of Solutions to Nonlinear Minimax Problems" > If anyone knows of similar work in the literature, I would appreciate > knowing of it. Shen Zuhe, A. Neumaier and M.C. Eiermann, Solving minimax problems by interval methods, BIT 30 (1990), 742-751. L. Jaulin (2001). Reliable minimax parameter estimation. Reliable Computing, issue 3, volume 7. 231-246. preprint from: http://www.istia.univ-angers.fr/~jaulin/publications.html Jaulin L., M. Kieffer, O. Didrit and E. Walter, Applied Interval Analysis, Springer, 2001 also contains some material on minimax problems. From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 2 11:48:06 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g42Gm6812837 for reliable_computing-outgoing; Thu, 2 May 2002 11:48:06 -0500 (CDT) Received: from kathmandu.sun.com (kathmandu.sun.com [192.18.98.36]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g42Glxh12833 for ; Thu, 2 May 2002 11:48:00 -0500 (CDT) Received: from engmail1.Eng.Sun.COM ([129.146.1.13]) by kathmandu.sun.com (8.9.3+Sun/8.9.3) with ESMTP id KAA11381; Thu, 2 May 2002 10:47:31 -0600 (MDT) Received: from phys-mpkmaila (phys-mpkmaila.Eng.Sun.COM [129.146.18.131]) by engmail1.Eng.Sun.COM (8.9.3+Sun/8.9.3/ENSMAIL,v2.1p1) with ESMTP id JAA15268; Thu, 2 May 2002 09:47:30 -0700 (PDT) Received: from gww (gww.Eng.Sun.COM [129.146.78.116]) by mpkmail.eng.sun.com (iPlanet Messaging Server 5.2 Patch 1 (built Apr 2 2002)) with SMTP id <0GVH00AM9TDQXL [at] mpkmail [dot] eng.sun.com>; Thu, 02 May 2002 09:49:02 -0700 (PDT) Date: Thu, 02 May 2002 09:47:30 -0700 (PDT) From: William Walster Subject: Re: Manifolds To: tupper [at] peda [dot] com, Arnold.Neumaier [at] univie [dot] ac.at Cc: reliable_computing [at] interval [dot] louisiana.edu Reply-to: William Walster Message-id: <0GVH00AMATDQXL [at] mpkmail [dot] eng.sun.com> MIME-version: 1.0 X-Mailer: dtmail 1.3.0 @(#)CDE Version 1.4.2 SunOS 5.8 sun4u sparc Content-type: TEXT/plain; charset=ISO-8859-1 Content-MD5: B2WdSITYG3vRkRflGFU4Yw== Content-Transfer-Encoding: 8bit X-MIME-Autoconverted: from QUOTED-PRINTABLE to 8bit by interval.louisiana.edu id g42Gm1h12834 Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Friends, According to James & James Mathematics Dictionary, Fourth Edition, there appears to be room for the word "manifold" to be used in a variety of ways. They say in part: Manifold: In general, manifold may mean any collection or set of objects. E.g., a Riemannian space is also called a Riemannian manifold; a subset of a vector space is said to be a *linear set* or a *linear manifold* if it contains all linear combinations of its members. However, manifold frequently has technical meaning beyond being a mere set, as illustrated by the following definitions: They go on to list a number, which I will not copy, but just list: topological manifold; pseudomanifold; orientable manifold non-orientable manifold. My interpretation is that as long as we are clear how we are using the term, we should not feel inhibited from using it in the way George and others have done in the past. Cheers, Bill >Date: Thu, 18 Apr 2002 10:21:24 +0200 >From: Arnold Neumaier >Subject: Re: Manifolds >To: Jeff Tupper / Pedagoguery Software Inc >Cc: reliable_computing [at] interval [dot] louisiana.edu >MIME-version: 1.0 >Content-transfer-encoding: 8BIT >X-Accept-Language: en, de > >Jeff Tupper wrote: > >> >When verifying solutions, one can verify only the nonsingular part >> >of the solution set (its relative interior), and this always forms a >> >manifold. But sometimes (outside theorems), there is a sloppy use of >> >manifold as a synonym for solution set. >> >> It seems to me that if, for example, Bolzano's theorem, or the >> (Poincaré-)Miranda theorem, is used by an interval method, that what >> is known is that the portion of the solution set just verified is >> known to be closed (not necessarily a manifold). > >To be 100% precise, what is verified is always a closed subset of a >manifold (without boundary). > >> More abstractly, I don't see why a special-purpose verifier couldn't >> be built for one particular (non-manifold) solution set. (I.e. look >> up the formula in a table and give back the exact solution set.) > >Unless the singular points are exactly representable and represented >(which is the case only in *very* few problems from applications), >and no rounding error occurs in the verification phase (which excludes >practically all interval techniques and necessitates symbolic methods) >they cannot be verified, since they are unstable under perturbations. >(Look at double roots or bifurcation points...) >One can produce valid enclosures for the region, but no associated >existence statement. > >Only generic situations can be verified by general purpose >interval techniques. > >Arnold From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 2 12:01:21 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g42H1JL12953 for reliable_computing-outgoing; Thu, 2 May 2002 12:01:19 -0500 (CDT) Received: from cs.utep.edu (mail.cs.utep.edu [129.108.5.3]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g42H1Fh12949 for ; Thu, 2 May 2002 12:01:15 -0500 (CDT) Received: from aragorn (aragorn [129.108.5.35]) by cs.utep.edu (8.11.3/8.11.3) with SMTP id g42H0IB07940; Thu, 2 May 2002 11:00:18 -0600 (MDT) Message-Id: <200205021700.g42H0IB07940 [at] cs [dot] utep.edu> Date: Thu, 2 May 2002 11:00:18 -0600 (MDT) From: Vladik Kreinovich Reply-To: Vladik Kreinovich Subject: Re: Manifolds To: tupper [at] peda [dot] com, Arnold.Neumaier [at] univie [dot] ac.at, Bill.Walster [at] Eng [dot] Sun.COM Cc: reliable_computing [at] interval [dot] louisiana.edu MIME-Version: 1.0 Content-Type: TEXT/plain; charset=us-ascii Content-MD5: 83IQX07C2toyl03OAhTLVw== X-Mailer: dtmail 1.3.0 @(#)CDE Version 1.4 SunOS 5.8 sun4u sparc Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk > Date: Thu, 02 May 2002 09:47:30 -0700 (PDT) > From: William Walster > Subject: Re: Manifolds > To: tupper [at] peda [dot] com, Arnold.Neumaier [at] univie [dot] ac.at > Cc: reliable_computing [at] interval [dot] louisiana.edu > MIME-version: 1.0 > Content-MD5: B2WdSITYG3vRkRflGFU4Yw== > Content-Transfer-Encoding: 8bit > X-MIME-Autoconverted: from QUOTED-PRINTABLE to 8bit by interval.louisiana.edu id g42Gm1h12834 > Dear Bill, Thanks for your thoughtful comments. > Manifold: In general, manifold may mean any collection or set of > objects. ... > My interpretation is that as long as we are clear how we are using the > term, we should not feel inhibited from using it in the way George and > others have done in the past. I think there is no disagreement between what you are saying here and what Jeff was saying. It is very important to spell out what exactly is meant by a manifold. Jeff's point was that is a paper says "manifold" without any further explanation of what is meant, then readers will assume that it is a topological or smooth manifold (because these are the most frequent usages of this word in math), and become confused if this is not what was meant. >From this viewpoint, if the intent is to say a set of solutions, it is better not to use the word manifold without explaining this non-standard meaning because otherwise, readers may erroneously assume that the word is taken in its technical meaning as opposed to commonsense meaning. My understanding is that you are not objecting to making such clarifications, so I really do not see any disagreement. Vladik From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 2 22:41:19 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g433fJW14087 for reliable_computing-outgoing; Thu, 2 May 2002 22:41:19 -0500 (CDT) Received: from espresso.cafe.net (espresso.cafe.net [204.244.119.1]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g433fCh14083 for ; Thu, 2 May 2002 22:41:12 -0500 (CDT) Received: from [192.168.1.2] (CPE00045a0ac3bf.cpe.net.cable.rogers.com [24.156.35.65]) by espresso.cafe.net (8.9.3/8.9.3) with ESMTP id WAA28997 for ; Thu, 2 May 2002 22:51:19 -0700 Mime-Version: 1.0 X-Sender: tupper [at] www [dot] peda.com Message-Id: In-Reply-To: <0GVH00AMATDQXL [at] mpkmail [dot] eng.sun.com> References: <0GVH00AMATDQXL [at] mpkmail [dot] eng.sun.com> Date: Thu, 2 May 2002 23:40:13 -0400 To: reliable_computing [at] interval [dot] louisiana.edu From: Jeff Tupper / Pedagoguery Software Inc Subject: Re: Manifolds Content-Type: text/plain; charset="us-ascii" ; format="flowed" Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk This topic seems to be a bit of a sleeper. At 9:47 AM -0700 5/2/02, William Walster wrote: >My interpretation is that as long as we are clear how we are using the >term, we should not feel inhibited from using it in the way George and >others have done in the past. I'm not trying to make others feel inhibited! At 11:00 AM -0600 5/2/02, Vladik Kreinovich wrote: >Jeff's point was that is a paper says "manifold" without any further >explanation of >what is meant, then readers will assume that it is a topological or smooth >manifold (because these are the most frequent usages of this word in >math), and >become confused if this is not what was meant. I'm not sure if my communiques had much of a point other than to distract readers from the mundane details of day-to-day life. At 11:00 AM -0600 5/2/02, Vladik Kreinovich wrote: >I think there is no disagreement between what you are saying here >and what Jeff >was saying. I tend to agree with this. At 9:47 AM -0700 5/2/02, William Walster wrote: >According to James & James Mathematics Dictionary, Fourth Edition, there >appears to be room for the word "manifold" to be used in a variety of >ways. Out of interest, is there a reason you ascribe such authority to this text? I ask since I've noticed that you've referenced it in at least one of your papers. At 9:47 AM -0700 5/2/02, William Walster wrote: >They say in part: > > Manifold: In general, manifold may mean any collection or set of > objects. I would argue that this shows the age of the dictionary. Manifolds are now, and have been for some time, a standard part of proper upper-level undergraduate university mathematics curricula around the world (courses designed for students going on to further study / work in mathematics) and are not the somewhat new idea they were many decades ago. I think you'd have a hard time finding a mathematician comfortable with using the term in this very general use. At 9:47 AM -0700 5/2/02, William Walster wrote: >They go on to list a number, which I will not copy, but just list: > > topological manifold; > pseudomanifold; > orientable manifold > non-orientable manifold. Solution sets need not fit into any of these types of manifold. I am not aware of any technical definition for any type of manifold that is loose enough to include general solutions sets. Keep in mind that solution sets do not even have to have a uniform dimension: for example, consider the solution set of x(y-|y|)=0. By the way, there are some other terms available as well. One example: the set of solutions to a set of real polynomials over R^n is a real algebraic variety. (And these do admit sets of non-uniform dimension.) Using these terms correctly may help bridge the gap between various disciplines. Misusing them may help keep this web site more lively. At 4:21 AM -0400 4/18/02, Arnold Neumaier wrote: >Unless the singular points are exactly representable and represented >(which is the case only in *very* few problems from applications), >and no rounding error occurs in the verification phase (which excludes >practically all interval techniques and necessitates symbolic methods) >they cannot be verified, since they are unstable under perturbations. >(Look at double roots or bifurcation points...) >One can produce valid enclosures for the region, but no associated >existence statement. It seems that we have a somewhat different view on this. I see more and more symbolic reasoning being folded in to interval methods. (Or I may just become more and more of an outcast of the interval community.) Furthermore, one of the applications that I am interested in is mathematical visualization, where the structure of solutions and their representation can be somewhat different than with scientific applications. Please keep in mind that I'm but a lowly Ph.D. student that doesn't know much about anything! But I'm hoping to become a Ph.D. candidate in the not-too-distant future (i.e. get past my depth exam). No disrespect is meant to anyone. :-) Best wishes, Jeff -- From owner-reliable_computing [at] interval [dot] louisiana.edu Fri May 3 09:53:30 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g43ErUZ15539 for reliable_computing-outgoing; Fri, 3 May 2002 09:53:30 -0500 (CDT) Received: from nwkea-mail-1.sun.com (nwkea-mail-1.sun.com [192.18.42.13]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g43ErMh15535 for ; Fri, 3 May 2002 09:53:23 -0500 (CDT) Received: from engmail1.Eng.Sun.COM ([129.146.1.13]) by nwkea-mail-1.sun.com (8.9.3+Sun/8.9.3) with ESMTP id HAA27507; Fri, 3 May 2002 07:52:44 -0700 (PDT) Received: from phys-mpkmaila (phys-mpkmaila.Eng.Sun.COM [129.146.18.131]) by engmail1.Eng.Sun.COM (8.9.3+Sun/8.9.3/ENSMAIL,v2.1p1) with ESMTP id HAA26919; Fri, 3 May 2002 07:52:44 -0700 (PDT) Received: from gww (gww.Eng.Sun.COM [129.146.78.116]) by mpkmail.eng.sun.com (iPlanet Messaging Server 5.2 Patch 1 (built Apr 2 2002)) with SMTP id <0GVJ0089LIQHNU [at] mpkmail [dot] eng.sun.com>; Fri, 03 May 2002 07:54:17 -0700 (PDT) Date: Fri, 03 May 2002 07:52:44 -0700 (PDT) From: William Walster Subject: Re: Manifolds To: reliable_computing [at] interval [dot] louisiana.edu, tupper [at] peda [dot] com Reply-to: William Walster Message-id: <0GVJ0089MIQHNU [at] mpkmail [dot] eng.sun.com> MIME-version: 1.0 X-Mailer: dtmail 1.3.0 @(#)CDE Version 1.4.2 SunOS 5.8 sun4u sparc Content-type: TEXT/plain; charset=us-ascii Content-transfer-encoding: 7BIT Content-MD5: F/OaBVLaU/aEgIUUpElOLg== Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk > > >At 9:47 AM -0700 5/2/02, William Walster wrote: >>According to James & James Mathematics Dictionary, Fourth Edition, there >>appears to be room for the word "manifold" to be used in a variety of >>ways. Jeff Tupper responded: > >Out of interest, is there a reason you ascribe such authority to this >text? I ask since I've noticed that you've referenced it in at least >one of your papers. > > > Jeff, No particular reason. It is handy. I should probably purchase a more recent edition. I will be surprised if general usage (the way interval researchers have used the word) has changed. My only point was to to reinforce the idea that the word manifold has in the past, at least, been used in a general way. BTW, The Merrian Webster's Collegiate Dictionay, Tenth Edition 1997, lists for the mathematical definition of manifold: A topological space in which every point has a neighborhood that is homemorphic to the interior of a sphere in Euclidean space of the same number of dimensions. My other citation to the James and James dictionary was regarding the definition of an indeterminate form, which I found to be very helpful and provided me with a metaphor to use in thinking about what their containment sets should be. I have always been very impressed by your work, and look forward to you finishing you Ph.D. As I am sure is the case with other members of the interval research community, if I can do anything to help, do not hesitate to contact me. Best regards, Bill G. William (Bill) Walster, Ph.D. Interval Technology Engineering Manager Sun Microsystems, Inc. 16 Network Circle, MS UMPK16-304 Menlo Park, CA 94025 (650) 786-9004 Direct (650) 786-9551 Fax bill.walster [at] eng [dot] sun.com From owner-reliable_computing [at] interval [dot] louisiana.edu Sat May 4 08:51:33 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g44DpXv18103 for reliable_computing-outgoing; Sat, 4 May 2002 08:51:33 -0500 (CDT) Received: from merlin.nerim.net (aboukir-101-1-11-gdr.adsl.nerim.net [62.212.99.186]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g44DpPh18099 for ; Sat, 4 May 2002 08:51:26 -0500 (CDT) Received: (from gdr@localhost) by merlin.nerim.net (8.11.6/8.11.6/SuSE Linux 0.5) id g44Dp8p04794; Sat, 4 May 2002 15:51:08 +0200 To: reliable_computing [at] interval [dot] louisiana.edu Subject: Re: Manifolds References: <0GVH00AMATDQXL [at] mpkmail [dot] eng.sun.com> From: Gabriel Dos Reis In-Reply-To: Jeff Tupper / Pedagoguery Software Inc's message of "Thu, 2 May 2002 23:40:13 -0400" Organization: CodeSourcery, LLC Mime-Version: 1.0 (generated by tm-edit 7.106) Content-Type: text/plain; charset=ISO-8859-1 Content-Transfer-Encoding: 8bit Date: 04 May 2002 15:51:08 +0200 Message-ID: Lines: 30 X-Mailer: Gnus v5.7/Emacs 20.7 Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Jeff Tupper / Pedagoguery Software Inc writes: [...] | At 9:47 AM -0700 5/2/02, William Walster wrote: | >They say in part: | > | > Manifold: In general, manifold may mean any collection or set of | > objects. | | I would argue that this shows the age of the dictionary. Manifolds | are now, and have been for some time, a standard part of proper | upper-level undergraduate university mathematics curricula around the | world (courses designed for students going on to further study / work | in mathematics) and are not the somewhat new idea they were many | decades ago. I think you'd have a hard time finding a mathematician | comfortable with using the term in this very general use. I would like to second Jeff here. I don't think that liberal (if not uncommon) usage of "manifold" is appropriate in such a technical paper. I must confess that as a Differential Geometer, I'm uncomfortable with it :-) Best, -- Gaby Gabriel Dos Reis | Centre de Mathématiques et de dosreis [at] cmla [dot] ens-cachan.fr | Leurs Applications Fax : (33) 01 47 40 59 01 | ENS de Cachan -- CNRS (URA 1611) 61, Avenue du Pdt Wilson, 94235 Cachan - Cedex From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 9 13:07:24 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g49I7Oc01115 for reliable_computing-outgoing; Thu, 9 May 2002 13:07:24 -0500 (CDT) Received: from cs.utep.edu (mail.cs.utep.edu [129.108.5.3]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g49I7Ji01111 for ; Thu, 9 May 2002 13:07:19 -0500 (CDT) Received: from aragorn (aragorn [129.108.5.35]) by cs.utep.edu (8.11.3/8.11.3) with SMTP id g49I7DM24723; Thu, 9 May 2002 12:07:13 -0600 (MDT) Message-Id: <200205091807.g49I7DM24723 [at] cs [dot] utep.edu> Date: Thu, 9 May 2002 12:07:13 -0600 (MDT) From: Vladik Kreinovich Reply-To: Vladik Kreinovich Subject: is there ready software for exact range estimation of quadratic functions? To: reliable_computing [at] interval [dot] louisiana.edu Cc: scott [at] ramas [dot] com, vladik [at] cs [dot] utep.edu MIME-Version: 1.0 Content-Type: TEXT/plain; charset=us-ascii Content-MD5: s8HwKL9cMd0psph5EYQfvQ== X-Mailer: dtmail 1.3.0 @(#)CDE Version 1.4 SunOS 5.8 sun4u sparc Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Dear Friends, It is well known that for a quadratic function f of a reasonably small number of variables (e.g., 5) defined on a box [a1,b1] x ... x [an,bn] with rational endpoints ai and bi, it is possible to compute its exact range. For example, we can use the fact that if the maximum of f on the box is attained at a point x1,...,xn, then for each variable xi, either xi=ai or xi=bi, or xi is inside the interval in which case the partial derivative df/dxi is equal to 0. Since f is quadratic, the equation df/dxi=0 is a linear equation. Thus, for each variable i, one of the three linear equation holds. Overall, we have a quite tractable number of 3^5=243 systems of 5 x 5 linear equations. We can solve all these systems, dismiss solutions which are outside the box, and find the larget value for all solutions that are inside. There are also probably better algorithms as well. The question that we have is as follows: *****is there available software that can compute this exact range - whether using the above algorithm or some better one?********* Such software will help us - and others - in testing different techniques and algorithms for computing the range of quadratic functions. Please reply to both of us, to vladik [at] cs [dot] utep.edu and scott [at] ramas [dot] com Thanks a lot. Scott and Vladik From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 9 14:10:11 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g49JAAK01290 for reliable_computing-outgoing; Thu, 9 May 2002 14:10:10 -0500 (CDT) Received: from mailbox.univie.ac.at (mailbox.univie.ac.at [131.130.1.27]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g49JA5i01286 for ; Thu, 9 May 2002 14:10:05 -0500 (CDT) Received: from univie.ac.at (hektor.mat.univie.ac.at [131.130.16.21]) by mailbox.univie.ac.at (8.12.2/8.12.2) with ESMTP id g49J9dx31249516; Thu, 9 May 2002 21:09:42 +0200 Message-ID: <3CDAC973.A5D150BD [at] univie [dot] ac.at> Date: Thu, 09 May 2002 21:09:39 +0200 From: Arnold Neumaier Organization: University of Vienna X-Mailer: Mozilla 4.78 [en] (X11; U; Linux 2.4.9-31 i686) X-Accept-Language: en, de MIME-Version: 1.0 To: Vladik Kreinovich CC: reliable_computing [at] interval [dot] louisiana.edu, scott [at] ramas [dot] com Subject: Re: is there ready software for exact range estimation of quadratic functions? References: <200205091807.g49I7DM24723 [at] cs [dot] utep.edu> Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Vladik Kreinovich wrote: > > Dear Friends, > > It is well known that for a quadratic function f of a reasonably small number > of variables (e.g., 5) defined on a box [a1,b1] x ... x [an,bn] with rational > endpoints ai and bi, it is possible to compute its exact range. > > For example, we can use the fact that if the maximum of f on the box is > attained at a point x1,...,xn, then for each variable xi, either xi=ai or > xi=bi, or xi is inside the interval in which case the partial derivative df/dxi > is equal to 0. Since f is quadratic, the equation df/dxi=0 is a linear > equation. > > Thus, for each variable i, one of the three linear equation holds. Overall, we > have a quite tractable number of 3^5=243 systems of 5 x 5 linear equations. We > can solve all these systems, dismiss solutions which are outside the box, and > find the larget value for all solutions that are inside. This is the peeling process described in Section 5.2.3 of Baker Kearfott's book, and rediscovered in Makino's thesis. Probably some version of it is implemented in GLOBSOL. Arnold From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 9 14:25:15 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g49JPE401408 for reliable_computing-outgoing; Thu, 9 May 2002 14:25:14 -0500 (CDT) Received: from cs.utep.edu (mail.cs.utep.edu [129.108.5.3]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g49JPAi01404 for ; Thu, 9 May 2002 14:25:10 -0500 (CDT) Received: from aragorn (aragorn [129.108.5.35]) by cs.utep.edu (8.11.3/8.11.3) with SMTP id g49JOup25253; Thu, 9 May 2002 13:24:56 -0600 (MDT) Message-Id: <200205091924.g49JOup25253 [at] cs [dot] utep.edu> Date: Thu, 9 May 2002 13:24:56 -0600 (MDT) From: Vladik Kreinovich Reply-To: Vladik Kreinovich Subject: Re: is there ready software for exact range estimation of quadratic functions? To: vladik [at] cs [dot] utep.edu, Arnold.Neumaier [at] univie [dot] ac.at Cc: reliable_computing [at] interval [dot] louisiana.edu, scott [at] ramas [dot] com MIME-Version: 1.0 Content-Type: TEXT/plain; charset=us-ascii Content-MD5: icY7dl0eSRoPYjpZgWTFCw== X-Mailer: dtmail 1.3.0 @(#)CDE Version 1.4 SunOS 5.8 sun4u sparc Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Dear Arnold, Thanks for the prompt reply. Let us see what Globsol people answer. Vladik We described this algorithm (as a "textbook method") on pp. 9-10 of our complexity book Vladik Kreinovich, Anatoly Lakeyev, Jiri Rohn, and Patrick Kahl, "Computational complexity and feasibility of data processing and interval computations", Kluwer, Dordrecht, 1997. (actually, we did cite Baker Kearfott's 1996 book). I apologize for not citing it in my query, I definitely did not mean to indicate that we invented this idea. > Date: Thu, 09 May 2002 21:09:39 +0200 > From: Arnold Neumaier > X-Accept-Language: en, de > MIME-Version: 1.0 > To: Vladik Kreinovich > CC: reliable_computing [at] interval [dot] louisiana.edu, scott [at] ramas [dot] com > Subject: Re: is there ready software for exact range estimation of quadratic functions? > Content-Transfer-Encoding: 7bit > > Vladik Kreinovich wrote: > > > > Dear Friends, > > > > It is well known that for a quadratic function f of a reasonably small number > > of variables (e.g., 5) defined on a box [a1,b1] x ... x [an,bn] with rational > > endpoints ai and bi, it is possible to compute its exact range. > > > > For example, we can use the fact that if the maximum of f on the box is > > attained at a point x1,...,xn, then for each variable xi, either xi=ai or > > xi=bi, or xi is inside the interval in which case the partial derivative df/dxi > > is equal to 0. Since f is quadratic, the equation df/dxi=0 is a linear > > equation. > > > > Thus, for each variable i, one of the three linear equation holds. Overall, we > > have a quite tractable number of 3^5=243 systems of 5 x 5 linear equations. We > > can solve all these systems, dismiss solutions which are outside the box, and > > find the larget value for all solutions that are inside. > > > This is the peeling process described in Section 5.2.3 of > Baker Kearfott's book, and rediscovered in Makino's thesis. > Probably some version of it is implemented in GLOBSOL. > > Arnold From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 9 15:50:54 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g49KorA01619 for reliable_computing-outgoing; Thu, 9 May 2002 15:50:53 -0500 (CDT) Received: from relay2.EECS.Berkeley.EDU (relay2.EECS.Berkeley.EDU [169.229.60.28]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g49Koni01614 for ; Thu, 9 May 2002 15:50:49 -0500 (CDT) Received: from relay3.EECS.Berkeley.EDU (localhost.Berkeley.EDU [127.0.0.1]) by relay2.EECS.Berkeley.EDU (8.9.3/8.9.3) with ESMTP id NAA02152 for ; Thu, 9 May 2002 13:50:47 -0700 (PDT) Received: from cs.berkeley.edu (dhcp-131-246.CS.Berkeley.EDU [128.32.131.246]) by relay3.EECS.Berkeley.EDU (8.9.3/8.9.3) with ESMTP id NAA16033 for ; Thu, 9 May 2002 13:50:46 -0700 (PDT) Message-ID: <3CDAE127.1E126035 [at] cs [dot] berkeley.edu> Date: Thu, 09 May 2002 13:50:47 -0700 From: James Demmel X-Mailer: Mozilla 4.76 [en] (Windows NT 5.0; U) X-Accept-Language: en MIME-Version: 1.0 To: interval Subject: Accurate Floating Point Summation Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk You may be interested in a paper by Yozo Hida and myself on how to exploit the arithmetic formats in the IEEE 754 floating point standard to get guaranteed accurate sums and dot products in an efficient manner. You can find the paper on my web page www.cs.berkeley.edu/~demmel under Research Projects and Books/BeBOP/Selected Papers and Reports or directly at www.cs.berkeley.edu/~demmel/AccurateSummation.{ps,pdf} Regards, Jim Demmel From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 9 17:53:16 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g49MrFl01916 for reliable_computing-outgoing; Thu, 9 May 2002 17:53:15 -0500 (CDT) Received: from bologna.vision.caltech.edu (bologna.vision.caltech.edu [131.215.134.19]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g49MrAi01912 for ; Thu, 9 May 2002 17:53:10 -0500 (CDT) Received: from peking.vision.caltech.edu (peking [131.215.134.18]) by bologna.vision.caltech.edu (8.9.3/8.8.7) with ESMTP id PAA06572 for ; Thu, 9 May 2002 15:53:09 -0700 Received: (from arrigo@localhost) by peking.vision.caltech.edu (8.9.3+Sun/8.9.1) id PAA05494; Thu, 9 May 2002 15:53:05 -0700 (PDT) X-Authentication-Warning: peking.vision.caltech.edu: arrigo set sender to arrigo [at] vision [dot] caltech.edu using -f To: reliable_computing [at] interval [dot] louisiana.edu Subject: Re: is there ready software for exact range estimation of quadratic functions? References: <200205091807.g49I7DM24723 [at] cs [dot] utep.edu> <3CDAC973.A5D150BD [at] univie [dot] ac.at> From: Arrigo Benedetti Date: 09 May 2002 15:53:05 -0700 In-Reply-To: <3CDAC973.A5D150BD [at] univie [dot] ac.at> Message-ID: Lines: 39 User-Agent: Gnus/5.0808 (Gnus v5.8.8) Emacs/20.4 MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk I am actually working at a polynomial (2nd order) range bounder that is based on this "peeling" process described by Kearfott. The method that I am using for generating the lower dimension boundary elements is a bit different from what Kearfott described in that it is not a recursive algorithm. I should have results very soon, and I am looking forward to post them to the list. Best -Arrigo Arnold Neumaier writes: > Vladik Kreinovich wrote: > > > > Dear Friends, > > > > It is well known that for a quadratic function f of a reasonably small number > > of variables (e.g., 5) defined on a box [a1,b1] x ... x [an,bn] with rational > > endpoints ai and bi, it is possible to compute its exact range. > > > > For example, we can use the fact that if the maximum of f on the box is > > attained at a point x1,...,xn, then for each variable xi, either xi=ai or > > xi=bi, or xi is inside the interval in which case the partial derivative df/dxi > > is equal to 0. Since f is quadratic, the equation df/dxi=0 is a linear > > equation. > > > > Thus, for each variable i, one of the three linear equation holds. Overall, we > > have a quite tractable number of 3^5=243 systems of 5 x 5 linear equations. We > > can solve all these systems, dismiss solutions which are outside the box, and > > find the larget value for all solutions that are inside. > > > This is the peeling process described in Section 5.2.3 of > Baker Kearfott's book, and rediscovered in Makino's thesis. > Probably some version of it is implemented in GLOBSOL. > > Arnold > From owner-reliable_computing [at] interval [dot] louisiana.edu Thu May 9 18:20:09 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g49NK9M02053 for reliable_computing-outgoing; Thu, 9 May 2002 18:20:09 -0500 (CDT) Received: from imf03bis.bellsouth.net (mail203.mail.bellsouth.net [205.152.58.143]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g49NK4i02048 for ; Thu, 9 May 2002 18:20:04 -0500 (CDT) Received: from u8174 ([66.20.81.15]) by imf03bis.bellsouth.net (InterMail vM.5.01.04.05 201-253-122-122-105-20011231) with SMTP id <20020509232124.IUKC28927.imf03bis.bellsouth.net@u8174>; Thu, 9 May 2002 19:21:24 -0400 Message-Id: <2.2.32.20020509231959.00befd10 [at] pop [dot] louisiana.edu> X-Sender: rbk5287 [at] pop [dot] louisiana.edu X-Mailer: Windows Eudora Pro Version 2.2 (32) Mime-Version: 1.0 Content-Type: text/plain; charset="us-ascii" Date: Thu, 09 May 2002 18:19:59 -0500 To: Arnold Neumaier , Vladik Kreinovich From: "R. Baker Kearfott" Subject: Re: is there ready software for exact range estimation of quadratic functions? Cc: reliable_computing [at] interval [dot] louisiana.edu, scott [at] ramas [dot] com Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Arnold, Yes, the peeling process is implemented in GlobSol. GlobSol finds the min of a function. We can call GlobSol twice to find the range of a quadratic over a box, provided the box is sufficiently low-dimensional. Best regards, Baker At 09:09 PM 5/9/02 +0200, Arnold Neumaier wrote: >Vladik Kreinovich wrote: >> >> Dear Friends, >> >> It is well known that for a quadratic function f of a reasonably small number >> of variables (e.g., 5) defined on a box [a1,b1] x ... x [an,bn] with rational >> endpoints ai and bi, it is possible to compute its exact range. >> >> For example, we can use the fact that if the maximum of f on the box is >> attained at a point x1,...,xn, then for each variable xi, either xi=ai or >> xi=bi, or xi is inside the interval in which case the partial derivative df/dxi >> is equal to 0. Since f is quadratic, the equation df/dxi=0 is a linear >> equation. >> >> Thus, for each variable i, one of the three linear equation holds. Overall, we >> have a quite tractable number of 3^5=243 systems of 5 x 5 linear equations. We >> can solve all these systems, dismiss solutions which are outside the box, and >> find the larget value for all solutions that are inside. > > >This is the peeling process described in Section 5.2.3 of >Baker Kearfott's book, and rediscovered in Makino's thesis. >Probably some version of it is implemented in GLOBSOL. > >Arnold > > --------------------------------------------------------------- R. Baker Kearfott, rbk [at] louisiana [dot] edu (337) 482-5346 (fax) (337) 482-5270 (work) (337) 981-9744 (home) URL: http://interval.louisiana.edu/kearfott.html Department of Mathematics, University of Louisiana at Lafayette Box 4-1010, Lafayette, LA 70504-1010, USA --------------------------------------------------------------- From owner-reliable_computing [at] interval [dot] louisiana.edu Fri May 10 10:45:17 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4AFjGL04035 for reliable_computing-outgoing; Fri, 10 May 2002 10:45:16 -0500 (CDT) Received: from cs.utep.edu (mail.cs.utep.edu [129.108.5.3]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4AFj9i04031 for ; Fri, 10 May 2002 10:45:10 -0500 (CDT) Received: from aragorn (aragorn [129.108.5.35]) by cs.utep.edu (8.11.3/8.11.3) with SMTP id g4AFj4s00547 for ; Fri, 10 May 2002 09:45:04 -0600 (MDT) Message-Id: <200205101545.g4AFj4s00547 [at] cs [dot] utep.edu> Date: Fri, 10 May 2002 09:45:02 -0600 (MDT) From: Vladik Kreinovich Reply-To: Vladik Kreinovich Subject: Issac'2002 Call For Posters To: reliable_computing [at] interval [dot] louisiana.edu MIME-Version: 1.0 Content-Type: TEXT/plain; charset=us-ascii Content-MD5: lLf0REowCcV5mdWssEvKVA== X-Mailer: dtmail 1.3.0 @(#)CDE Version 1.4 SunOS 5.8 sun4u sparc Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk deadline for posters is June 29. Issac stands for International Symposium on Symbolic and Algebraic Computation, they would welcome interval-related papers. I am adding this link to the Forthcoming Conferences part of the interval computations webpage http://www.cs.utep.edu/interval-comp ------------- Begin Forwarded Message ------------- From: Conference.Issac2002 [at] lifl [dot] fr Date: Fri, 10 May 2002 02:44:39 +0200 (MET DST) TO: vladik [at] cs [dot] utep.edu Subject: Issac'2002 Call For Posters Mime-Version: 1.0 The Issac'02 Call for Posters has started. You will find information on the Issac'02 website at the following address: http://www.lifl.fr/issac2002/callposters.html Best regards, Fran\c{c}ois Lemaire. ------------- End Forwarded Message ------------- From owner-reliable_computing [at] interval [dot] louisiana.edu Sat May 11 11:34:36 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4BGYaS06871 for reliable_computing-outgoing; Sat, 11 May 2002 11:34:36 -0500 (CDT) Received: from lcyoung.math.wisc.edu (lcyoung.math.wisc.edu [144.92.166.90]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4BGYUi06866 for ; Sat, 11 May 2002 11:34:30 -0500 (CDT) Received: from ultra10.math.wisc.edu (ultra10.math.wisc.edu [144.92.166.180]) by lcyoung.math.wisc.edu (8.11.6/8.11.6) with ESMTP id g4BGYEM08762; Sat, 11 May 2002 11:34:14 -0500 (CDT) Received: from localhost (hans@localhost) by ultra10.math.wisc.edu (8.11.6+Sun/8.11.6) with ESMTP id g4BGYBE10504; Sat, 11 May 2002 11:34:11 -0500 (CDT) Date: Sat, 11 May 2002 11:34:11 -0500 (CDT) From: Hans Schneider To: NETS -- at-net , "Hershkowitz, Danny -- Hershkowitz Daniel" , Danny Hershkowitz , E-LETTER , "na.digest" , ipnet-digest [at] math [dot] msu.edu, wim@bell-labs.com, hjt [at] eos [dot] ncsu.edu, vkm [at] eedsp [dot] gatech.edu, reliable_computing [at] interval [dot] louisiana.edu Subject: LAA contents Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Dear Net Organizer: Please circulate the attached LAA contents over your net. Thanks hans ------------------------------------------------------------------------------ Hans Schneider hans [at] math [dot] wisc.edu. Department of Mathematics 608-262-1402 (Work) Van Vleck Hall 608-271-7252 (Home) 480 Lincoln Drive 608-263-8891 (Work FAX) University of Wisconsin-Madison No Home FAX at present Madison WI 53706 USA http://www.math.wisc.edu/~hans (URL) ------------------------------------------------------------------------------ Journal: Linear Algebra and its Applications ISSN : 0024-3795 Volume : 350 Issue : 1-3 Date : 15-Jul-2002 Visit the journal at http://www.elsevier.nl/locate/jnlnr/07738 The full text of articles is available via ScienceDirect : http://www.sciencedirect.com/ The numerical range of a nonnegative matrix C.-K. Li, B.-S. Tam, P. Yuan Wu pp 1-23 A weighted Drazin inverse and applications V. Rakocevic, Y. Wei pp 25-39 Commutativity preserving linear maps and Lie automorphisms of strictly triangular matrix space Y. Cao, Z. Chen, C. Huang pp 41-66 State space calculations for discrete probability densities B. Hanzon, R.J. Ober pp 67-87 A priori error bounds on invariant subspace approximations by block Krylov subspaces M. Robbe, M. Sadkane pp 89-103 Covering numbers under small perturbations R. Fernandes pp 105-124 Containment regions for zeros of polynomials from numerical ranges of companion matrices H. Linden pp 125-145 The polytope of degree sequences of hypergraphs N.L. Bhanu Murthy, M.K. Srinivasan pp 147-170 Some explicit formulas for the polynomial decomposition of the matrix exponential and applications R. Ben Taher, M. Rachidi pp 171-184 Intersections of nest algebras in finite dimensions P.A. Fillmore, W.E. Longstaff, G.W.M. H. Radjavi, Y. Zhong pp 185-197 On the reconstruction of Toeplitz matrix inverses from columns G. Heinig pp 199-212 The Lie algebra structure and controllability of spin systems F. Albertini, D. D'Alessandro pp 213-235 Condition numbers of algebraic Riccati equations in the Frobenius norm J.-g. Sun pp 237-261 Full text via ScienceDirect : http://www.sciencedirect.com/science?_ob=GatewayURL&_origin=CONTENTS&_method=citationSearch&_piikey=S002437950200294X&_version=1&md5=b26f613add034e5972d5a4c203f57324 Linear operators that preserve pairs of matrices which satisfy extreme rank properties L.B. Beasley, S.-G. Lee, S.-Z. Song pp 263-272 A new proof of Mayer's theorem C.-T. Pang, Y.-Y. Lur, S.-M. Guu pp 273-278 Total positivity and Toda flow G.M.L. Gladwell pp 279-284 A note on Krylov subspace methods for singular systems Z.-H. Cao, M. Wang pp 285-288 Convexoid and generalized derivations M. Barraa pp 289-292 Author index pp 293 From owner-reliable_computing [at] interval [dot] louisiana.edu Mon May 13 02:40:42 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4D7efw11195 for reliable_computing-outgoing; Mon, 13 May 2002 02:40:41 -0500 (CDT) Received: from sun1.udg.es (sun1.udg.es [130.206.45.89]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4D7eYi11191 for ; Mon, 13 May 2002 02:40:35 -0500 (CDT) Received: from eia.udg.es (silver.udg.es [130.206.129.16]) by sun1.udg.es (8.9.3/8.9.3/otb) with ESMTP id JAA13093 for ; Mon, 13 May 2002 09:39:24 +0100 (WET DST) Received: from PCEIA000 (pceia000.udg.es [130.206.129.66]) by eia.udg.es (8.8.5/8.8.5) with ESMTP id JAA22424 for ; Mon, 13 May 2002 09:41:46 +0100 (WET DST) From: "Josep Vehi" To: Subject: Promoting interval computations Date: Mon, 13 May 2002 09:40:30 +0200 Organization: Universitat de Girona Message-ID: <000001c1fa51$752cf4a0$4281ce82@PCEIA000> MIME-Version: 1.0 Content-Type: multipart/alternative; boundary="----=_NextPart_000_0001_01C1FA62.38B5C4A0" X-Priority: 3 (Normal) X-MSMail-Priority: Normal X-Mailer: Microsoft Outlook, Build 10.0.2627 Importance: Normal X-MimeOLE: Produced By Microsoft MimeOLE V6.00.2600.0000 Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk This is a multi-part message in MIME format. ------=_NextPart_000_0001_01C1FA62.38B5C4A0 Content-Type: text/plain; charset="us-ascii" Content-Transfer-Encoding: 7bit Dear all, The technical committee on Smart Adaptive Systems of the European network of excellence EUNITE (http://www.eunite.org ) is preparing a bibliography on methods ans techniques for smart adaptive systems. It will not be a huge database but just a collection of sellected (key) papers for each technique. As intervals is a way to deal with uncertainty and to manage such uncertainty is crucial in chanching environements and I'm sure that the are many applications of intervals to adaptive systems, I proposed to include a section on intervals in the bibliography and it was accepted. Please, would you be so kind to send me references on applications of intervals with the stress on adaptivity? You will find bellow the details of the bibliographic database organisation. Josep -------------- Scope: Soft computing techniques that can be seen as adaptive. Adaptation Levels: 1. Changing environment. 2. Similar Setting without explicit porting (the system should adapt itself, for a similar application). 3. Adaptation to a new/unknown application. Other definitions: from the task force on taxonomy on the EUNITE web site. Learning levels (Dourado): 1. Parameter learning 2. Structure learning 3. Trajectory learning 4. Task learning 5. Goal learning 6. Learning organisations Methods: * Fuzzy * Neural * Evolutionary * Hybrid * Multi-modular * Multi-agent * Universal learning * Symbolic Machine learning * Probabilistic * Rough sets * Heuristic methods * Chaotic systems * Interval computations * . ------=_NextPart_000_0001_01C1FA62.38B5C4A0 Content-Type: text/html; charset="us-ascii" Content-Transfer-Encoding: quoted-printable

Dear = all,

 

The technical = committee on Smart Adaptive Systems of the European network of excellence EUNITE (http://www.eunite.org) is preparing = a bibliography on methods ans techniques for smart adaptive systems. It will not be a = huge database but just a collection of sellected (key) papers for each = technique.

 

As intervals is a way = to deal with uncertainty and to manage such uncertainty is crucial in chanching environements and I’m sure that the are many applications of = intervals to adaptive systems, I proposed to include a section on intervals in the bibliography and it was accepted.

 

Please, would you be = so kind to send me references on applications of intervals with the stress on adaptivity?

 

You will find bellow = the details of the bibliographic database organisation.

 

Josep

 

--------------

 

Scope= : Soft computing = techniques that can be seen as adaptive.

 

Adaptation = Levels:

  1. Changing = environment.
  2. Similar Setting without explicit porting = (the system should adapt itself, for a similar = application).
  3. Adaptation to a new/unknown = application.

 

 

Other definitions: = from the task force on taxonomy on the EUNITE web = site.

 

Learn= ing levels (Dourado):

  1. Parameter = learning
  2. Structure = learning
  3. Trajectory = learning
  4. Task = learning
  5. Goal = learning
  6. Learning = organisations

 

Metho= ds:

  • Fuzzy
  • Neural
  • Evolutionary
  • Hybrid
  • Multi-modular
  • Multi-agent
  • Universal = learning
  • Symbolic Machine = learning
  • Probabilistic
  • Rough sets
  • Heuristic = methods
  • Chaotic = systems
  • Interval = computations

 

------=_NextPart_000_0001_01C1FA62.38B5C4A0-- From owner-reliable_computing [at] interval [dot] louisiana.edu Tue May 14 06:06:33 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4EB6W714487 for reliable_computing-outgoing; Tue, 14 May 2002 06:06:32 -0500 (CDT) Received: from bxgs.fj.cn ([211.97.147.91]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4EB6Pi14483 for ; Tue, 14 May 2002 06:06:26 -0500 (CDT) Received: from ceti.com.pl [62.194.194.5] by bxgs.fj.cn with ESMTP (SMTPD32-7.05) id AE8D110150; Tue, 14 May 2002 19:01:33 +0800 Message-ID: <00003ed02ffd$00002c7f$00003202 [at] chello [dot] se> To: From: "Customer Care Center" Subject: Conference calls/best quality/$.18 per minute! Date: Tue, 14 May 2002 04:03:13 -1900 MIME-Version: 1.0 Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable X-Mailer: Microsoft Outlook Express Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk New Web Technology

UNLIMITED WEB CON= FERENCING
Subscribe to the Web Conference C= enter for only $40.00 per month!
(Connects up to 15 participants at a time p= lus audio charges)

<= /TABLE>

  • Manage your meetings virtually on line
  • Application sharing
  • Multi-platform compatible (no software needed)
  • Call anytime, from anywhere, to anywhere
  • Unlimited usage for up to 15 participants (larger groups ava= ilabale)
  • Lowest rate $.18 cents per minunte for audio (toll charge= s included)
  • Quality, easy to use service
  • Numerous interactive features
  • FREE DEMO

    Eliminate or Reduce Travel Expense
    Try this on for savings... and turn those unnecessary trips into problem solving and fact finding con= ferences, saving time and wear and tear on your overworked business staff.

    To find out more abo= ut this revolutionary concept, fill out the form below.

    Required Input Field*

    Name*
    Web Address
    Company Name*
    State*
    Business Phone*
    Home Phone
    Email Address*



    From owner-reliable_computing [at] interval [dot] louisiana.edu Wed May 15 09:42:10 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4FEgAP17572 for reliable_computing-outgoing; Wed, 15 May 2002 09:42:10 -0500 (CDT) Received: from clmboh1-smtp3.columbus.rr.com (clmboh1-smtp3.columbus.rr.com [65.24.0.112]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4FEg3i17568 for ; Wed, 15 May 2002 09:42:04 -0500 (CDT) Received: from oemcomputer (dhcp065-024-174-102.columbus.rr.com [65.24.174.102]) by clmboh1-smtp3.columbus.rr.com (8.11.2/8.11.2) with SMTP id g4FBx5608417 for ; Wed, 15 May 2002 07:59:06 -0400 (EDT) Message-ID: <002201c1fc07$d0891560$66ae1841 [at] columbus [dot] rr.com> From: "Ramon Moore" To: "interval" Subject: Fw: Difficulty in your paper Date: Wed, 15 May 2002 07:58:21 -0400 MIME-Version: 1.0 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: 7bit X-Priority: 3 X-MSMail-Priority: Normal X-Mailer: Microsoft Outlook Express 5.50.4133.2400 X-MimeOLE: Produced By Microsoft MimeOLE V5.50.4133.2400 Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Dear Colleagues, Is there a kind soul among you who will help the gentleman who sent me the plea below? I suspect he has a bug in his program, but I do not have the energy to help him find it. Ray (Ramon Moore) ----- Original Message ----- From: "N. S. Kubal _00423002_" To: "Ramon Moore" Sent: Wednesday, May 15, 2002 3:44 AM Subject: Re: Difficulty in your paper > On Wed, 15 May 2002, Ramon Moore wrote: > > > Unfortunately I am unable to open the document you attached. Perhaps you can > > just explain your question in text in an email. > > > > R. Moore > > Dear Sir, > Sorry for your inconvienence. Please go through the following > matter. > > ------------------------------------------------------------------------ > Recently, I have implemented an optimization algorithm in FORTRAN 95 to > find global minimum based on your paper "On Computing the Range of Values" > published in Computing, 28:225-237, 1982 with Asaithambi, N.S. and Shen > Zuhe. > When I ran the algorithm on one of the test function in [1] it > did not give me the correct global minimum. > > Test function in [1] > > f(x,y) = (x^2 + y^2)/200 - cos(x)*cos(y/sqrt(2)) + 1 > X=[-100,100], Y=[-100,100] > > Results obtained using the algorithm given in your paper > Global minimum = 0.145618691261803 > Global minimum box = [-3.1104241497814656,-3.1104241032153367] > [-4.3556719087064267,-4.3556718155741691] > > But the results given in [1] are > Global minimum = 0 > Global minimum point = (0,0) > > When I went through each steps of the algorithm, I had come to > know that the algorithm deletes the boxes, which contain the global > minimum, which is "zero", due to steps 8 and step 13. Now in these steps > the variable 'S' decides to keep/delete the reduced boxes (Z). The values > of 'S' are based on the formulae given in your paper followed by equation > (3.2). > I used new variable 'S1' in steps 8 and 13, which is based on the > following conditions, > > S1 .neq. 0 if (Infimum (DiF(X)) <= 0 and Supremum (DiF(X)) >= 0 in any one > direction) > S1 .eq. 0 if (Infimum (DiF(X)) > 0 or Supremum (DiF(X)) < 0 in all the > directions) > > (i.e. delete the box if the function is strictly monotonic in all the > directions). > > The modified steps 8 and 13 are, > Step 8: if S1(1) is 0 and Z(1) does not intersect the boundary of Z(0), > discard Z(1) and go to step 12. > Step 13: if S1(2) is 0 and Z(2) does not intersect the boundary of Z(0), > discard Z(2) and go to step 17. > > Results obtained after modification > Global minimum = 0.0E0 > Global minimum box = [0.0E0, 0.0E0] [0.0E0, 0.0E0] > > I would highly appreciate if you kindly spend some of your > precious time to go through the algorithm and tell me whether the changes > I have made are correct or not. > > References: > [1] C. Jansson, O. Knuppel. "A Global Minimization Method: The > Multi-Dimensional Case." TUHH (Technische Universitat Hamburg-Harburg), > Bericht 92.1, Januar, 1992, pp.47-48. > > > > From owner-reliable_computing [at] interval [dot] louisiana.edu Wed May 15 12:09:04 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4FH93417964 for reliable_computing-outgoing; Wed, 15 May 2002 12:09:04 -0500 (CDT) Received: from mx1.su.se (mx1.su.se [130.237.162.73]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4FH8wi17960 for ; Wed, 15 May 2002 12:08:59 -0500 (CDT) Received: from joeb.matematik.su.se (joeb.matematik.su.se [130.237.198.14]) by mx1.su.se (8.11.6/8.11.6) with ESMTP id g4FH8nE00348; Wed, 15 May 2002 19:08:49 +0200 Received: (from joeb@localhost) by joeb.matematik.su.se (8.9.3/8.9.1) id TAA26653; Wed, 15 May 2002 19:08:49 +0200 Date: Wed, 15 May 2002 19:08:49 +0200 From: Message-Id: <200205151708.TAA26653 [at] joeb [dot] matematik.su.se> To: reliable_computing [at] interval [dot] louisiana.edu, rmoore17 [at] columbus [dot] rr.com Subject: Re: Fw: Difficulty in your paper Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk I've just read through your letter, including the 'plea'; and I think I interpret it a bit differently. As far as I can see, Mr. Kubal tries to inform you (in a polite and somewhat involved manner) that he thinks he found an error in a paper published by you, and also a remedy for the error. He writes that when he implemented your algorithm right off he got an error, but that when he modified your algorithm in a couple of steps, and implemented the modified algorithm, then he got the correct result. At least, this is what I read into the text. Joergen Backelin From owner-reliable_computing [at] interval [dot] louisiana.edu Wed May 15 13:50:29 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4FIoT618223 for reliable_computing-outgoing; Wed, 15 May 2002 13:50:29 -0500 (CDT) Received: from marnier.ucs.louisiana.edu (root [at] marnier [dot] ucs.louisiana.edu [130.70.132.233]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4FIoNi18219 for ; Wed, 15 May 2002 13:50:24 -0500 (CDT) Received: from Inspiron-8200 (h143215.louisiana.edu [130.70.143.215]) by marnier.ucs.louisiana.edu (8.11.3/8.11.3/ull-ucs-mx-host_1.6) with SMTP id g4FIoLw20845; Wed, 15 May 2002 13:50:21 -0500 (CDT) Message-Id: <2.2.32.20020515185021.00c6b7ac [at] pop [dot] louisiana.edu> X-Sender: rbk5287 [at] pop [dot] louisiana.edu X-Mailer: Windows Eudora Pro Version 2.2 (32) Mime-Version: 1.0 Content-Type: multipart/mixed; boundary="=====================_1021506621==_" Date: Wed, 15 May 2002 13:50:21 -0500 To: "Ramon Moore" , "interval" From: "R. Baker Kearfott" Subject: Re: Fw: Difficulty in your paper X-Attachments: C:\GlobSol\examples\Moore_2002_05_15\Moore.OT1; C:\GlobSol\examples\Moore_2002_05_15\Moore.DT1; Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk --=====================_1021506621==_ Content-Type: text/plain; charset="us-ascii" Friends, For what it's worth, I tried this problem in GlobSol. It found the optimum (0, at (0,0)) without bisection; the output is attached :-) I note that the gradient expressions are SUE's. Best regards, Baker At 07:58 AM 5/15/2002 -0400, Ramon Moore wrote: >Dear Colleagues, > >Is there a kind soul among you who will help the gentleman who sent me the >plea below? >I suspect he has a bug in his program, but I do not have the energy to help >him find it. > >Ray (Ramon Moore) > >----- Original Message ----- >From: "N. S. Kubal _00423002_" >To: "Ramon Moore" >Sent: Wednesday, May 15, 2002 3:44 AM >Subject: Re: Difficulty in your paper > > >> On Wed, 15 May 2002, Ramon Moore wrote: >> >> > Unfortunately I am unable to open the document you attached. Perhaps you >can >> > just explain your question in text in an email. >> > >> > R. Moore >> >> Dear Sir, >> Sorry for your inconvienence. Please go through the following >> matter. >> >> ------------------------------------------------------------------------ >> Recently, I have implemented an optimization algorithm in FORTRAN 95 to >> find global minimum based on your paper "On Computing the Range of Values" >> published in Computing, 28:225-237, 1982 with Asaithambi, N.S. and Shen >> Zuhe. >> When I ran the algorithm on one of the test function in [1] it >> did not give me the correct global minimum. >> >> Test function in [1] >> >> f(x,y) = (x^2 + y^2)/200 - cos(x)*cos(y/sqrt(2)) + 1 >> X=[-100,100], Y=[-100,100] >> >> Results obtained using the algorithm given in your paper >> Global minimum = 0.145618691261803 >> Global minimum box = [-3.1104241497814656,-3.1104241032153367] >> [-4.3556719087064267,-4.3556718155741691] >> >> But the results given in [1] are >> Global minimum = 0 >> Global minimum point = (0,0) >> >> When I went through each steps of the algorithm, I had come to >> know that the algorithm deletes the boxes, which contain the global >> minimum, which is "zero", due to steps 8 and step 13. Now in these steps >> the variable 'S' decides to keep/delete the reduced boxes (Z). The values >> of 'S' are based on the formulae given in your paper followed by equation >> (3.2). >> I used new variable 'S1' in steps 8 and 13, which is based on the >> following conditions, >> >> S1 .neq. 0 if (Infimum (DiF(X)) <= 0 and Supremum (DiF(X)) >= 0 in any one >> direction) >> S1 .eq. 0 if (Infimum (DiF(X)) > 0 or Supremum (DiF(X)) < 0 in all the >> directions) >> >> (i.e. delete the box if the function is strictly monotonic in all the >> directions). >> >> The modified steps 8 and 13 are, >> Step 8: if S1(1) is 0 and Z(1) does not intersect the boundary of Z(0), >> discard Z(1) and go to step 12. >> Step 13: if S1(2) is 0 and Z(2) does not intersect the boundary of Z(0), >> discard Z(2) and go to step 17. >> >> Results obtained after modification >> Global minimum = 0.0E0 >> Global minimum box = [0.0E0, 0.0E0] [0.0E0, 0.0E0] >> >> I would highly appreciate if you kindly spend some of your >> precious time to go through the algorithm and tell me whether the changes >> I have made are correct or not. >> >> References: >> [1] C. Jansson, O. Knuppel. "A Global Minimization Method: The >> Multi-Dimensional Case." TUHH (Technische Universitat Hamburg-Harburg), >> Bericht 92.1, Januar, 1992, pp.47-48. >> >> >> >> > > --=====================_1021506621==_ Content-Type: text/plain; charset="us-ascii" Content-Disposition: attachment; filename="Moore.OT1" Output from FIND_GLOBAL_MIN on 05/15/2002 at 13:43:08. Version for the system is: October 10, 2000 Codelist file name is: MooreG.CDL Box data file name is: Moore.DT1 Initial box: [ -0.1000D+03, 0.1000D+03 ] [ -0.1000D+03, 0.1000D+03 ] BOUND_CONSTRAINT: F F F F --------------------------------------- CONFIGURATION VALUES: EPS_DOMAIN: 0.1000D-07 MAXITR: 20000 MAX_CPU_SECONDS: 0.1080D+05 DO_INTERVAL_NEWTON: T QUADRATIC: T FULL_SPACE: F VERY_GOOD_INITIAL_GUESS: F USE_SUBSIT: F OUTPUT UNIT: 7 PRINT_LENGTH: 1 PHI_MUST_CONVERGE: T EQ_CNS_MUST_CONVERGE: T INEQ_CNS_MUST_CONVERGE: T PHI_THICKNESS_FACTOR: 0.5000D+00 EQ_CNS_THICKNESS_FACTOR: 0.5000D+00 INEQ_CNS_THICKNESS_FACTOR: 0.5000D+00 PHI_MUST_CONVERGE: T EQ_CNS_MUST_CONVERGE: T INEQ_CNS_MUST_CONVERGE: T PHI_CONVERGENCE_FACTOR: 0.1000D-13 EQ_CNS_CONVERGENCE_FACTOR: 0.1000D-13 INEQ_CNS_CONVERGENCE_FACTOR: 0.1000D-13 CONTINUITY_ACROSS_BRANCHES: F SINGULAR_EXPANSION_FACTOR: 0.1000D+02 HEURISTIC PARAMETER ALPHA: 0.5000D+00 APPROX_OPT_BEFORE_BISECTION: T USE_TAYLOR_EQUALITY_CONSTRAINTS F USE_TAYLOR_INEQ_CONSTRAINTS F USE_TAYLOR_OBJECTIVE F USE_TAYLOR_EQ_CNS_GRD F USE_TAYLOR_GRAD F USE_TAYLOR_INEQ_CNS_GRD F USE_TAYLOR_REDUCED_INEWTON F COSY_POLYNOMIAL_ORDER 5 LEAST_SQUARES_FUNCTIONS: F NONLINEAR_SYSTEM: F Default point optimizer was used. THERE WERE NO BOXES IN THE LIST OF SMALL BOXES. LIST OF BOXES CONTAINING VERIFIED FEASIBLE POINTS: Box no.: 1 Box coordinates: [ -0.4950D-09, 0.4950D-09 ] [ -0.1095-306, 0.1095-306 ] PHI: [ -0.1110D-14, 0.2221D-14 ] B%LIUI(1,*): F F B%LIUI(2,*): F F B%SIDE(*): F F B%PEEL(*): F F Level: 0 Box contains the following approximate root: 0.0000D+00 0.0000D+00 OBJECTIVE ENCLOSURE AT APPROXIMATE ROOT: [ -0.1332D-14, 0.2220D-14 ] Unknown = T Contains_root = T Changed coordinates: F F ------------------------------------------------- ALGORITHM COMPLETED WITH LESS THAN THE MAXIMUM NUMBER, 20000 OF BOXES. Number of evaluations of a row of the point expanded Jacobi matrix: 27 No. dense interval residual evaluations -- gradient code list: 14 Number of orig. system C-LP preconditioner rows: 2 Total number of forward_substitutions: 48 Number of Gauss--Seidel steps on the dense system: 2 Number point dense residual evaluations, gradient codelist: 2 No. point dense Hessian matrix evaluations -- gradient codelist: 1 Total number of iterations of one of the point Newton methods: 1 Total number of dense slope matrix evaluations: 8 Total number second-order interval evaluations of the original function: 6 Total number dense interval reduced gradient evaluations: 3 Number of times the interval Newton method made a coordinate interval smaller: 2 Number of times the approximate solver was called: 2 Number Fritz-John matrix evaluations: 1 BEST_ESTIMATE: 0.2220D-14 Overall CPU time: 0.1001D-01 CPU time in PEEL_BOUNDARY: 0.0000D+00 CPU time in REDUCED_INTERVAL_NEWTON: 0.0000D+00 --=====================_1021506621==_ Content-Type: text/plain; charset="us-ascii" Content-Disposition: attachment; filename="Moore.DT1" 1D-8 -100 100 -100 100 F F F F --=====================_1021506621==_ Content-Type: text/plain; charset="us-ascii" --------------------------------------------------------------- R. Baker Kearfott, rbk [at] louisiana [dot] edu (337) 482-5346 (fax) (337) 482-5270 (work) (337) 981-9744 (home) URL: http://interval.louisiana.edu/kearfott.html Department of Mathematics, University of Louisiana at Lafayette Box 4-1010, Lafayette, LA 70504-1010, USA --------------------------------------------------------------- --=====================_1021506621==_-- From owner-reliable_computing [at] interval [dot] louisiana.edu Wed May 15 14:00:54 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4FJ0rB18353 for reliable_computing-outgoing; Wed, 15 May 2002 14:00:53 -0500 (CDT) Received: from clmboh1-smtp3.columbus.rr.com (clmboh1-smtp3.columbus.rr.com [65.24.0.112]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4FJ0ni18349 for ; Wed, 15 May 2002 14:00:49 -0500 (CDT) Received: from oemcomputer (dhcp065-024-174-102.columbus.rr.com [65.24.174.102]) by clmboh1-smtp3.columbus.rr.com (8.11.2/8.11.2) with SMTP id g4FJ0bv18448; Wed, 15 May 2002 15:00:42 -0400 (EDT) Message-ID: <002d01c1fc42$b6513160$66ae1841 [at] columbus [dot] rr.com> From: "Ramon Moore" To: "N. S. Kubal _00423002_" Cc: "interval" Subject: Fw: Fw: Difficulty in your paper Date: Wed, 15 May 2002 14:59:53 -0400 MIME-Version: 1.0 Content-Type: multipart/mixed; boundary="----=_NextPart_000_002A_01C1FC21.2BA5EFE0" X-Priority: 3 X-MSMail-Priority: Normal X-Mailer: Microsoft Outlook Express 5.50.4133.2400 X-MimeOLE: Produced By Microsoft MimeOLE V5.50.4133.2400 Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk This is a multi-part message in MIME format. ------=_NextPart_000_002A_01C1FC21.2BA5EFE0 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: 7bit Dear N. S. Kubal, I am forwarding to you a response to your question from Prof. Kearfott, whom I thank for his help. As you can see, the algorithm suggested in the paper of 1982 has been supplanted by others by now. In any case, I am sorry to say that I have no energy at this time for digging back into the details of that 20 year-old paper. Perhaps your suggested changes are necessary. I do not know. I no longer have a working program implementing the algorithm of the 1982 paper, and I do not wish to write another one at this point in time. I am hoping that others may help with a clear answer, if you are really determined to know. Best wishes, Ramon Moore ----- Original Message ----- From: "R. Baker Kearfott" To: "Ramon Moore" ; "interval" Sent: Wednesday, May 15, 2002 2:50 PM Subject: Re: Fw: Difficulty in your paper > Friends, > > For what it's worth, I tried this problem in GlobSol. It found the > optimum (0, at (0,0)) without bisection; the output is attached :-) > > I note that the gradient expressions are SUE's. > > Best regards, > > Baker > > At 07:58 AM 5/15/2002 -0400, Ramon Moore wrote: > >Dear Colleagues, > > > >Is there a kind soul among you who will help the gentleman who sent me the > >plea below? > >I suspect he has a bug in his program, but I do not have the energy to help > >him find it. > > > >Ray (Ramon Moore) > > > >----- Original Message ----- > >From: "N. S. Kubal _00423002_" > >To: "Ramon Moore" > >Sent: Wednesday, May 15, 2002 3:44 AM > >Subject: Re: Difficulty in your paper > > > > > >> On Wed, 15 May 2002, Ramon Moore wrote: > >> > >> > Unfortunately I am unable to open the document you attached. Perhaps you > >can > >> > just explain your question in text in an email. > >> > > >> > R. Moore > >> > >> Dear Sir, > >> Sorry for your inconvienence. Please go through the following > >> matter. > >> > >> ------------------------------------------------------------------------ > >> Recently, I have implemented an optimization algorithm in FORTRAN 95 to > >> find global minimum based on your paper "On Computing the Range of Values" > >> published in Computing, 28:225-237, 1982 with Asaithambi, N.S. and Shen > >> Zuhe. > >> When I ran the algorithm on one of the test function in [1] it > >> did not give me the correct global minimum. > >> > >> Test function in [1] > >> > >> f(x,y) = (x^2 + y^2)/200 - cos(x)*cos(y/sqrt(2)) + 1 > >> X=[-100,100], Y=[-100,100] > >> > >> Results obtained using the algorithm given in your paper > >> Global minimum = 0.145618691261803 > >> Global minimum box = [-3.1104241497814656,-3.1104241032153367] > >> [-4.3556719087064267,-4.3556718155741691] > >> > >> But the results given in [1] are > >> Global minimum = 0 > >> Global minimum point = (0,0) > >> > >> When I went through each steps of the algorithm, I had come to > >> know that the algorithm deletes the boxes, which contain the global > >> minimum, which is "zero", due to steps 8 and step 13. Now in these steps > >> the variable 'S' decides to keep/delete the reduced boxes (Z). The values > >> of 'S' are based on the formulae given in your paper followed by equation > >> (3.2). > >> I used new variable 'S1' in steps 8 and 13, which is based on the > >> following conditions, > >> > >> S1 .neq. 0 if (Infimum (DiF(X)) <= 0 and Supremum (DiF(X)) >= 0 in any one > >> direction) > >> S1 .eq. 0 if (Infimum (DiF(X)) > 0 or Supremum (DiF(X)) < 0 in all the > >> directions) > >> > >> (i.e. delete the box if the function is strictly monotonic in all the > >> directions). > >> > >> The modified steps 8 and 13 are, > >> Step 8: if S1(1) is 0 and Z(1) does not intersect the boundary of Z(0), > >> discard Z(1) and go to step 12. > >> Step 13: if S1(2) is 0 and Z(2) does not intersect the boundary of Z(0), > >> discard Z(2) and go to step 17. > >> > >> Results obtained after modification > >> Global minimum = 0.0E0 > >> Global minimum box = [0.0E0, 0.0E0] [0.0E0, 0.0E0] > >> > >> I would highly appreciate if you kindly spend some of your > >> precious time to go through the algorithm and tell me whether the changes > >> I have made are correct or not. > >> > >> References: > >> [1] C. Jansson, O. Knuppel. "A Global Minimization Method: The > >> Multi-Dimensional Case." TUHH (Technische Universitat Hamburg-Harburg), > >> Bericht 92.1, Januar, 1992, pp.47-48. > >> > >> > >> > >> > > > > > ---------------------------------------------------------------------------- ---- > > --------------------------------------------------------------- > R. Baker Kearfott, rbk [at] louisiana [dot] edu (337) 482-5346 (fax) > (337) 482-5270 (work) (337) 981-9744 (home) > URL: http://interval.louisiana.edu/kearfott.html > Department of Mathematics, University of Louisiana at Lafayette > Box 4-1010, Lafayette, LA 70504-1010, USA > --------------------------------------------------------------- > ------=_NextPart_000_002A_01C1FC21.2BA5EFE0 Content-Type: application/octet-stream; name="Moore.OT1" Content-Transfer-Encoding: quoted-printable Content-Disposition: attachment; filename="Moore.OT1" Output from FIND_GLOBAL_MIN on 05/15/2002 at 13:43:08. Version for the system is: October 10, 2000 =20 Codelist file name is: MooreG.CDL Box data file name is: Moore.DT1 =20 =20 =20 Initial box: [ -0.1000D+03, 0.1000D+03 ] [ -0.1000D+03, 0.1000D+03 ] =20 BOUND_CONSTRAINT: F F F F =20 --------------------------------------- CONFIGURATION VALUES: =20 EPS_DOMAIN: 0.1000D-07 MAXITR: 20000 MAX_CPU_SECONDS: 0.1080D+05 DO_INTERVAL_NEWTON: T QUADRATIC: T FULL_SPACE: F VERY_GOOD_INITIAL_GUESS: F USE_SUBSIT: F OUTPUT UNIT: 7 PRINT_LENGTH: 1 PHI_MUST_CONVERGE: T EQ_CNS_MUST_CONVERGE: T INEQ_CNS_MUST_CONVERGE: T PHI_THICKNESS_FACTOR: 0.5000D+00 EQ_CNS_THICKNESS_FACTOR: 0.5000D+00 INEQ_CNS_THICKNESS_FACTOR: 0.5000D+00 PHI_MUST_CONVERGE: T EQ_CNS_MUST_CONVERGE: T INEQ_CNS_MUST_CONVERGE: T PHI_CONVERGENCE_FACTOR: 0.1000D-13 EQ_CNS_CONVERGENCE_FACTOR: 0.1000D-13 INEQ_CNS_CONVERGENCE_FACTOR: 0.1000D-13 CONTINUITY_ACROSS_BRANCHES: F SINGULAR_EXPANSION_FACTOR: 0.1000D+02 HEURISTIC PARAMETER ALPHA: 0.5000D+00 APPROX_OPT_BEFORE_BISECTION: T USE_TAYLOR_EQUALITY_CONSTRAINTS F USE_TAYLOR_INEQ_CONSTRAINTS F USE_TAYLOR_OBJECTIVE F USE_TAYLOR_EQ_CNS_GRD F USE_TAYLOR_GRAD F USE_TAYLOR_INEQ_CNS_GRD F USE_TAYLOR_REDUCED_INEWTON F COSY_POLYNOMIAL_ORDER 5 LEAST_SQUARES_FUNCTIONS: F NONLINEAR_SYSTEM: F Default point optimizer was used.=20 =20 THERE WERE NO BOXES IN THE LIST OF SMALL BOXES. =20 LIST OF BOXES CONTAINING VERIFIED FEASIBLE POINTS: =20 Box no.: 1 Box coordinates: [ -0.4950D-09, 0.4950D-09 ] [ -0.1095-306, 0.1095-306 ] PHI: [ -0.1110D-14, 0.2221D-14 ] B%LIUI(1,*): F F B%LIUI(2,*): F F B%SIDE(*): F F B%PEEL(*): F F Level: 0 Box contains the following approximate root: 0.0000D+00 0.0000D+00 OBJECTIVE ENCLOSURE AT APPROXIMATE ROOT: [ -0.1332D-14, 0.2220D-14 ] Unknown =3D T Contains_root =3D T Changed coordinates: F F ------------------------------------------------- ALGORITHM COMPLETED WITH LESS THAN THE MAXIMUM NUMBER, 20000 OF BOXES. Number of evaluations of a row of the point expanded Jacobi matrix: 27 No. dense interval residual evaluations -- gradient code list: = 14 Number of orig. system C-LP preconditioner rows: 2 Total number of forward_substitutions: 48 Number of Gauss--Seidel steps on the dense system: 2 Number point dense residual evaluations, gradient codelist: = 2 No. point dense Hessian matrix evaluations -- gradient codelist: 1 Total number of iterations of one of the point Newton methods: 1 Total number of dense slope matrix evaluations: 8 Total number second-order interval evaluations of the original function: 6 Total number dense interval reduced gradient evaluations: 3 Number of times the interval Newton method made a coordinate interval smaller: 2 Number of times the approximate solver was called: 2 Number Fritz-John matrix evaluations: 1 BEST_ESTIMATE: 0.2220D-14 Overall CPU time: 0.1001D-01 CPU time in PEEL_BOUNDARY: 0.0000D+00 CPU time in REDUCED_INTERVAL_NEWTON: 0.0000D+00 ------=_NextPart_000_002A_01C1FC21.2BA5EFE0 Content-Type: application/octet-stream; name="Moore.DT1" Content-Transfer-Encoding: 7bit Content-Disposition: attachment; filename="Moore.DT1" 1D-8 -100 100 -100 100 F F F F ------=_NextPart_000_002A_01C1FC21.2BA5EFE0-- From owner-reliable_computing [at] interval [dot] louisiana.edu Wed May 15 15:51:55 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4FKptI18619 for reliable_computing-outgoing; Wed, 15 May 2002 15:51:55 -0500 (CDT) Received: from mailbox.univie.ac.at (mailbox.univie.ac.at [131.130.1.27]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4FKpki18615 for ; Wed, 15 May 2002 15:51:50 -0500 (CDT) Received: from univie.ac.at (hektor.mat.univie.ac.at [131.130.16.21]) by mailbox.univie.ac.at (8.12.2/8.12.2) with ESMTP id g4FKpX8S068992; Wed, 15 May 2002 22:51:35 +0200 Message-ID: <3CE2CA55.643F2BC [at] univie [dot] ac.at> Date: Wed, 15 May 2002 22:51:33 +0200 From: Arnold Neumaier Organization: University of Vienna X-Mailer: Mozilla 4.78 [en] (X11; U; Linux 2.4.9-31 i686) X-Accept-Language: en, de MIME-Version: 1.0 To: RC-mailing-list Subject: Taylor forms - use and limits References: <3CAC7B04.C43A9306 [at] univie [dot] ac.at> Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Following last months discussion on the mailing list about the Taylor models of Martin Berz, I wrote a survey article detailing on 37 pages my current knowledge about Taylor arithmetic in its various forms. Below you can find the abstract and the conclusions. It can be downloaded from the following address: Arnold Neumaier, Taylor forms - use and limits http://www.mat.univie.ac.at/~neum/papers.html#taylor Arnold Neumaier Abstract ======== This review is a response to recent discussions on the reliable computing mailing list, and to continuing uncertainties about the properties and merits of Taylor forms, multivariate higher degree generalizations of centered forms. They were invented around 1980 by Lanford, documented in detail in 1984 by Eckmann, Koch and Wittwer, and independently studied and popularized since 1996 by Berz, Makino and Hoefkens. A highlight is their application to the verified integration of asteroid dynamics in the solar system in 2001, although the details given are not sufficient to check the validity of their claims. Apart from summarizing what Taylor forms are and do, this review puts them into the perspective of more traditional methods, in particular centered forms, discusses the major applications, and analyzes some of their elementary properties. Particular emphasis is given to overestimation properties and the wrapping effect. A deliberate attempt has been made to offer value statements with appropriate justifications; but all opinions given are my own and might be controversial. Conclusions =========== 1. The various versions of Taylor arithmetic constitute a significant enhancement of the toolkit of interval analysis techniques. Indeed, interval coefficient versions were used by mathematical physicists to prove estimates important for computer-assisted proofs, and floating-point coefficient versions were used by Berz and his group to verify solutions of celestial mechanics problems that so far defied interval techniques. Berz and his group also used Taylor models for applications to multivariate integration over a box, differential algebraic equations, and other problems that have not been discussed before in the interval literature. 2. Claims that Taylor arithmetic approximates ranges to higher than second order are unfounded, at least for the currently implemented versions that do not attempt to compute the exact range of the quadratic term. What is true is that in problems where in real arithmetic substantial cancellation of leading digits occurs, Taylor arithmetic may drastically reduce overestimation over intervals narrow enough such that the remainder intervals remain sufficiently small. These are sometimes quite small intervals and sometimes quite wide ones, depending on the problem. If dependent intervals are present but without the cancellation structure, gains from a Taylor arithmetic are slight, and the much higher cost and storage does not justify their use. 3. Claims that Taylor arithmetic eliminates the wrapping effect are unfounded. In the asteroid integrations by Berz et al., the wrapping effect is reduced by means of a technique called shrink wrapping which is poorly documented, but without which rounding errors are quickly magnified to huge intervals, ruining the gain from cancellation. What is true is that Taylor arithmetic may drastically reduce the part of the wrapping effect due to nonlinearities, since the nonlinear dependence on initial conditions is taken into account to a large extent, while methods like Lohner's AWA incur here additional wrapping from the need to enclose bended boxes into parallelepipeds. However, the wrapping incurred by the error term is not addressed at all by Taylor arithmetic, and shrink wrapping seems to be a good remedy only for problems whose Jacobians have no eigenvalues with significantly differing real parts, or where drastic phase space volume reductions compensate the wrapping. This is apparently the case for the asteroid problem and (for some time) for the Lorenz attractor, but it excludes already simple 2-dimensional Hamiltonian dynamical systems with a hyperbolic fixed point, which are very benign from an analytic point of view. 4. In the COSY implementation of Taylor models used by Berz et al. for the claimed verified integration of asteroid dynamics in the solar system, rounding error control is completely undocumented, so that it is difficult to assess whether the implementation is mathematically rigorous. Moreover, the (almost lacking) details given about shrink wrapping are not sufficient to check the validity of their claims. From owner-reliable_computing [at] interval [dot] louisiana.edu Wed May 15 15:55:22 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4FKtLN18710 for reliable_computing-outgoing; Wed, 15 May 2002 15:55:21 -0500 (CDT) Received: from clmboh1-smtp3.columbus.rr.com (clmboh1-smtp3.columbus.rr.com [65.24.0.112]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4FKtGi18706 for ; Wed, 15 May 2002 15:55:16 -0500 (CDT) Received: from oemcomputer (dhcp065-024-174-102.columbus.rr.com [65.24.174.102]) by clmboh1-smtp3.columbus.rr.com (8.11.2/8.11.2) with SMTP id g4FKtBv03959; Wed, 15 May 2002 16:55:11 -0400 (EDT) Message-ID: <00cf01c1fc52$b388f660$66ae1841 [at] columbus [dot] rr.com> From: "Ramon Moore" To: "N. S. Kubal _00423002_" Cc: "interval" Subject: question Date: Wed, 15 May 2002 16:54:26 -0400 MIME-Version: 1.0 Content-Type: multipart/alternative; boundary="----=_NextPart_000_00CC_01C1FC31.2C3D5AA0" X-Priority: 3 X-MSMail-Priority: Normal X-Mailer: Microsoft Outlook Express 5.50.4133.2400 X-MimeOLE: Produced By Microsoft MimeOLE V5.50.4133.2400 Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk This is a multi-part message in MIME format. ------=_NextPart_000_00CC_01C1FC31.2C3D5AA0 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable Dear N. S. Kubal, I have tried to understand the difference between what you propose: I used new variable 'S1' in steps 8 and 13, which is based on the following conditions, S1 .neq. 0 if (Infimum (DiF(X)) <=3D 0 and Supremum (DiF(X)) >=3D 0 = in any one direction) S1 .eq. 0 if (Infimum (DiF(X)) > 0 or Supremum (DiF(X)) < 0 in all = the directions) (i.e. delete the box if the function is strictly monotonic in all = the directions). The modified steps 8 and 13 are, Step 8: if S1(1) is 0 and Z(1) does not intersect the boundary of = Z(0), discard Z(1) and go to step 12. Step 13: if S1(2) is 0 and Z(2) does not intersect the boundary of = Z(0), discard Z(2) and go to step 17. and what was in the 1982 paper at steps 8 and 13 Step 8: if S(1) is empty [meaning that the function is strictly = monotonic in all directions] and Z(1) does not intersect the = boundary of Z(0), discard Z(1) and go to step 12. . . . Step 13: etc. I do not see any difference whatever. I do not understand how you got = the wrong answer other than by having a bug in the implementation. Ramon Moore ------=_NextPart_000_00CC_01C1FC31.2C3D5AA0 Content-Type: text/html; charset="iso-8859-1" Content-Transfer-Encoding: quoted-printable
    Dear N. S. Kubal,
     
    I have tried to understand the difference = between what you=20 propose:
     
    I used new variable 'S1' in steps 8 and 13, which is based = on
    the=20 following conditions,
        S1 .neq. 0 if (Infimum = (DiF(X))=20 <=3D 0 and Supremum (DiF(X)) >=3D 0 in any one=20 direction)
        S1 .eq.  0 if (Infimum (DiF(X)) = > 0 or=20 Supremum (DiF(X)) < 0 in all the directions)
        = (i.e.=20 delete the box if the function is strictly monotonic in all the=20 directions).

        The = modified=20 steps 8 and 13 are,
        Step 8:  if S1(1) is 0 = and Z(1)=20 does not intersect the boundary of Z(0),
        discard = Z(1) and=20 go to step 12.
        Step 13: if S1(2) is 0 and Z(2) = does not=20 intersect the boundary of Z(0),
        discard Z(2) and = go to=20 step 17.
     
    and what was in the 1982 paper at steps 8 and=20 13
     
        Step 8: if S(1) is empty = [meaning that=20 the function is strictly monotonic in all directions] and Z(1) does not=20         intersect the boundary of=20 Z(0),
        discard Z(1) and go to step=20 12.
        . . .
        Step 13: etc.
     
    I do not see any difference whatever. I do not = understand=20 how you got the wrong answer other than by having a bug in the=20 implementation.
     
    Ramon Moore
    ------=_NextPart_000_00CC_01C1FC31.2C3D5AA0-- From owner-reliable_computing [at] interval [dot] louisiana.edu Fri May 17 09:46:35 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4HEkYV23661 for reliable_computing-outgoing; Fri, 17 May 2002 09:46:34 -0500 (CDT) Received: from clmboh1-smtp3.columbus.rr.com (clmboh1-smtp3.columbus.rr.com [65.24.0.112]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4HEkSi23657 for ; Fri, 17 May 2002 09:46:29 -0500 (CDT) Received: from oemcomputer (dhcp065-024-174-102.columbus.rr.com [65.24.174.102]) by clmboh1-smtp3.columbus.rr.com (8.11.2/8.11.2) with SMTP id g4HEkFK21521; Fri, 17 May 2002 10:46:15 -0400 (EDT) Message-ID: <007601c1fdb1$7b3082a0$66ae1841 [at] columbus [dot] rr.com> From: "Ramon Moore" To: "N. S. Kubal _00423002_" Cc: "interval" References: Subject: Re: question Date: Fri, 17 May 2002 10:45:23 -0400 MIME-Version: 1.0 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: 7bit X-Priority: 3 X-MSMail-Priority: Normal X-Mailer: Microsoft Outlook Express 5.50.4133.2400 X-MimeOLE: Produced By Microsoft MimeOLE V5.50.4133.2400 Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Dear Nandkishor, I apologize for taking so long to read carefully the problem you raised. Yes, you are correct. The algorithm given in the 1982 paper throws the baby out with the bath water in the example you showed. In fact, on page 2 of your attachment (pdf version), at step 7 in Iteration 2, we already have the exact solution as b_1= 0 and Z(1)=([0,0],[0,0]). Clearly we do NOT want to delete that. Your way of fixing the problem is fine, testing for strict monotonicity, as you have pointed out in your example. In the 1982 paper, on page 232, after (3.2) is: "where S is the set of integers i such that D_i F(X) properly contains zero and . . ." but the following definition of (u_i,v_i) should have been, as you point out, using strict inequalities to correspond with what was intended in the definition of S. Thank you very much for pointing this out. Best Wishes:- Ramon Moore ----- Original Message ----- From: "N. S. Kubal _00423002_" To: "Ramon Moore" Sent: Friday, May 17, 2002 9:39 AM Subject: Re: question > > > Dear Sir, > I am bothering you again with the same doubt I have raised in one of my > mails to you. By strictly monotonic I understand D_F(X)<0 or D_F(X)>0. > But in your paper 'S' is empty when sup(D_F(X))<=0 or inf(D_F(X))>=0. When > 'S' is defined as given in the paper it deletes the solution box. This is > the reason why I felt there is a need of new variable in place of 'S'. I > attached a document (.ps as well as .pdf) which gives results of each > steps of the algorithm as given in the paper for a particular function and > domain. > > With Regards > > Nandkishor Kubal > > > > > > On Thu, 16 May 2002, Ramon Moore wrote: > > > Dear N. S. Kubal, > > > > I have tried to understand the difference between what you propose: > > > > I used new variable 'S1' in steps 8 and 13, which is based on > > the following conditions, > > S1 .neq. 0 if (Infimum (DiF(X)) <= 0 and Supremum (DiF(X)) >= 0 in any one direction) > > S1 .eq. 0 if (Infimum (DiF(X)) > 0 or Supremum (DiF(X)) < 0 in all the directions) > > (i.e. delete the box if the function is strictly monotonic in all the directions). > > > > The modified steps 8 and 13 are, > > Step 8: if S1(1) is 0 and Z(1) does not intersect the boundary of Z(0), > > discard Z(1) and go to step 12. > > Step 13: if S1(2) is 0 and Z(2) does not intersect the boundary of Z(0), > > discard Z(2) and go to step 17. > > > > and what was in the 1982 paper at steps 8 and 13 > > > > Step 8: if S(1) is empty [meaning that the function is strictly monotonic in all directions] and Z(1) does not intersect the boundary of Z(0), > > discard Z(1) and go to step 12. > > . . . > > Step 13: etc. > > > > I do not see any difference whatever. I do not understand how you got the wrong answer other than by having a bug in the implementation. > > > > Ramon Moore > > > > -- > **************************************************************************** * > contact me on > > nskubal [at] ee [dot] iitb.ac.in > nanduskubal [at] yahoo [dot] com > > Department Residence > > NANDKISHOR S. KUBAL ROOM NO 287 > SYSTEM AND CONTROL ENGG. HOSTEL NO 1 > ELECTRICAL DEPT. > IIT,BOMBAY PHONE NO > POWAI,MUMBAI. > (internal) 5601 > (external) 5721017 / 5720093 > PHONE NO > 5722545 EXT 5741 > From owner-reliable_computing [at] interval [dot] louisiana.edu Fri May 17 11:23:57 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4HGNud23918 for reliable_computing-outgoing; Fri, 17 May 2002 11:23:56 -0500 (CDT) Received: from nwkea-mail-1.sun.com (nwkea-mail-1.sun.com [192.18.42.13]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4HGNqi23914 for ; Fri, 17 May 2002 11:23:52 -0500 (CDT) Received: from engmail2.Eng.Sun.COM ([129.146.1.25]) by nwkea-mail-1.sun.com (8.9.3+Sun/8.9.3) with ESMTP id JAA12976; Fri, 17 May 2002 09:23:45 -0700 (PDT) Received: from phys-mpkmaila (phys-mpkmaila.Eng.Sun.COM [129.146.18.131]) by engmail2.Eng.Sun.COM (8.9.3+Sun/8.9.3/ENSMAIL,v2.1p1) with ESMTP id JAA02332; Fri, 17 May 2002 09:23:45 -0700 (PDT) Received: from gww (gww.Eng.Sun.COM [129.146.78.116]) by mpkmail.eng.sun.com (iPlanet Messaging Server 5.2 Patch 1 (built Apr 2 2002)) with SMTP id <0GW900B8UKANV0 [at] mpkmail [dot] eng.sun.com>; Fri, 17 May 2002 09:25:35 -0700 (PDT) Date: Fri, 17 May 2002 09:23:45 -0700 (PDT) From: William Walster Subject: Re: Fw: Difficulty in your paper To: nskubal [at] ee [dot] iitb.ac.in, rmoore17 [at] columbus [dot] rr.com Cc: reliable_computing [at] interval [dot] louisiana.edu Reply-to: William Walster Message-id: <0GW900B8VKANV0 [at] mpkmail [dot] eng.sun.com> MIME-version: 1.0 X-Mailer: dtmail 1.3.0 @(#)CDE Version 1.4.2 SunOS 5.8 sun4u sparc Content-type: TEXT/plain; charset=us-ascii Content-transfer-encoding: 7BIT Content-MD5: Q3o1dA60lXyUHvw6JBGZyw== Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Dear Ray and Nandkishor, The following URL contains a small test program that our interval team prepared some time ago. It gets the correct answer to Griewank's problem. http://wwws.sun.com/software/Developer-products/fortran/interval/general/Minimum .html#griewank2.f90 While it is important to find and fix problems, whatever their source, I think it is important to point out that there is no doubt about the correctness of the interval approach to solving this and other nonlinear problems. It is *this* fact for which I will always be grateful to you, Ray. Best regards, Bill G. William (Bill) Walster, Ph.D. Interval Technology Engineering Manager Sun Microsystems, Inc. 16 Network Circle, MS UMPK16-304 Menlo Park, CA 94025 (650) 786-9004 Direct (650) 786-9551 Fax bill.walster [at] eng [dot] sun.com From owner-reliable_computing [at] interval [dot] louisiana.edu Sun May 19 17:42:09 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4JMg9T29869 for reliable_computing-outgoing; Sun, 19 May 2002 17:42:09 -0500 (CDT) Received: from cs.utep.edu (mail.cs.utep.edu [129.108.5.3]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4JMg4i29865 for ; Sun, 19 May 2002 17:42:04 -0500 (CDT) Received: from aragorn (aragorn [129.108.5.35]) by cs.utep.edu (8.11.3/8.11.3) with SMTP id g4JMfc704281 for ; Sun, 19 May 2002 16:41:38 -0600 (MDT) Message-Id: <200205192241.g4JMfc704281 [at] cs [dot] utep.edu> Date: Sun, 19 May 2002 16:41:38 -0600 (MDT) From: Vladik Kreinovich Reply-To: Vladik Kreinovich Subject: from NA Digest To: reliable_computing [at] interval [dot] louisiana.edu MIME-Version: 1.0 Content-Type: TEXT/plain; charset=us-ascii Content-MD5: stMq/b3He8XstoBzQ3ymBw== X-Mailer: dtmail 1.3.0 @(#)CDE Version 1.4 SunOS 5.8 sun4u sparc Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk From: Arnold Neumaier Date: Fri, 10 May 2002 12:26:58 +0200 Subject: The Mathematics Genealogy Project I came across a very interesting WWW-site: The Mathematics Genealogy Project http://hcoonce.math.mankato.msus.edu/index.html The site contains names, advisors, and titles of dissertations of mathematicians with a Ph.D., throughout the history Currently (May 6th, 2002) there are records of 55249 mathematicians. While fairly complete on US mathematicians, the directory is still quite incomplete for the rest of the world. The authors invite submissions updating incomplete or missing information. Please contribute and make the directory more complete and hence more useful. Arnold Neumaier From owner-reliable_computing [at] interval [dot] louisiana.edu Sun May 19 19:22:02 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4K0M2h00131 for reliable_computing-outgoing; Sun, 19 May 2002 19:22:02 -0500 (CDT) Received: from clmboh1-smtp3.columbus.rr.com (clmboh1-smtp3.columbus.rr.com [65.24.0.112]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4K0Lvi00125 for ; Sun, 19 May 2002 19:21:57 -0500 (CDT) Received: from oemcomputer (dhcp065-024-174-102.columbus.rr.com [65.24.174.102]) by clmboh1-smtp3.columbus.rr.com (8.11.2/8.11.2) with SMTP id g4K0LrK15172; Sun, 19 May 2002 20:21:53 -0400 (EDT) Message-ID: <000801c1ff94$1cb8cb80$66ae1841 [at] columbus [dot] rr.com> From: "Ramon Moore" To: "Vladik Kreinovich" Cc: "interval" References: <200205192241.g4JMfc704281 [at] cs [dot] utep.edu> Subject: Re: from NA Digest Date: Sun, 19 May 2002 20:20:13 -0400 MIME-Version: 1.0 Content-Type: text/plain; charset="iso-8859-1" Content-Transfer-Encoding: 7bit X-Priority: 3 X-MSMail-Priority: Normal X-Mailer: Microsoft Outlook Express 5.50.4133.2400 X-MimeOLE: Produced By Microsoft MimeOLE V5.50.4133.2400 Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk Fascinating! Thanks. Just for fun I looked up a few people I know or knew, and there they surely were! I looked up Charles Noble, who was an 85 year old emeritus professor (and former chairman) of the mathematics dept at Berkeley when I was a student there. I lived as a guest in his house for more than a year during 1949-50. I recall many fascinating stories he told about his days as a student at Go"ttingen at the end of the 19th century. Indeed, if you look him up in the Math Geneology webpages, you will find he got his Ph.D. from, David Hilbert in 1901. If you look up Hilbert, you will find that he was advisor to 71 Ph.D.students. When you look up Henri Poincare', however, you will find no students. I was bit surprised by that, because Charles Noble told stories about hikes together with both Hilbert and Poincare', when Henri was visiting at Go"ttingen. I do not know of any Ph.D. advisees of Poincare', but I know he was around students, at least at Go"ttingen. Anyway, thanks Arnold and Vladik for pointing out this fun website. Ramon Moore ----- Original Message ----- From: "Vladik Kreinovich" To: Sent: Sunday, May 19, 2002 6:41 PM Subject: from NA Digest > From: Arnold Neumaier > Date: Fri, 10 May 2002 12:26:58 +0200 > Subject: The Mathematics Genealogy Project > > I came across a very interesting WWW-site: > > The Mathematics Genealogy Project > http://hcoonce.math.mankato.msus.edu/index.html > > The site contains names, advisors, and titles of dissertations > of mathematicians with a Ph.D., throughout the history > > Currently (May 6th, 2002) there are records of 55249 mathematicians. > While fairly complete on US mathematicians, the directory > is still quite incomplete for the rest of the world. > The authors invite submissions updating incomplete or missing > information. > > Please contribute and make the directory more complete and > hence more useful. > > Arnold Neumaier > > > From owner-reliable_computing [at] interval [dot] louisiana.edu Mon May 20 13:30:39 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4KIUdL02290 for reliable_computing-outgoing; Mon, 20 May 2002 13:30:39 -0500 (CDT) Received: from mail.comset.net (mail.comset.net [213.172.0.3]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4KIUWi02286 for ; Mon, 20 May 2002 13:30:32 -0500 (CDT) Received: from 8-026.dialup.comset.net ([213.172.8.26] helo=e0gumi46) by mail.comset.net with smtp (Exim 3.33 #1) id 179rwt-0003QK-00 for reliable_computing [at] interval [dot] louisiana.edu; Mon, 20 May 2002 22:31:59 +0400 Message-ID: <005a01c2002b$a2f42c20$710bfea9 [at] wplus [dot] net> From: "Vyacheslav Nesterov" To: "RC mailing list" Subject: Reliable Computing, Volume 8, issue 4, 2002 Date: Mon, 20 May 2002 22:15:27 +0400 MIME-Version: 1.0 Content-Type: multipart/alternative; boundary="----=_NextPart_000_0023_01C2004B.D8C91840" X-Priority: 3 X-MSMail-Priority: Normal X-Mailer: Microsoft Outlook Express 5.00.2615.200 X-MimeOLE: Produced By Microsoft MimeOLE V5.00.2615.200 Sender: owner-reliable_computing [at] interval [dot] louisiana.edu Precedence: bulk This is a multi-part message in MIME format. ------=_NextPart_000_0023_01C2004B.D8C91840 Content-Type: text/plain; charset="koi8-r" Content-Transfer-Encoding: quoted-printable 4-2002 Reliable Computing Volume 8, issue 4, 2002 Mathematical Research Sparse Systems in Fixed Point Form Ramon E. Moore 249-265 On Existence and Uniqueness Verification for Non-Smooth Functions R. Baker Kearfott 267-282 Interval Estimations of Solution Sets to Real-Valued Systems of Linear or Non-Linear Equations Miguel A. Sainz, Ernest Gardenyes, Lambert Jorba 283-305 Short Communication Interval Arithmetic on Multimedia Architectures Juergen Wolff von Gudenberg 307-312 Letter to the Editor Grand Challenges and Scientific Standards in Interval Analysis Arnold Neumaier 313-320 ------=_NextPart_000_0023_01C2004B.D8C91840 Content-Type: text/html; charset="koi8-r" Content-Transfer-Encoding: quoted-printable
              &nbs= p;         =20 4-2002
              &nb= sp;    =20 Reliable=20 Computing
              =   =20 Volume 8, issue 4, 2002
     
              &nbs= p; =20 Mathematical Research
     
    Sparse Systems in Fixed Point = Form
    Ramon=20 E. Moore
    249-265
     
    On Existence and Uniqueness = Verification=20 for  Non-Smooth Functions
    R. Baker = Kearfott
    267-282
     
    Interval Estimations of Solution = Sets=20 to  Real-Valued Systems of Linear
    or Non-Linear = Equations
    Miguel A.=20 Sainz, Ernest Gardenyes, Lambert Jorba
    283-305
     
             Short=20 Communication
     
    Interval Arithmetic on Multimedia = Architectures
    Juergen Wolff von Gudenberg
    307-312
     
             Letter to the=20 Editor
     
    Grand Challenges and Scientific = Standards in=20 Interval Analysis
    Arnold Neumaier
    313-320
     
     
    ------=_NextPart_000_0023_01C2004B.D8C91840-- From owner-reliable_computing [at] interval [dot] louisiana.edu Tue May 21 09:56:21 2002 Received: (from daemon@localhost) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) id g4LEuLt04701 for reliable_computing-outgoing; Tue, 21 May 2002 09:56:21 -0500 (CDT) Received: from kathmandu.sun.com (kathmandu.sun.com [192.18.98.36]) by interval.louisiana.edu (8.11.3/8.11.3/ull-interval-math-majordomo-1.3) with ESMTP id g4LEuEi04697 for ; Tue, 21 May 2002 09:56:15 -0500 (CDT) Received: from engmail1.Eng.Sun.COM ([129.146.1.13]) by kathmandu.sun.com (8.9.3+Sun/8.9.3) with ESMTP id IAA25043 for ; Tue, 21 May 2002 08:56:13 -0600 (MDT) Received: from phys-mpkmaila (phys-mpkmaila.Eng.Sun.COM [129.146.18.131]) by engmail1.Eng.Sun.COM (8.9.3+Sun/8.9.3/ENSMAIL,v2.1p1) with ESMTP id HAA22812 for ; Tue, 21 May 2002 07:56:13 -0700 (PDT) Received: from gww (gww.Eng.Sun.COM [129.146.78.116]) by mpkmail.eng.sun.com (iPlanet Messaging Server 5.2 Patch 1 (built Apr 2 2002)) with SMTP id <0GWG0076GUWWN8 [at] mpkmail [dot] eng.sun.com> for reliable_computing [at] interval [dot] louisiana.edu; Tue, 21 May 2002 07:58:08 -0700 (P
    All information g= iven herein is strictly confidential and will not be re-distributed for any reason other than for the specific use intended.<= br> To be removed from our distribution lists, please Cli= ck here..