dersonebasisfun.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dersonebasisfun (idxi, pxi, xi, kxi, ders1, ders2)

Function/Subroutine Documentation

dersonebasisfun()

subroutine dersonebasisfun	(	integer	idxi,
		integer	pxi,
		intent(in)	xi,
		intent(in)	kxi,
			ders1,
			ders2 )

Definition at line 28 of file dersonebasisfun.F.

C--------------------------------------------------------------------------------------------------------
C
C This subroutine calculates B-spline function and derivate of B-spline function
C Assembling B spline functions for NURBS is outside this subroutine 
C
C--------------------------------------------------------------------------------------------------------
C VAR      | SIZE      | TYP  |  RW  |  DEFINITION
C--------------------------------------------------------------------------------------------------------
C IDXI     | 1         | I    |  R   | ELEMENT INDEX IN KNOT VECTOR IN XI DIRECTION
C PXI      | 1         | I    |  R   | POLYNOMIAL INTERPOLATION DEGREE IN XI DIRECTION
C XI       | 1         | F    |  R   | COUNTER PARAMETER VALUE (WHERE THE FUNCTION AND DERIVATE ARE CALCULATED)
C--------------------------------------------------------------------------------------------------------
C KXI      | NKXI      | F    |  R   | (FULL) KNOT VECTOR IN XI DIRECTION FOR THE CURRENT PATCH (GROUP)
C--------------------------------------------------------------------------------------------------------
C DERS1    | PXI+1     | F    |  W   | INTERPOLATION FUNCTION
C DERS2    | PXI+1     | F    |  W   | DERIVATE OF INTERPOLATION FUNCTION
C--------------------------------------------------------------------------------------------------------
C Implicite Types
C--------------------------------------------------------------------------------------------------------
#include     "implicit_f.inc"
C--------------------------------------------------------------------------------------------------------
C Dummy Arguments 
C--------------------------------------------------------------------------------------------------------
      INTEGER PXI, IDXI
      my_real, 
     .  INTENT(IN) :: xi
      my_real, 
     .  DIMENSION(*), INTENT(IN) :: kxi
      my_real ders1, ders2
C--------------------------------------------------------------------------------------------------------
C Local variables 
C--------------------------------------------------------------------------------------------------------
      INTEGER J, JJ, K, NDERS
      my_real saved, temp, aleft, right
      my_real, DIMENSION(2) :: ders
      my_real, DIMENSION(PXI+1,PXI+1) :: andu
      my_real, DIMENSION(PXI+1) :: nd
 
      nders=1
      andu(:,:)=zero
 
      DO j=0,pxi
        IF ((xi>=kxi(idxi+j)).AND.(xi<kxi(idxi+j+1))) THEN
          andu(j+1,1) = one
        ELSE
          andu(j+1,1) = zero
        ENDIF
      ENDDO
 
      DO k=1,pxi
        IF (andu(1,k) == 0) THEN
          saved = zero
        ELSE 
          saved = ((xi-kxi(idxi))*andu(1,k))/(kxi(idxi+k)-kxi(idxi))
        ENDIF
        DO j=0,pxi-k
          aleft = kxi(idxi+j+1)
          right = kxi(idxi+j+k+1)
          IF (andu(j+2,k) == 0) THEN
            andu(j+1,k+1) = saved
            saved = zero
          ELSE
            temp = andu(j+2,k)/(right-aleft)
            andu(j+1,k+1) = saved+(right-xi)*temp
            saved = (xi-aleft)*temp
          ENDIF  
        ENDDO
      ENDDO
      
      ders(1) = andu(1,pxi+1) 
 
CCC ERREUR A PARTIR DE LA
 
      DO k=1,nders
        DO j=1,k+1
          nd(j) = andu(j,pxi-k+1)
        ENDDO
        DO jj=1,k
          IF (nd(1) == 0) THEN
            saved = zero
          ELSE 
            saved = nd(1)/(kxi(idxi+pxi-k+jj)-kxi(idxi))
          ENDIF
          DO j=1,k-jj+1
            aleft = kxi(idxi+j)
            right = kxi(idxi+j+pxi+jj-1)
c            RIGHT = KXI(IDXI+J+PXI+JJ+1)
            IF (nd(j+1) == 0) THEN
              nd(j) = (pxi-k+jj)*saved
              saved = zero
            ELSE
              temp = nd(j+1)/(right-aleft)
              nd(j) = (pxi-k+jj)*(saved-temp)
              saved = temp
            ENDIF
          ENDDO
        ENDDO
      ders(2) = nd(1)
      ENDDO
 
      ders1 = ders(1)
      ders2 = ders(2)
 
      RETURN