ninterp.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/.
!||====================================================================
!||    ninterp   ../engine/source/interfaces/int14/ninterp.F
!||--- called by ------------------------------------------------------
!||    i14dmp    ../engine/source/interfaces/int14/i14dmp.F
!||    i14frt    ../engine/source/interfaces/int14/i14frt.F
!||    i6damp    ../engine/source/interfaces/inter3d/i6damp.F
!||    i6fri3    ../engine/source/interfaces/inter3d/i6fri3.F
!||====================================================================
      SUBROUTINE ninterp(IFUNC,NPC,PLD,NPOINT,XV,YV)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IFUNC,NPC(*),NPOINT
C     REAL
      my_real
     .  pld(*),xv(*),yv(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .   pente, y
      INTEGER RANG(MVSIZ)
      INTEGER NDEB, NREST, I,
     .        IPT, NPT, II
C-----------------------------------------------
      ndeb=0
      nrest=npoint
 100  CONTINUE
C---------------------------------
      DO i=1,min(mvsiz,nrest)
        rang(i)=0
      ENDDO
C---------------------------------
      npt = (npc(ifunc+1)-npc(ifunc))/2
      DO ipt=1,npt-1
        ii=npc(ifunc)+2*(ipt-1)
        DO i=1,min(mvsiz,nrest)
          IF (pld(ii)<=xv(ndeb+i)
     .    .AND.xv(ndeb+i)<=pld(ii+2)) rang(i)=ipt
        ENDDO
      ENDDO
      ii=npc(ifunc)+2*(npt-1)
      DO i=1,min(mvsiz,nrest)
          IF (xv(ndeb+i)>pld(ii)) rang(i)=npt
      ENDDO
C---------------------------------
      DO i=1,min(mvsiz,nrest)
        IF (rang(i)==0) THEN
          ii=npc(ifunc)
          pente=(pld(ii+3)-pld(ii+1))/(pld(ii+2)-pld(ii))
          yv(ndeb+i) =pld(ii+1)+pente*(xv(ndeb+i)-pld(ii))
        ELSEIF (rang(i)==npt) THEN
          ii=npc(ifunc)+2*(npt-1)
          pente=(pld(ii+1)-pld(ii-1))/(pld(ii)-pld(ii-2))
          yv(ndeb+i) =pld(ii+1)
     .       +max(-pld(ii+1),pente*(xv(ndeb+i)-pld(ii)))
        ELSE
          ii=npc(ifunc)+2*(rang(i)-1)
          pente=(pld(ii+3)-pld(ii+1))/(pld(ii+2)-pld(ii))
          yv(ndeb+i) =pld(ii+1)
     .       +pente*(xv(ndeb+i)-pld(ii))
        ENDIF
      ENDDO
C---------------------------------
C     Groupe suivant.
C---------------------------------
      IF (nrest-mvsiz>0) THEN
        nrest=nrest-mvsiz
        ndeb =ndeb +mvsiz
        GOTO 100
      END IF
C-----------------------------------------------
      RETURN
      END