home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #30 / NN_1992_30.iso / spool / sci / math / numanal / 3645 < prev    next >
Encoding:
Internet Message Format  |  1992-12-21  |  5.1 KB
  1. Path: sparky!uunet!dtix!mimsy!ingber
  2. From: ingber@umiacs.umd.edu (Lester Ingber)
  3. Newsgroups: sci.math.num-analysis
  4. Subject: Re: A question about VFSR
  5. Message-ID: <62948@mimsy.umd.edu>
  6. Date: 19 Dec 92 18:54:33 GMT
  7. References: <1992Dec18.193249.22059@murdoch.acc.Virginia.EDU>
  8. Sender: news@mimsy.umd.edu
  9. Organization: UMIACS, University of Maryland, College Park, MD 20742
  10. Lines: 122
  11.  
  12. Sorry, somehow my reply got posted to the wrong query.
  13. _____________________
  14.  
  15. Perhaps I can be of some help since I developed the code.  It might be
  16. of use to the general community if you would post all your replies.
  17.  
  18. In article <1992Dec18.193249.22059@murdoch.acc.Virginia.EDU> you write:
  19. >Hi,
  20. >
  21. >  I am thinking about using VFSR (A Simulated Reannealing code from
  22. >  NETLIB) on my problem which is to minimize a nonlinear function
  23. >  subject to nonlinear inequality constraints. In addition, the
  24. >  constraints are nondifferentiable wrt the parameters although
  25. >  the objective function is differentiable.
  26. >
  27. >  I dont know anything about the theory behind VFSR and I wanted to
  28. >  use it as a black box. My question is:
  29. >  Can I use VFSR for this purpose? 
  30.  
  31. I've used similar constraints.  One case that was published was in
  32. %A L. Ingber
  33. %T Statistical mechanics of neocortical interactions:
  34. A scaling paradigm applied to electroencephalography
  35. %J Phys. Rev. A
  36. %N 6
  37. %V 44
  38. %P 4017-4060
  39. %D 1991
  40.  
  41. >  It is my understanding that VFSR requires the user to supply a
  42. >  subroutine which returns with the cost fuction and a flag which
  43. >  determines whether the (parameter) state at that point violates some
  44. >  constraints or not. Apparently, if the constraints are violated, a
  45. >  new set of trial parameters are generated. Is there a guarantee that
  46. >  the new set of trial parameters are valid? If not, then it would seem
  47. >  to me that you could get stuck at such points.
  48.  
  49. Yes, you could get stuck if the space being sampled had no other
  50. good points.  I.e., on return from the cost function with an invalid
  51. flag, a new random point is selected.  If you think your space might
  52. be flooded with invalid points, you might set up some adaptive
  53. schedule to change the ranges during the run.  I have done this
  54. in other problems, e.g., by passing the parameter maximum and
  55. minimum values along with the cost function.  On lines 133-135
  56. in vfsr.h I have placed some comment for this possibility.
  57.  
  58. >  I would appreciate any information on the issues above.
  59. >  Thanks.
  60. >
  61. >
  62. >  -sriram srinivasan
  63. >   (ss4i@kelvin.seas.virginia.edu)
  64.  
  65. Article 6507 of sci.math.num-analysis:
  66. Path: mimsy!ingber
  67. From: ingber@umiacs.umd.edu (Lester Ingber)
  68. Newsgroups: sci.math.num-analysis
  69. Subject: Re: Nonlinear systems of equations
  70. Message-ID: <62947@mimsy.umd.edu>
  71. Date: 19 Dec 92 18:49:22 GMT
  72. References: <92Dec17.174419.27750@acs.ucalgary.ca> <168C1EBFA.MTURELLA@cmsa.gmr.com>
  73. Sender: news@mimsy.umd.edu
  74. Organization: UMIACS, University of Maryland, College Park, MD 20742
  75. Lines: 57
  76.  
  77.  
  78. Perhaps I can be of some help since I developed the code.  It might be
  79. of use to the general community if you would post all your replies.
  80.  
  81. In article <1992Dec18.193249.22059@murdoch.acc.Virginia.EDU> you write:
  82. >Hi,
  83. >
  84. >  I am thinking about using VFSR (A Simulated Reannealing code from
  85. >  NETLIB) on my problem which is to minimize a nonlinear function
  86. >  subject to nonlinear inequality constraints. In addition, the
  87. >  constraints are nondifferentiable wrt the parameters although
  88. >  the objective function is differentiable.
  89. >
  90. >  I dont know anything about the theory behind VFSR and I wanted to
  91. >  use it as a black box. My question is:
  92. >  Can I use VFSR for this purpose? 
  93.  
  94. I've used similar constraints.  One case that was published was in
  95. %A L. Ingber
  96. %T Statistical mechanics of neocortical interactions:
  97. A scaling paradigm applied to electroencephalography
  98. %J Phys. Rev. A
  99. %N 6
  100. %V 44
  101. %P 4017-4060
  102. %D 1991
  103.  
  104. >  It is my understanding that VFSR requires the user to supply a
  105. >  subroutine which returns with the cost fuction and a flag which
  106. >  determines whether the (parameter) state at that point violates some
  107. >  constraints or not. Apparently, if the constraints are violated, a
  108. >  new set of trial parameters are generated. Is there a guarantee that
  109. >  the new set of trial parameters are valid? If not, then it would seem
  110. >  to me that you could get stuck at such points.
  111.  
  112. Yes, you could get stuck if the space being sampled had no other
  113. good points.  I.e., on return from the cost function with an invalid
  114. flag, a new random point is selected.  If you think your space might be
  115. flooded with invalid points, you might set up some adaptive schedule to
  116. change the ranges during the run.  I have done this in other problems,
  117. e.g., by passing the parameter maximum and minimum values along with
  118. the cost function.  On lines 133-135 in vfsr.h I have placed some
  119. comment for this possibility.
  120.  
  121. >  I would appreciate any information on the issues above.
  122. >  Thanks.
  123. >
  124. >
  125. >  -sriram srinivasan
  126. >   (ss4i@kelvin.seas.virginia.edu)
  127.  
  128. Lester
  129. -- 
  130. || Prof. Lester Ingber                         [10ATT]0-700-L-INGBER ||
  131. || Lester Ingber Research                        Fax: 0-700-4-INGBER ||
  132. || P.O. Box 857                           Voice Mail: 1-800-VMAIL-LI ||
  133. || McLean, VA  22101                EMail: ingber@alumni.caltech.edu ||
  134.