finter.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
#include "my_real.inc"
!||====================================================================
!||    finter                    ../starter/source/tools/curve/finter.F
!||--- called by ------------------------------------------------------
!||    cm58_refsta               ../starter/source/materials/mat/mat058/cm58_refsta.F
!||    fct_fiber_2               ../starter/source/materials/mat/mat098/lossfun_98.F
!||    hm_read_eos_compaction2   ../starter/source/materials/eos/hm_read_eos_compaction2.F90
!||    inimom_fvm                ../starter/source/elements/solid/solide/inimom_fvm.F
!||    law81_upd                 ../starter/source/materials/mat/mat081/law81_upd.F90
!||    law87_upd                 ../starter/source/materials/mat/mat087/law87_upd.F90
!||    lossfun_98                ../starter/source/materials/mat/mat098/lossfun_98.F
!||    m36init                   ../starter/source/materials/mat/mat036/m36init.F90
!||    m38init                   ../starter/source/materials/mat/mat038/m38init.F
!||    m75init                   ../starter/source/materials/mat/mat075/m75init.F
!||    m77init                   ../starter/source/materials/mat/mat077/m77init.F
!||    sigeps35                  ../starter/source/materials/mat/mat035/sigeps35.F
!||    sigeps37                  ../starter/source/materials/mat/mat037/sigeps37.F
!||    sigeps38                  ../starter/source/materials/mat/mat038/sigeps38.F
!||    sigeps42                  ../starter/source/materials/mat/mat042/sigeps42.F
!||    sigeps70                  ../starter/source/materials/mat/mat070/sigeps70.F
!||    sigeps90                  ../starter/source/materials/mat/mat090/sigeps90.F
!||====================================================================
      my_real FUNCTION finter(IFUNC,XX,NPF,TF,DERI)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
      INTEGER ifunc,npf(*),i
      my_real tf(*), deri, xx, dx1, dx2
 
      finter = 0
C
      dx2 = tf(npf(ifunc)) - xx
C
C     Constant Function
      IF ((npf(ifunc+1)-npf(ifunc))==2) THEN
        finter = tf(npf(ifunc)+1)
        RETURN
      ELSE
        DO 100 i=npf(ifunc)+2,npf(ifunc+1)-2,2
        dx1 = -dx2
        dx2 = tf(i) - xx
        IF(dx2>=zero.OR.i==npf(ifunc+1)-2)THEN
          deri = (tf(i+1) - tf(i-1)) / (tf(i) - tf(i-2))
          IF(dx1<=dx2)THEN
            finter = tf(i-1) + dx1 * deri
          ELSE
            finter = tf(i+1) - dx2 * deri
          ENDIF
          RETURN
        ENDIF
 100    CONTINUE
      END IF
C
      RETURN
      END