Skip to Content.
Sympa Menu

opal - Re: [Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t)

opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

Re: [Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t)


Chronological Thread  
  • From: Can Davut <can.davut AT manchester.ac.uk>
  • To: Adelmann Andreas <andreas.adelmann AT psi.ch>
  • Cc: "opal AT lists.psi.ch" <Opal AT lists.psi.ch>
  • Subject: Re: [Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t)
  • Date: Fri, 2 May 2025 14:49:25 +0000
  • Accept-language: tr-TR, en-GB, en-US
  • Arc-authentication-results: i=3; mx.microsoft.com 1; spf=pass (sender ip is 143.55.144.62) smtp.rcpttodomain=lists.psi.ch smtp.mailfrom=manchester.ac.uk; dmarc=pass (p=none sp=none pct=100) action=none header.from=manchester.ac.uk; dkim=pass (signature was verified) header.d=manchester.ac.uk; arc=pass (0 oda=1 ltdi=1 spf=[1,1,smtp.mailfrom=manchester.ac.uk] dkim=[1,1,header.d=manchester.ac.uk] dmarc=[1,1,header.from=manchester.ac.uk])
  • Arc-authentication-results: i=2; mx.microsoft.com 1; spf=pass (sender ip is 143.55.144.62) smtp.rcpttodomain=lists.psi.ch smtp.mailfrom=manchester.ac.uk; dmarc=pass (p=none sp=none pct=100) action=none header.from=manchester.ac.uk; dkim=pass (signature was verified) header.d=manchester.ac.uk; arc=pass (0 oda=1 ltdi=1 spf=[1,1,smtp.mailfrom=manchester.ac.uk] dkim=[1,1,header.d=manchester.ac.uk] dmarc=[1,1,header.from=manchester.ac.uk])
  • Arc-authentication-results: i=1; mx.microsoft.com 1; spf=pass smtp.mailfrom=manchester.ac.uk; dmarc=pass action=none header.from=manchester.ac.uk; dkim=pass header.d=manchester.ac.uk; arc=none
  • Arc-message-signature: i=3; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector10001; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-AntiSpam-MessageData-ChunkCount:X-MS-Exchange-AntiSpam-MessageData-0:X-MS-Exchange-AntiSpam-MessageData-1; bh=Aa5pnH/RPYcsthSdHBHU9VONP0yDtCN7t9bYlDQcsuc=; b=W+jAkWQ3kuN6XnGeOOSct1lEvnffQLtltHwfXfN6y6H3fXJF9aUUEpqov3/SVbNr/8aV3SfGWqKuXcC2zO4nchp1ZE5omsNMQdvAlnrMrKn0SeyBbh8yk5Q3fRdjMVZV8qGeVsa5ycgXVycS1WsSH/xs0lMFrWp+2xROQdHuPLAkGFu9qv25zgq19DZy/xwQrpx2EZSzHVdxHasgxNwafAIW8LilVskgznl3WrGqudX6AcEsVpE9K5HBUPosudnVlNOqSvl2+PYntA8OmHuji8Q0tLm+yFFwkevPITL36F5UhR08wUs/zm48YcO0gte0DbtvKECsVY2wE4JaaeFOpA==
  • Arc-message-signature: i=2; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector10001; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-AntiSpam-MessageData-ChunkCount:X-MS-Exchange-AntiSpam-MessageData-0:X-MS-Exchange-AntiSpam-MessageData-1; bh=Aa5pnH/RPYcsthSdHBHU9VONP0yDtCN7t9bYlDQcsuc=; b=dDLVuRF9h8xeSVyxTSd29yu8+uUbLWivc5CDU7DEVLSsIabvcrUyVQqbKXJF4N+iVFZcTN+8WcbwYMGbzd/UROkKYAAkHrGXTCUqiY57BHhjiep5w7RvxPpdfQ6vQnRVQjJw2kWrhnRIuZbSYqf3woa/mxOYadHEDQ6YTWnwBr9lT0jD4Ig0b6Ov3SDsyoxYsyNdvFzooLtFKy3uHtKPPITpsK3NL30NbLErhBQeWuFGGSDIWbE7hvW7VSGBem8S1fGSv/LZ9U3YIU9bEfz3CC2N0hsiwDBqWim/3S4iccdEvJW616j+JpLRRyfn4ed/Yx3Ij4K9cE8hufDTqS8i/g==
  • Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector10001; h=From:Date:Subject:Message-ID:Content-Type:MIME-Version:X-MS-Exchange-AntiSpam-MessageData-ChunkCount:X-MS-Exchange-AntiSpam-MessageData-0:X-MS-Exchange-AntiSpam-MessageData-1; bh=Aa5pnH/RPYcsthSdHBHU9VONP0yDtCN7t9bYlDQcsuc=; b=vtrnl5hv9Hf5XNrih+vhOjbR27dttXMCd8TpiPJ4VgsqJ/jfzMNvboF7/gGQyEceHukE26Zpd0RtcVIZewPappLp0Wyr2/RcnZEpXKD/JaslWLNhcghlNj3EcDV9wrzbWyExNQ1uvEkCeCrEGzhX8Lzlt13lH8+rV1oOHLQ8oDM0F3A2cJI6CI8UWGhbfd+8X/8BxMRhklG6H4NXyjNhRn5voaighwcRa9N/mzLqI1nhcdzrH072ClfKjGPox5mCUtkSFNbImpUDrax/E3bq9toqvf9oBx4GdR/Tz11eJAuCgCLHVA6OgPlsxlj5bo0+HyBL4GxmFkZee77VbsQYCA==
  • Arc-seal: i=3; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=pass; b=tK24bJbTlri5D6cIzjME9fu8Lgmm2uri0SoQijVPF1LB2dYK8UtF/PDINOdv3SnDMN9oghv2/u219FoFiqQxwJlpx2vXYFffKFKRRuNRwBgLtnjlmT5RusRGaB7ladxhFQBNe3aHPhUX7iJVBnNrulLOthROBbyHPe6faIYLMzRyK5+uw7TSA7i7W+Ggo8NXoFThqAAh0cXp1bDaI0h0lSCX6DrpMyyDIIuIPCSVuLS58bOBaiMEC9Cxs6XBuNV+nti1oCJapj7j+0fAyLi7qUP8JcWxZO26sLi/XveMtK5NLoG3hdgN7Jsy99XcbN7/pe07y6O6Bof/ZdgUNj0f2Q==
  • Arc-seal: i=2; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=pass; b=GnPmZH/132LeZKbYLtB8bc+duAqZVzp2c3hlTTaM8uJq2LOCmxCTnszYVSKLjBzFkeDiYiXBcehhptN3P7XfAGsJ1BuiHk9/fwvm5Ye//UByO+Wxc6bEy9eX8ihcyrpSd/+dlNE/5PeAI2q2hQp5EKbk2IGdm1F+o74Z8RJhNSPVmMofpR+//i81uklDAGoQWCw55BbBRqFBS8VaRj6nkuqjwVVc1IGD/1L/WA9hCClR8FpngLttY79bn15NkyiHG+Pmpl4jsK3rKJQt6g2zCjReHwjp4DZwntyj1sJgaoQVHd8aF3lKoH5c31qlsD38gw0j2GVEp+VwPh187gdr9Q==
  • Arc-seal: i=1; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=none; b=sqwrtWK6dHYyhxMopVkNdFIOZ5EwnrsUzqI+yhzw16ETddjVSSV4b2JvKOgomEAhetTyYXgUInBoyRZWCHL9DpRD6Abx3mtbab7e0ZusPvCHTu7WI85IaUf8kun9Ez+Y6D0f2V9XYvhxrryMr48U4faUugpfe51qIFwuuyLMkUMaL9Y7OKaGSCDSh3J9YOrW+L3M16AXPfeI6Ur5ILdgkNR2NNYqH6nhFW7s9ErblKg50IoWm14JXZeW9QEIL4VeDdWapF158Aky2SPz5RurnVFevFgTDjNmi2WtXbGyRkoUTfHGD1hc0x+qDqpLBRITbHC6SuYw4ZDM2rOki15Evg==
  • Authentication-results: spf=pass (sender IP is 143.55.144.62) smtp.mailfrom=manchester.ac.uk; dkim=pass (signature was verified) header.d=manchester.ac.uk;dmarc=pass action=none header.from=manchester.ac.uk;compauth=pass reason=100

Hi Andy,

I use the previous version (2022.1) as it was the latest when I started using OPAL.

Best Wishes,
Can


From: opal-request AT lists.psi.ch <opal-request AT lists.psi.ch> on behalf of Adelmann Andreas <andreas.adelmann AT psi.ch>
Sent: Friday, May 2, 2025 3:45:16 pm
To: Can Davut <can.davut AT manchester.ac.uk>
Cc: opal AT lists.psi.ch <Opal AT lists.psi.ch>
Subject: Re: [Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t)

Hi Can thanks for the additional information. Will have a look at the problem now.
Last think, I assume you use the latest release of OPAL, correct?



On 2 May 2025, at 05:58, Can Davut <can.davut AT manchester.ac.uk> wrote:

Hi Andy,
 
I initially thought the error was caused by the time step I used for the chicane—since the integration depends on it—so I reduced it to 30 μm, but the error persisted. I also tweaked the filter and wake-function parameters, but without success.
 
Because I use the same field map for all four dipoles and can save the CSR data for the first dipole and the drift without issue, I’ve ruled out the field map as the source of the problem.
 
It seems that increasing the time step before the chicane causes the beam to arrive in under 30 s. If it would help, I can share an HDF5 file so you can use the restart command to diagnose what’s going wrong in the chicane.
 
I’d greatly appreciate any insights you might have.
 
Best regards,
Can
 
Kimden: Adelmann Andreas <andreas.adelmann AT psi.ch>
Tarih: 01 May 2025 17:57
Kime: Can Davut <can.davut AT manchester.ac.uk>
Konu: Re: [Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t)
 
Hi Can, yes I got the reproducer and had a quick look. I can not see an immediate 
error on your side. Can you make the problem smaller i.e. just the chicane? This would help in debugging
the problem. Maybe it is just an unfortunate combination of parameters.
 
We need to solve a nonlinear set of equations and the Newton solver is not converging
within the 100 iterations. Maybe your case needs more than 100 but I do not yet know. 
 
Thanks Andy  


On 1 May 2025, at 11:18, Can Davut <can.davut AT manchester.ac.uk> wrote:
 
Hi Andy,
 
I sent the reproducer in private yesterday. I hope you received it. Just in case I forward it.
 
Thanks in advance for your help.
 
Best Regards,
Can
 
Kimden: Can Davut
Tarih: 30 April 2025 13:35
Kime: Adelmann Andreas <andreas.adelmann AT psi.ch>
Konu: YNT: [Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t)
 
Hi Andy,
 
Please find attached reproducer file including the field maps and input script.
 
Best Regards,
Can
 
Kimden: Adelmann Andreas <andreas.adelmann AT psi.ch>
Tarih: 30 April 2025 13:29
Kime: Can Davut <can.davut AT manchester.ac.uk>
Bilgi: opal AT lists.psi.ch
Konu: Re: [Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t)
 
Hi Can can you send us a reproducer?
 
Thanks Andy 
 

 

On 30 Apr 2025, at 04:13, Can Davut <can.davut AT manchester.ac.uk> wrote:
 
Dear OPAL developers and users,
 
I am experiencing issues with the CSR calculations in my OPAL-t simulations.
 
Following the examples provided in the manual, I set up a symmetric C-shaped chicane using RBEND dipoles and the default field map ("1DPROFILE1-DEFAULT.T7") for 25cm magnets. The bunch compresses by the desired factor of 2, and the tracking completes successfully when CSR effects are not monitored.
 
However, when I attempt to monitor CSR effects, I define a filter as:
CSR_FILTER: FILTER, TYPE = "Savitzky-Golay", NPOINTS = 20, NLEFT = 4, NRIGHT = 4, POLYORDER = 4;
 
and then define the wake function as:
FS_CSR_WAKE: WAKE, TYPE = "1D-CSR-IGF", NBIN = 200, FILTERS = "CSR_FILTER";
 
This wake function is assigned to all dipoles and subsequent drifts. When the simulation runs with this setup, it crashes at the entrance of the second dipole. I also attempted to enable CSRDUMP = TRUE to monitor the raw CSR wake data but encountered the same issue.
The relevant error message is:


<image001.png>


Could you please help me understand what might be causing this issue, so that I can investigate the CSR impact on the bunch more effectively?
 
Thank you very much in advance for your support.
 
Best Regards,
Can
 
Dr. Can Davut | Research Associate
 
Accelerator Physics Group | Department of Physics & Astronomy | The University of Manchester
 
Office : 7.25 Schuster Building / A24B Cockcroft Institute





Archive powered by MHonArc 2.6.24.

Top of Page