Skip to Content.
Sympa Menu

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

opal AT lists.psi.ch

Subject: The OPAL Discussion Forum

List archive

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


Chronological Thread  
  • 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: Tue, 4 Apr 2023 10:44:08 +0000
  • Accept-language: en-US
  • Arc-authentication-results: i=1; mx.microsoft.com 1; spf=pass (sender ip is 104.47.22.104) 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=ie0rOi0mN/UqImIm/XISHgFQVw0H63Ey1BzaPZFXh/Y=; b=WMdWCB9zh7G6HZxgpaPAbfNg7HsjRnpLr13dKCZFlUUu0WWzIx9nfYE0Kky44BkjlHz94D1n4Byt62pdPc38yFo6lCQ1r91dMJu0IH4KBc+VgauJRCh5hMgE1I0FTjbTOZOPbxhousd5F8B3IqC5KUUcqhynVpFFth1vHcLyWyOSdmHIP5diqTnmV6mFT67gdnnz+2fECknNxJOxb31tE9n1oXG85nZK7sg2HpOSSUGzqsCFG+1MF/Ox4UegFBTMtxpz7I5yMsHD2zATE4yco2+zIqkXZh2lqS621Ay/BYFKexuQ1iaPKhDQJJJM/iwJYS/r7w4Yyw3KnlL4iXok6Q==
  • Arc-seal: i=1; a=rsa-sha256; s=arcselector9901; d=microsoft.com; cv=none; b=h5kvFs1X5tZOEx7rN2FCsBZBXBXfT/qrrNVaLJtRbYAHQ7qSh0h6iFgEewjXkcz7rTqTDpbNC9X0j5MXAue1+VfMUNY6mw2MSq2KtdlPsThCsCRMAz58ilAd8qF9XQrZgcwsTynUYav8oKt/VCyrVGXAwpm0OW6Eqa6bPiR7Ck7t9TJ3Sw2FQr0Z2KCU38fBB+01MNoxPKvGLTt8gisnl5N3SqsQtjN6Q/mJS3tHVlYWGAjfJ5kI224CwGxQCJToXKEBX0dPeGhoJwfv0qY7175L4FAPxHkj4S0N6HicOkUcsyejV2+XdhjETHowTZmXQYbszbyvnmSlpgh0660ueg==
  • Authentication-results: mc4.ethz.ch; iprev=pass (mail-zr0che01on2041.outbound.protection.outlook.com) smtp.remote-ip=40.107.24.41; 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 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>




Archive powered by MHonArc 2.6.24.

Top of Page