diff options
Diffstat (limited to '2.3-1/tests/unit_tests/Moebius.dia.ref')
| -rw-r--r-- | 2.3-1/tests/unit_tests/Moebius.dia.ref | 3200 |
1 files changed, 3200 insertions, 0 deletions
diff --git a/2.3-1/tests/unit_tests/Moebius.dia.ref b/2.3-1/tests/unit_tests/Moebius.dia.ref new file mode 100644 index 00000000..d331e043 --- /dev/null +++ b/2.3-1/tests/unit_tests/Moebius.dia.ref @@ -0,0 +1,3200 @@ +// +// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab +// Copyright (C) 2010-2010 - DIGITEO - Bruno JOFRET +// +// This file must be used under the terms of the CeCILL. +// This source file is licensed as described in the file COPYING, which +// you should have received as part of this distribution. The terms +// are also available at +// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt +// +// +// <-- ENGLISH IMPOSED --> +// <-- JVM NOT MANDATORY --> +global SCI2CHOME; +scilab2c(SCI2CHOME+"/tests/unit_tests/test305_Moebius/Moebius.sci",TMPDIR, SCI2CHOME+"/tests/unit_tests/test305_Moebius"); + + ----------------------------------- + + --- Executing your SCILAB code. --- + + ----------------------------------- + + x(1:10) + + 19.6 + 18.645391 + 15.872572 + 11.546217 + 6.0803399 + 1.207D-15 + - 6.1076856 + - 11.648965 + - 16.08034 + - 18.961619 + + y(1:10) + + 0. + 6.0582549 + 11.532098 + 15.892004 + 18.713362 + 19.717157 + 18.797524 + 16.033425 + 11.683051 + 6.1610035 + + z(1:10) + + 0. + - 0.0625738 + - 0.1236068 + - 0.1815962 + - 0.2351141 + - 0.2828427 + - 0.3236068 + - 0.3564026 + - 0.3804226 + - 0.3950753 + + x1(1:10) + + 2.4142136 + 2.4063283 + 2.3827967 + 2.34399 + 2.2905202 + 2.2232306 + 2.1431822 + 2.0516376 + 1.9500404 + 1.8399929 + + y1(1:10) + + 0. + 0. + 0. + 0. + 0. + 0. + 0. + 0. + 0. + 0. + + z1(1:10) + + 0. + 0.1243449 + 0.2408768 + 0.3422736 + 0.4221640 + 0.4755283 + 0.4990134 + 0.4911436 + 0.4524135 + 0.3852566 + + ------------------------------------------ + + --- End Execution of your SCILAB code. --- + + ------------------------------------------ + + ==> SCI2C hArtes/POLIBA Tool!!! + + Removing directory: TMPDIR/SCI2CTmpResultsReports + + Removing directory: TMPDIR + + ==> Initialize SCI2C and USER2C Libraries. + + ==> Load SCI2C and USER2C Libraries. + + ==> Start translation of function "Moebius" +Warning: file 'TMPDIR/SCI2CTmpResultsReports/Moebius/Moebius_copy.sci' already opened in Scilab. +Warning: file 'TMPDIR/SCI2CTmpResultsReports/Moebius/Moebius_copy.sci' already opened in Scilab. + + ==> Generate the AST. + + ==> Generate C code in TMPDIR/main.c + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: //the program construct moebius + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // and klein bottle surface from + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: //http://www.metlin.org/mathematics/ + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // topological-constructs/non-orientable-surfaces/ + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // --------------------- + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // --- Moebius strip --- + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // --------------------- + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: num_points = 20; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: R = 20; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: a = -0.4; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: b = 0.4; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: step=(b-a)./(num_points-1); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: step1=(2*%pi)./num_points; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: s=a:step:b; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: t=0:step1:2*%pi; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: [S, T] = XYgrid(s, t); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: x = [R + (S.*cos(T./2))].*cos(T); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: y = [R + (S.*cos(T./2))].*sin(T); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: z = S.*sin(T./2); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp('x(1:10)') + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp(x(1:10)) + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp('y(1:10)') + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp(y(1:10)) + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp('z(1:10)') + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp(z(1:10)) + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // ENABLE/DISABLE PLOT + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // scf(0); surf(x, y, z); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // ----------------------- + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // --- Klein Bottle II --- + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // ----------------------- + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: num_points = 50; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: //u = linspace (0, 2*(%pi), num_points); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: //v = linspace (0, 2*(%pi), num_points); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: step_u = (2*%pi)./num_points + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: u = 0:step_u:2*%pi; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: v = u; + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: [U,V] = XYgrid(u,v); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: x1 = cos(U).*(cos(U./2).*(sqrt(2)+cos(V))+(sin(U./2).*sin( + V).*cos(V))); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: y1 = sin(U).*(cos(U./2).*(sqrt(2)+cos(V))+(sin(U./2).*sin( + V).*cos(V))); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: z1 = -1.*sin(U./2).*(sqrt(2)+cos(V))+cos(U./2).*sin(V).*co + s(V); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp('x1(1:10)'); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp(x1(1:10)); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp('y1(1:10)') + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp(y1(1:10)) + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp('z1(1:10)') + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: disp(z1(1:10)) + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // ENABLE/DISABLE PLOT + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: // scf(1) ;surf(x1, y1, z1) + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: endfunction + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + ==> Start translation of function "XYgrid" +Warning: file 'TMPDIR/SCI2CTmpResultsReports/XYgrid/XYgrid_copy.sci' already opened in Scilab. +Warning: file 'TMPDIR/SCI2CTmpResultsReports/XYgrid/XYgrid_copy.sci' already opened in Scilab. + + ==> Generate the AST. + + ==> Generate C code in TMPDIR/d2d2XYgridd2d2.c + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: n = length(x); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: m = length(y); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: X=zeros(m,n); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: Y=zeros(m,n); + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: for k=1:m + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ####################################################### + ########### + + SCI2C: X(k,1:n)=x; + + SCI2C: ####################################################### + ###########*/ + + /*SCI2C: ####################################################### + ########### + + SCI2C: end + + SCI2C: ####################################################### + ###########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: for h=1:n + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ####################################################### + ########### + + SCI2C: Y(1:m,h)=y'; + + SCI2C: ####################################################### + ###########*/ + + /*SCI2C: ####################################################### + ########### + + SCI2C: end + + SCI2C: ####################################################### + ###########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: endfunction + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + /*SCI2C: ########################################################## + ######## + + SCI2C: + + SCI2C: ########################################################## + ########*/ + + ==> Copying sources + + ==> Copying headers + + ==> Copying interfaces + + ==> Generating Builder TMPDIR/Makefile + + ==> Translation Successfully Completed!!! +fd = mopen(TMPDIR+"/main.c"); +mgetl(fd) + ans = + +!/* ! +! ! +!** ************************************************ ! +! ! +!** This file has been generated using ! +! ! +!** Scilab2C (Version 2.0) ! +! ! +!** ! +! ! +!** Please visit following links for more informations: ! +! ! +!** Atoms Module: http://atoms.scilab.org/toolboxes/scilab2c ! +! ! +!** Scilab2C Forge: http://forge.scilab.org/index.php/p/scilab2c/ ! +! ! +!** Scilab2C ML: http://forge.scilab.org/index.php/p/scilab2c/ ! +! ! +!** ************************************************ ! +! ! +!*/ ! +! ! +! ! +! ! +! ! +! ! +!/* ! +! ! +!** ----------------- ! +! ! +!** --- Includes. --- ! +! ! +!** ----------------- ! +! ! +!*/ ! +! ! +!#include "main.h" ! +! ! +!/* ! +! ! +!** --------------------- ! +! ! +!** --- End Includes. --- ! +! ! +!** --------------------- ! +! ! +!*/ ! +! ! +! ! +! ! +! ! +! ! +! ! +! ! +!/* ! +! ! +!** ------------------------------------- ! +! ! +!** --- Global Variables Declaration. --- ! +! ! +!** ------------------------------------- ! +! ! +!*/ ! +! ! +! ! +! ! +! ! +! ! +!/* ! +! ! +!** ----------------------------------------- ! +! ! +!** --- End Global Variables Declaration. --- ! +! ! +!** ----------------------------------------- ! +! ! +!*/ ! +! ! +! ! +! ! +!/* ! +! ! +! SCI2C: -------------------------------------------------------------! +! ----- ! +! ! +! SCI2C: //SCI2C: DEFAULT_PRECISION= DOUBLE ! +! ! +! SCI2C: ! +! ! +! SCI2C: function Moebius() ! +! ! +! SCI2C: -------------------------------------------------------------! +! ----- ! +! ! +!*/ ! +! ! +!int main() ! +! ! +!{ ! +! ! +!/* ! +! ! +!** ----------------------------- ! +! ! +!** --- Variable Declaration. --- ! +! ! +!** ----------------------------- ! +! ! +!*/ ! +! ! +! ! +! ! +! double num_points = 20; ! +! ! +! ! +! ! +! double R = 20; ! +! ! +! ! +! ! +! double a = -0.4; ! +! ! +! ! +! ! +! double b = 0.4; ! +! ! +! ! +! ! +! double __temp1; ! +! ! +! ! +! ! +! double __temp2; ! +! ! +! ! +! ! +! double step; ! +! ! +! ! +! ! +! double __temp3; ! +! ! +! ! +! ! +! double step1; ! +! ! +! ! +! ! +! double* s = NULL; ! +! ! +! int __sSize[2]; ! +! ! +! ! +! ! +! double __temp4; ! +! ! +! ! +! ! +! double* t = NULL; ! +! ! +! int __tSize[2]; ! +! ! +! ! +! ! +! double* S = NULL; ! +! ! +! int __SSize[2]; ! +! ! +! ! +! ! +! double* T = NULL; ! +! ! +! int __TSize[2]; ! +! ! +! ! +! ! +! double* __temp5 = NULL; ! +! ! +! int ____temp5Size[2]; ! +! ! +! ! +! ! +! double* __temp6 = NULL; ! +! ! +! int ____temp6Size[2]; ! +! ! +! ! +! ! +! double* __temp7 = NULL; ! +! ! +! int ____temp7Size[2]; ! +! ! +! ! +! ! +! double* __temp8 = NULL; ! +! ! +! int ____temp8Size[2]; ! +! ! +! ! +! ! +! double* __temp9 = NULL; ! +! ! +! int ____temp9Size[2]; ! +! ! +! ! +! ! +! double* x = NULL; ! +! ! +! int __xSize[2]; ! +! ! +! ! +! ! +! double* __temp10 = NULL; ! +! ! +! int ____temp10Size[2]; ! +! ! +! ! +! ! +! double* __temp11 = NULL; ! +! ! +! int ____temp11Size[2]; ! +! ! +! ! +! ! +! double* __temp12 = NULL; ! +! ! +! int ____temp12Size[2]; ! +! ! +! ! +! ! +! double* __temp13 = NULL; ! +! ! +! int ____temp13Size[2]; ! +! ! +! ! +! ! +! double* __temp14 = NULL; ! +! ! +! int ____temp14Size[2]; ! +! ! +! ! +! ! +! double* y = NULL; ! +! ! +! int __ySize[2]; ! +! ! +! ! +! ! +! double* __temp15 = NULL; ! +! ! +! int ____temp15Size[2]; ! +! ! +! ! +! ! +! double* __temp16 = NULL; ! +! ! +! int ____temp16Size[2]; ! +! ! +! ! +! ! +! double* z = NULL; ! +! ! +! int __zSize[2]; ! +! ! +! ! +! ! +! char __temp17[1*8] = {"x(1:10)"}; ! +! ! +! int ____temp17Size[2] = {1,8}; ! +! ! +! ! +! ! +! double ans1; ! +! ! +! ! +! ! +! double __temp18[1 * 10]; ! +! ! +! int ____temp18Size[2] = {1, 10};; ! +! ! +! ! +! ! +! double __temp19[10 * 1]; ! +! ! +! int ____temp19Size[2] = {10, 1};; ! +! ! +! ! +! ! +! double ans2; ! +! ! +! ! +! ! +! char __temp20[1*8] = {"y(1:10)"}; ! +! ! +! int ____temp20Size[2] = {1,8}; ! +! ! +! ! +! ! +! double ans3; ! +! ! +! ! +! ! +! double __temp21[1 * 10]; ! +! ! +! int ____temp21Size[2] = {1, 10};; ! +! ! +! ! +! ! +! double __temp22[10 * 1]; ! +! ! +! int ____temp22Size[2] = {10, 1};; ! +! ! +! ! +! ! +! double ans4; ! +! ! +! ! +! ! +! char __temp23[1*8] = {"z(1:10)"}; ! +! ! +! int ____temp23Size[2] = {1,8}; ! +! ! +! ! +! ! +! double ans5; ! +! ! +! ! +! ! +! double __temp24[1 * 10]; ! +! ! +! int ____temp24Size[2] = {1, 10};; ! +! ! +! ! +! ! +! double __temp25[10 * 1]; ! +! ! +! int ____temp25Size[2] = {10, 1};; ! +! ! +! ! +! ! +! double ans6; ! +! ! +! ! +! ! +! double __temp26; ! +! ! +! ! +! ! +! double step_u; ! +! ! +! ! +! ! +! double __temp27; ! +! ! +! ! +! ! +! double* u = NULL; ! +! ! +! int __uSize[2]; ! +! ! +! ! +! ! +! double* v = NULL; ! +! ! +! int __vSize[2]; ! +! ! +! ! +! ! +! double* U = NULL; ! +! ! +! int __USize[2]; ! +! ! +! ! +! ! +! double* V = NULL; ! +! ! +! int __VSize[2]; ! +! ! +! ! +! ! +! double* __temp28 = NULL; ! +! ! +! int ____temp28Size[2]; ! +! ! +! ! +! ! +! double* __temp29 = NULL; ! +! ! +! int ____temp29Size[2]; ! +! ! +! ! +! ! +! double* __temp30 = NULL; ! +! ! +! int ____temp30Size[2]; ! +! ! +! ! +! ! +! double __temp31; ! +! ! +! ! +! ! +! double* __temp32 = NULL; ! +! ! +! int ____temp32Size[2]; ! +! ! +! ! +! ! +! double* __temp33 = NULL; ! +! ! +! int ____temp33Size[2]; ! +! ! +! ! +! ! +! double* __temp34 = NULL; ! +! ! +! int ____temp34Size[2]; ! +! ! +! ! +! ! +! double* __temp35 = NULL; ! +! ! +! int ____temp35Size[2]; ! +! ! +! ! +! ! +! double* __temp36 = NULL; ! +! ! +! int ____temp36Size[2]; ! +! ! +! ! +! ! +! double* __temp37 = NULL; ! +! ! +! int ____temp37Size[2]; ! +! ! +! ! +! ! +! double* __temp38 = NULL; ! +! ! +! int ____temp38Size[2]; ! +! ! +! ! +! ! +! double* __temp39 = NULL; ! +! ! +! int ____temp39Size[2]; ! +! ! +! ! +! ! +! double* __temp40 = NULL; ! +! ! +! int ____temp40Size[2]; ! +! ! +! ! +! ! +! double* __temp41 = NULL; ! +! ! +! int ____temp41Size[2]; ! +! ! +! ! +! ! +! double* x1 = NULL; ! +! ! +! int __x1Size[2]; ! +! ! +! ! +! ! +! double* __temp42 = NULL; ! +! ! +! int ____temp42Size[2]; ! +! ! +! ! +! ! +! double* __temp43 = NULL; ! +! ! +! int ____temp43Size[2]; ! +! ! +! ! +! ! +! double* __temp44 = NULL; ! +! ! +! int ____temp44Size[2]; ! +! ! +! ! +! ! +! double __temp45; ! +! ! +! ! +! ! +! double* __temp46 = NULL; ! +! ! +! int ____temp46Size[2]; ! +! ! +! ! +! ! +! double* __temp47 = NULL; ! +! ! +! int ____temp47Size[2]; ! +! ! +! ! +! ! +! double* __temp48 = NULL; ! +! ! +! int ____temp48Size[2]; ! +! ! +! ! +! ! +! double* __temp49 = NULL; ! +! ! +! int ____temp49Size[2]; ! +! ! +! ! +! ! +! double* __temp50 = NULL; ! +! ! +! int ____temp50Size[2]; ! +! ! +! ! +! ! +! double* __temp51 = NULL; ! +! ! +! int ____temp51Size[2]; ! +! ! +! ! +! ! +! double* __temp52 = NULL; ! +! ! +! int ____temp52Size[2]; ! +! ! +! ! +! ! +! double* __temp53 = NULL; ! +! ! +! int ____temp53Size[2]; ! +! ! +! ! +! ! +! double* __temp54 = NULL; ! +! ! +! int ____temp54Size[2]; ! +! ! +! ! +! ! +! double* __temp55 = NULL; ! +! ! +! int ____temp55Size[2]; ! +! ! +! ! +! ! +! double* y1 = NULL; ! +! ! +! int __y1Size[2]; ! +! ! +! ! +! ! +! double* __temp56 = NULL; ! +! ! +! int ____temp56Size[2]; ! +! ! +! ! +! ! +! double* __temp57 = NULL; ! +! ! +! int ____temp57Size[2]; ! +! ! +! ! +! ! +! double* __temp58 = NULL; ! +! ! +! int ____temp58Size[2]; ! +! ! +! ! +! ! +! double __temp59; ! +! ! +! ! +! ! +! double* __temp60 = NULL; ! +! ! +! int ____temp60Size[2]; ! +! ! +! ! +! ! +! double* __temp61 = NULL; ! +! ! +! int ____temp61Size[2]; ! +! ! +! ! +! ! +! double* __temp62 = NULL; ! +! ! +! int ____temp62Size[2]; ! +! ! +! ! +! ! +! double* __temp63 = NULL; ! +! ! +! int ____temp63Size[2]; ! +! ! +! ! +! ! +! double* __temp64 = NULL; ! +! ! +! int ____temp64Size[2]; ! +! ! +! ! +! ! +! double* __temp65 = NULL; ! +! ! +! int ____temp65Size[2]; ! +! ! +! ! +! ! +! double* __temp66 = NULL; ! +! ! +! int ____temp66Size[2]; ! +! ! +! ! +! ! +! double* __temp67 = NULL; ! +! ! +! int ____temp67Size[2]; ! +! ! +! ! +! ! +! double* __temp68 = NULL; ! +! ! +! int ____temp68Size[2]; ! +! ! +! ! +! ! +! double* __temp69 = NULL; ! +! ! +! int ____temp69Size[2]; ! +! ! +! ! +! ! +! double* z1 = NULL; ! +! ! +! int __z1Size[2]; ! +! ! +! ! +! ! +! char __temp70[1*9] = {"x1(1:10)"}; ! +! ! +! int ____temp70Size[2] = {1,9}; ! +! ! +! ! +! ! +! double ans7; ! +! ! +! ! +! ! +! double __temp71[1 * 10]; ! +! ! +! int ____temp71Size[2] = {1, 10};; ! +! ! +! ! +! ! +! double __temp72[10 * 1]; ! +! ! +! int ____temp72Size[2] = {10, 1};; ! +! ! +! ! +! ! +! double ans8; ! +! ! +! ! +! ! +! char __temp73[1*9] = {"y1(1:10)"}; ! +! ! +! int ____temp73Size[2] = {1,9}; ! +! ! +! ! +! ! +! double ans9; ! +! ! +! ! +! ! +! double __temp74[1 * 10]; ! +! ! +! int ____temp74Size[2] = {1, 10};; ! +! ! +! ! +! ! +! double __temp75[10 * 1]; ! +! ! +! int ____temp75Size[2] = {10, 1};; ! +! ! +! ! +! ! +! double ans10; ! +! ! +! ! +! ! +! char __temp76[1*9] = {"z1(1:10)"}; ! +! ! +! int ____temp76Size[2] = {1,9}; ! +! ! +! ! +! ! +! double ans11; ! +! ! +! ! +! ! +! double __temp77[1 * 10]; ! +! ! +! int ____temp77Size[2] = {1, 10};; ! +! ! +! ! +! ! +! double __temp78[10 * 1]; ! +! ! +! int ____temp78Size[2] = {10, 1};; ! +! ! +! ! +! ! +! double ans12; ! +! ! +! ! +! ! +! ! +! ! +! ! +! ! +! ! +! ! +!/* ! +! ! +!** --------------------------------- ! +! ! +!** --- End Variable Declaration. --- ! +! ! +!** --------------------------------- ! +! ! +!*/ ! +! ! +!/* ! +! ! +!** --------------- ! +! ! +!** --- C code. --- ! +! ! +!** --------------- ! +! ! +!*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: //the program construct moebius ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // and klein bottle surface from ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: //http://www.metlin.org/mathematics/ ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // topological-constructs/non-orientable-surfaces/ ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // --------------------- ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // --- Moebius strip --- ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // --------------------- ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: num_points = 20; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! num_points = d0OpEquald0(20); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: R = 20; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! R = d0OpEquald0(20); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: a = -0.4; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: b = 0.4; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! b = d0OpEquald0(0.4); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: step=(b-a)./(num_points-1); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! __temp1 = d0d0OpMinusd0(b,a); ! +! ! +! __temp2 = d0d0OpMinusd0(num_points,1); ! +! ! +! step = d0d0OpDotSlashd0(__temp1,__temp2); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: step1=(2*%pi)./num_points; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! __temp3 = d0d0OpStard0(2,SCI2C_PI); ! +! ! +! step1 = d0d0OpDotSlashd0(__temp3,num_points); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: s=a:step:b; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! __sSize[0]=1; ! +! ! +! __sSize[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! s = (double*) realloc(s,(1)*(d0absd0(d0floord0(((0.8/d0reald0(step)! +! )+1))))*sizeof(double)); ! +! ! +! d0d0d0OpColond2(a,step,b,s); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: t=0:step1:2*%pi; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! __temp4 = d0d0OpStard0(2,SCI2C_PI); ! +! ! +! __tSize[0]=1; ! +! ! +! __tSize[1]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0(step1! +! ))+1))); ! +! ! +! t = (double*) realloc(t,(1)*(d0absd0(d0floord0((((d0reald0(__temp4)! +! -0)/d0reald0(step1))+1))))*sizeof(double)); ! +! ! +! d0d0d0OpColond2(0,step1,__temp4,t); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: [S, T] = XYgrid(s, t); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! __SSize[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0(step1! +! ))+1))); ! +! ! +! __SSize[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! S = (double*) realloc(S,(d0absd0(d0floord0((((d0reald0(__temp4)-0)/! +! d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0reald0(step))+1! +! ))))*sizeof(double)); ! +! ! +! __TSize[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0(step1! +! ))+1))); ! +! ! +! __TSize[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! T = (double*) realloc(T,(d0absd0(d0floord0((((d0reald0(__temp4)-0)/! +! d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0reald0(step))+1! +! ))))*sizeof(double)); ! +! ! +! d2d2XYgridd2d2(s, __sSize,t, __tSize,S,T); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: x = [R + (S.*cos(T./2))].*cos(T); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ____temp5Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0! +! (step1))+1))); ! +! ! +! ____temp5Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp5 = (double*) realloc(__temp5,(d0absd0(d0floord0((((d0reald0(! +! __temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0rea! +! ld0(step))+1))))*sizeof(double)); ! +! ! +! d2d0OpDotSlashd2(T, __TSize,2,__temp5); ! +! ! +! ____temp6Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0! +! (step1))+1))); ! +! ! +! ____temp6Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp6 = (double*) realloc(__temp6,(d0absd0(d0floord0((((d0reald0(! +! __temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0rea! +! ld0(step))+1))))*sizeof(double)); ! +! ! +! d2cosd2(__temp5, ____temp5Size,__temp6); ! +! ! +! ____temp7Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0! +! (step1))+1))); ! +! ! +! ____temp7Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp7 = (double*) realloc(__temp7,(d0absd0(d0floord0((((d0reald0(! +! __temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0rea! +! ld0(step))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(S, __SSize,__temp6, ____temp6Size,__temp7); ! +! ! +! ____temp8Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0! +! (step1))+1))); ! +! ! +! ____temp8Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp8 = (double*) realloc(__temp8,(d0absd0(d0floord0((((d0reald0(! +! __temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0rea! +! ld0(step))+1))))*sizeof(double)); ! +! ! +! d0d2OpPlusd2(R,__temp7, ____temp7Size,__temp8); ! +! ! +! ____temp9Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0! +! (step1))+1))); ! +! ! +! ____temp9Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp9 = (double*) realloc(__temp9,(d0absd0(d0floord0((((d0reald0(! +! __temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0rea! +! ld0(step))+1))))*sizeof(double)); ! +! ! +! d2cosd2(T, __TSize,__temp9); ! +! ! +! __xSize[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0(step1! +! ))+1))); ! +! ! +! __xSize[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! x = (double*) realloc(x,(d0absd0(d0floord0((((d0reald0(__temp4)-0)/! +! d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0reald0(step))+1! +! ))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp8, ____temp8Size,__temp9, ____temp9Size,x);! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: y = [R + (S.*cos(T./2))].*sin(T); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ____temp10Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald! +! 0(step1))+1))); ! +! ! +! ____temp10Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp10 = (double*) realloc(__temp10,(d0absd0(d0floord0((((d0reald! +! 0(__temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0r! +! eald0(step))+1))))*sizeof(double)); ! +! ! +! d2d0OpDotSlashd2(T, __TSize,2,__temp10); ! +! ! +! ____temp11Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald! +! 0(step1))+1))); ! +! ! +! ____temp11Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp11 = (double*) realloc(__temp11,(d0absd0(d0floord0((((d0reald! +! 0(__temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0r! +! eald0(step))+1))))*sizeof(double)); ! +! ! +! d2cosd2(__temp10, ____temp10Size,__temp11); ! +! ! +! ____temp12Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald! +! 0(step1))+1))); ! +! ! +! ____temp12Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp12 = (double*) realloc(__temp12,(d0absd0(d0floord0((((d0reald! +! 0(__temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0r! +! eald0(step))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(S, __SSize,__temp11, ____temp11Size,__temp12); ! +! ! +! ____temp13Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald! +! 0(step1))+1))); ! +! ! +! ____temp13Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp13 = (double*) realloc(__temp13,(d0absd0(d0floord0((((d0reald! +! 0(__temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0r! +! eald0(step))+1))))*sizeof(double)); ! +! ! +! d0d2OpPlusd2(R,__temp12, ____temp12Size,__temp13); ! +! ! +! ____temp14Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald! +! 0(step1))+1))); ! +! ! +! ____temp14Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp14 = (double*) realloc(__temp14,(d0absd0(d0floord0((((d0reald! +! 0(__temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0r! +! eald0(step))+1))))*sizeof(double)); ! +! ! +! d2sind2(T, __TSize,__temp14); ! +! ! +! __ySize[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0(step1! +! ))+1))); ! +! ! +! __ySize[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! y = (double*) realloc(y,(d0absd0(d0floord0((((d0reald0(__temp4)-0)/! +! d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0reald0(step))+1! +! ))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp13, ____temp13Size,__temp14, ____temp14Size! +! ,y); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: z = S.*sin(T./2); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ____temp15Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald! +! 0(step1))+1))); ! +! ! +! ____temp15Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp15 = (double*) realloc(__temp15,(d0absd0(d0floord0((((d0reald! +! 0(__temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0r! +! eald0(step))+1))))*sizeof(double)); ! +! ! +! d2d0OpDotSlashd2(T, __TSize,2,__temp15); ! +! ! +! ____temp16Size[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald! +! 0(step1))+1))); ! +! ! +! ____temp16Size[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! __temp16 = (double*) realloc(__temp16,(d0absd0(d0floord0((((d0reald! +! 0(__temp4)-0)/d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0r! +! eald0(step))+1))))*sizeof(double)); ! +! ! +! d2sind2(__temp15, ____temp15Size,__temp16); ! +! ! +! __zSize[0]=d0absd0(d0floord0((((d0reald0(__temp4)-0)/d0reald0(step1! +! ))+1))); ! +! ! +! __zSize[1]=d0absd0(d0floord0(((0.8/d0reald0(step))+1))); ! +! ! +! z = (double*) realloc(z,(d0absd0(d0floord0((((d0reald0(__temp4)-0)/! +! d0reald0(step1))+1))))*(d0absd0(d0floord0(((0.8/d0reald0(step))+1! +! ))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(S, __SSize,__temp16, ____temp16Size,z); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp('x(1:10)') ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ans1 = g2dispd0(__temp17, ____temp17Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp(x(1:10)) ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! d0d0OpColond2(1,10,__temp18); ! +! ! +! d2d2OpExtd2(x, __xSize,__temp18, ____temp18Size,__temp19); ! +! ! +! ans2 = d2dispd0(__temp19, ____temp19Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp('y(1:10)') ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ans3 = g2dispd0(__temp20, ____temp20Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp(y(1:10)) ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! d0d0OpColond2(1,10,__temp21); ! +! ! +! d2d2OpExtd2(y, __ySize,__temp21, ____temp21Size,__temp22); ! +! ! +! ans4 = d2dispd0(__temp22, ____temp22Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp('z(1:10)') ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ans5 = g2dispd0(__temp23, ____temp23Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp(z(1:10)) ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! d0d0OpColond2(1,10,__temp24); ! +! ! +! d2d2OpExtd2(z, __zSize,__temp24, ____temp24Size,__temp25); ! +! ! +! ans6 = d2dispd0(__temp25, ____temp25Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // ENABLE/DISABLE PLOT ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // scf(0); surf(x, y, z); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // ----------------------- ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // --- Klein Bottle II --- ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // ----------------------- ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: num_points = 50; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! num_points = d0OpEquald0(50); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: //u = linspace (0, 2*(%pi), num_points); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: //v = linspace (0, 2*(%pi), num_points); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: step_u = (2*%pi)./num_points ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! __temp26 = d0d0OpStard0(2,SCI2C_PI); ! +! ! +! step_u = d0d0OpDotSlashd0(__temp26,num_points); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: u = 0:step_u:2*%pi; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! __temp27 = d0d0OpStard0(2,SCI2C_PI); ! +! ! +! __uSize[0]=1; ! +! ! +! __uSize[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(step! +! _u))+1))); ! +! ! +! u = (double*) realloc(u,(1)*(d0absd0(d0floord0((((d0reald0(__temp27! +! )-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d0d0d0OpColond2(0,step_u,__temp27,u); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: v = u; ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! __vSize[0]=1; ! +! ! +! __vSize[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(step! +! _u))+1))); ! +! ! +! v = (double*) realloc(v,(1)*(d0absd0(d0floord0((((d0reald0(__temp27! +! )-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2OpEquald2(u, __uSize,v); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: [U,V] = XYgrid(u,v); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! __USize[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(step! +! _u))+1))); ! +! ! +! __USize[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(step! +! _u))+1))); ! +! ! +! U = (double*) realloc(U,(d0absd0(d0floord0((((d0reald0(__temp27)-0)! +! /d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0reald0(__temp27)! +! -0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! __VSize[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(step! +! _u))+1))); ! +! ! +! __VSize[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(step! +! _u))+1))); ! +! ! +! V = (double*) realloc(V,(d0absd0(d0floord0((((d0reald0(__temp27)-0)! +! /d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0reald0(__temp27)! +! -0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2XYgridd2d2(u, __uSize,v, __vSize,U,V); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: x1 = cos(U).*(cos(U./2).*(sqrt(2)+cos(V))+(sin(U./2).*sin(! +! V).*cos(V))); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ____temp28Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp28Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp28 = (double*) realloc(__temp28,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(U, __USize,__temp28); ! +! ! +! ____temp29Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp29Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp29 = (double*) realloc(__temp29,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d0OpDotSlashd2(U, __USize,2,__temp29); ! +! ! +! ____temp30Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp30Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp30 = (double*) realloc(__temp30,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(__temp29, ____temp29Size,__temp30); ! +! ! +! __temp31 = d0sqrtd0(2); ! +! ! +! ____temp32Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp32Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp32 = (double*) realloc(__temp32,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(V, __VSize,__temp32); ! +! ! +! ____temp33Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp33Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp33 = (double*) realloc(__temp33,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d0d2OpPlusd2(__temp31,__temp32, ____temp32Size,__temp33); ! +! ! +! ____temp34Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp34Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp34 = (double*) realloc(__temp34,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp30, ____temp30Size,__temp33, ____temp33Size! +! ,__temp34); ! +! ! +! ____temp35Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp35Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp35 = (double*) realloc(__temp35,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d0OpDotSlashd2(U, __USize,2,__temp35); ! +! ! +! ____temp36Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp36Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp36 = (double*) realloc(__temp36,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2sind2(__temp35, ____temp35Size,__temp36); ! +! ! +! ____temp37Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp37Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp37 = (double*) realloc(__temp37,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2sind2(V, __VSize,__temp37); ! +! ! +! ____temp38Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp38Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp38 = (double*) realloc(__temp38,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp36, ____temp36Size,__temp37, ____temp37Size! +! ,__temp38); ! +! ! +! ____temp39Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp39Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp39 = (double*) realloc(__temp39,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(V, __VSize,__temp39); ! +! ! +! ____temp40Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp40Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp40 = (double*) realloc(__temp40,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp38, ____temp38Size,__temp39, ____temp39Size! +! ,__temp40); ! +! ! +! ____temp41Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp41Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp41 = (double*) realloc(__temp41,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpPlusd2(__temp34, ____temp34Size,__temp40, ____temp40Size,__! +! temp41); ! +! ! +! __x1Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(ste! +! p_u))+1))); ! +! ! +! __x1Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(ste! +! p_u))+1))); ! +! ! +! x1 = (double*) realloc(x1,(d0absd0(d0floord0((((d0reald0(__temp27)-! +! 0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0reald0(__temp2! +! 7)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp28, ____temp28Size,__temp41, ____temp41Size! +! ,x1); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: y1 = sin(U).*(cos(U./2).*(sqrt(2)+cos(V))+(sin(U./2).*sin(! +! V).*cos(V))); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ____temp42Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp42Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp42 = (double*) realloc(__temp42,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2sind2(U, __USize,__temp42); ! +! ! +! ____temp43Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp43Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp43 = (double*) realloc(__temp43,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d0OpDotSlashd2(U, __USize,2,__temp43); ! +! ! +! ____temp44Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp44Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp44 = (double*) realloc(__temp44,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(__temp43, ____temp43Size,__temp44); ! +! ! +! __temp45 = d0sqrtd0(2); ! +! ! +! ____temp46Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp46Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp46 = (double*) realloc(__temp46,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(V, __VSize,__temp46); ! +! ! +! ____temp47Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp47Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp47 = (double*) realloc(__temp47,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d0d2OpPlusd2(__temp45,__temp46, ____temp46Size,__temp47); ! +! ! +! ____temp48Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp48Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp48 = (double*) realloc(__temp48,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp44, ____temp44Size,__temp47, ____temp47Size! +! ,__temp48); ! +! ! +! ____temp49Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp49Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp49 = (double*) realloc(__temp49,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d0OpDotSlashd2(U, __USize,2,__temp49); ! +! ! +! ____temp50Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp50Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp50 = (double*) realloc(__temp50,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2sind2(__temp49, ____temp49Size,__temp50); ! +! ! +! ____temp51Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp51Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp51 = (double*) realloc(__temp51,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2sind2(V, __VSize,__temp51); ! +! ! +! ____temp52Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp52Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp52 = (double*) realloc(__temp52,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp50, ____temp50Size,__temp51, ____temp51Size! +! ,__temp52); ! +! ! +! ____temp53Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp53Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp53 = (double*) realloc(__temp53,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(V, __VSize,__temp53); ! +! ! +! ____temp54Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp54Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp54 = (double*) realloc(__temp54,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp52, ____temp52Size,__temp53, ____temp53Size! +! ,__temp54); ! +! ! +! ____temp55Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp55Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp55 = (double*) realloc(__temp55,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpPlusd2(__temp48, ____temp48Size,__temp54, ____temp54Size,__! +! temp55); ! +! ! +! __y1Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(ste! +! p_u))+1))); ! +! ! +! __y1Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(ste! +! p_u))+1))); ! +! ! +! y1 = (double*) realloc(y1,(d0absd0(d0floord0((((d0reald0(__temp27)-! +! 0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0reald0(__temp2! +! 7)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp42, ____temp42Size,__temp55, ____temp55Size! +! ,y1); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: z1 = -1.*sin(U./2).*(sqrt(2)+cos(V))+cos(U./2).*sin(V).*co! +! s(V); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ____temp56Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp56Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp56 = (double*) realloc(__temp56,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d0OpDotSlashd2(U, __USize,2,__temp56); ! +! ! +! ____temp57Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp57Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp57 = (double*) realloc(__temp57,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2sind2(__temp56, ____temp56Size,__temp57); ! +! ! +! ____temp58Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp58Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp58 = (double*) realloc(__temp58,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d0d2OpStard2(1,__temp57, ____temp57Size,__temp58); ! +! ! +! __temp59 = d0sqrtd0(2); ! +! ! +! ____temp60Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp60Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp60 = (double*) realloc(__temp60,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(V, __VSize,__temp60); ! +! ! +! ____temp61Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp61Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp61 = (double*) realloc(__temp61,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d0d2OpPlusd2(__temp59,__temp60, ____temp60Size,__temp61); ! +! ! +! ____temp62Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp62Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp62 = (double*) realloc(__temp62,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp58, ____temp58Size,__temp61, ____temp61Size! +! ,__temp62); ! +! ! +! ____temp63Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp63Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp63 = (double*) realloc(__temp63,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2OpMinusd2(__temp62, ____temp62Size,__temp63); ! +! ! +! ____temp64Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp64Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp64 = (double*) realloc(__temp64,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d0OpDotSlashd2(U, __USize,2,__temp64); ! +! ! +! ____temp65Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp65Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp65 = (double*) realloc(__temp65,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(__temp64, ____temp64Size,__temp65); ! +! ! +! ____temp66Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp66Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp66 = (double*) realloc(__temp66,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2sind2(V, __VSize,__temp66); ! +! ! +! ____temp67Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp67Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp67 = (double*) realloc(__temp67,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp65, ____temp65Size,__temp66, ____temp66Size! +! ,__temp67); ! +! ! +! ____temp68Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp68Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp68 = (double*) realloc(__temp68,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2cosd2(V, __VSize,__temp68); ! +! ! +! ____temp69Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! ____temp69Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0real! +! d0(step_u))+1))); ! +! ! +! __temp69 = (double*) realloc(__temp69,(d0absd0(d0floord0((((d0reald! +! 0(__temp27)-0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0re! +! ald0(__temp27)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpDotStard2(__temp67, ____temp67Size,__temp68, ____temp68Size! +! ,__temp69); ! +! ! +! __z1Size[0]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(ste! +! p_u))+1))); ! +! ! +! __z1Size[1]=d0absd0(d0floord0((((d0reald0(__temp27)-0)/d0reald0(ste! +! p_u))+1))); ! +! ! +! z1 = (double*) realloc(z1,(d0absd0(d0floord0((((d0reald0(__temp27)-! +! 0)/d0reald0(step_u))+1))))*(d0absd0(d0floord0((((d0reald0(__temp2! +! 7)-0)/d0reald0(step_u))+1))))*sizeof(double)); ! +! ! +! d2d2OpPlusd2(__temp63, ____temp63Size,__temp69, ____temp69Size,z1! +! ); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp('x1(1:10)'); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ans7 = g2dispd0(__temp70, ____temp70Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp(x1(1:10)); ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! d0d0OpColond2(1,10,__temp71); ! +! ! +! d2d2OpExtd2(x1, __x1Size,__temp71, ____temp71Size,__temp72); ! +! ! +! ans8 = d2dispd0(__temp72, ____temp72Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp('y1(1:10)') ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ans9 = g2dispd0(__temp73, ____temp73Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp(y1(1:10)) ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! d0d0OpColond2(1,10,__temp74); ! +! ! +! d2d2OpExtd2(y1, __y1Size,__temp74, ____temp74Size,__temp75); ! +! ! +! ans10 = d2dispd0(__temp75, ____temp75Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp('z1(1:10)') ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ans11 = g2dispd0(__temp76, ____temp76Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: disp(z1(1:10)) ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! d0d0OpColond2(1,10,__temp77); ! +! ! +! d2d2OpExtd2(z1, __z1Size,__temp77, ____temp77Size,__temp78); ! +! ! +! ans12 = d2dispd0(__temp78, ____temp78Size); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // ENABLE/DISABLE PLOT ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: // scf(1) ;surf(x1, y1, z1) ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: endfunction ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +! ! +! ! +! /* ! +! ! +! ** --------------------- ! +! ! +! ** --- Free Section. --- ! +! ! +! ** --------------------- ! +! ! +! */ ! +! ! +! free(s); ! +! ! +! free(t); ! +! ! +! free(S); ! +! ! +! free(T); ! +! ! +! free(__temp5); ! +! ! +! free(__temp6); ! +! ! +! free(__temp7); ! +! ! +! free(__temp8); ! +! ! +! free(__temp9); ! +! ! +! free(x); ! +! ! +! free(__temp10); ! +! ! +! free(__temp11); ! +! ! +! free(__temp12); ! +! ! +! free(__temp13); ! +! ! +! free(__temp14); ! +! ! +! free(y); ! +! ! +! free(__temp15); ! +! ! +! free(__temp16); ! +! ! +! free(z); ! +! ! +! free(u); ! +! ! +! free(v); ! +! ! +! free(U); ! +! ! +! free(V); ! +! ! +! free(__temp28); ! +! ! +! free(__temp29); ! +! ! +! free(__temp30); ! +! ! +! free(__temp32); ! +! ! +! free(__temp33); ! +! ! +! free(__temp34); ! +! ! +! free(__temp35); ! +! ! +! free(__temp36); ! +! ! +! free(__temp37); ! +! ! +! free(__temp38); ! +! ! +! free(__temp39); ! +! ! +! free(__temp40); ! +! ! +! free(__temp41); ! +! ! +! free(x1); ! +! ! +! free(__temp42); ! +! ! +! free(__temp43); ! +! ! +! free(__temp44); ! +! ! +! free(__temp46); ! +! ! +! free(__temp47); ! +! ! +! free(__temp48); ! +! ! +! free(__temp49); ! +! ! +! free(__temp50); ! +! ! +! free(__temp51); ! +! ! +! free(__temp52); ! +! ! +! free(__temp53); ! +! ! +! free(__temp54); ! +! ! +! free(__temp55); ! +! ! +! free(y1); ! +! ! +! free(__temp56); ! +! ! +! free(__temp57); ! +! ! +! free(__temp58); ! +! ! +! free(__temp60); ! +! ! +! free(__temp61); ! +! ! +! free(__temp62); ! +! ! +! free(__temp63); ! +! ! +! free(__temp64); ! +! ! +! free(__temp65); ! +! ! +! free(__temp66); ! +! ! +! free(__temp67); ! +! ! +! free(__temp68); ! +! ! +! free(__temp69); ! +! ! +! free(z1); ! +! ! +! /* ! +! ! +! ** ------------------------- ! +! ! +! ** --- End Free Section. --- ! +! ! +! ** ------------------------- ! +! ! +! */ ! +! ! +! ! +! ! +! return(0); ! +! ! +! ! +! ! +! /*SCI2C: ##########################################################! +! ######## ! +! ! +! SCI2C: ! +! ! +! SCI2C: ##########################################################! +! ########*/ ! +! ! +!} ! +! ! +! ! +mclose(fd); |