raaupdate.jl

########################################################################################################
### GCMMA-MMA-Julia                                                                                  ### 
###                                                                                                  ###
### This file is part of GCMMA-MMA-Julia. GCMMA-MMA-Julia is licensed under the terms of GNU         ###
### General Public License as published by the Free Software Foundation. For more information and    ###
### the LICENSE file, see <https://github.com/pollinico/GCMMA-MMA-Julia/blob/main/LICENSE>.          ###
###                                                                                                  ###
### The orginal work is written by Krister Svanberg in MATLAB.                                       ###
### This is the Julia version of the code written by Nicolò Pollini.                                 ###
### version 18-05-2023                                                                               ###
########################################################################################################

#-------------------------------------------------------------
#
#    Copyright (C) 2007 Krister Svanberg
#
#    This file, raaupdate.m, is part of GCMMA-MMA-code.
#    
#    GCMMA-MMA-code is free software; you can redistribute it and/or
#    modify it under the terms of the GNU General Public License as 
#    published by the Free Software Foundation; either version 3 of 
#    the License, or (at your option) any later version.
#    
#    This code is distributed in the hope that it will be useful,
#    but WITHOUT ANY WARRANTY; without even the implied warranty of
#    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
#    GNU General Public License for more details.
#    
#    You should have received a copy of the GNU General Public License
#    (file COPYING) along with this file.  If not, see 
#    <http://www.gnu.org/licenses/>.
#    
#    You should have received a file README along with this file,
#    containing contact information.  If not, see
#    <http://www.smoptit.se/> or e-mail mmainfo@smoptit.se or krille@math.kth.se.
#
#------
#
#  Version April 2007.
#
#  Values of the parameters raa0 and raa are updated
#  during an inner iteration.
#
function raaupdate(xmma,xval,xmin,xmax,low,upp,f0valnew,fvalnew,f0app,fapp,raa0,raa,raa0eps,raaeps,epsimin)
  #
  raacofmin = 10^(-12)
  eeem = ones(size(raa))
  eeen = ones(size(xmma))
  xmami = xmax-xmin
  xmamieps = 0.00001*eeen
  xmami = max(xmami,xmamieps)
  xxux = (xmma-xval)./(upp-xmma)
  xxxl = (xmma-xval)./(xmma-low)
  xxul = xxux.*xxxl
  ulxx = (upp-low)./xmami
  raacof = xxul'*ulxx
  raacof = max(raacof,raacofmin)
  #
  f0appe = f0app + 0.5*epsimin
  if f0valnew > f0appe
    deltaraa0 = (1/raacof)*(f0valnew-f0app)
    zz0 = 1.1*(raa0 + deltaraa0)
    zz0 = min(zz0,10*raa0)
  #  zz0 = min(zz0,1000*raa0)
    raa0 = copy(zz0)
  end
  #
  fappe = fapp + 0.5*epsimin*eeem
  fdelta = fvalnew-fappe
  deltaraa = (1/raacof)*(fvalnew-fapp)
  zzz = 1.1*(raa + deltaraa)
  zzz = min(zzz,10*raa)
  #zzz = min(zzz,1000*raa)
  raa[findall(x -> x>0,fdelta)] = zzz[findall(x -> x>0, fdelta)]
#---------------------------------------------------------------------
return raa0,raa
end