opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

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

From: Adelmann Andreas <andreas.adelmann AT psi.ch>
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)
Date: Wed, 30 Apr 2025 12:28:35 +0000
Accept-language: en-US
Arc-authentication-results: i=1; mx.microsoft.com 1; spf=pass smtp.mailfrom=psi.ch; dmarc=pass action=none header.from=psi.ch; dkim=pass header.d=psi.ch; arc=none
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=pfYeaePoYTGuUFmAK1Z+hjPkMe+EodUUg6VCrWjyTiI=; b=O3GR6MGWbODgbm+UiZryGEph4t7eMIINOyu3fqKTAJ7au+teciTZc53gT58LK0241sUBaAMtx7gazdP5YKEnYTHxXLaE3DLfffFKr6u+rOq5LntCSCCjEjsA4kBbGUSPJuGPHpWq4y6zHa6Zverno44ispNYq6B14HF9xTVoOtOK0iDAI20+WVtDiY5Cyh0Qm7lWIMD6z7qoOMRCj6SKlG9MkV6L8uCmhOBLIq4xR9jd2VTSCWTKc0/azCPZDbdf6G9OJW5AmUH4AxHnrUllcS9Mv0RPsGkjxIXkea+T1bwXIFQJ7KYpT/xbctRF3SFl7rH/d4DPFkPSe3IbQ7bvdQ==
Arc-seal: i=1; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=none; b=AXF8RSUZ0z+l3skK/ogm9a16P7ErYTrTMc37bAKwnzFUgDFdPuJGDdMq6yA7IsUhJahMPixQdDfwF/9s7+UsrL5GaqeepRmbOSL1WiJumdrsvDxmH36wFHP/M0h7afS3+wRmuSGmxnDIDM2suu1axaDeJrrjgLsmCooIlJWw7Sji/yG8pob2HqtIKySvAXHPC7EsS6SCI+jNgPjJKHgUWu6cuLqks1Wq4KPq2uOMqxqiaIaK8ezGUQweXHILI0zzvVF+lA7B5LcyNbbHTLyk3tXW2cSlV1jXBTFvirzmNAFmc3O6RWZSz9qCMEbOCmV4M03dZY2ArJvdEjMWIIPU9Q==
Authentication-results: dkim=none (message not signed) header.d=none;dmarc=none action=none header.from=psi.ch;

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 25 cm 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

Attachment: smime.p7s
Description: S/MIME cryptographic signature

[Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t), Can Davut, 04/30/2025
- Re: [Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t), Adelmann Andreas, 04/30/2025