Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- CURV_KVEC%KNOTS(knot)=(/0.d0/)
- The shapes of the array expressions do not conform.
- cdegree=2
- count=5
- do knot=0,cdegree+count
- if (knot.le.cdegree) then
- CURV_KVEC%KNOTS(knot)=0.d0
- test(knot+1)=0.d0
- elseif (knot.ge.count) then
- test(knot+1)=1.d0
- CURV_KVEC%KNOTS(knot)=1.d0
- elseif(nvec.eq.0) then
- num=num+1
- do kn=1,cdegree
- CURV_KVEC%KNOTS(knot+kn-1)=0.5d0 !dble(num)/(dble(count-cdegree-1)/dble(cdegree)+1.d0)
- test(knot+kn)=0.5d0 !dble(num)/(dble(count-cdegree-1)/dble(cdegree)+1.d0)
- nvec=cdegree-1
- enddo
- elseif(nvec.ne.0) then
- nvec=nvec-1
- endif
- write(14,*)knot,num,nvec,test(knot+1)
- enddo
- ! CURV_KVEC%KNOTS(0:CURV_KVEC%LENGTH) = (/0.D0,0.D0,0.D0,0.5d0,0.5d0,1.d0,1.d0,1.d0/)
- CURV_KVEC%KNOTS(knot) = 0.d0
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement