Skip to content

GitLab

  • Menu
Projects Groups Snippets
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
  • Sign in
  • S src
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
    • Locked Files
  • Issues 59
    • Issues 59
    • List
    • Boards
    • Service Desk
    • Milestones
    • Iterations
  • Merge requests 4
    • Merge requests 4
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Code review
    • Issue
    • Repository
  • Wiki
    • Wiki
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • OPAL
  • src
  • Issues
  • #277

Closed
Open
Created Feb 20, 2019 by frey_m@frey_mMaintainer

OPAL-Cycl - Particle out of bounds

Close to the septum of the PSI Ring cyclotron we lose a lot of particles and particles might end up out of the field map. Thus, we get magnetic field values of \mathcal{O}(10^{177}) which is not physical. This ends up of producing NAN values of the particle positions what causes problems in case of AMR (adaptive mesh refinement) simulations.

I propose to destroy the particles directly after applying the field line 3494 of ParallelCyclotronTracker::computeExternalFields_m by

if ( outOfBound ) {
    itsBunch_m->destroy(1, i, true);
}

With this fix, it works -- that is: a small test worked.

Are there any concerns or better suggestions?

CC: @snuverink_j

Edited Feb 20, 2019 by frey_m
Assignee
Assign to
Time tracking