opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

Re: [Opal] quaternion isn't unit quaternion. Norm: -nan

From: Christopher Pierce <cmpierce AT uchicago.edu>
To: Adelmann Andreas <andreas.adelmann AT psi.ch>
Cc: "opal AT lists.psi.ch" <opal AT lists.psi.ch>
Subject: Re: [Opal] quaternion isn't unit quaternion. Norm: -nan
Date: Tue, 11 Apr 2023 20:31:27 +0000
Accept-language: en-US
Arc-authentication-results: i=1; mx.microsoft.com 1; spf=pass smtp.mailfrom=uchicago.edu; dmarc=pass action=none header.from=uchicago.edu; dkim=pass header.d=uchicago.edu; arc=none
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=gW2vdjyRdbJ6nNX6xr/9eFEImL8RBGJJaftbUesi/BA=; b=gYMMKyW6hluK8v2qC7uv1z+D6zIP14IcspJXU8oSjVr3ZSFrsEYtyO5evd9UhX+Czn32vf9i0AcuwJnuoFqYxWXDTo8XmuOkvX6b6Oah4dRAbxdQO9YeEdF8iaguVtpbeL7ODjA2ngPYEgp18e79YESd5rWltlDq/7x3UaMvhv3MqbOl89QF1lDXa/Avf1V39IM1Kc6tTqirb1b5JlNQoake71QLgqNrZ2MzuOVMgKDibxYgnlWiwUaC8bcwEroRwu3JB0wOZgp4W7uYQaXplK9C5xiRpayETkVzKyKEumF6YCiYEZ8gYACFVCAPU9K5fEFxiBnG8OgPLlEA9C3VoQ==
Arc-seal: i=1; a=rsa-sha256; s=arcselector9901; d=microsoft.com; cv=none; b=WAD5fFH5IRxA7/q5XOE7KwmOTqwi52gcvrRqUdpd/iJgJJvWIZhDmtq0AWWkGQb3eyCYmL2HZzHYWRE/tvgckkPZi5N6XU6aZXN0fyjsXsh+y01OUJwGkszr1RXa0Uauuv5y93ylQ6TcoUYnWcwBnn5/lGlLYTpf1UjKMin1DiVdANn3iqgrrcu/TUtz9KzFfJEqPhhaB35iwrs3xw4WEYOHYap/pjrl1Y/071VmEUumuEOEr3IQoBBV3R6z/wQGTzPPkk7Fr0QTKmUAsUJtcOYx8eP2z8hNLHOJGAxum745zbkUG3yEOShtSsbObRVcB5Rh1n22ZMfmisEc4BQWxQ==
Authentication-results: mc2.ethz.ch; iprev=pass (mail-dm6nam10on2095.outbound.protection.outlook.com) smtp.remote-ip=40.107.93.95; spf=pass smtp.mailfrom=uchicago.edu; dkim=pass header.d=uchicagoedu.onmicrosoft.com header.s=selector2-uchicagoedu-onmicrosoft-com header.a=rsa-sha256; dmarc=pass header.from=uchicago.edu
Authentication-results: dkim=none (message not signed) header.d=none;dmarc=none action=none header.from=uchicago.edu;
Msip_labels:

Yep, sorry for the late response, but feel free to include the files I sent you.

Cheers,

Chris

From: Adelmann Andreas <andreas.adelmann AT psi.ch>
Sent: Friday, April 7, 2023 4:28 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: -nan

Chris can I use your field map and input file for the enhancement request, i.e. store the

information in our gitlab?

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 6 Apr 2023, at 21:29, Christopher Pierce <cmpierce AT uchicago.edu> wrote:

Ahh, Ok. Yes, the fieldmap really is a strange one.

Yeah, if it would be difficult to make things work correctly when Ez=0 through the whole file, then some sort of error message mentioning the problem would probably be a good solution. Setting Ez to like 0.00000001 really shouldn't affect the results.

Cheers,

Chris

From: Adelmann Andreas <andreas.adelmann AT psi.ch>
Sent: Thursday, April 6, 2023 9:24 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: -nan

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
============================================

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: -nan

Hi 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 is

running well. I do not see however, the error you encountering. What version of OPAL

are 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