opal AT lists.psi.ch
Subject: The OPAL Discussion Forum
List archive
- From: Adelmann Andreas <andreas.adelmann AT psi.ch>
- To: Christopher Pierce <cmpierce AT uchicago.edu>
- Cc: "opal AT lists.psi.ch" <opal AT lists.psi.ch>
- Subject: Re: [Opal] quaternion isn't unit quaternion. Norm: -nan
- Date: Thu, 6 Apr 2023 14:24:46 +0000
- Accept-language: en-US
- Arc-authentication-results: i=1; mx.microsoft.com 1; spf=pass (sender ip is 104.47.22.49) smtp.rcpttodomain=lists.psi.ch smtp.mailfrom=psi.ch; dmarc=pass (p=none sp=none pct=100) action=none header.from=psi.ch; dkim=none (message not signed); arc=none (0)
- Arc-message-signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=microsoft.com; s=arcselector9901; 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=AupcgpEkvJFOaUEmW6Yi/wp/oVo9GWKqcg0ZyLx/Bho=; b=hYo2vUKTfh1h8npkaXg8/UtA85BOpQs8Oe1BOr/lS6h4Mpm/k5YTf5ih5139hbYkUpiX5s3/AfKyIkZNtSIFlEJ8lLbG9Kht35ABrdS4HzQGKSMso6aAdwmpAoWS0vCTYLTVXD32NltjsoTcwkYMtnMP1/8Prg7nX+KSZX/jXOd4nzbW1FKXWX44BKh5vOR2ZAYH7tpEHHJCA6y6aCcVuzwufxjTUi75U1CJNvojXE6nsTi2ylHdV0XDkxnD+N4nrBpekuvfh9RWGU6x8mang3V75evzphaYiUP2oYJp1OCoGWBNMb6vV7Rs5DO56G0Kq6WWD2BfTCTYXrZIa73erQ==
- Arc-seal: i=1; a=rsa-sha256; s=arcselector9901; d=microsoft.com; cv=none; b=AyIJWZVYgOQPocUHZDhhz5f9c7El2/I/gKxaPIoVOiHHCH4WTvaLzrXjF5Pxl4gZNcxAlOlfRR7MEkvzsNPUGCX75pg9WsLYtG6K4DbAxu+pJ+pSjPPPE7uAcvYiXtBAYQyJH5bDtNf8S8s9c653TaE8Aq+8H4qUZXTtLUcsl94yuSSBHZDAS45GDLjJmfIqeZKP91okQS+UyH6B2CA7DNpPzvSEE/ehOBTm0vGblMOFGmOhLB+G4lFZZcybPHEHoEnda9GRRpn1Bueyz+p3KAuBqB3f5SmAdPTfrAfadeeMqwzmEdkGXKoxkSgJdyblAIx8SrOvbfZtHI4OAMvKhg==
- Authentication-results: mc4.ethz.ch; iprev=pass (mail-gv0che01on2060.outbound.protection.outlook.com) smtp.remote-ip=40.107.23.60; spf=pass smtp.mailfrom=psi.ch; dkim=pass header.d=psi.ch header.s=selector2 header.a=rsa-sha256; dmarc=pass header.from=psi.ch
- Authentication-results-original: dkim=none (message not signed) header.d=none;dmarc=none action=none header.from=psi.ch;
Hi Chris, with the latest developer version I get the same error.
All field components zero, I view as a pathological case and also understand
why the error (exception) occurs. The way you “fixed” this is to me a good
path forwards.
In case you disagree, I can open an issue, otherwise I would try to create
a better message to the use, in such a situation.
What do you think?
Cheers A
------
Dr. sc. math. Andreas (Andy) Adelmann
Head a.i. Labor for Scientific Computing and Modelling
Paul Scherrer Institut OHSA/D09 CH-5232 Villigen PSI
Phone Office: xx41 56 310 42 33 Fax: xx41 56 310 31 91
Zoom ID: 470-582-4086 Password: AdA
Zoom Link: https://ethz.zoom.us/j/4705824086?pwd=dFcvT1pMMGY0bHg0dTNncUNZZTJkZz09
-------------------------------------------------------
Friday: ETH HPK G 28 +41 44 633 3076
============================================
The more exotic, the more abstract the knowledge,
the more profound will be its consequences.
Leon Lederman
============================================
------
Dr. sc. math. Andreas (Andy) Adelmann
Head a.i. Labor for Scientific Computing and Modelling
Paul Scherrer Institut OHSA/D09 CH-5232 Villigen PSI
Phone Office: xx41 56 310 42 33 Fax: xx41 56 310 31 91
Zoom ID: 470-582-4086 Password: AdA
Zoom Link: https://ethz.zoom.us/j/4705824086?pwd=dFcvT1pMMGY0bHg0dTNncUNZZTJkZz09
-------------------------------------------------------
Friday: ETH HPK G 28 +41 44 633 3076
============================================
The more exotic, the more abstract the knowledge,
the more profound will be its consequences.
Leon Lederman
============================================
On 4 Apr 2023, at 19:58, Christopher Pierce <cmpierce AT uchicago.edu> wrote:
Huh, interesting. I am using version 2022.1.0 (git rev. a4e0e14466557289b4f7a843fb9eb72700ff1992). Is that the development version that doesn't show the error?
Cheers,
Chris
From: Adelmann Andreas <andreas.adelmann AT psi.ch>
Sent: Tuesday, April 4, 2023 5:44 AM
To: Christopher Pierce <cmpierce AT uchicago.edu>
Cc: opal AT lists.psi.ch <opal AT lists.psi.ch>
Subject: Re: [Opal] quaternion isn't unit quaternion. Norm: -nanHi Chris great that OPAL is of help to your research!
And yes your original field map is special :), indeed.
I can reproduce and confirm that with a small Ez component OPAL isrunning well. I do not see however, the error you encountering. What version of OPALare you using?
Cheers A
------
Dr. sc. math. Andreas (Andy) Adelmann
Head a.i. Labor for Scientific Computing and Modelling
Paul Scherrer Institut OHSA/D09 CH-5232 Villigen PSI
Phone Office: xx41 56 310 42 33 Fax: xx41 56 310 31 91
Zoom ID: 470-582-4086 Password: AdA
Zoom Link: https://ethz.zoom.us/j/4705824086?pwd=dFcvT1pMMGY0bHg0dTNncUNZZTJkZz09
-------------------------------------------------------
Friday: ETH HPK G 28 +41 44 633 3076
============================================
The more exotic, the more abstract the knowledge,
the more profound will be its consequences.
Leon Lederman
============================================
On 4 Apr 2023, at 05:53, Christopher Pierce <cmpierce AT uchicago.edu> wrote:
Hey folks,
First, I have recently started using opal and wanted to thank you for it. It looks like it will be very useful for my research.
I am emailing because I am running into a problem with the attached input file. It is supposed to represent a beam in a artificial electrostatic focusing channel. I do this by loading a fieldmap into the simulation with a constant negative radially pointing E-field and no Ez component.
The problem I encounter is that after loading the input file, the program returns the error "quaternion isn't unit quaternion. Norm: -nan".
I found out that I can make the error go away by setting a single element of the fieldmap to have finite Ez (also attaching the file as "e-lens-fixed.T7". This fieldmap seems to give me the results I expect. Is there something that I am doing wrong here?
Cheers,
Chris Pierce<20230403-opal-focusing-channel.txt><e-lens.T7><e-lens-fixed.T7>
-
[Opal] quaternion isn't unit quaternion. Norm: -nan,
Christopher Pierce, 04/04/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Adelmann Andreas, 04/04/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Christopher Pierce, 04/04/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Adelmann Andreas, 04/06/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Christopher Pierce, 04/06/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Adelmann Andreas, 04/07/2023
- Re: [Opal] quaternion isn't unit quaternion. Norm: -nan, Christopher Pierce, 04/11/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Adelmann Andreas, 04/07/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Christopher Pierce, 04/06/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Adelmann Andreas, 04/06/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Christopher Pierce, 04/04/2023
-
Re: [Opal] quaternion isn't unit quaternion. Norm: -nan,
Adelmann Andreas, 04/04/2023
Archive powered by MHonArc 2.6.24.