opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

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

From: Can Davut <can.davut AT manchester.ac.uk>
To: opal <opal AT lists.psi.ch>
Subject: [Opal] CSR Calculation Crash in Symmetric C-Chicane Simulation (OPAL-t)
Date: Wed, 30 Apr 2025 09:13:47 +0000
Accept-language: tr-TR, en-GB, en-US
Arc-authentication-results: i=3; mx.microsoft.com 1; spf=pass (sender ip is 205.220.185.114) 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 205.220.185.114) 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=oY6lyl9s7kuLhgODOtkiJwU6GbS+oxPZviyMk86I/NU=; b=J6LI5gsOxXDiarzB5d+Z/nztCWz4UNAS5uwQmKXYyjf5alvZng+bnC0RIkjFZUPVOAAUmvpLU6/ubfoaIZbU7tZyWafEbyTXCa4mtpIv6f1DP1SPm0VafqNT5cy55iZz5dcOLafLTpKzhecAwCk9ZVM77b7ecX6IntYsgv+eIMIJ9Wn5jbXwzS3wob/IvimcByT2Bm2aujzscTx1EGhXFZ9aUpo+Hxsw+mHS1VFzgTlLXpq20GV0J/YoYZw2OGPXfmOuh1lH6WCUWM2XUMC0yvnVKZEINpZTLNfQ8IFyxSm9R6ULmt3YtWjTs4iEW4XMbC996/VHGEJUKQWQH1xV/A==
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=oY6lyl9s7kuLhgODOtkiJwU6GbS+oxPZviyMk86I/NU=; b=kO0ti+jj6dQz2tgChFls8SMX7NRpfvMBHtH/NjllDeMmIPmSZvM/8VAiK8A5Pg85m2M8+4pP/thpQJAgShPf5iv1lW/jy1YjyomYqSHBA4dtlc4ngq7+ccjd/VymR+v3bUKHiu6/2uDHEkaQQ4TSeG/zqxdrfDLTeLhA1W+21aMHrSByU5UH5x4athTdFds+G6CzbAkNSQ6VUSAlBbVtfAWnDKVkHA3DuLPad96aR04AVj9QKQQ+dvRqx2NVX6QWyK5SLIRVmkUT42YUrvY7k7/cEf+i9QCCR1Mha3Yp+yf+FDftiKMlO3CpeqvSRjvmMHcgl2RHgMEVK+iOBsQtkQ==
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=oY6lyl9s7kuLhgODOtkiJwU6GbS+oxPZviyMk86I/NU=; b=ZIJ/26NgSrHHxy7zAclNVxP+r3ygTsE25AfdLNid8MddCh6XvWDOZ3WrdCSKxQxZKdQT4/amVB6gRWjL+vbXQVCL4jvoCo3RUEfzmrQG8WnyHNQVAEKnRe72IPpNx0xBHL+Ng516q8Cm0+DEO8YWSbL6fds0GaRzqsG4ViwVPWO3mXNiXocAX8RAD01X9dV55W4ySOwERRH7Dvtmc64WB/R4yENu9JF5MCG5RvL+eHRRGS3qK9TPh5rTdri2DS4OJEaHjy+koGk2gBib2afRMYaHNuNy8eTOeMC1yPk3NedZmM1lMToRvjXyNmCqjo6esrbaqrKnvHndeWgwZ+7ONw==
Arc-seal: i=3; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=pass; b=WEbmsu8O254Au6FqwCuAzXpYLoDuaPgWg+x1WzE3bXSxwppSoxEsHM+SF54TDm9azx/XpbZ6eN83/BRGuE9jRYBCpncJlNoU9E7o6gbTXzWwyU/KPMxFWs6O3TVA16uQ/ko2DCxsKegyvUXIYkYxTIhBo8VALByUPDoWpqi4U+uugMduuvotMIZMDJ7ATL8oKrBSot4yxYmNEVmNYK4Ql8XyK+roIdqPLUj68ai9izEaZQXchHbJqK4pJzvKKK5L7sQAuA5wUpueHc3ZTvK1lH1dxS8c+CZRH9Lh2RByXJrzngtaczltBBiKklxJbFXF2J/xi2+i728EyXwx2UBAcg==
Arc-seal: i=2; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=pass; b=t8AuEHCweUFflmBFBb7ovRxuSfXExRj93w8wpL6w9kMYAMRpXrJZZv9MKlMGHdlMeXEFw931NWGwvy34+jo4wxgpwqvZnqHOoYlTnfMzyuE94KYdG0rCtQhN8ockdQ8BefOpmRFrBK92Rb4en9Z2ztLimTcBVAYU8LztBcnPtE8aa4uOcFHnmj6ufAOd5HT+P27q0a78z2dYD4BUeVkWRxg3CRoqDryotl+sRNTGOPYCTk6RQ4UprL5A9PHxehJ+lHBtcfUCMqv05buA2pTZRrBXpVAl/ml/QAbuVvVV9YlksO98vkLyJU3fFF4fD5rWTOpn5xwYt9/QDZ+Eb+D5cg==
Arc-seal: i=1; a=rsa-sha256; s=arcselector10001; d=microsoft.com; cv=none; b=w28MTra9OUEJzWtkmJA/i1JzxvFa0BuFHGiB5tTaobnQbe5rI0fiEN4oy0hMeZwUaiUaCfTEtThgP5CCCxrGfEqpSD7pCMJQYDE7PW/neCcWK7UrWMi+WoBvNUobhzrq9Bl039GRNfNqhlpNc0zFx1Nax63T+RkeeG0TWYvzWYK/Dbt2pB+d7k3AyeXW9mf1KnmV48KlRVveQaca4ruUPof7lO3khPE5HjlqTe77HzAif44Rs+avtFSxePNOYvoTROF3Wkol94x7g7Lkc574zAErLpxWSNJNFUK6gCM9uxzKPHKIRoIJC7yqPdKCdyr0bFaRytmXJC7JkdW0YicJoA==
Authentication-results: spf=pass (sender IP is 205.220.185.114) 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

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:

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

[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