Geant4 Cross Reference |
1 # Category field History 2 3 See `CONTRIBUTING.rst` for details of **requir 4 which **must** added in reverse chronological 5 It must **not** be used as a substitute for wr 6 7 ---------------------------------------------- 8 9 ## 2024-11-26 Ivana Hrivnacova (field-V11-02-0 10 - Better names for the G4EquationType enum con 11 12 ## 2024-11-04 Philippe Canal (field-V11-02-05) 13 - Improve stability of static initialization f 14 StepperDescription() in G4DormandPrince745. 15 16 ## 2024-09-10 Ivana Hrivnacova (field-V11-02-0 17 - Coverity fixes: 18 - Use std::move, const auto& instead of auto 19 20 ## 2024-06-20 Ivana Hrivnacova (field-V11-02-0 21 - Fix in G4FieldBuilder for local fields: 22 apply the local fields in the order of user 23 24 ## 2024-04-29 Ivana Hrivnacova (field-V11-02-0 25 - Added classes for automatic field constructi 26 - G4FieldBuilder[Messenger], G4FieldParamete 27 - No changes in other Geant4 classes 28 29 ## 2024-04-17 Ben Morgan (field-V11-02-01) 30 - Fix maybe-uninitialized warnings found in AT 31 - Reported in internal ATLASSIM-6058 ticket. 32 33 ## 2024-04-02 Gabriele Cosmo (field-V11-02-00) 34 - Applied trivial clang-tidy fixes to classes, 35 constructor/destructor, auto and using state 36 - Minor cleanup in QSS classes. 37 38 ## 2023-11-03 Gabriele Cosmo (field-V11-01-06) 39 - Reinstated default DormandPrince745 stepper. 40 41 ## 2023-10-27 Lucio Santi (field-V11-01-05) 42 - Improvements on the QSS Stepper 43 The main QSS integration loop is interrupted 44 exceeds a predefined threshold (set to 1000 45 46 ## 2023-06-14 Gabriele Cosmo (field-V11-01-04) 47 - Applied clang-tidy fixes (readability, moder 48 49 ## 2023-06-12 John Apostolakis (field-V11-01-0 50 - Inclusion of QSS integration method. 51 Most recent changes: 52 * Fixes for issue with SetPrecision (used un 53 * Updated values of parameters in G4Propagat 54 * Withdrew interim changes that made QSS the 55 56 ## 2023-05-26 Lucio Santi 57 - Fixes for open issues: 58 * default values for QSS made stricter dQRel 59 to improve behaviour in full experiment te 60 * implemented missing DistChord method -- ne 61 62 ## 2023-01-31 John Apostolakis 63 - Initial inclusion of QSS integration method. 64 Renamed classes to include 'G4' prefix, and 65 Added Geant4 disclaimer / license to files. 66 Removed QSSmessenger from RunManager and G4P 67 - Temporarily made QSS the default method for 68 69 ## 2023-01-30 Rodrigo Castro, Lucio Santi & Le 70 - Included Quantum State Simulation (QSS) inte 71 Alternative method of integration which crea 72 functions for each coordinate (x,p) - changi 73 polynomials when the deviation from exact so 74 a threshold value. 75 Offers built-in interpolation capability - e 76 of the intersection of the trajectory with s 77 Currently restricted to using only with a pu 78 Created by Rodrigo Castro and Lucio Santi, w 79 by Leandro A. Gómez Vidal (and others?) 80 81 ## 2023-03-16 Ben Morgan (field-V11-01-02) 82 - Fixed compilation warnings on clang-15 for m 83 in input parameters for steppers. 84 85 ## 2023-01-09 Ben Morgan (field-V11-01-01) 86 - Use PROJECT_ scoped _SOURCE/BINARY_DIR varia 87 scope version to ease use of Geant4 as a CMa 88 - Fixes [GitHub PR #52](https://github.com/G 89 90 ## 2022-12-09 Laurie Nevay (field- 91 - Reduced printout for setting any valid value 92 G4FieldManager. 93 94 ## 2022-11-28 Gabriele Cosmo (field- 95 - Fixed restore of stream precision in G4Field 96 97 ## 2022-11-14 John Apostolakis (field- 98 - Revised G4FieldManager to ensure that epsilo 99 are less than a 'maximum accepted' accuracy 100 behaviour of the integration. Improved their 101 - warnings if min > max, with corrective be 102 - a fatal exception in case of values outsi 103 104 To cope with needs of legacy applications or 105 the value of the 'ceiling' maximum accepted 106 the new static method 107 G4FieldManager::SetMaxAcceptedEpsilon( 108 but must remain under or equal to a final ce 109 fMaxFinalEpsilon=0.03 110 111 ## 2022-11-10 Gabriele Cosmo (field-V11-00-03) 112 - Fixed compilation warnings for implicit type 113 114 ## 2022-11-05 Divyansh Tiwari, John Apostolaki 115 116 - Introduced G4BorisScheme and G4BorisDriver, 117 integration method, created as part of GSoC 118 119 ## 2022-10-05 Gabriele Cosmo (field-V11-00-01) 120 - Fixed compilation warnings on Intel/icx comp 121 but not used. 122 123 ## 2021-12-10 Ben Morgan (field-V11-00-00) 124 - Change to new Markdown History format 125 126 --- 127 128 # History entries prior to 11.0 129 130 March 30, 2021 B.Morgan 131 ----------------------- 132 - Migrate sources.cmake to modular build API 133 134 November 26, 2020 J.Apostolakis 135 ------------------------------- 136 - Fix to enable G4TDormandPrince45 to be used 137 138 This fix was needed because G4Interpolation 139 copies of the stepper class. Due to this it 140 with the type of the templated stepper. 141 142 For a stepper to be used with 'G4Interpolati 143 implement the method 144 Equation_type GetSpecificEquation() 145 Templated steppers must ensure that the retu 146 is the specific equation type of the templat 147 148 This method was added to G4DormandPrince745 149 150 October 12, 2020 J.Apostolakis 151 ------------------------------ 152 - Fixes & refinements - mostly in templated cl 153 Compilation fixes, e.g. missing include <cas 154 Refinements, e.g. int -> unsigned int in tem 155 156 October 9, 2020 J.Apostolakis 157 ----------------------------- 158 - Configure G4ChordFinder to use templated G4T 159 default stepper for magnetic fields (when on 160 explicitly by the application.) 161 162 October 9, 2020 J.Apostolakis & J. Xie 163 -------------------------------------- 164 - Adapted classes that avoid virtual calls for 165 created by Josh Xie (Google Summer of Code 2 166 S. Wenzel & J.A. 167 168 The revised design and implementations inclu 169 * equation of motion templated on the field 170 * steppers templated on the type of equatio 171 variables. 172 In addition the key methods are marked as 'i 173 they judge possible, to embed them and avoid 174 175 These classes can be combined also with the 176 (G4IntegrationStepper, G4FSALIntegrationStep 177 to avoid virtual calls in all the levels up 178 the integration driver. ] 179 180 The types of stepper currently available are 181 i) the originally developed intermediate ord 182 183 - Name (templated) - original -- Order/ 184 / E 185 =========================================== 186 G4TClassicalRK4: G4ClassicalRK4 4th / 187 G4TCashKarpRKF45: G4CashKarpRKF45 5th / 188 189 ii) the original low order steppers (for spe 190 G4TSimpleHeum: G4SimpleHeum 3rd / 191 G4TSimpleRunge: G4SimpleRunge 2nd / 192 G4TExplicitEuler: G4ExplicitEuler 1st / 193 194 and 195 iii) some newly adapted stepper(s) / ( J.A. 196 G4TDormandPrince45 G4DormandPrince745 5th/ 197 =========================================== 198 199 Note: Additional lower-order steppers are no 200 'Vector' steppers (using blaze) are no 201 202 Experimental versions of fields with inline 203 G4TUniformMagField, G4TQuadrupoleMagFiel 204 are included to show how the field evaluatio 205 206 Note 2: 207 Steppers which are not embedded obtain their 208 breaking the step into two halves, then comp 209 As a result many more field / derivative eva 210 embedded steppers of the same order, but the 211 estimation has different properties; so it m 212 213 June 15, 2020 G.Cosmo 214 --------------------- 215 - Implemented move constructor and operator fo 216 217 March 24, 2020 J.Apostolakis 218 ---------------------------- 219 - New class G4DriverReporter to print progress 220 Used for debugging differences between Old 221 - Revisions in testing CMake config file and s 222 for magnetic field. 223 224 March 13, 2020 J.Apostolakis 225 ---------------------------- 226 - G4VIntegrationDriver & dependent driver clas 227 Added new virtual StreamInfo() method, 228 & used it to implement operator << for 229 Made (virtual) DoesReIntegrate() const. 230 - G4MagInt_Driver: Fixed max iterations & clar 231 ComputeNewStepSize does NOT (yet) respect 232 ( Done to enable comparisons with new G4Inte 233 - G4OldMagIntDriver maintains all old behaviou 234 235 236 January 22, 2020 G.Cosmo 237 ------------------------ 238 - Fixed compilation errors and configuration f 239 240 January 13, 2020 G.Cosmo 241 ------------------------ 242 - Turn off verbosity flags by default in G4Int 243 G4InterpolationDriver and G4MagIntegratorDri 244 245 December 10, 2019 B.Morgan 246 -------------------------- 247 - Cleanup CMake build, removing obsolete granu 248 explicit include_directories. 249 250 November 14, 2019 J.Apostolakis 251 ------------------------------- 252 - Added method / attribute to G4VIntegrationDr 253 to inform whether the driver recomputes seg 254 AccurateAdvance is called. 255 256 October 27, 2019 D.Sorokin 257 ---------------------------- 258 - Remove statistics printouts from destructor 259 260 October 24, 2019 H.Burkhardt 261 ---------------------------- 262 - Standard sextupole magnet first implementati 263 264 October 18, 2019 D. Sorokin 265 ----------------------------------------- 266 - Fix floating point exception in G4BFieldInte 267 268 October 1, 2019 G. Cosmo 269 ----------------------------------------- 270 - Fixed unitialised data in G4EqMagElectricFie 271 reported by Coverity. 272 273 September 12, 2019 G. Cosmo 274 ----------------------------------------- 275 - Minor C++11 revision. Code formatting and cl 276 277 July 8, 2019 D. Sorokin 278 ----------------------------------------- 279 - G4BFieldIntegrationDriver: an attempt to sol 280 281 July 3, 2019 D. Sorokin 282 ----------------------------------------- 283 - Make G4InterpolationDriver default for perfo 284 285 June 25, 2019 D. Sorokin 286 ----------------------------------------- 287 - Fix Dormand-Prince 5th order interpolation m 288 289 June 20, 2019 D. Sorokin 290 ----------------------------------------- 291 - New version of G4InterpolationDriver with mu 292 293 June 12, 2019 D. Sorokin & J. Apostolakis 294 ----------------------------------------- 295 - Revert to use G4IntegrationDriver as default 296 ( c-tor for magnetic field. ). 297 Change to use concrete type of Stepper as te 298 299 June 10, 2019 Dmitry Sorokin 300 ------------------------- 301 - Remove division in G4DormandPrince745::Inter 302 303 June 7, 2019 Dmitry Sorokin 304 ------------------------- 305 - Switch to original chord finding algorithm i 306 - Reset step estimate in G4MagInt_Driver only 307 - Change G4VERBOSE to G4FIELD_DEBUG in G4Inter 308 309 June 4, 2019 Dmitry Sorokin 310 ------------------------- 311 - Use FSAL property of G4DormandPrince745 in G 312 313 June 3, 2019 Dmitry Sorokin 314 ------------------------- 315 - Fix memory churn in G4InterpolationDriver. 316 317 May 20, 2019 Dmitry Sorokin 318 ------------------------- 319 - Make interpolation driver default for immedi 320 321 April 8, 2019 H.Burkhardt 322 ------------------------- 323 - Correcting transformation to local system in 324 325 March 20, 2019 G.Cosmo 326 ---------------------- 327 - Corrected 'RegularStepperName' in G4ChordFin 328 default stepper G4DormandPrince745 rather th 329 330 November 9, 2018 J.Apostolakis 331 ------------------------------ 332 - Small fixes for warnings in compilation. 333 - Revisions of descriptions of recent tags (fo 334 335 September 26, 2018 D.Sorokin 336 ---------------------------- 337 - Implemented G4InterpolationDriver. This driv 338 stepper with interpolation property to integ 339 The field evaluations used to calculate the 340 be reused to obtain one or multiple internal 341 So in steps in which an intersection is foun 342 evaluations are expected to be reduced subst 343 344 September 13, 2018 D.Sorokin 345 ---------------------------- 346 - G4NystromRK4 creates G4CachedMagneticField i 347 - Removed G4MagIntegratorStepper::ComputeRight 348 (all calls now use non-virtual RightHandSide 349 - Reversed declaration of final for the class 350 Returned to previous status: only GetFieldVa 351 ( This avoid need to change examples/tests a 352 users' code. ) 353 354 September 12, 2018 D.Sorokin 355 ---------------------------- 356 - Created new experimental driver G4BFieldInte 357 magnetic fields. In QuickAdvance it uses he 358 estimated 'helix' rotation angle is larger t 359 - G4NystromRK4 uses G4CachedMagneticField for 360 NOTE: It now the user's responsibility to cr 361 ( The alternative is to create G4CachedMagne 362 But this would have several disadvantages: 363 * It would change the field in the user's 364 * G4NystromRK4 is not able to check whethe 365 that changes the equation using SetEquat 366 367 July 10, 2018 D. Sorokin - field- 368 ------------------------ Tags by J. Apostolaki 369 - Corrections in calculation of chord distance 370 ( use 4th order interpolation from Shampine 371 - Revised G4FieldUtils: moved copy() method fr 372 added 2 new setValue: simple & variadic meth 373 - Fixes to copy input time in ModifiedMidpoint 374 Bogacki-Shampine 4/5 and 2/3 (note: B-S 2/3 375 Dormand Prince 4/5, 5/6 and 7/8, and Tsitou 376 - Trial revision in G4BulirschStoer to reduce 377 operations (J.A.) 378 - Fix for unitialised pIntStepper pointer in G 379 constructor. Revision of data member variab 380 381 June 15, 2018 G.Cosmo - fiel 382 --------------------- 383 - Fixed minor Coverity defects in G4BurlischSt 384 Code cleanup. 385 386 June 13, 2018 G.Cosmo - fiel 387 --------------------- 388 - Leave override final only for clone() method 389 not in destructor... 390 391 June 12, 2018 G.Cosmo - fiel 392 --------------------- 393 - Use G4Exception warnings instead of streamin 394 - Added override final for virtual methods in 395 - Some code cleanup. 396 397 Apr 30, 2018 J. Apostolakis - fiel 398 --------------------------- 399 - G4VIntegrationDriver and derived classes 400 401 Added new virtual method in G4VIntegrationDr 402 void RenewStepperAndAdjust(G4MagIntegrat 403 which used to exist in G4MagInt_Driver. 404 ( Its absence caused a problem in some exist 405 406 This method must be defined in those derived 407 a Runge Kutta stepper ( i.e. some or any typ 408 409 - G4IntegrationDriver, G4FSALIntegrationDriver 410 Defined default implementation (which throws 411 implementations in the templated classes G4I 412 G4FSALIntegrationDriver which check the type 413 change the stepper if it is appropriate. 414 415 Feb 13, 2018 J. Apostolakis - fiel 416 --------------------------- 417 - Change in G4FieldManager::CreateChordFinder 418 better with case that field argument is null 419 420 Feb 13, 2018 D. Sorokin - fiel 421 ----------------------- - fiel 422 - Introduced Bulirsch-Stoer method of integrat 423 to Runge-Kutta based on the mid-point method 424 Separate driver derived G4VIntegrationDriver 425 using specialisation of G4IntegrationDriver 426 Note: this method does not have a RK 'Steppe 427 428 Jan 31, 2018 J. Apostolakis 429 --------------------------- 430 - G4MagneticField no longer inherits from G4El 431 This was incorrect, and allowed a user to pr 432 field, but obtain the force only from the B- 433 the equation created in the G4ChordFinder co 434 => potential interface change. 435 436 Result: It is no longer possible to use the 437 pure magnetic field for the case of a mix 438 This change is desired, and is a *protect 439 use led to incorrect results. 440 Instead the full chain of classes for equati 441 steppers and driver must be used for the 442 443 Feb 9, 2017 G.Cosmo - fiel 444 ------------------- 445 - Removed static data from G4FSALIntegrationDr 446 G4MagIntegratorDriver and moved to base clas 447 Fix required for allowing proper symbol expo 448 449 Nov 17, 2017 J. Apostolakis - fiel 450 --------------------------- 451 - G4ChordFinder: Reinstated old driver G4MagIn 452 with G4DormandPrince745 (5th order) in G4C 453 2nd constructor (with magnetic field.) 454 ( Reverses change made in field-V10-03-27 ) 455 456 Nov 16, 2017 J. Apostolakis - fiel 457 --------------------------- 458 - Suppressed statistics printing in G4ChordFin 459 ( G4IntegrationDriver and G4FSALIntegration 460 protecting them by verbosity flags. 461 462 - Minor fix to initialise scalar data member i 463 G4RK547FEq1, G4RK547FEq2 and G4RK547FEq3. 464 Fixes coverity issues 97886, 97888 and 97889 465 466 467 Nov 14, 2017 J. Apostolakis - fiel 468 --------------------------- 469 - Suppressed extra verbosity in G4ChordFinder. 470 - Fix in G4FSALDormandPrince745 for unitialise 471 Initialise ak8[], ak9[] in constructor. This 472 memory allocation, without deletion - i.e. a 473 - Initialise fLastStepLength and/or fAuxSteppe 474 in constructors of G4DormandPrince745 and G4 475 Coverity issues: 67394 (G4FSALDormandPrince7 476 67391 (G4DormandPrince7 477 478 Nov 14, 2017 J. Apostolakis - fiel 479 --------------------------- - fiel 480 - G4ChordFinder: changes in default stepper & 481 in which only a magnetic field is provided ( 482 483 * Reinstated G4DormandPrince745 (5th order) 484 (corrected -23, which missed to initialise 485 486 * New templated driver G4IntegrationDriver i 487 the simple constuctors of G4ChordFinder an 488 which only require a (pure) magnetic field 489 490 Additional information and details: 491 * New optional parameter to constructor of G 492 allow to create an FSAL stepper and corres 493 This aims to reuse some derivatives, thus 494 of derivatives and thus field evaluations. 495 496 * The templated 'driver' classes, available 497 G4ChordFinder's 2nd constructor, are used 498 class: 499 - G4IntegrationDriver for non-FSAL 500 - G4FSALIntegrationDriver for FSAL-capab 501 These make non-virtual calls to the steppe 502 503 - First fix in G4FSALDormandPrince745 for memo 504 delete[] fInitialDyDx in destructor. Coverit 505 506 ---------------------------------------------- 507 Multiple tags (below) were to better test step 508 Date Tagger Tag name Stepper 509 ============================================== 510 Nov 3, D. Sorokin field-V10-03-25 G4FSALBo 511 Nov 3, J. Apostolakis field-V10-03-24 G4FSALDo 512 Nov 2, J. Apostolakis field-V10-03-23 G4Dorman 513 Nov 2, J. Apostolakis field-V10-03-22 G4RK547F 514 ---------------------------------------------- 515 Glossary of drivers: 516 the FSAL driver = using G4FSALIntegrationDr 517 the new driver = using G4IntegrationDriver 518 ============================================== 519 520 Nov 2, 2017 D. Sorokin - fiel 521 --------------------------------------- 522 - Added new steppers G4RK547FEq{1,2,3} which i 523 Runge-Kutta 4,5 pairs with: 524 - FSAL property (return derivative at end 525 - stable equilibrium states ( improved ha 526 integration: better stability in choo 527 Each one can be used with all existing Runge 528 - G4FSALIntegrationDriver, using its FSAL 529 - G4IntegrationDriver without its FSAL 530 - the original G4MagInt_Driver, also igno 531 532 - Reinstated G4DormandPrince745 ( known widel 533 the default RK method in G4ChordFinder for p 534 ( Withdrew G4BogackiShampine45. ) J. Apos 535 536 Nov 1, 2017 J. Apostolakis - fiel 537 --------------------------- 538 - Made G4BogackiShampine45 the default RK meth 539 for pure magnetic fields. 540 541 Oct 30, 2017 D. Sorokin/J.Apostolakis - fiel 542 --------------------------------------- 543 - add GetStepper method to G4VIntegrationDrive 544 545 Oct 20, 2017 D. Sorokin - fiel 546 --------------------------------------- 547 Introduced base class for any driver algorithm 548 New concrete classes added: 549 - G4IntegrationDriver: implementation of templ 550 for non-fsal explicit Runge-Kutta steppers 551 - G4FSALIntegrationDriver: implementation of d 552 FSAL explicit Runge-Kutta steppers. 553 Existing class G4MagInt_Drv remains, deriving 554 555 - Expanded testing in testProElectroMagField ( 556 ( J. Apostolakis) 557 558 Oct 13, 2017 D. Sorokin - fiel 559 --------------------------------------- 560 - add G4FieldUtils.hh, G4FieldUtils.icc, G4Fie 561 utilities to extract data from arrays in con 562 magneticfield integration. 563 564 Oct 13, 2017 G. Cosmo - fiel 565 --------------------------------------- 566 - G4NystromRK4: use consistenly G4 types for d 567 568 Aug 31, 2017 J. Apostolakis - fiel 569 --------------------------------------- 570 - Made G4DormandPrince457 the default Stepper 571 G4ChordFinder - when a driver/stepper is not 572 573 Aug 31, 2017 J. Apostolakis - fiel 574 --------------------------------------- 575 - G4NystromRK4: used caching of momentum magni 576 to reduce overhead of fix in field-V10-03-10 577 578 Aug 25, 2017 J. Apostolakis - fiel 579 --------------------------------------- 580 - G4MagInt_Driver: added GetEquationOfMotion m 581 code with future develoment. 582 583 Aug 24, 2017 D.Sorokin & J. Apostolakis - fiel 584 --------------------------------------- 585 - G4NystromRK4: correction in call to getField 586 ( Spotted by D. Sorokin ) 587 Also small cleanup of testing code into 2 sm 588 ( J. Apostolakis ) 589 590 Aug 23, 2017 J.Apostolakis - fiel 591 -------------------------- 592 - G4NystromRK4 - Fix the caching of momentum, 593 Fixes issue reported by ALICE ( problem repo 594 Thanks to Ruben Shahoyan for identifying and 595 - Temporarly made G4NystromRK4 the default Ste 596 constructor (to ensure adequate testing). 597 598 Jul 04, 2017 J.Apostolakis - field 599 -------------------------- 600 - G4DormandPrinceRK78: initialise fLastStepLen 601 members in constructor. Coverity issue: 6739 602 Improved description and reference to origin 603 - G4MagIntegratorDriver: initialised fNoCalls 604 605 Jul 03, 2017 J.Apostolakis - field 606 -------------------------- 607 - G4FSALBogackiShampline45: initialise 'bi[][] 608 in method PrepareConstants() called by const 609 ( Similar to G4BogackiShampine => seek con 610 Small fix: initialise fAuxStepper to nullptr 611 ( Part of Coverity issue 67394 ) 612 - G4DormandPrinceRK56: fix for significant mem 613 space in SetupInterpolate_low(), apparently 614 because it was initialised in each use, not 615 Improved initialisation of ak10, 11, 12 (mov 616 and fixed small memory leak (added deletion 617 - G4MagIntegratorStepper: initialise fNoRHSCal 618 ( fIntegrationOrder is not yet used widely - 619 Coverity issue 67392. 620 - G4VFSALIntegrationStepper: initialise counte 621 - G4BogackiShampine23: initialised fIsFSAL, fI 622 Suppressed obsolete virtual method isFSAL, 623 Deleted duplicate virtual method 'isFSAL()'. 624 now exists in base class G4MagIntegratorStep 625 Made copy constructor and equal operator pri 626 Added explanatory comments in source file; c 627 (it calculates derivative @ endpoint. ) 628 - Several steppers ( G4DoLoMcPriRK34, G4FSALBo 629 G4DormandPrince745, G4BogackiShampine45 ): i 630 fAuxStepper data members in constructor. 631 Coverity issues: 67393, 67394, 67391, 67390, 632 633 Jul 03, 2017 J.Apostolakis - field 634 -------------------------- 635 - fix for memory leak in Stepper() in two step 636 G4FSALBogackiShampline45 ( Coverity issue 637 G4FSALDormandPrince745 ( Coverity issue 638 - G4BogackiShampline23: fix for memory leak in 639 ( Coverity issue 67355 ) 640 641 May 17, 2017 J.Apostolakis - field 642 -------------------------- 643 - G4FieldManager::SetDetectorField 644 Fixed default of 'fail mode' parameter in 645 to 0 (silent) - not 1 which was used for 646 In case of existing G4ChordFinder, use 'max' 647 648 649 May 11, 2017 J.Apostolakis - field 650 -------------------------- 651 - Reintroduced changes of field-V10-01-05, 06 652 issue recalled by ATLAS. 653 - New fixes/refinements (since field-V10-01-07 654 o Extra check for missing equation. (Fix) 655 o Renamed G4FieldManager::SetFieldChangesEne 656 InitialiseFieldChangesEnergy() to clarify 657 - Original functionality change (copied from f 658 o G4FieldManager.cc : Revised behaviour of S 659 passes the Field to the Equation, as is th 660 The previous behaviour is possible, using 661 Values: 0 = silent ; 662 1 (default) = warning; 663 2 or other = fatal except 664 665 This change is to ensure that user expecta 666 if the field is changed for an existing, w 667 be propagated to the equation of motion. 668 Note: this causes potential conflict betwe 669 and the field, which will need to be 670 Old behaviour (just set our own pointer) i 671 ProposeDetectorField(G4Field *detectorFiel 672 Note: Error reporting behaviour of SetDete 673 complain by default if ChordFinder i 674 ( The default values of failMode = 0 675 If a chord-finder exists, it will still wa 676 To obtain have noisy extra checks call Set 677 or by using the new method for an existing 678 ChangeDetectorField(G4Field *detectorField 679 OLD: Co-works and requires new tag in geomet 680 681 Mar 22, 2017 G.Cosmo - fiel 682 --------------------- 683 - Revert change made in previous tag and add s 684 singleton pointer. 685 686 Mar 17, 2017 G.Cosmo - fiel 687 --------------------- 688 - Define G4FieldManagerStore as normal singlet 689 690 Oct 7, 2016 J.Apostolakis - fiel 691 --------------------------- 692 - Checked loops for termination for infinite l 693 them with "Loop checking" comment 694 695 Oct 7, 2016 D.Sorokin - fiel 696 ----------------------- Commit: J. Apostolaki 697 - Protection for multiple inclusion of G4Class 698 (thanks to Dmitrii Sorokin, GSoC 2016 stude 699 700 Aug 16, 2016 J.Apostolakis - fiel 701 --------------------------- 702 - Added initialization of fAuxStepper in const 703 G4DoLoMcPriRK34, G4DormandPrince745 & G4Dorm 704 705 June 11, 2016 J.Apostolakis - fiel 706 --------------------------- 707 - Restored G4ClassicalRK4 as default stepper f 708 (due to performance issues reported in bench 709 710 June 10, 2016 G.Cosmo - fiel 711 --------------------- 712 - Fixed code indentation in G4HelixMixedSteppe 713 fixes compilation warnings on gcc-6.1. 714 715 June 1, 2016 J.Apostolakis - fiel 716 -------------------------- 717 - Fix in multiple steppers for the size of two 718 intermediate values; they must be at least 719 'fNumberOfStateVariables' long. 720 Changed: G4CashKarp, 721 G4BogackiShampine 23 & 34 722 G4DormandPrince 745, RK56 & RK78 723 G4DoLoMcPriRK34 724 G4FSALBogackiShampine45, G4FSALDorm 725 - Made G4CashKarp the default stepper in G4Cho 726 - Deleted trial class G4NewChordFinder. 727 728 June 1, 2016 J.Apostolakis - fiel 729 -------------------------- 730 - Deleted test code (RkTestRun) due to cmake i 731 732 June 1, 2016 J.Apostolakis - fiel 733 -------------------------- 734 - Made G4ClassicalRK4 the default stepper agai 735 - Fixes in test code (RkTestRun) for new names 736 737 May 31, 2016 J.Apostolakis - fiel 738 -------------------------- 739 - New class G4NewChordFinder, that uses old G4 740 the newer G4FSALIntegratorDriver -- first 741 Created new test testNewChordFinder.cc to ch 742 743 - Renamed new Stepper classes to use 'G4' pref 744 DormandPrince{745,RK56,RK67} TsitourasRK45 745 746 - Renamed FSAL base classes to use 'G4' prefix 747 FSALMagIntegratorStepper to G4VFSALIntegr 748 FSALIntegratorDriver to G4FSALIntegr 749 and the experimental FSAL versions of steppe 750 FBogackiShampine45.hh to G4FSALBogackiS 751 FDormandPrince745.hh to G4FSALDormandP 752 ( Note: FSAL steppers remain 'behind' the impr 753 754 - Made DormandPrince745 (DP45) the default ste 755 756 May 26, 2016 J.Apostolakis & Som. Banerjee - 757 ------------------------------------------ 758 - First version of FSAL classes (as additional 759 FSAL Integrator Driver : New *concrete*, 760 FSAL Integrator Stepper : New base class ( 761 FBogackiShampine45: FSAL-version of Bogac 762 FDormandPrince745: FSAL-version of FDorm 763 Developed by Somnath Banerjee in GSoC 2015 764 765 Context: 766 * FSAL Stepper's key 'Stepper()' method for St 767 argument to return the derivative at the e 768 If the step was successful, this estimate 769 accuracy as the stepper (i.e. 5th order fo 770 771 - FSAL Driver method uses the revised Stepper 772 the final (endpoint) derivative for succes 773 as the start derivative for the next step, 774 continuing of integration. 775 776 - New test (RkTestRun.cc) covering FSAL method 777 Developed by Somnath Banerjee in GSoC 2015 778 779 May 26, 2016 J.Apostolakis - 780 -------------------------- 781 Preparation for addition of FSAL steppers (& d 782 - G4MagIntegratorStepper: Added counter for ca 783 with Get/Reset methods. 784 - Cosmetic fixes in sources.cmake (tabs to spa 785 786 May 25, 2016 J.Apostolakis & Som. Banerjee - 787 ------------------------------------------ 788 - Added new stepper classes embedded RK method 789 * DoLoMcPriRK34, 6-stage 3/4 RK (interp 790 * DormandPrinceRK56, 9-stage 5/6 RK (interp 791 * DormandPrinceRK78, 13-stage 7/8 RK (interp 792 created by Somnath Banerjee as part of GSoC 793 794 May 25, 2016 J.Apostolakis - field 795 -------------------------- 796 - Revised testPropagateMagField, to reduce err 797 steps. An error message is still written 798 but test does not stop. 799 800 - G4HelixMixedStepper: Fixes & added new (5th 801 Included (G4)Bogacki-Shampine 23/45 , Dor 802 and newest TsitourasRK45. 803 804 May 24, 2016 J.Apostolakis - field 805 -------------------------- 806 - Added TsitourasRK45 stepper. 807 - Made Cash-Karp the default stepper (in G4Cho 808 809 May 23, 2016 J.Apostolakis - field 810 -------------------------- 811 - Made DormandPrince745 (DP45) the default ste 812 813 May 20, 2016 J.Apostolakis - field 814 ------------------------- 815 - Added steppers created by Somnath Banerjee a 816 G4BogackiShampine23 (BS23 below), G4Bogack 817 DormandPrince745 (DP45). 818 These implement third order (BS23) and fifth 819 embedded RK tableaus. 820 - DP45 includes interpolation method, to evalu 821 ('y') at intermediate points of a successf 822 BS45's interpolation method exists, but is n 823 824 References: 825 - BS45 is embedded pair from paper by P. Bog 826 “An efficient Runge-Kutta (4,5) pair,â 827 Comput. Math. with Appl., vol. 32, no. 6 828 - BS23 is four stage RK method from the pape 829 “A 3(2) pair of Runge - Kutta formulas 830 Appl. Math. Lett., vol 2, no 4, p.321 (1 831 - DP45 are from Table 2 of J. R. Dormand an 832 “A family of embedded Runge-Kutta form 833 Journal of comput. and applied Math., vo 834 835 May 4, 2016 G.Cosmo - field- 836 ------------------- 837 - Fixed Clone() signatures in derived Field cl 838 More corrections on implementations of opera 839 840 Apr 29, 2016 G.Cosmo - field- 841 -------------------- 842 - Corrected implementation of operator=() and 843 G4ElectroMagneticField, G4ElectricField and 844 845 Feb 25, 2016 P.Gumplinger - field- 846 ------------------------- 847 - address Problem Report: 1831 848 change: G4Mag_SpinEqRhs, G4EqEMFieldWithSpin 849 G4RepleteEofM.cc from g_BMT = (magMo 850 to g_BMT = (std:: 851 852 853 Dec 11, 2015 A.Dotti - field-V10-02-0 854 -------------------------- 855 Fixing typo in tests/CMakeList.txt file: subdi 856 instead of add_subdirectory() 857 858 Nov 2, 2015 J.Apostolakis - field 859 Oct 22, 2015 J.Apostolakis - field 860 -------------------------- 861 - Reverted changes in G4FieldManager from tag 862 inadvertently reintroduced in tag field-V10 863 (fixed in tag -15: committed new copy from 864 865 - Added explanation comment in G4MagHelicalSte 866 867 Oct 13, 2015 J.Apostolakis - field 868 -------------------------- 869 - G4ChordFinder::FindNextChord: Added check ag 870 to identify potential issue of lack of conve 871 - test: fixed testProElectricMagField.cc 872 873 - Note: tag-13 inadvertently reintroduced chan 874 875 Jun 26, 2015 P.Gumplinger - field 876 ------------------------- 877 - Include 'NumberOfVariables' in G4RepleteEofM 878 879 Jun 22, 2015 G.Cosmo - field 880 -------------------- 881 - Withdrawn also changes introduced in tag "fi 882 responsible for memory violation errors in V 883 Now corresponding to tag "field-V10-01-04". 884 885 Jun 18, 2015 G.Cosmo - field 886 -------------------- 887 - Withdrawn changes introduced from tag "field 888 "field-V10-01-07" included, responsible for 889 run-time errors on Valgrind checks. 890 891 Jun 10, 2015 P.Gumplinger - field 892 ------------------------- 893 - initialize all elements of dydx[] in G4Reple 894 Thanks to Tatsuya Kikawa (Osaka Univ.) 895 896 Jun 10, 2015 J.Apostolakis - field 897 -------------------------- 898 - G4FieldManager.cc : Made behaviour of SetDet 899 longer complains by default if ChordFinder 900 ( The default values of failMode = 0 i.e. 901 If a chord-finder exists, it will still war 902 903 It is still possible to have louder/extra c 904 SetDetectorField with failMode > 0 905 or by using a new method, which expects an 906 ChangeDetectorField(G4Field *detectorFiel 907 908 Jun 8, 2015 J.Apostolakis - field 909 -------------------------- 910 - G4FieldManager.cc : Revised behaviour of Set 911 passes the Field to the Equation, as is the 912 The previous behaviour is possible, using e 913 flag. Values: 0 = silent ; 914 1 (default) = warning; 915 2 or other = fatal except 916 Or can use alternative (new) suggestion met 917 918 This change is to ensure that user expectat 919 if the field is changed for an existing, 920 be propagated to the equation of motion. 921 Note: this causes potential conflict betwee 922 and the field, which will need to be 923 924 Old behaviour (just set our own pointer) is 925 ProposeDetectorField(G4Field *detectorFie 926 927 Co-works and requires new tag in geometry/n 928 929 Jun 4, 2015 P. Gumplinger - field 930 ------------------------- 931 - change assignment to field instead of Field 932 933 Jan 28, 2015 J.Apostolakis - field 934 -------------------------- 935 - Revert unintended inclusion of changes in G4 936 937 Jan 15, 2015 J.Apostolakis - field 938 -------------------------- 939 - Small refinements in G4FieldTrack - access r 940 Avoid division by zero for Unit direction 941 Added method to get Rest Mass. 942 Needed for fixes related to Problem report 1 943 Thanks to Nathan S. Froemming <nfroemm@uw.ed 944 and created dedicated test to demonstrate i 945 946 Nov 03, 2014 G.Cosmo - field 947 -------------------- 948 - Moved constructors and simple methods to in 949 G4LineSection for CPU speedup (courtesy of S 950 951 May 21, 2014 J.Apostolakis - field 952 -------------------------- 953 - G4FieldTrack: 954 * Fixed SetChargeAndMoments() to use the new 955 G4ChargeState::SetChargesAndMoments(); it 956 erroneous change in G4ChargeState::SetChar 957 is now renamed SetChargeSpinMoments() to s 958 * Added pdg-spin as argument to constructors 959 * Renamed 'spin' to 'Polarization' in method 960 avoid confusion; 961 * Revised constructors to add pdg-spin as la 962 - G4ChargeState: 963 - Addressed Coverity issue #56281 'Arguments 964 G4ChargeState::SetChargeMoments() when use 965 G4FieldTrack::SetChargeMoments() from addi 966 without revising name of method to reflect 967 - Changed order of arguments in constructor, 968 (was 2nd) and making all arguments mandato 969 - Renamed SetChargeAndMoments() method with 970 to SetChargeSpinMoments(); 971 - Added new methods including SetChargesAndM 972 SetChargeDipoleMoments() (no spin), and al 973 SetChargeMdmSpin() (which includes spin) t 974 combinations. 975 - G4MagIntegratorDriver::OneGoodStep(): 976 * Trial changes to address Valgrind issue, a 977 978 Apr 01, 2014 J.Apostolakis - fi 979 -------------------------- 980 - Moved G4FieldTrack::LoadFromArray() to .cc a 981 thread-local static data local to the functi 982 983 Mar 14, 2014 J.Apostolakis - fi 984 -------------------------- 985 - Fix assignment operator for G4FieldTrack. (C 986 - Improve printing of G4FieldTrack 987 988 Jan 17, 2014 P.Gumplinger - fi 989 ------------------------- 990 - add G4Mag_EqRhs::SetChargeMomentumMass( part 991 to G4Mag_SpinEqRhs::SetChargeMomentumMass, s 992 http://hypernews.slac.stanford.edu/HyperNews 993 994 Nov 18, 2013 P.Gumplinger - fi 995 ------------------------- 996 - fix Coverity warning in G4RepleteEofM 997 998 Nov 17, 2013 J.Apostolakis 999 -------------------------- 1000 - Fixed and improved HelixMixed Stepper: 1001 Corrected the default stepper for small s 1002 Added Angle-threshold parameter. (Default 1003 Added more Steppers, including G4NystromR 1004 Revised Id numbers of Steppers, to make t 1005 More comments in Stepper(). 1006 1007 Nov 05, 2013 P.Gumplinger - f 1008 ------------------------- 1009 - remove writing out-of-bounds dydx[8] = 0. i 1010 1011 Nov 04, 2013 P.Gumplinger - f 1012 ------------------------- 1013 - add G4RepleteEofM to sources.cmake 1014 f 1015 Nov 01st, 2013 P.Gumplinger - f 1016 --------------------------- 1017 - add G4RepleteEofM: EofM that includes force 1018 as well as mu*grad(B) and spin tracking. 1019 1020 May 10th, 2013 P.Gumplinger/J.Apostolakis - f 1021 ----------------------------------------- 1022 - Revised number of field components (to allo 1023 18 are required for new applications (B, E 1024 - Fixes for initialisation errors found by Co 1025 1026 May 10th, 2013 P.Gumplinger - field-V09-06-06 1027 --------------------------- 1028 - prevent divided by zero and shortcut spin t 1029 particles: G4Mag_SpinEqRhs.cc, G4EqEMFieldW 1030 G4EqEMFieldWithEDM.cc 1031 f 1032 May 8th, 2013 P.Gumplinger/J.Apostolakis - f 1033 ----------------------------------------- 1034 - Introduced new class G4ChargeState to hold 1035 (promoted from internal class of G4Field 1036 - Changed the signature of EquationOfMotion:: 1037 take G4ChargeState in place of G4double (cha 1038 - Deleted SetChargeMomentumMass from G4ChordF 1039 G4MagInt_Driver 1040 1041 - Migrated testPropagateMagField and testProp 1042 1043 March 19th, 2013 G.Cosmo - field-V09-06-03, f 1044 ------------------------- 1045 - Cleanup of MT changes in G4FieldManagerStor 1046 1047 February 26th, 2013 A.Dotti - field-V09-06-01 1048 --------------------------- 1049 - First implementation of Cloning mechanism ( 1050 1051 February 22nd, 2013 G.Cosmo - field-V09-06-00 1052 --------------------------- 1053 - Cleanup of unnecessary changes for multi-th 1054 - Removed unused private data-members in G4Ma 1055 1056 November 14th, 2012 P.Mato - field-V09-05-07 1057 -------------------------- 1058 - G4QuadrupoleMagField: allow displaced origi 1059 - Added G4SystemOfUnits.hh in unit tests 1060 1061 August 30th, 2012 P.Mato - field-V09-05-06 1062 ------------------------ 1063 - Adapted CMakeLists.txt files to new unit te 1064 1065 June 11th, 2012 G.Cosmo - field-V09-05-05 1066 ----------------------- 1067 - Explicitly use inclusion of headers for sys 1068 constants, in plan to remove implicit inclu 1069 1070 May 14th, 2012 G.Cosmo - field-V09-05-04 1071 ---------------------- 1072 - Fixed spurious case of hidden variable visi 1073 detected with '-Wshadow' compilation option 1074 1075 February 3rd, 2012 J.Apostolakis - field-V09 1076 --------------------------------- 1077 - G4FieldTrack: 1078 Added new (inline) method to set the Rest 1079 reset this - and keep reuse the same G 1080 (across tracks.) 1081 Made SetCurvePnt() method private. 1082 Simplified implementation of constructors 1083 1 division per call. 1084 1085 February 3rd, 2012 G.Cosmo - field-V09-05-0 1086 ---------------------------- 1087 - Removed tests in CMakeLists.txt causing cat 1088 testing! 1089 1090 February 2nd, 2012 J.Apostolakis - field-V0 1091 ---------------------------------- 1092 - G4Mag_UsualEqRhs: in SetChargeMomentumMass( 1093 division by zero; stopped updating 'fInvCur 1094 is obsolete). 1095 - Added tests in CMakeLists.txt 1096 1097 December 13th, 2011 G.Cosmo - field-V09-05-0 1098 ---------------------------- 1099 - Fixed Coverity defects in G4SimpleRunge (re 1100 data member) and G4FieldTrack (missing G4Ch 1101 - Synchronised unit tests to trunk (new test 1102 1103 November 11th, 2011 J.Apostolakis - field-V0 1104 ---------------------------------- 1105 - Added G4Field.cc to list of files in source 1106 1107 November 4th, 2011 J.Apostolakis - field-V0 1108 ---------------------------------- 1109 - G4Field: revised to improve ability to dete 1110 Added new data member to identify gravity f 1111 Added new method IsGravityActive() to acces 1112 Created new source file: moving constructor 1113 - Changed G4UniformGravityField and G4Electro 1114 to accomodate this refinement. 1115 1116 November 2nd, 2011 P.Gumplinger - field-V09- 1117 -------------------------------- 1118 - update list of files in sources.cmake 1119 1120 October 28th, 2011 P.Gumplinger - field-V09- 1121 -------------------------------- 1122 - Added (uniform) gravity field capability. 1123 New classes G4EqGravityField, G4UniformGrav 1124 1125 *Required* for co-working tag in Transporta 1126 1127 October 12th, 2011 P.Gumplinger - field-V09- 1128 -------------------------------- 1129 - replaced: omegac = 0.105658387*GeV/mass * 2 1130 with: omegac = (eplus/particleMass)*c_l 1131 in G4EqEMFieldWithEDM.cc, G4Mag_SpinEqRhs.c 1132 changed calculation: 1133 std::sqrt(kineticEnergy*kineticEnergy + 1134 to: 1135 std::sqrt(kineticEnergy*(kineticEnergy+2 1136 in G4FieldTrack.cc 1137 1138 August 11th, 2011 G.Cosmo - field-V09-04-01 1139 -------------------------- 1140 - Migrated to new G4Exception scheme for warn 1141 - Defined scheme for warnings/errors coding: 1142 + 1xxx for warnings, 0xxx for errors. 1143 Code numbering: 1144 + 0001: not supported or implemented featur 1145 + 0002: invalid setup/construct; 1146 + 0003: invalid run condition; 1147 + 1001: incorrect setup; automatically adap 1148 + 1002: notification. 1149 1150 May 13th, 2011 J.Apostolakis 1151 ---------------------------- 1152 - Added new class G4TrialsCounter for keeping 1153 - Avoid printing out pointer in G4CachedMagne 1154 - Updated unit tests. 1155 1156 May 12th, 2011 G.Cosmo - field-V09-04-00 1157 ----------------------- 1158 - Fixed compilation warnings from gcc-4.6.0 f 1159 1160 May 5th, 2011 J.Apostolakis 1161 ---------------------------- 1162 - Reordered data members to group invariants, 1163 1164 Sep 10th, 2010 J.Apostolakis - field-V09-03- 1165 ----------------------------- 1166 - Revised constructor of G4MagErrorStepper to 1167 - Corrected MagErrorStepper to copy State Var 1168 - Enable G4ConstRK4 to copy remaining State V 1169 1170 Sep 10th, 2010 J.Apostolakis - field-V09-03- 1171 ----------------------------- 1172 - New class G4MonopoleEq for rhs of d(pc)/ds 1173 in SI units (V.Grichine). 1174 - Fixed passing of time in G4NystromRK4 1175 1176 Jul 21st, 2010 T.Nikitina - field-V09-03-01 1177 -------------------------- 1178 - Fixed cases of memory corruption in G4RKG3_ 1179 - Fixed case of unused array data member in G 1180 - Removed useless code never executed in G4Co 1181 - Fixed initialization in G4NystromRK4 constr 1182 1183 Jul 14th, 2010 G.Cosmo - field-V09-03-00 1184 ----------------------- 1185 - Added dummy initialization of members in co 1186 G4ConstRK4, G4EqEMFieldWithEDM, G4EqEMField 1187 G4FieldTrack, G4MagHelicalStepper, G4MagInt 1188 G4Mag_SpinEqRhs, G4Mag_UsualEqRhs, G4Nystro 1189 G4UniformElectricField, G4UniformElectricFi 1190 1191 Nov 12th, 2009 J.Apostolakis - field-V09-02- 1192 ----------------------------- 1193 - G4MagIntegratorDriver: activate check on i 1194 1195 Nov 12th, 2009 J.Apostolakis - field-V09-02- 1196 ----------------------------- (fix only) 1197 - G4Nystrom: Corrected interface method getF 1198 (Problem found by gcc 4.3 - i 1199 1200 Nov 6th, 2009 P.Gumplinger - field-V09-02-07 1201 --------------------------- 1202 - bug fix in G4EqEMFieldWithSpin and G4EqEMFi 1203 thanks to Hiromi Iinuma (KEK) see: 1204 http://hypernews.slac.stanford.edu/HyperNew 1205 1206 Nov 5th, 2009 J.Apostolakis - field-V09-02-0 1207 ---------------------------- 1208 - G4MagIntegratorDriver.cc : Enabled call to 1209 - G4NystromRK4.cc : Disabled auxilia 1210 ComputeRightHan 1211 1212 Nov 5th, 2009 J.Apostolakis - field-V09-02-0 1213 ---------------------------- 1214 - Added new virtual method CalculateRightHand 1215 in caching momentum (and field value) b 1216 Default implementation in G4MagIntegrat 1217 - G4Nystrom: New Set/Get method for cache di 1218 Changed private data members in 1219 - G4MagIntegratorDriver: alternative call to 1220 As a result G4NystromRK4 operat 1221 1222 Nov 5th, 2009 J.Apostolakis - field-V09-02-0 1223 ---------------------------- 1224 - G4CachedMagneticField: New Simple class to 1225 Speeds up code when calculation of field 1226 - G4NystromRK4 : New Stepper with Nys 1227 with analytic estimation of integration 1228 Greatly reduces number of field value pe 1229 - Revised testPropagateMagField to use Cached 1230 and to cover G4NystromRK4 and G4ConstRK4 1231 1232 Nov 4th, 2009 P.Gumplinger - field-V09-02-03 1233 --------------------------- 1234 - (minor change) remove comment from G4EqEMFi 1235 add G4EqEMFieldWithEDM class: this is the R 1236 electric and magnetic field, with spin trac 1237 Thanks to Kevin Lynch, Phys. Dept. at Bosto 1238 1239 May 18th, 2009 T.Nikitina - field-V09-02-02 1240 ------------------------- 1241 - Enhanced algorithm G4ChordFinder::ApproxCur 1242 BrentLocator. 1243 1244 March 25th, 2009 G.Cosmo - field-V09-02-01 1245 ------------------------ 1246 - Some code cleanup and formatting... 1247 1248 March 6th, 2009 P.Gumplinger - field-V09-02-0 1249 ---------------------------- 1250 - Added 3rd term of BMT equation (Spin x Beta 1251 addresses emfields forum posting #155 (bug 1252 Phys. Dept. at Boston University. 1253 - Moved renormalization of spin from G4EqEMFi 1254 G4ClassicalRK4 and G4SimpleHeum. 1255 - Added Spin propagation errors to the criter 1256 G4MagIntegratorDriver but actually don't ad 1257 1258 November, 19th, 2008 P.Gumplinger - field-V0 1259 ---------------------------------- 1260 - Renormalized the spin to 1 in G4EqEMFieldWi 1261 1262 November, 7th, 2008 P.Gumplinger - field-V0 1263 ---------------------------------- 1264 - Added SetAnomaly() and GetAnomaly() accesso 1265 1266 October, 28th, 2008 T.Nikitina - field-V09- 1267 -------------------------------- 1268 - Added 'ApproxCurveV' to method G4ChordFinde 1269 order to better calculate the curve length. 1270 - Some code cleanup. 1271 1272 September, 18th, 2008 T.Nikitina 1273 ---------------------------------- 1274 - Added new stepper class G4ConstRK4, perform 1275 with error calculation in constant magnetic 1276 from G4ClassicalRK4. 1277 1278 July, 15th, 2008 J.Apostolakis - field-V09- 1279 -------------------------------- 1280 - Revised signature of G4ChordFinder::FinderN 1281 argument passed a reference. 1282 1283 May, 28th, 2008 T.Nikitina 1284 ---------------------------- 1285 - G4ChordFinder: added new method ApproxCurve 1286 location mechanism. 1287 1288 April 24th, 2008 T.Nikitina - field-V09-01- 1289 ----------------------------- 1290 - Fix 'unused' non-initialized values of dydx 1291 and G4EqEMFieldWithSpin. 1292 1293 January 17th, 2008 G.Cosmo - field-V09-01-0 1294 ---------------------------- 1295 - G4CashKarpRKF45: deleted un-needed duplicat 1296 'fNumberOfVariables', already defined in G4M 1297 - Some code cleanup. 1298 1299 December 7th, 2007 J.Apostolakis & T.Nikiti 1300 --------------------------------------------- 1301 - New class G4FieldManagerStore. It enables r 1302 finders' state. Needed to address problem r 1303 1304 September 9th, 2007 G.Cosmo - field-V09-00-0 1305 ---------------------------- 1306 - Some code cleanup in G4MagHelicalStepper.cc 1307 1308 August 30th, 2007 P.Gumplinger - field-V09 1309 --------------------------------- 1310 - Add G4EqEMFieldWithSpin class to allow trac 1311 G4ElectroMagneticFields 1312 1313 August 21st, 2007 T.Nikitina - field-V09-0 1314 ------------------------------- 1315 - Fix to DistChord() method for Helical Stepp 1316 also for Ang<2pi not only for Ang<pi; use c 1317 improve CPU performance. 1318 - Improved AdvanceHelix() to be able to make 1319 in one call. 1320 - In G4HelixMixedStepper added possibility to 1321 small steps and add some statistics how off 1322 - G4RKG3_Stepper: improved internal algorithm 1323 the stepper now gives comparable CPU perfor 1324 - Updated unit tests. 1325 1326 August 17th, 2007 G.Cosmo - field-V09-00-0 1327 ---------------------------- 1328 - Changed FatalException to EventMustBeAborte 1329 if proposed integration step is negative. 1330 1331 June 8th, 2007 T.Nikitina - field-V08-03-0 1332 ---------------------------- 1333 - Relaxed FatalException to Warning in G4MagI 1334 for case of step zero and return previous v 1335 observed by ATLAS in release 8.3 when requi 1336 -- ie using very small values for acceptabl 1337 1338 May 18th, 2007 T.Nikitina - field-V08-03-0 1339 ---------------------------- 1340 - New class G4HelixMixedStepper, mixing helic 1341 G4MagHelicalStepper 1342 - Keep values of inverse Radius, .. 1343 - New calculation of DistChord for SteppingAn 1344 G4RKG3Stepper 1345 - New calculation of DistChord for SteppingAn 1346 1347 May 16th, 2007 G.Cosmo - field-V08-03-01 1348 ------------------------- 1349 - Introduced new class G4ErrorMag_UsualEqRhs 1350 field when propagation is backwards for err 1351 - Coworks with tag "global-V08-03-02". 1352 1353 May 11th, 2007 G.Cosmo - field-V08-03-00 1354 ------------------------- 1355 - Use call to G4GeometryTolerance instead of 1356 G4MagIntegratorDriver. 1357 - Requires tag "global-V08-03-00" and related 1358 1359 May 10th, 2007 J.Apostolakis - field-V08-0 1360 ------------------------------- 1361 - G4MagIntegratorDriver: 1362 * added protections against zero initial 1363 * revised (re)calculation of h, and quit 1364 1365 April 26th, 2007 T.Nikitina - field-V08-02-0 1366 ---------------------------- 1367 - Corrected integration constants in G4MagHel 1368 - Added implementation of G4RKG3_Stepper, ste 1369 - Updated NTST unit test. New unit test "fiel 1370 1371 November 14th, 2006 J.Apostolakis - field-V0 1372 ---------------------------------- 1373 - G4MagIntegratorDriver 1374 * Made first argument of GetDerivatives a r 1375 - G4MagIntegratorStepper 1376 * Include the SetEquationOfMotion() method 1377 (for Geant4e, as in field-V08-01-00) 1378 * Update of disclaimers. 1379 1380 November 13th, 2006 G.Cosmo - field-V08-01-0 1381 ---------------------------- 1382 - G4FieldTrack: 1383 o Implemented missing definition for GetCha 1384 o Corrected 'const' qualifier for accessors 1385 o Cleanup ... 1386 - NOTE: this tag is NOT including the change 1387 1388 June 21st, 2006 J.Apostolakis - field-V08-01 1389 ------------------------------ 1390 - Including the SetEquationOfMotion() method 1391 1392 June 22nd, 2006 J.Apostolakis - field-V08-00 1393 ----------------------------- 1394 - G4ExactHelixStepper: reintroduced corrected 1395 --> testPropagateMagField now tests it al 1396 - G4MagHelicalStepper: made Stepper virtual ( 1397 to enable G4ExactHelixStepper to 1398 1399 June 21st, 2006 J.Apostolakis - field-V08-00 1400 ----------------------------- 1401 - G4ExactHelixStepper: taking out class 1402 (provisional tag in case of difficulties w 1403 1404 --- The field-V08-00-0x , x > 6, tags do not 1405 1406 May 27th, 2006 J.Apostolakis - field-V08-00- 1407 ----------------------------- 1408 - G4FieldTrack: Fix for simple compilation er 1409 1410 May 23rd, 2006 J.Apostolakis - field-V08-00- 1411 ----------------------------- 1412 - G4FieldTrack: Fix for memory leak for Charg 1413 Made G4ChargeState a data member (has an i 1414 1415 May 18th, 2006 J.Apostolakis - field-V08-00- 1416 ----------------------------- 1417 - G4FieldTrack: Fix for issue with gcc 4.0 ( 1418 1419 May 16th, 2006 J.Apostolakis - field-V08-00- 1420 ----------------------------- 1421 - G4ExactHelixStepper 1422 Taken out alpha version - it was included 1423 1424 May 11th, 2006 J.Apostolakis - field-V08-00- 1425 ----------------------------- 1426 - G4FieldTrack 1427 Significant, provisional change: 1428 * Added charge and moments to FieldTrack 1429 (SetMethods, and pointer to embedde 1430 Technical revisions: 1431 * New constructor, with charge and witho 1432 * New almost-default constructor (taking 1433 Also: 1434 * New constructor to take const referenc 1435 * New implementation method UpdateFourMo 1436 * Streaming operator takes const referen 1437 and outputs momentum magnitude & dir 1438 1439 - G4ExactHelixStepper 1440 Included in error 1441 1442 Apr 6th, 2006 J.Apostolakis - field-V08-00-0 1443 ---------------------------- 1444 - In G4FieldManager 1445 Added missing const GetChordFinder() method 1446 Corrected description (at top of include fil 1447 - It is possible to set the field manager o 1448 1449 1450 Nov 9th, 2005 G.Cosmo - field-V07-01-00 1451 ---------------------- 1452 - Archived source code of field01 unit test, 1453 - Migrated unit test NTST to use <sstream> an 1454 with G4VParticleChange new interface. 1455 1456 Dec 2nd, 2004 G.Cosmo - field-V06-02-02 1457 ---------------------- 1458 - Implemented migration to <cmath> for standa 1459 1460 Dec 2nd, 2004 J.Apostolakis - field-V06-02-0 1461 ---------------------------- 1462 - Improvements to comments in key Field class 1463 - Changes to test/NTST, renaming 'main()' fil 1464 - Several improvements to test files 1465 testProElectroMagField.cc 1466 * updating for new Electric Field inheri 1467 different ChordFinder creation ) 1468 * adding columns to print-outs in spin c 1469 testPropagateMagField.cc 1470 * Added verbosity level 1 to field Propa 1471 testPropagateSpin.cc 1472 * Added more printout of parameters (sta 1473 and G4FieldManager before running each 1474 test/NTST 1475 * Renamed main() file to testNTST.cc and 1476 1477 1478 Oct 20th, 2004 P.Gumplinger - field-V06-02-0 1479 ---------------------------- 1480 - Added Set/GetAnomaly() methods in G4Mag_Spi 1481 1482 Feb 9th, 2004 G.Cosmo - field-V06-00-00 1483 ---------------------- 1484 - G4UniformMagField.cc: added missing impleme 1485 GetConstantFieldValue(). Fixes bug report # 1486 1487 Nov 13th, 2003 J.Apostolakis - field-V05-02- 1488 ----------------------------- 1489 - G4ChordFinder: 1490 o extended the arguments of AdvanceChordLim 1491 adding the current safety sphere (origin, 1492 o fixed set functions which had side effect 1493 - Added new class G4ChordFinderSaf (using abo 1494 o which uses safety to create longer chords 1495 propagation in fewer steps --> performa 1496 - G4FieldManager: 1497 Changed default accuracy parameters to mu 1498 Epsilon Min Default to 5.0e-5 from 1499 Epsilon Max Default to 0.001 from 1500 1501 Nov 7th, 2003 J.Apostolakis - field-V05-02- 1502 ----------------------------- 1503 - G4ChordFinder 1504 o Made this a base class 1505 * to enable specialised implementations 1506 * virtual methods: FindNextChord, PrintSt 1507 * new method AccumulateStatistics to aid 1508 1509 o New default value of 1510 Delta Chord to 0.25 mm (from 1511 1512 - G4FieldManager 1513 o Changed default values of 1514 Delta One Step to 0.010 mm (from 1515 Delta Intersection to 0.001 mm (from 1516 for increased default accuracy 1517 o Use field's DoesFieldChangesEnergy method 1518 1519 Nov 7th, 2003 J.Apostolakis - field-V05-02-0 1520 ---------------------------- 1521 - G4ChordFinder 1522 Created new class members for parameters u 1523 fFirstFraction for first st 1524 fFractionLast fraction of 1525 fFractionNextEstimate confidence i 1526 fMultipleRadius for small ra 1527 New method to Print Statistics. 1528 Verbosity to enable optional printing of s 1529 1530 Nov 7th, 2003 J.Apostolakis - field-V05-02-0 1531 ---------------------------- 1532 - G4ChordFinder 1533 o Added simple statistics (with access meth 1534 - test/NTST 1535 o Added new test, taken from D. Williams' e 1536 modified and maintained by J. Apostolaki 1537 1538 Nov 5th, 2003 J.Apostolakis - field-V05-02-0 1539 ---------------------------- 1540 - G4FieldManager 1541 o Added Get/Set Min/Max imumEpsilonStep 1542 o Added fEpsilonMin fEpsilonMax & default t 1543 >> Now these can be tailored for different 1544 1545 - G4Field G4ElectroMagneticField G4MagneticFi 1546 o Adding the virtual method DoesFieldChan 1547 o So that G4ElectroMagneticField derives fr 1548 and G4MagneticField derives from G4Electr 1549 1550 - G4MagErrorStepper and derived classes 1551 (G4ClassicalRK4 G4SimpleRunge G4ExplicitE 1552 o the constructors take G4EquationOfMotion 1553 as it these applicable to all types of eq 1554 (not just ones for Magnetic field) 1555 1556 - G4EqMagElectricField 1557 o derives from G4EquationOfMotion 1558 1559 - G4ChordFinder 1560 o Moved simple constructor to .cc file (tog 1561 o Changed fDefaultDeltaChord to a constant 1562 1563 Nov 4th, 2003 G.Cosmo - field-V05-02-03 1564 ---------------------- 1565 Undocumented. 1566 1567 Nov 4th, 2003 G.Cosmo - field-V05-02-02 1568 ---------------------- 1569 - G4ChordFinder.cc: 1570 o Fixed compilation problem occuring when G 1571 1572 Oct 31st, 2003 G.Cosmo - field-V05-02-01 1573 ----------------------- 1574 - Code cleanup and removal of redundant heade 1575 where necessary by forward declarations. 1576 1577 Oct 24th, 2003 G.Cosmo - field-V05-02-00 1578 ----------------------- 1579 - Migrated to use new G4Exception scheme. 1580 - Updated outputs of unit tests. 1581 1582 Jun 25th, 2003 J.Apostolakis - field-V05-01- 1583 ----------------------------- 1584 - Suppressed printing of statistics (by defau 1585 - Fix for unused variables warning in ChordFi 1586 1587 Jun 24th, 2003 J.Apostolakis - field-V05-01- 1588 ----------------------------- 1589 - Fix for problems introduced in field-V05-01 1590 o G4ChordFinder: 1591 - Fixed unitialised variable for next ste 1592 - Protected assignment of 'fLastStepEstim 1593 unitialised value. 1594 1595 Jun 21st, 2003 J.Apostolakis - field-V05-01- 1596 ----------------------------- 1597 - G4FieldManager: 1598 o Revised into 'slightly' abstract class wi 1599 ConfigureForTrack( const G4Track *) 1600 to enable the user to configure it with d 1601 depending on energy, particle id, volume, 1602 Currently by default, no such configurati 1603 1604 Jun 21st, 2003 J.Apostolakis - field-V05-01- 1605 ----------------------------- 1606 - Revised G4ChordFinder: 1607 * FindNextChord now provides next Step esti 1608 * AdvanceChordLimited to use 1609 - above extra argument to FindNextChord 1610 - the next Step argument as initial-step 1611 - G4MagIntegratorDriver: 1612 * Added extra argument to Accurate Advance: 1613 hinitial: suggested initial step size 1614 - Note: 1615 * Purpose: performance boost with small/li 1616 including minor revisions of interfaces. 1617 1618 Jun 20th, 2003 J.Apostolakis - field-V05-01- 1619 ----------------------------- 1620 - G4MagIntegratorDriver Quick Advance: 1621 o Took out two unneccessary square roots an 1622 - Notes: 1623 1) Pure magnetic, electric, spin field test 1624 2) Purpose: performance boost with small c 1625 1626 Jun 20th, 2003 J.Apostolakis - field-V05-01- 1627 ----------------------------- 1628 - G4MagIntegratorDriver: 1629 o Corrected neglected conversion to <std> f 1630 - Migrated unit tests to <std>. 1631 1632 Jun 19th, 2003 J.Apostolakis - field-V05-01- 1633 ----------------------------- 1634 - G4ChordFinder: 1635 o Corrected simple constructor to initialis 1636 including LastStepEstimate_Unconstrained. 1637 - G4MagIntegratorDriver: 1638 o Keep simple statistics on number of steps 1639 o Took out an unneccessary square root (in 1640 - NOTE: tests field0X are to be migrated. 1641 Purpose: fixes, small performance boost, si 1642 1643 Apr 2nd, 2003 G. Cosmo - field-V05-00-01 1644 --------------------------- 1645 - Created .cc files for G4MagneticField and G 1646 - Fixed definition of copy constructors in fi 1647 - Moved obsolete method StepWithEst() to priv 1648 G4CashKarpRKF45 and G4ClassicalRK4. 1649 - Cleared strict pedantic warnings. 1650 1651 Feb 5th, 2003 G. Cosmo - field-V05-00-00 1652 --------------------------- 1653 - Removed redundant inclusions of PhysicalCon 1654 Now using directly CLHEP's headers. 1655 1656 Dec 12th, 2002 G. Cosmo - field-V04-01-06 1657 --------------------------- 1658 - Moved unused variable 'nStpPr' to debug sec 1659 method AccurateAdvance(). 1660 1661 Dec 06th, 2002 G. Cosmo - field-V04-01-05 1662 --------------------------- 1663 - Merged changes by Makoto for migration to n 1664 messengers of unit tests. 1665 1666 Nov 29th, 2002 J. Apostolakis - field-V04 1667 --------------------------------- 1668 Further fixes for unitialised variable - in o 1669 G4ExplicitEuler, G4ImplicitEuler, G4ImplicitE 1670 1671 Large memory leak fix for G4ImplicitEuler, sm 1672 Test accuracy requirements for spin adjusted. 1673 1674 Nov 20th, 2002 J. Apostolakis - field-V04 1675 --------------------------------- 1676 Further fixes for unitialised variables: 1677 1678 * G4MagIntegratorStepper 1679 - Created instance variables for number of 1680 1681 * G4MagErrorStepper 1682 - Utilise number of State variables to ini 1683 1684 * G4Mag_SpinEqRhs 1685 - Initialise the values of dydx that we do 1686 1687 Note: Undoing these changes and the ones belo 1688 if undefined values are used in varia 1689 have integrated. 1690 BUT the fix assists us to find truly un 1691 -- if they exist -- using valgrind. 1692 1693 Nov 9th, 2002 J. Apostolakis - field-V04 1694 --------------------------------- 1695 Fixes to address use of uninitialised variabl 1696 (this was done by design, but was neither ide 1697 Changed: 1698 1699 * G4MagIntegratorDriver 1700 - Added instance variables for number of In 1701 - Initialise all State variables of Output 1702 - Only the Integration variables will be in 1703 1704 * G4ClassicalRK4 and G4MagErrorStepper 1705 - Copy the initial time to intermediate and 1706 in order to provide a zero-eth estimate 1707 Note1: Should clean-up this fix, making it 1708 Note2: Should *not* be used for time varyin 1709 with at least 8 integration variab 1710 --------------------------------- 1711 G4ChordFinder.cc 1712 - Correction from zero dChord 1713 which is does not update the Unconstrai 1714 - Added instance variables to store number 1715 and new print method for 'debugging' 1716 1717 G4Mag_UsualEqRhs.cc 1718 - Attempted approximations for speedup *ov 1719 1720 1721 Oct 29th, 2002 J. Apostolakis - test-V04-0 1722 -------------------------------- 1723 - Changed output of tests for small output ch 1724 1725 Oct 29th, 2002 J. Apostolakis - test-V04-0 1726 -------------------------------- 1727 - (Only changed and tagged test subdirectory) 1728 Changed tests to use the new method SetChar 1729 PropagatorInField in order to be ready for 1730 G4PropagatorInField. 1731 (A key problem at end-July was this non-pre 1732 1733 Jul 24th, 2002 G. Cosmo - field-V04-01-01 1734 -------------------------- 1735 - G4FieldManager[.hh.cc]: made fDefault_Delta 1736 fDefault_Delta_Intersection_Val normal cons 1737 resolve initialisation problems. Coworks wi 1738 Fixes problem report #383. 1739 1740 Jul 23rd, 2002 G. Cosmo - field-V04-01-00 1741 -------------------------- 1742 - G4MagIntegratorDriver.cc: cleared warning o 1743 Protected variable 'nStpPr' with G4DEBUG_FI 1744 - Updated reference outputs and some other up 1745 1746 Jun 25th, 2002 J. Apostolakis - field-V04- 1747 -------------------------------- 1748 - G4ChordFinder.cc 1749 In NewStep: fixed case of stepTrial >> ste 1750 (it was treated the same as stepTrial 1751 Added guard for stepTrial==0.0 in FindNextC 1752 - G4MagIntegratorDriver.cc 1753 Added protection against error of 0.0 in Co 1754 1755 - G4EquationOfMotion 1756 Made RightHandSide an inline method 1757 In GetFieldValue signature of Point is doub 1758 Defined G4maximum_number_of_field_component 1759 1760 Jun 7th, 2002 J. Apostolakis - field-V04- 1761 -------------------------------- 1762 - G4Mag_SpinEqRhs.cc 1763 o Fix to momentum component of Rhs. 1764 1765 May 31st, 2002 J. Apostolakis - field-V04- 1766 -------------------------------- 1767 - G4MagHelicalStepper.icc 1768 o Fix to LinearStep (correcting previous 1769 - G4ChordFinder.cc 1770 o Corrective behaviour for very small or 1771 [observed in complex 3-d fields (purgi 1772 approximate equations of 1773 - G4MagIntegratorDriver 1774 o More information in debug printing (tan 1775 1776 May 7th, 2002 J. Apostolakis 1777 -------------------------------- 1778 - G4MagIntegratorDriver - ( on trunk ) 1779 * Fix to update the field track for steps s 1780 minstep. 1781 1782 1783 April 19th, 2002 J. Apostolakis (not inclu 1784 -------------------------------- 1785 - G4MagIntegratorDriver.hh,cc (on trunk) 1786 * Bug fix: last step can now only be trigg 1787 * Debugging methods improved. 1788 1789 May 7th, 2002 J. Apostolakis - field-V04- 1790 -------------------------------- 1791 Tag including only two changes below: 1792 - Putting aside the other fixes/changes (as 1793 G4MagIntegratorDriver, tagging only the fix t 1794 for steps smaller than the chord finder's min 1795 1796 May 7th, 2002 J. Apostolakis 1797 -------------------------------- 1798 - G4MagIntegratorDriver - ** branch ** on 1799 * Fix to update the field track for steps s 1800 minstep. 1801 1802 May 7th, 2002 J. Apostolakis 1803 -------------------------------- 1804 - G4MagIntegratorStepper 1805 * Changed RightHandSide from a virtual meth 1806 method -- it is not and should not be c 1807 1808 Mar 28th, 2002 J. Apostolakis - field-V04-00 1809 ----------------------- 1810 - Fixed LinearStep in G4MagHelicalStepper (ex 1811 that had not been updated for new integra 1812 --> found in testing helical steppers in 1813 in region of negligible field. 1814 1815 - Changed maximum number of steps in driver, 1816 halving from 500/order to 250/order. 1817 1818 Mar 28th, 2002 J. Apostolakis 1819 ----------------------- 1820 - Committed directory "OtherField" with examp 1821 Delphi, Quadropole etc. 1822 1823 - Updated outputs of testPropagateMagField fo 1824 that has 3rd test case. 1825 1826 - Added new runs and outputs for testProElect 1827 several steppers, and modified test scrip 1828 1829 1830 Feb 2nd, 2002 V. Grichine 1831 ------------------------ 1832 Update of test/field02 to utilise new feature 1833 - for Electric Field, to signal that a fiel 1834 - updating of time by integration, for time 1835 1836 Jan 10th, 2002 G. Cosmo - field-V04-00-00 1837 ------------------------ 1838 - Activated G4DEBUG_FIELD in GNUmakefile. 1839 1840 Jan 7th, 2002 J. Apostolakis 1841 ----------------------------- 1842 - G4MagHelicalStepper.cc: renamed local varia 1843 1844 Dec 4th, 2001 G. Cosmo - field-V03-02-08 1845 ----------------------- 1846 - Fixed interface to virtual functions in the 1847 Added 'virtual' qualifier to destructor and 1848 of G4Field which may act as base class for 1849 Modified files: G4MagneticField.hh, G4Unifo 1850 G4UniformMagField.hh. 1851 - G4ChordFinder.cc: the unphysical curve leng 1852 inside #ifdef G4DEBUG_FIELD. Other #ifdef r 1853 - G4LineSection[.hh.cc]: correction for A=B w 1854 is returned (V.Grichine). 1855 1856 Nov 28th, 2001 J. Apostolakis - field-V03-02 1857 ------------------------------ 1858 G4FieldManager 1859 - Added fFieldChangesEnergy data member and g 1860 ( required to allow transportation to know h 1861 enabling the integration of time only for 1862 1863 Nov 28th, 2001 J. Apostolakis - field-V03-02 1864 ------------------------------ 1865 Small fixes only: 1866 1867 G4CashKarpRKF45 1868 - fixed deletion of unassigned pointer/array 1869 1870 G4EqMagElectricField 1871 - bug fixed in the formula for inverse veloci 1872 1873 Nov 21st, 2001 G. Cosmo - field-V03-02-05 1874 ------------------------ 1875 - Protected verbose warnings with G4FIELD_DEB 1876 (V.Grichine). 1877 - Replaced G4DEBUG with G4DEBUG_FIELD in G4Ma 1878 - Updates to unit test field02. 1879 1880 Nov 13th, 2001 G. Cosmo - field-V03-02-04 1881 ------------------------ 1882 - Removed warnings on gcc-2.95.2. 1883 - Included unit tests field01, field02, field 1884 1885 Nov 9th, 2001 J. Apostolakis - field-V03-02- 1886 ----------------------------- 1887 Motivated by the use case of an Electric Fiel 1888 - in order to achieve reasonable accuracy in 1889 - and to enable fields that vary with the gl 1890 we (John Ap. and V. Grichine) changed the fol 1891 1892 - Modified G4Field to allow the time in the G 1893 - as 4th component of the "position" arr 1894 1895 - Modified EquationOfMotion to provide the ti 1896 of G4Field 1897 * as the 4th component of a position / 1898 * in place of only position as a 3-arry 1899 1900 - Added the global time as variable to be int 1901 case of an Electric Field 1902 1903 Support for these modifications is required i 1904 1905 (Better support will require either 1906 - a modification in G4FieldTrack to indicate 1907 integrated or not. 1908 - or additionally making the time calculatio 1909 bring in code to do it in from G4Transpor 1910 1911 1912 Nov 8th, 2001 J. Apostolakis - field-V03-02 1913 ------------------------ 1914 - Introduced several changes in G4MagIntegrat 1915 1916 It now respects its minimum driver step si 1917 * whereas before it did not. 1918 1919 Modified maximum number of steps, lowering 1920 -> it should reduce the time spent integ 1921 in some cases where it is unneccessary. 1922 1923 Protected printing of (nearly?) all warnin 1924 1925 1926 Nov 05th, 2001 G. Cosmo - field-V03-02-01 1927 ------------------------ 1928 - Fixed trivial warnings on Linux-g++. 1929 1930 Aug 27th, 2001 J. Apostolakis - field-V03-02 1931 ------------------------------ 1932 Tagged last change. 1933 This tag is utilised in transport-V03-02-01 ( 1934 1935 Aug 16th, 2001 J. Apostolakis 1936 ------------------------------ 1937 G4ChordFinder: 1938 Added method to erase/reset the state of the 1939 This is unwanted internal state that was car 1940 events as a side effect. 1941 1942 Calling this at the start of each track, it 1943 eliminate a non-repeatability that was possi 1944 ( It has been observed in particular in v 1945 1946 May 23rd, 2001 J. Apostolakis - field-V03-01 1947 ------------------------------ 1948 - Added Exception to G4RKG3Stepper constructo 1949 1950 May 23rd, 2001 J. Apostolakis - field-V03-01 1951 ------------------------------ 1952 - Fixed G4FieldManager constructors to give d 1953 parameters. 1954 1955 - Updated testPropagateSpin (it now compiles) 1956 and some outputs of testPropagateMagField 1957 1958 1959 Apr 12th, 2001 P. Gumplinger - field-V03-01- 1960 ----------------------------- 1961 - Updated G4Mag_SpinEqRhs.hh and G4Mag_SpinEq 1962 to spin tracking. 1963 1964 Apr 04th, 2001 G. Cosmo - field-V03-00-05 1965 ------------------------ 1966 - G4FieldTrack.icc: fixed pedantic warnings o 1967 1968 Mar 23rd, 2001 J. Apostolakis - field-V03-00 1969 ------------------------------ 1970 - Electric Field unit test now works: test ru 1971 Test outputs saved (on Linux). 1972 1973 - Fixed src/G4UniformElectricField.cc to set 1974 (compat 1975 G4EqMa 1976 1977 - Tidied G4FieldTrack 1978 * deleted obsolete methods Set/Get Momentu 1979 * in constructor to set the value of the s 1980 if pSpin is null. (It was undefine 1981 1982 Mar 23rd, 2001 J. Apostolakis - field-V03-00 1983 ------------------------------ 1984 - Fixed (hacked) CashKarp to work using Auxil 1985 1986 - Modified G4MagIntegratorStepper 1987 to use pointer to G4EquationOfMotion 1988 instead of G4Mag_EqRhs 1989 1990 - Design iteration appears neccessary 1991 1992 Mar 23rd, 2001 J. Apostolakis - field-V03-00 1993 ------------------------------ 1994 - Fixed G4MagHelicalStepper::AdvanceHelix for 1995 the possibility of (correct) use of all H 1996 1997 - test/testProElectroMagField.cc updated for 1998 (Not yet updated: Spin tests) 1999 2000 2001 Mar 19th, 2001 J. Apostolakis - field-V03-00 2002 ------------------------------ 2003 - Tagging the modification of several classes 2004 instead of velocity. 2005 Reason for change: Velocity integration ca 2006 velocities, due to arithmetic inaccuraci 2007 Classes modified: 2008 G4FieldTrack, (Changed components 3,4,5 2009 G4Mag_EqRhs, G4Mag_UsualEqRhs, G4EqMagEle 2010 2011 Feb 20th, 2001 J. Apostolakis 2012 2013 - Modified several classes to integrate in mo 2014 2015 G4FieldTrack: 2016 Changed components 3,4,5 of SixVector to 2017 Modified constructor 2018 Added energy calculation to "LoadFromArr 2019 GetEnergy method is now implemented. 2020 Deleted obsolete (old) method names for 2021 2022 G4ChordFinder 2023 Renamed old-method names of Field Track 2024 2025 Jan 29th, 2001 G. Cosmo - field-V03-00-00 2026 ------------------------- 2027 2028 - Minor fix to G4MagneticField.hh: added expl 2029 copy constructor. Fixes report #205. 2030 2031 Nov 20th, 2000 G. Cosmo - field-V02-00-02 2032 2033 - Fixes to remove warnings from "-Wall -ansi 2034 o commented out variables declared and not 2035 o fixed declaration of variables used befor 2036 o fixed order of initialisation of member d 2037 o fixed usage of unsigned-int (size_t) for 2038 2039 Nov 9th, 2000 G. Cosmo - field-V02-00-01 2040 2041 - G4FieldManager: added check on existence of 2042 in CreateChordFinder(G4MagneticField*) meth 2043 2044 Nov 1st, 2000 G. Cosmo - field-V02-00-00 2045 2046 - QA code revision and cleanup (+fixes from C 2047 o Added (private) declarations of copy cons 2048 where needed. 2049 o Added "const" qualifier to accessor metho 2050 o Added equality check on operator= where n 2051 o Made destructor "virtual" for classes hav 2052 wherever missing. 2053 o Changed virtual functions to non-inline w 2054 o Changed return value to be "const" pointe 2055 - G4EquationOfMotion::GetFieldObj() 2056 - G4FieldManager::GetDetectorField() 2057 - G4MagIntegratorDriver::GetStepper() 2058 o Fixed and corrected signature to G4MagInt 2059 o Added "inline" qualifier to methods' decl 2060 o General cosmetics 2061 2062 Jun 2nd, 2000 J. Apostolakis 2063 - G4FieldManager[.hh.icc] 2064 o Added DeltaIntersection and DeltaOneStep 2065 Get/Set methods and SetAccuraciesWithDelt 2066 2067 May 31st, 2000 J. Apostolakis - field-V01-01 2068 2069 include/G4FieldManager.hh,icc 2070 o Moved Delta Intersection and Delta One St 2071 (from Propagator in Field), to allow th 2072 2073 May 11th, 2000 J. Apostolakis - field-V01-01 2074 2075 - G4ChordFinder[.hh.cc] 2076 o Modified method for finding new trial ste 2077 New method assumes a quadratic relationsh 2078 d_chord proportional to step-l 2079 o Added two small refinement in order to av 2080 small numerical inaccuracies in estimatio 2081 FindNextChord 2082 to first step use (1-0.001) 2083 New Step: 2084 to subsequent step use trial_st 2085 2086 May 10th, 2000 G. Cosmo 2087 2088 - G4MagneticField.hh: minor fix to dummy impl 2089 operator. 2090 2091 May 9th, 2000 J. Apostolakis - field-V01 2092 2093 - G4MagErrorStepper: 2094 Modified DistChord to check whether the c 2095 2096 - G4RKG3_Stepper: 2097 Deleted error message from StepNoErr, as 2098 2099 - G4MagIntegratorDriver: 2100 Changed the condition for the debug outpu 2101 2102 2103 April 27th, 2000 G. Cosmo 2104 2105 - Reorganised comments for the Software Refer 2106 - Moved inlined definition to .icc files, whe 2107 Created: G4ChordFinder.icc and G4FieldManag 2108 2109 April 13th, 2000 J. Apostolakis - field-V01 2110 2111 - Modified signature to function DumbStepper( 2112 G4HelixExplicitEuler, G4HelixHeum, G4HelixI 2113 G4HelixSimpleRunge and G4MagHelicalStepper. 2114 - Modified signature and implementation of fu 2115 in G4MagHelicalStepper, to use Field (as Th 2116 - Made Richardson extrapolation optional (it 2117 it offers benefit) in G4MagHelicalStepper. 2118 - Updated unit test testPropagateMagField and 2119 - G4LineSection: fixes to prevent divisions b 2120 (V.Grichine). 2121 2122 January 19th, 2000 G. Cosmo - field-V01-00 2123 2124 Syncronized versioning of files with HEAD aft 2125 (essentially CVS header changes). 2126 2127 December 7th, 1999 J. Apostolakis - field- 2128 2129 For negative or very small curve lengths (h) 2130 (about the distance of endpoints being furthe 2131 as it is meaningless. 2132 2133 November 23rd, 1999 J. Apostolakis - field- 2134 2135 Tagged this sub-category's state in geometry- 2136 2137 2138 July 6th, 1999 J. Apostolakis - field- 2139 2140 These fixes and those in field-01-00-06 and f 2141 from the need to cope with the imprecision of 2142 interaction with paths that are rather straig 2143 2144 G4MagIntegrationDriver: 2145 i) AccurateAdvance now that the chord distan 2146 gth (times 1+epsilon). 2147 2148 ii) OneGoodStep and QuickAdvance considers th 2149 the error of the position - and ensure that b 2150 the "eps" maximum given. 2151 2152 G4ChordFinder: 2153 Having found that the curve_length vs ABdis 2154 the imprecision of the integration: 2155 i) Modified the exception in this case to h 2156 relative disrepancies ( 10 * eps ). 2157 ii) Undertaken a default corrective action i 2158 Note that this disrepancy only happens for ve 2159 ie step lengths much smaller than the curvatu 2160 2161 July 1st, 1999 J. Apostolakis - field- 2162 2163 In order to cope with the anomalous condition 2164 in G4PropagatorInField, we ensure that the f 2165 0 and 1.0. (giving it a default value of 0.5 2166 ( in G4ChordFinder::ApproxCurvePointV ) 2167 2168 June 29th, 1999 J. Apostolakis - field- 2169 2170 I re-instated a check on whether the input 2171 the two input curve points is truly greater t 2172 ( in G4ChordFinder::ApproxCurvePointV ) 2173 2174 The condition that prints out a warning messa 2175 take into account the case of a very small fi 2176 ( in G4MagInt_Driver::AccurateAdvance ) 2177 2178 June ??th, 1999 J. Apostolakis - field- 2179 2180 The condition that prints out a warning messa 2181 an incorrect comparison. 2182 ( in G4MagInt_Driver::AccurateAdvance ) 2183 2184 April 19th, 1999 J. Apostolakis - field- 2185 2186 G4MagErrorStepper: Cosmetic changes to name o 2187 2188 March 4th, 1999 J. Apostolakis - field- 2189 2190 Corrected delete statement of arrays to use ' 2191 in destructors of 2 steppers (G4CashKarpRKF45 2192 2193 February 17th, 1999 J. Apostolakis - field 2194 2195 Moved G4PropagatorInField files to geometry/ 2196 Updated GNUmakefile not to depend on geometr 2197 -> the known circular dependency is remove 2198 2199 2200 February 12th, 1999 J. Apostolakis - field 2201 2202 A tag that contains the first version that 2203 directory). G4PropagateInSpin is still in t 2204 in a branch of the header file was needed. 2205 2206 2207 February 10th, 1999 J. Apostolakis 2208 2209 All the following changes are needed for ener 2210 2211 G4ChordFinder now passes to the constructed d 2212 Because of the need for the stepper to know a 2213 2214 -> Moved theNumberOfVariables from derived c 2215 to base class 2216 2217 As a consequense the classes that derived fro 2218 were changed: G4CashKarpRKF45, G4MagHelicalSt 2219 2220 2221 February 9th, 1999 J. Apostolakis 2222 2223 Equation of Motion: 2224 G4MagEqRhs.hh Made SetChargeMomentumMass vir 2225 can initialise its constants 2226 2227 2228 February 9th, 1999 P. Gumplinger 2229 2230 Added Equation of Motion for Spin in Magnetic 2231 2232 2233 February 17th, 1999 J. Apostolakis - fiel 2234 ---> date is correct, made sure that tag d 2235 2236 Moved G4PropagatorInField files to geometry/ 2237 Updated GNUmakefile not to depend on geometr 2238 -> the known circular dependency is remove 2239 2240 2241 --------------------------------------------- 2242 ------ Geant 4 Production Release - 2243 --------------------------------------------- 2244 Nov 19th, 1998 J. Apostolakis - field-0 2245 2246 G4PropagatorInField files modified: 2247 2248 - to store the last safety computed and it 2249 - to give them to a requester 2250 --------------------------------------------- 2251 2252 Nov 19th, 1998 J. Apostolakis - field-0 2253 2254 Subtle fix to interface of G4UniformMagField 2255 ( that tripped Visual C++ ) 2256 2257 G4PropagatorInField: added printing method 2258 2259 Nov 19th, 1998 J. Apostolakis - field-0 2260 2261 Steppers: Fixed instance variables constr 2262 and clash of names in CashKarp: 2263 2264 Added post-const to all Field classes for 2265 void GetFieldValue(const G4double 2266 G4double B[] 2267 2268 Added SetFieldValue methods to UniformMagF 2269 2270 Added GetConstantFieldValue() methods to U 2271 2272 2273 Nov 19th, 1998 J. Apostolakis - field-0 2274 2275 Fixes to Propagate and Steppers 2276 2277 Nov 19th, 1998 J. Apostolakis - field-0 2278 2279 Minor corrections (???) 2280 2281 Nov 19th, 1998 J. Apostolakis - field-0 2282 2283 Major development to allow Electric Field 2284 2285 2286 2287 June 17th, 1998 J. Apostolakis - magfiel 2288 2289 Modified Files: 2290 include/G4PropagatorInField.hh 2291 src/G4PropagatorInField.cc 2292 test/testPropagateInField.cc 2293 --------------------------------------------- 2294 2295 1) Added an assertion after each call to Lo 2296 It checks that the located volume is equ 2297 If this check fails, the Propagator does 2298 2299 2) Added an argument to ComputeStep method 2300 The new volume argument is used to check th 2301 2302 3) When G4Navigator's ComputeStep returned 2303 used to be treated as a failure to intersect. 2304 2305 stepLength==stepRequested is an intersect 2306 so I am now treating it as such in each 2307 Ignoring it causes problems, including 2308 --> but treating it as an intersection 2309 2310 2311 May 7th, 1998 J. Apostolakis - magfie 2312 2313 - Tag includes new version of G4ChordFinder ( 2314 2315 2316 May 7th, 1998 W. Wander (recorded/c 2317 2318 - New "helical" RK steppers that move in heli 2319 * Several changes to add steppers 2320 2321 Modified Files: 2322 G4MagErrorStepper.hh 2323 G4MagErrorStepper.cc 2324 added Helical stepping methods to 2325 2326 G4MagIntegratorStepper.hh 2327 made RightHandSide virtual (this 2328 2329 testPropagateMagField.cc 2330 Test of Helical stepping functiona 2331 2332 Added Files for new Helical RK steppers, der 2333 G4HelixExplicitEuler.cc ( and .hh ) 2334 G4HelixHeum.cc 2335 G4HelixImplicitEuler.cc 2336 G4HelixSimpleRunge.cc 2337 2338 Added File for implementation: 2339 G4MagErrorStepper.icc 2340 2341 2342 May 7th, 1998 J. Apostolakis - magfi 2343 2344 - added small changes to 2345 * G4ChordFinder: methods to get 2346 * G4MagIntegratorDriver: made OneGoodSt 2347 - created the file.