R/optimbase.isfeasible.R

In sbihorel/optimbase: R Port of the 'Scilab' Optimbase Module

# Copyright (C) 2008-2009 - INRIA - Michael Baudin
# Copyright (C) 2009-2010 - DIGITEO - Michael Baudin
# Copyright (C) 2010-$year$ - Sebastien Bihorel
#
# This file must be used under the terms of the CeCILL.
# This source file is licensed as described in the file COPYING, which
# you should have received as part of this distribution. The terms
# are also available at
# http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
#
# This source code is a R port of the optimbase component
# originally written by Michael Baudin for Scilab.

optimbase.isfeasible <- function(this=NULL,x=NULL){

  isfeasible = 1
  #
  # Check if the point is in the bounds.
  #
  if (isfeasible==1){
    hasbounds <- optimbase.hasbounds(this=this)
    if (hasbounds){
      for (ix in 1:this$numberofvariables){
        xmin <- this$boundsmin[ix]
        xmax <- this$boundsmax[ix]
        xix <- x[ix]
        if (xix < xmin){
          isfeasible <- 0
          this <- optimbase.log(this=this,msg=sprintf('Component #%d/%d of x is lower than min bound %e',
                                                      ix,this$numberofvariables,xmin))
          break
        }
        if (xix > xmax){
          isfeasible <- 0
          this <- optimbase.log(this=this,msg=sprintf('Component #%d/%d of x is greater than max bound %e',
                                                      ix,this$numberofvariables,xmax))
          break
        }
      }
    }
  }
  #
  # Check inequality constraints
  #
  if (isfeasible==1){
    if (this$nbineqconst>0){
      if (this$withderivatives){ 
        tmp <- optimbase.function(this=this,x=x,index=5)
          this <- tmp$this
          f <- tmp$f
          g <- tmp$g
          c <- tmp$c
          gc <- tmp$gc
          index <- tmp$index
        rm(tmp)
      } else {
        tmp <- optimbase.function(this=this,x=x,index=5)
          this <- tmp$this
          f <- tmp$f
          c <- tmp$c
          index <- tmp$index
        rm(tmp)
      }
      for (ic in 1:this$nbineqconst){
        if (c[ic] < 0.0 ){
          this <- optimbase.log(this=this,msg=sprintf('Inequality constraint #%d/%d is not satisfied for x',
                                                      ic,this$nbineqconst))
          isfeasible <- -1
          break
        }
      }
    }
  }

  varargout <- list(this=this,isfeasible=isfeasible)

  return(varargout)

}