Advertisement
Guest User

Untitled

a guest
Apr 16th, 2014
52
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.77 KB | None | 0 0
  1. CURV_KVEC%KNOTS(knot)=(/0.d0/)
  2.  
  3. The shapes of the array expressions do not conform.
  4.  
  5. cdegree=2
  6. count=5
  7.  
  8. do knot=0,cdegree+count
  9. if (knot.le.cdegree) then
  10. CURV_KVEC%KNOTS(knot)=0.d0
  11. test(knot+1)=0.d0
  12. elseif (knot.ge.count) then
  13. test(knot+1)=1.d0
  14. CURV_KVEC%KNOTS(knot)=1.d0
  15. elseif(nvec.eq.0) then
  16. num=num+1
  17. do kn=1,cdegree
  18. CURV_KVEC%KNOTS(knot+kn-1)=0.5d0 !dble(num)/(dble(count-cdegree-1)/dble(cdegree)+1.d0)
  19. test(knot+kn)=0.5d0 !dble(num)/(dble(count-cdegree-1)/dble(cdegree)+1.d0)
  20. nvec=cdegree-1
  21. enddo
  22. elseif(nvec.ne.0) then
  23. nvec=nvec-1
  24. endif
  25. write(14,*)knot,num,nvec,test(knot+1)
  26. enddo
  27.  
  28. ! CURV_KVEC%KNOTS(0:CURV_KVEC%LENGTH) = (/0.D0,0.D0,0.D0,0.5d0,0.5d0,1.d0,1.d0,1.d0/)
  29.  
  30. CURV_KVEC%KNOTS(knot) = 0.d0
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement