// // 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);