All Forums  Origin Forum  Origin Forum  Error incompatible variables

Note: You must be registered in order to post a reply.
To register, click here. Registration is FREE!

Screensize: 640 x 480 800 x 600 1024 x 768 1280 x 1024 UserName: Password: Anti-Spam Code: Format Mode: Basic  Help  Prompt  Format: Font Andale Mono Arial Arial Black Book Antiqua Century Gothic Comic Sans MS Courier New Georgia Impact Lucida Console Script MT Bold Stencil Tahoma Times New Roman Trebuchet MS Verdana   Size 1 2 3 4 5 6   Color Black Red Yellow Pink Green Orange Purple Blue Beige Brown Teal Navy Maroon LimeGreen Message: * HTML is OFF * Forum Code is ON Smilies [quote][i]Originally posted by DB-L[/i] [br]Origin Ver. and Service Release (Select Help-->About Origin): 9.1 Operating System: windows vista on VMware os X Hello, for fitting data, I need to find the root of a polynomial (vx). For this I include Newton-Rhapson method in the fit. I determine a function vfx and its derivate vdf. after fews iterations I obtain the value of vx. This works well with "double" but with my data "vector" . I can't compile, the message : "error: incompatibles variables" appears... My code is pasted just below. Is someone could help me please! Best regards. #pragma numlittype(push, TRUE) #pragma numlittype(push, TRUE) #pragma warning(error : 15618) #include <origin.h> #include <math.h> #include <stdio.h> // Add your special include files here. // For example, if you want to fit with functions from the NAG library, // add the header file for the NAG functions here. #include <ONLSF.H> #include <OC_nag.h> // Add code here for other Origin C functions that you want to define in this file, // and access in your fitting function. // You can access C functions defined in other files, if those files are loaded and compiled // in your workspace, and the functions have been prototyped in a header file that you have // included above. // You can access NLSF object methods and properties directly in your function code. // You should follow C-language syntax in defining your function. // For instance, if your parameter name is P1, you cannot use p1 in your function code. // When using fractions, remember that integer division such as 1/2 is equal to 0, and not 0.5 // Use 0.5 or 1/2.0 to get the correct value. // For more information and examples, please refer to the "User-Defined Fitting Function" // section of the Origin Help file. //---------------------------------------------------------- // void _nlsfmultiSites( // Fit Parameter(s): double n, double L, double K12, double H12, double H1, double K1, double H2, double K2, // Independent Variable(s): double x, // Dependent Variable(s): double& y) { // Beginning of editable part Worksheet wks = Project.ActiveLayer(); NLFitContext *pCtxt = Project.GetNLFitContext(); if ( pCtxt ) { static vector vX, vY; static int nSize; BOOL bIsNewParamValues = pCtxt->IsNewParamValues(); if ( bIsNewParamValues ) { Dataset dsMt( wks, 2 ), dsXMt( wks, 0 ); int nn = dsXMt.GetSize(); vector vTheta( nn ), vC12, vC1, vC2; vector vdH1dx, vdH2dx, vdHdx, vdrdx; vC12 = n*K12*dsMt; vC1 = n*K1*dsMt; vC2 = n*n*K2*dsMt*dsMt; vector vx, vx1, vx2, vx3, vfx, vdf, vdx, vrtn; int J; vx1 = 0; vx2 = L*dsXMt/n; vx3 = dsXMt/1000000000; vrtn = 0.5*(vx1 + vx2); for (J=1; J<=100; J++) { vfx = vrtn^2; // (-(1+vC1*vrtn+vC2*vrtn^2)+sqrt(8*vC12^2*vrtn+(1+vC1*vrtn+vC2*vrtn^2)^2))*(1-vC1*vrtn-3*vC2*vrtn^2)-((4*vC12^2*vrtn)*(1+2*vrtn-(2*L*dsXMt/n))); // function vdf = 2*vrtn; //((-6*vC2*vrtn-vC1)*(sqrt((vC2*vrtn^2+vC1*vrtn+1)+8*vC12^2*vrtn)-vC2*vrtn^2-vC1*vrtn-1) + (-3*vC2*vrtn^2-vC1*vrtn+1)*((2*(2*vC2*vrtn+vC1)*(vC2*vrtn^2+vC1*vrtn+1)+8*vC12^2)/(2*sqrt((vC2*vrtn^2+vC1*vrtn+1)^2+8*vC12^2*vrtn))-2*vC2*vrtn-vC1) - 4*vC12^2*(2*vrtn-(2*L*dsXMt/n)-8*vC12^2*vrtn) ; // derivate of function vdx = - vfx/vdf; vrtn += vdx; [red]if ((vdx*vdx) < vx3) { vx = vrtn; // maybe an exit(0)? }[/red] } vdrdx = -(-vC1-2*vC2*vx+(((1+vC1*vx+vC2*vx^2)*(vC1+2*vC2*vx)+4*vC12^2)/sqrt((1+vC1*vx+vC2*vx^2)^2+8*vC12*vC12*vx)))*((1/vx)-vC1-3*vC2*vx)*(n/(8*vC12*vC12*L))+(-(1+vC1*vx+vC2*vx^2)+sqrt((1+vC1*vx+vC2*vx^2)^2+8*vC12*vC12*vx))*((1/(vx^2)+3*vC2))*(n/(8*vC12^2*L))+(n/L); vdH1dx = (n/(4*vC12^2))*(-vC1-2*vC2*vx+(((1+vC1*vx+vC2*vx^2)*(vC1+2*vC2*vx)+4*vC12^2)/sqrt((1+vC1*vx+vC2*vx^2)^2+8*vC12^2*vx)))*(0,5*H12*((-1/vx)-vC1-vC2*vx)+vC1*H1+vC2*H2*vx); vdH2dx = (n/(4*vC12^2))*(-(1+vC1*vx+vC2*vx^2)+ sqrt((1+vC1*vx+vC2*vx^2)^2+8*vC12^2*vx))*(0,5*H12*((1/(vx^2))-vC2)+H2*vC2); vdHdx = vdH1dx + vdH2dx; vY = vdHdx/vdrdx; } ocmath_interpolate( &x, &y, 1, vX, vY, nSize ); } // End of editable part } [/quote]   Check here to include your profile signature. Check here to subscribe to this topic.