फोरट्रान

Question

में सबराउटिन के लिए आवंटित सरणी कैसे पास करें निम्न कोड एक सेगमेंटेशन फॉल्ट लौटा रहा है क्योंकि आवंटित सरणी जिसे मैं पास करने का प्रयास कर रहा हूं उसे ठीक से पहचाना नहीं जा रहा है (आकार रिटर्न 1, जब यह 3 होना चाहिए)। इस पृष्ठ में (http://www.eng-tips.com/viewthread.cfm?qid=170599) एक समान उदाहरण यह इंगित करता है कि इसे F95 में ठीक काम करना चाहिए; मेरी कोड फ़ाइल में एक एफ 0 9 0 एक्सटेंशन है, लेकिन मैंने इसे F95 में बदलने की कोशिश की, और मैं संकलन करने के लिए gfortran का उपयोग कर रहा हूँ।

मेरा अनुमान है कि समस्या इस तरह से होनी चाहिए कि मैं subroutine को आवंटित सरणी को पार कर रहा हूं; मैं क्या गलत कर रहा हूं?

!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%! 
PROGRAM test 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%! 
IMPLICIT NONE 
DOUBLE PRECISION,ALLOCATABLE :: Array(:,:) 
INTEGER      :: iii,jjj 

ALLOCATE(Array(3,3)) 
DO iii=1,3 
DO jjj=1,3 
    Array(iii,jjj)=iii+jjj 
    PRINT*,Array(iii,jjj) 
ENDDO 
ENDDO 
CALL Subtest(Array) 

END PROGRAM 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%! 
SUBROUTINE Subtest(Array) 
DOUBLE PRECISION,ALLOCATABLE,INTENT(IN) :: Array(:,:) 
INTEGER         :: iii,jjj 

PRINT*,SIZE(Array,1),SIZE(Array,2) 
DO iii=1,SIZE(Array,1) 
DO jjj=1,SIZE(Array,2) 
    PRINT*,Array(iii,jjj) 
ENDDO 
ENDDO 

END SUBROUTINE 
!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%! 

Answer 1

एक प्रक्रिया एक डमी तर्क एक allocatable, तो एक स्पष्ट इंटरफेस किसी भी बुला दायरे में आवश्यक है कि है।

आपको लगता है कि स्पष्ट इंटरफेस प्रदान कर सकते हैं अपने आप को मुख्य कार्यक्रम के अंदर अपने सबरूटीन के लिए एक इंटरफेस ब्लॉक डालकर (ऐसे कई बातें एक स्पष्ट इंटरफेस की आवश्यकता है कि कर रहे हैं, एक allocatable डमी लेकिन एक है।)। एक वैकल्पिक और दूर, दूर, कहीं बेहतर विकल्प मॉड्यूल के अंदर सबराउटिन डालना है और फिर मुख्य प्रोग्राम में उस मॉड्यूल का उपयोग करना है - स्पष्ट इंटरफ़ेस स्वचालित रूप से बनाया जाता है। इस बारे में एक उदाहरण है जिस पर आपने एक लिंक प्रदान किया है - xwb द्वारा पोस्ट देखें। भी अपनी स्थिति को क्वेरी, यह पुनः निर्दिष्ट करें, यह पुनःआवंटन, आदि

Answer 2

कृपया -

ध्यान दें कि यह केवल एक डमी तर्क के लिए समझ में आता है allocatable विशेषता अगर आप अपनी आवंटन स्थिति से संबंधित कुछ करने के लिए जा रहे हैं के लिए ध्यान दें कि आपके आवंटनीय डमी तर्क array को intent(in) के साथ घोषित किया गया है, जिसका अर्थ है कि इसकी आवंटन स्थिति संबंधित वास्तविक तर्क की होगी (और प्रक्रिया के दौरान इसे बदला नहीं जा सकता है)। आपके सबराउटिन को पारित वास्तविक तर्क को एक स्पष्ट इंटरफ़ेस के साथ भी संदर्भित करने के लिए अवैध और इसलिए गैरकानूनी हो सकता है। कंपाइलर इसे नहीं जानता और size जैसे पूछताछ का व्यवहार ऐसे मामलों में अपरिभाषित है।

इसलिए, आपको पहले इसकी सामग्री को संदर्भित करने से पहले के allocated(array) के आवंटन स्थिति की जांच करनी होगी।

subroutine subtest(array) 
    double precision, allocatable, intent(in) :: array(:,:) 
    integer         :: iii, jjj 

    if(allocated(array)) then 
    print*, size(array, 1), size(array, 2) 
    do iii = lbound(array, 1), ubound(array, 1) 
     do jjj = lbound(array, 2), ubound(array, 2) 
     print*, array(iii,jjj) 
     enddo 
    enddo 
    endif 
end subroutine 

Answer 3

यह एक सरल उदाहरण एक साथ allocatable डमी तर्क का उपयोग करता है: मैं आगे, lbound और ubound के साथ पूर्ण सरणी पर छोरों को लागू करने के बाद से सामान्य रूप में आप array की सीमा के बारे में निश्चित नहीं किया जा सकता सुझाव है मॉड्यूल।

module arrayMod 
    real,dimension(:,:),allocatable :: theArray  
end module arrayMod 

program test 
    use arrayMod 
    implicit none 

    interface 
     subroutine arraySub 
     end subroutine arraySub 
    end interface 

    write(*,*) allocated(theArray) 
    call arraySub 
    write(*,*) allocated(theArray) 
end program test 

subroutine arraySub 
    use arrayMod 

    write(*,*) 'Inside arraySub()' 
    allocate(theArray(3,2)) 
end subroutine arraySub