Rodrigues’ rotation formula: Difference between revisions

 

=={{header|Ada}}==

use Ada.Text_Io;

with Ada.Numerics.Elementary_Functions;

Theta => Angle(Source, Target))));

end Rodrigues;

{{out}}

<pre>

=={{header|ALGOL 68}}==

{{Trans|JavaScript}}

MODE VECTOR = [ 3 ]REAL;

MODE MATRIX = [ 3 ]VECTOR;

)

)

{{out}}

<pre>

=={{header|AutoHotkey}}==

{{Trans|JavaScript}}

v2 := [8,5,-30]

a := getAngle(v1, v2)

, [z*x*t - y*sa, z*y*t + x*sa, ca + z*z*t]]

return matrixMultiply(r, p)

{{out}}

<pre>[2.232221, 1.395138, -8.370829]</pre>

=={{header|C}}==

{{trans|JavaScript}}

#include <math.h>

 

printf("[%.13f, %.13f, %.13f]\n", np.x, np.y, np.z);

return 0;

 

{{out}}

{{trans|JavaScript}}

{{works with|Factor|0.99 2021-06-02}}

math.vectors prettyprint sequences sequences.generalizations ;

 

 

{ 5 -6 4 } { 8 5 -30 }

{{out}}

<pre>

=={{header|FreeBASIC}}==

This example rotates the vector [-1, 2, -0.4] around the axis [-1, 2, 1] in increments of 18 degrees.

 

type vector

r = rodrigues( v, k, theta )

print using "##.### [##.### ##.### ##.###]"; theta; r.x; r.y; r.z

{{out}}<pre>

Theta rotated vector

=={{header|Go}}==

{{trans|JavaScript}}

 

import (

np := aRotate(v1, ncp, a)

fmt.Println(np)

 

{{out}}

=={{header|JavaScript}}==

===JavaScript: ES5===

return Math.sqrt(v[0]*v[0] + v[1]*v[1] + v[2]*v[2]);

}

var ncp = normalize(cp);

var np = aRotate(v1, ncp, a);

 

===JavaScript: ES6===

and is not available to all JavaScript interpreters)

 

"use strict";

 

null, 2

);

{{Out}}

<pre>[

of numbers. Some of the functions have been generalized to work with vectors

of arbitrary length.

# v1 and v2 should be vectors of the same length.

def dotProduct(v1; v2): [v1, v2] | transpose | map(.[0] * .[1]) | add;

;

 

{{out}}

<pre>

=={{header|Julia}}==

{{trans|Perl}}

using JSON3

 

np = rodrotate(v1, ncp, a)

JSON3.write(np) # "[2.2322210733082284,1.3951381708176411,-8.370829024905854]"

 

=={{header|Nim}}==

{{trans|Wren}}

Only changed most function names.

 

type

nvp = normalized(vp)

np = v1.rotate(nvp, a)

 

{{out}}

=={{header|Perl}}==

===Task-specific===

use strict;

use Math::Trig; # acos

my $json=JSON->new->canonical;

 

{{out}}

<pre>[2.23222107330823,1.39513817081764,-8.37082902490585]</pre>

===Generalized===

use warnings;

use feature <say signatures>;

 

my($v1,$v2) = ([5, -6, 4], [8, 5, -30]);

{{out}}

<pre>2.23222107330823 1.39513817081764 -8.37082902490585</pre>

=={{header|Phix}}==

{{trans|JavaScript}}

<span style="color: #008080;">with</span> <span style="color: #008080;">javascript_semantics</span>

<span style="color: #008080;">function</span> <span style="color: #000000;">norm</span><span style="color: #0000FF;">(</span><span style="color: #004080;">sequence</span> <span style="color: #000000;">v</span><span style="color: #0000FF;">)</span>

<span style="color: #000000;">np</span> <span style="color: #0000FF;">=</span> <span style="color: #000000;">aRotate</span><span style="color: #0000FF;">(</span><span style="color: #000000;">v1</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">ncp</span><span style="color: #0000FF;">,</span> <span style="color: #000000;">a</span><span style="color: #0000FF;">);</span>

<span style="color: #0000FF;">?</span><span style="color: #000000;">np</span>

{{out}}

<pre>

=={{header|Processing}}==

{{trans|C}}

//Aamrun, 30th June 2022

 

}

 

 

=={{header|Raku}}==

sub norm (@v) { sqrt @v⋅@v }

sub normalize (@v) { @v X/ @v.&norm }

my @v1 = [5,-6, 4];

my @v2 = [8, 5,-30];

{{out}}

<pre>2.232221073308229 1.3951381708176411 -8.370829024905852</pre>

Alternately, differing mostly in style:

 

 

sub infix:<❌> (@v1, @v2) { # cross product

}).join: "\n"

}

{{out}}

<pre>Task example - Point and composite axis / angle:

=={{header|Wren}}==

{{trans|JavaScript}}

 

var normalize = Fn.new { |v|

var ncp = normalize.call(cp)

var np = aRotate.call(v1, ncp, a)

 

{{out}}