Error while working with: Affine, quaternions and cos/sin
Subject: Error while working with: Affine, quaternions and cos/sin
From: Marcel Tünnissen
Date: Sun, 17 Apr 2005 07:23:23 +0000
Thanks!
>dotexptsimp:false;
>dotscrules:true;
>matrix_element_mult:".";
>current_variables:[i,j,k];
>declare_weights(i,1,j,1,k,1);
>order_function:polynomial_monomial_alphalessp;
>/* Contains redundant rules */
>set_up_dot_simplifications([i.i+1,j.j+1,k.k+1,i.j-k,j.i+k,j.k-i,k.j+i,k.i-j,i.k+j]);
>
>declare_scalar_list([cos_alpha,sin_alpha]);
>ya:(sqrt(5)+3)*k/(sqrt(3)*sqrt(5)+sqrt(3))+2*i/(sqrt(3)*sqrt(5)+sqrt(3));
>Sa:(sqrt(2)*sqrt(5)+3*sqrt(2))*k/(sqrt(3)*sqrt(5)+sqrt(3))
> -j+2*sqrt(2)*i/(sqrt(3)*sqrt(5)+sqrt(3));
>dotsimp(expand(Sa . (cos_alpha-ya*sin_alpha)));
With this setup I get the expected result, indeed.
> > FUNCALL: undefined function FERROR
I still get this FERROR if I use usual cos and sin functions:
declare_scalar_list([alpha]);
ya:(sqrt(5)+3)*k/(sqrt(3)*sqrt(5)+sqrt(3))+2*i/(sqrt(3)*sqrt(5)+sqrt(3));
Sa:(sqrt(2)*sqrt(5)+3*sqrt(2))*k/(sqrt(3)*sqrt(5)+sqrt(3))
-j+2*sqrt(2)*i/(sqrt(3)*sqrt(5)+sqrt(3));
dotsimp(expand(Sa . (cos(alpha)-ya*sin(alpha))));
Perhaps it has to do with the fact that load("affine") doesn't work for me
either?
(%i6) load("affine");
0 errors, 0 warnings
Load failed for /usr/local/share/maxima/5.9.1/share/affine/affine.lisp
-- an error. Quitting. To debug this try DEBUGMODE(TRUE);
I followed on the other hand the Affine README file, which states:
COMPILING AFFINE:
(load "sysdef.lisp")
(make:make :affine :compile t)
LOADING AFFINE:
(load "sysdef.lisp")
(make:make :affine)
This way of loading only works for me when I start affine in the affine
directory, otherwise maxima cannot find the sysdef.lisp file:
(%i1) to_lisp();
Type (to-maxima) to restart
MAXIMA> (load "sysdef.lisp")
Maxima encountered a Lisp error:
LOAD: A file with name sysdef.lisp does not exist
I guess I need to setup some kind of path somewhere.
/Teun
_________________________________________________________________
Express yourself instantly with MSN Messenger! Download today it's FREE!
http://messenger.msn.click-url.com/go/onm00200471ave/direct/01/