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