onebasisfun.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	onebasisfun (idxii, idxi, pxi, xi, kxi, ders1)

Function/Subroutine Documentation

onebasisfun()

subroutine onebasisfun	(	integer	idxii,
		integer	idxi,
		integer	pxi,
		intent(in)	xi,
		intent(in)	kxi,
			ders1 )

Definition at line 28 of file onebasisfun.F.

C--------------------------------------------------------------------------------------------------------
C
C This subroutine calculates 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 DERIVATIVE 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   | DERIVATIVE OF INTERPOLATION FUNCTION
C--------------------------------------------------------------------------------------------------------
C Implicite Types
C--------------------------------------------------------------------------------------------------------
#include     "implicit_f.inc"
C--------------------------------------------------------------------------------------------------------
C Dummy Arguments 
C--------------------------------------------------------------------------------------------------------
      INTEGER :: PXI, IDXI, IDXII
      my_real, INTENT(IN) :: xi
      my_real, DIMENSION(*), INTENT(IN) :: kxi
      my_real :: ders1
C--------------------------------------------------------------------------------------------------------
C Local variables 
C--------------------------------------------------------------------------------------------------------
      INTEGER :: J, K
      my_real :: saved, temp, aleft, right
      my_real, DIMENSION(PXI+1,PXI+1) :: andu
 
      andu(:,:)=zero
      andu(idxii,1)=one
 
c      DO J=0,PXI
c        IF ((XI>=KXI(IDXI+J)).AND.(XI<KXI(IDXI+J+1))) THEN
c          ANDU(J+1,1) = ONE
c        ELSE
c          ANDU(J+1,1) = ZERO
c        ENDIF
c      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
      
      ders1 = andu(1,pxi+1) 
 
      RETURN