Geant4 का निरंतर सत्यापन

मैंने Geant4 परियोजना के सत्यापन के बारे में सही लेख लिखा था। मुझे आपको पृष्ठभूमि याद दिलाने दें। हाल ही में, Geant4 लाइब्रेरी के एक पुराने संस्करण की जाँच की गई और " कॉपी-पेस्ट एंड मून्स" लेख लिखा गया। पुराने संस्करण की जाँच क्यों की गई? लोग परिपूर्ण नहीं हैं। निरीक्षण का सार पिछले लेख में पाया जा सकता है। अब Geant4 संस्करण 10.0-बीटा के सत्यापन पर एक छोटी रिपोर्ट आपके ध्यान में पेश की गई है।

पिछले लेख का सारांश

लेख "कॉपी-पेस्ट एंड मून्स" स्थैतिक विश्लेषण पद्धति और पीवीएस-स्टूडियो की क्षमताओं का उपयोग करने के लाभों का वर्णन करता है। दाढ़ी वाले Geant4 प्रोजेक्ट संस्करण 4.9.4 की जाँच की गई थी। कई संदिग्ध कोड टुकड़े पाए गए थे, जिन्हें लेख में वर्णित किया गया था।

Geant4 मोंटे कार्लो विधियों का उपयोग करके पदार्थ के माध्यम से प्राथमिक कणों के पारित होने का अनुकरण करने के लिए एक कार्यक्रम है। CERN में ऑब्जेक्ट ओरिएंटेड C ++ प्रोग्रामिंग लैंग्वेज में विकसित किया गया। यह GEANT के पिछले संस्करणों का एक और विकास है, काफी संशोधित और पूरक है।

परियोजना की वेबसाइट: http://geant4.org

पिछले लेख में, 16 अत्यधिक संदिग्ध कोड अंशों का वर्णन किया गया था। एक नए चेक से उनमें से केवल 10 का पता चला। बाकी स्पष्ट रूप से या तो तय हैं या यह कोड अब मौजूद नहीं है। मैं फिर से एक नए लेख में इन कोड अंशों का वर्णन नहीं करूंगा। उन्हें पिछले लेख में देखा जा सकता है। यह इंगित करता है कि नए संस्करण में विश्लेषक चेतावनी गायब हो गई या नहीं।

मैं इस तरह के एक अजीब शोध प्रारूप के लिए माफी माँगता हूँ। लेकिन मुझे उम्मीद है कि यह Geant4 परियोजना में कई कमियों को ठीक करने और पीवीएस-स्टूडियो टूल पर लोगों का ध्यान आकर्षित करने के लिए चोट नहीं पहुंचाती है।

अगर मैं गलत नहीं हूं, तो लाइब्रेरी का पिछला संस्करण 2011 को संदर्भित करता है। इस समय के दौरान, बहुत कुछ बदल गया है, और यह आश्चर्य की बात नहीं है कि कोड में कुछ नए अजीब स्थान पाए गए थे। आइए देखें कि क्या नया है या क्या मैंने पहले नहीं देखा था।

नया संदिग्ध कोड स्निपेट

जिन संदिग्ध स्थानों पर मैंने ध्यान दिया, उनकी पूरी सूची geant4_new.txt फ़ाइल में है। मैं इस बात पर जोर देना चाहता हूं कि आपको इस सूची पर पूरी तरह भरोसा नहीं करना चाहिए। परियोजना को स्वयं जांचना और सभी चेतावनियों का विश्लेषण करना बेहतर है। हम कुछ समय के लिए Geant4 डेवलपर्स की कुंजी के लिए प्रदान करने के लिए तैयार हैं: एक प्रतिक्रिया फ़ॉर्म ।

वही सुविधाएँ

G4double G4CsvAnalysisManager::GetH2Xmin(G4int /*id*/) const { ExceptionForHistograms("GetH2Xmin"); return 0; } G4double G4CsvAnalysisManager::GetH2Xmax(G4int /*id*/) const { ExceptionForHistograms("GetH2Xmin"); return 0; } 

PVS-Studio चेतावनी: V524 यह विचित्र है कि 'GetH2Xmax' फ़ंक्शन का शरीर 'GetH2Xmin' फ़ंक्शन के शरीर के समतुल्य है। _G4analysis- संग्रह g4csvanalysismanager.cc 933

GetH2Xmax () फ़ंक्शन, सबसे अधिक संभावना है, एक अन्य पैरामीटर के साथ ExceptionForHistograms () फ़ंक्शन को कॉल करना चाहिए:

 ExceptionForHistograms("GetH2Xmax"); 

त्रुटि गंभीर नहीं लगती है। जैसा कि मैं इसे समझता हूं, यह एक अपवाद हैंडलिंग है। हालाँकि, मैंने इस कॉपी-पेस्ट त्रुटि का उल्लेख करने का निर्णय लिया।

ऊर्जा शून्य है

कार्य कैलकुलेट करें अचानक, यह चर रीसेट हो गया है। यह बहुत ही संदिग्ध है।

 G4double G4RKFieldIntegrator::CalculateTotalEnergy(const G4KineticTrackVector& Barions) { G4double Etot = 0; .... for(G4int c2 = c1 + 1; c2 < nBarion; c2++) { .... // Esk2 Etot += t1*std::pow(Alpha/pi, 3/2)* std::exp(-Alpha*r12*r12); // Eyuk Etot += ....; // Ecoul Etot += 1.44*p1->GetDefinition()->GetPDGCharge()* p2->GetDefinition()->GetPDGCharge()/r12* Erf(std::sqrt(Alpha)*r12); // Epaul Etot = 0; .... } .... } 

PVS-Studio चेतावनी: V519 'Etot' चर को क्रमिक रूप से दो बार मान दिया गया है। शायद यह एक गलती है। चेक लाइनें: 80, 83. _G4processes -र्कचर g4rkfieldintegrator.cc 83

संदिग्ध तर्क

 G4double G4EmBiasingManager::ApplySecondaryBiasing(....) { .... if(0 == nsplit) { if(safety > fSafetyMin) .... } if(1 == nsplit) { weight = ApplyRussianRoulette(vd, index); } else { G4double tmpEnergy = pPartChange->GetProposedKineticEnergy(); G4ThreeVector tmpMomDir = .... weight = ApplySplitting(vd, track, currentModel, index, tcut); pPartChange->SetProposedKineticEnergy(tmpEnergy); pPartChange->ProposeMomentumDirection(tmpMomDir); } .... } 

पीवीएस-स्टूडियो चेतावनी: V646 आवेदन के तर्क का निरीक्षण करने पर विचार करें। यह संभव है कि 'और' कीवर्ड गायब है। _G4processes- संग्रह g4embiasingmanager.cc 299

कोड को "अगर और" निर्माण के उपयोग के रूप में डिज़ाइन किया गया है। लेकिन ऐसा है नहीं। यदि आप कोड को सही ढंग से प्रारूपित करते हैं, तो हमें यह मिलता है:

 if(0 == nsplit) { if(safety > fSafetyMin) .... } if(1 == nsplit) { weight = ApplyRussianRoulette(vd, index); } else { G4double tmpEnergy = pPartChange->GetProposedKineticEnergy(); G4ThreeVector tmpMomDir = .... weight = ApplySplitting(vd, track, currentModel, index, tcut); pPartChange->SetProposedKineticEnergy(tmpEnergy); pPartChange->ProposeMomentumDirection(tmpMomDir); } 

ध्यान दें कि यदि किसी अन्य शाखा से संबंधित ब्लॉक को हमेशा निष्पादित किया जाता है यदि "1! = Nsplit"। इसमें संदेह है कि काम के एक अलग तर्क की कल्पना की गई थी।

इसी तरह की स्थिति यहाँ देखी जा सकती है:

V646 आवेदन के तर्क का निरीक्षण करने पर विचार करें। यह संभव है कि 'और' कीवर्ड गायब है। _G4processes- संग्रह g4embiasingmanager.cc 347

अधूरा कोड?

 void G4MolecularDecayTable::AddExcitedState(const G4String& label) { channelsMap::iterator channelsIter = fDecayChannelsMap.find(label); if(channelsIter != fDecayChannelsMap.end()) { G4String errMsg = "Excited state" + label + " already registered in the decay table."; G4Exception("G4MolecularDecayTable::AddExcitedState", "G4MolecularDecayTable003", FatalErrorInArgument, errMsg); return; } fDecayChannelsMap[label] ; } 

पीवीएस-स्टूडियो चेतावनी: V607 मालिकाना अभिव्यक्ति 'fDecayChannelsMap [लेबल]'। _G4processes- संग्रह g4moleculardecaytable.cc 140

फ़ंक्शन का अंत बहुत संदिग्ध लगता है:

 fDecayChannelsMap[label] ; 

यह क्या है क्या यहाँ कुछ याद आ रहा है? आप ऐरे सेल के साथ क्या करना चाहते थे?

कीनेमेटीक्स

निम्नलिखित उदाहरण काफी लंबा है। दुर्भाग्य से, आगे भी इसे काटना मेरे लिए कठिन है। देखें कि 'id' में कौन से मान हो सकते हैं।

 void G4QMDCollision::CalKinematicsOfBinaryCollisions( G4double dt) { .... G4int id = 0; .... if ( secs ) { .... id++; .... } if ( std::abs ( eini - efin ) < fepse*10 ) .... else { .... for ( G4int i0i = 0 ; i0i < id-1 ; i0i++ ) { theSystem->DeleteParticipant( i0i+n0 ); } .... } .... } 

चेतावनी पीवीएस-स्टूडियो: वी 621 'ऑपरेटर' के लिए निरीक्षण करने पर विचार करें। यह संभव है कि लूप को गलत तरीके से निष्पादित किया जाएगा या बिल्कुल भी निष्पादित नहीं किया जाएगा। _G4processes- संग्रह g4qmdcollision.cc 228

यदि शर्त "अगर (सेकंड)" पूरी नहीं हुई है, तो चर 'आईडी' शून्य के बराबर रहेगा। फिर, प्रपत्र का एक चक्र उत्पन्न हो सकता है:

 for ( G4int i0i = 0 ; i0i < -1 ; i0i++ ) 

यह एक बहुत ही अजीब चक्र होगा। मुझे लगता है कि यह CalKinematicsOfBinaryCollisions () फ़ंक्शन के तर्क की जांच करने के लायक है।

अन्य

ऐसी कई चेतावनियाँ हैं जिनका मैंने पिछले लेख में वर्णन नहीं किया था। मैं इसमें नहीं रहूंगा। एक उदाहरण के लिए, मैं एक मामला दिखाऊंगा:

 class G4HadronicException : public std::exception { .... }; inline G4double G4GeneralPhaseSpaceDecay::Pmx( G4double e, G4double p1, G4double p2) { if (e-p1-p2 < 0 ) { G4HadronicException(__FILE__, __LINE__, "G4GeneralPhaseSpaceDecay::Pmx " "energy in cms > mass1+mass2"); } .... } 

पीवीएस-स्टूडियो चेतावनी। V596 ऑब्जेक्ट बनाया गया था, लेकिन इसका उपयोग नहीं किया जा रहा है। 'थ्रो' कीवर्ड गायब हो सकता है: G4HadronicException (FOO) को फेंक दें; _G4processes- संग्रह g4generalphasespacedecay.hh 116

'थ्रो' स्टेटमेंट को भूल जाने पर कई त्रुटियां होती हैं। परिणामस्वरूप, 'G4HadronicException' प्रकार की एक वस्तु बनाई जाती है और तुरंत नष्ट हो जाती है। और फिर प्रोग्राम गलत डेटा के साथ काम करना जारी रखता है।

आलेख में वर्णित ये टाइपोस , फ़ाइल geant4_new.txt में पाए जा सकते हैं। माइक्रोपोटिमिटेशन से संबंधित चेतावनी भी हैं।

निष्कर्ष

एक पुरानी परियोजना की जाँच की व्यर्थता के विषय पर। आप देखें, और मैंने छुआ। :)

मुझे ट्रोल करने का अच्छा कारण।