了解知识
Fortran中的指针如何使用,功能怎样,下面的的5个例子足可以让你明白一切!
 
对于单个值,用起来很简单,但是无法体现指针的强大功能,
示例1:
 
 
  1. program test_pointer_1   
  2. implicit none   
  3. integer,pointer::p1,p2 !定义指针   
  4. integer::s   
  5. !------------------------------------   
  6. ALLOCATE(p1,p2) !分配内存空间   
  7. read(*,*)p1,p2   
  8. s=p1+p2   
  9. write(*,*)'s=',s   
  10. Deallocate(p1,p2) !释放内存空间   
  11. !------------------------------------   
  12. end program test_pointer_1   
 
用指针来操作数组,可以随意引用数组中的任意片段,是指针的强大功能之一,
请看下面的示例2:
 


  1. program test_pointer_2   
  2. implicit none   
  3. integer,pointer::a(:)   
  4. integer,target::b(5)=(/1,2,3,4,5/)   
  5. !-----------------------------------------   
  6. !---引用整个数组   
  7. a=>b          ! a(1)=1;a(2)=2;a(3)=3;a(4)=4;a(5)=5   
  8. write(*,*) 'a=>b ',a   
  9. !-----------------------------------------   
  10. !---引用数组连续的片段   
  11. a=>b(2:3)     ! a(1)=2;a(2)=3   
  12. write(*,*) 'a=>b(2:3) ',a   
  13. !-----------------------------------------   
  14. !---引用数组不连续的片段   
  15. a=>b(1:5:2)   ! a(1)=1;a(2)=3;a(3)=5   
  16. write(*,*) 'a=>b(1:5:2) ',a   
  17. !-----------------------------------------   
  18. !---引用数组的倒序(相当于将数组左右颠倒)   
  19. a=>b(5:1:-1)  ! a(1)=5;a(2)=4;a(3)=3;a(4)=2;a(5)=1   
  20. write(*,*) 'a=>b(5:1:-1) ',a   
  21. !------------------------------------   
  22. end program test_pointer_2   
 
如果是二维数组,想引用某行或者某列的数据,
 
请看示例3:
 
 
  1. program test_pointer_3   
  2. implicit none   
  3. Integer::i,j,k   
  4. real,pointer::a(:)   
  5. real,target::b(3,4)   
  6. !-----------------------------------------   
  7. !---二维数组赋初值   
  8. !data b/1,2,3,4,5,6,7,8,9,10,11,12/   
  9. data ((b(i,j),i=1,3),j=1,4)/1,2,3,4,5,6,7,8,9,10,11,12/   
  10. !-----------------------------------------   
  11. !---打印二维数组内容   
  12. do i=1,3   
  13.     write(*,*) b(i,:)   
  14. end do   
  15. !-----------------------------------------   
  16. !---引用二维数组第二行   
  17. a=>b(2,:)     
  18. write(*,*) 'a=>b(2,:) ',a   
  19. !-----------------------------------------   
  20. !---引用二维数组第三列   
  21. a=>b(:,3)     
  22. write(*,*) 'a=>b(:,3)  ',a   
  23. !------------------------------------   
  24. end program test_pointer_3   
 
 
程序运行结果如下:
 
 
 
上面的指针是一维的,对于用二维指针来引用二维数组中的二维片段也是同样方便,
 
请看示例4:
 
  1. program test_pointer_4   
  2. implicit none   
  3. Integer::i,j,k   
  4. real,pointer::p(:,:) !二维的指针   
  5. real,target::b(3,4)   
  6. !-----------------------------------------   
  7. !---二维数组赋初值   
  8. !data b/1,2,3,4,5,6,7,8,9,10,11,12/   
  9. data ((b(i,j),i=1,3),j=1,4)/1,2,3,4,5,6,7,8,9,10,11,12/   
  10. !-----------------------------------------   
  11. !---打印二维数组内容   
  12. write(*,*) 'b(3,4)'   
  13. do i=1,3   
  14.     write(*,*) b(i,:)   
  15. end do   
  16. !-----------------------------------------   
  17. !---引用二维数组中的二维片段   
  18. p=>b(1:2,2:3)     
  19. write(*,*) 'p=>b(1:2,2:3)'   
  20. do i=1,2   
  21.     write(*,*) p(i,:)   
  22. end do   
  23. !-----------------------------------------   
  24. p=0    !将b(1:2,2:3)中的元素置为0   
  25. !---打印二维数组内容   
  26. write(*,*) 'b(3,4)'   
  27. do i=1,3   
  28.     write(*,*) b(i,:)   
  29. end do   
  30. !------------------------------------   
  31. end program test_pointer_4   
 
 
运行如下:
 
 
把指针当做动态数组来使用,也是很方便的:
 
 
  1. program test_pointer_5   
  2. implicit none   
  3. Integer::i   
  4. real,pointer::p(:) !利用指针为数组动态分配内存   
  5. !-----------------------------------------   
  6. allocate(p(5))     !---为数组(指针)分配内存空间   
  7. !-----------------------------------------   
  8. !---数组赋值   
  9. do i=1,5   
  10.     p(i)=i*2   
  11. end do   
  12. !-----------------------------------------   
  13. write(*,*) 'p(5)=',p !p(5)=2,4,6,8,10   
  14. !-----------------------------------------   
  15. deallocate(p)   !---为数组(指针)释放内存空间   
  16. end program test_pointer_5   
 
 
上面的指针使用是最初级的,高级的用法还有实现指针链表和type类型的指针数组等
对于指针下面的语句代码常常用到:
 
 
  1. ! 常用语句   
  2. POINTER_C=>NULL();          ! 指针指向NULL   
  3. ! NULLIFY(PTR_A);             ! 指针指向NULL   
  4. STATUS = ASSOCIATED (C);    ! 如果C有所指向返回TRUE   
  5. STATUS = ASSOCIATED (C, E); ! 如果C指向E返回TRUE,如果C和E指向的目标相同或者二者均为NULL返回TRUE(二者都是指针)   
 
标签: Fortran
扩展知识