[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Newbiew need help optimizing code and understanding Octave
From: |
malhotrag |
Subject: |
Re: Newbiew need help optimizing code and understanding Octave |
Date: |
Tue, 1 Apr 2008 04:13:51 -0700 (PDT) |
Just a correction the loop should be
for widthstrip = 36:2:84
malhotrag wrote:
>
> Hi All,
>
> Below is the code I am writing as a part of a bigger project. My Issues
> are I am an octave newbie and I used to do programming in C/C++ so I write
> octave code pretty similar to C/C++ code.
>
> The code below runs really slow (it's not that complex though), so I would
> like tips from all on how to make it fast. I know I need to vectorize this
> but I don't really understand the concept of vectorization and so need
> expert advice.
>
> Any help would be really appreciated..
>
> Thanks
>
> h = 80;
> r = 3.5;
> s = 12;
> xmin = 0;
> xmax = 120;
> widthstrip = 36;
> for widthstrip = 82:2:82
> ymin = (h/2) - (widthstrip/2);
> ymax = (h/2) + (widthstrip/2);
> Fd12 = 0;
> sum = 0;
> count = 0;
>
> for i = xmin:xmax
> for j = ymin:ymax
> x = i;
> y = j;
> S = s / r;
> X = x / r;
> Y = y / r;
> H = h / r;
>
> A = X^2 + Y^2 + S^2;
> B = S^2 + X^2;
> C = (H - Y)^2;
>
> part1 = S / B;
> part2 = S / (2 * B * pi);
> part3 = acos( (Y^2 - B + 1) / (A - 1) );
> part4 = acos( (C - B + 1) / (C + B - 1) );
> part5 = (A + 1) / (sqrt( (A - 1)^2 + 4 * Y^2 ));
> part6 = acos( (Y^2 - B + 1) / (sqrt(B) * (A - 1)) );
> part7 = (C + B + 1) / (sqrt( (C + B - 1)^2 + 4 * C ));
> part8 = acos( (C - B + 1) / (sqrt(B) * (C + B - 1)) );
> part9 = H * acos( (1/sqrt(B)) );
>
> sum = sum + (part1 - (part2 * (part3 + part4 - (Y * part5 *
> part6) -
> (sqrt(C) * part7 * part8) + part9)));
> ++count;
> endfor
> endfor
>
> Fd12 = sum;
>
> printf("\nWidth of %d inch Fd1-2 = %f, Count is %d and Fd1-2 final = %f
> ", widthstrip, Fd12,count,Fd12/count);
>
> endfor
>
--
View this message in context:
http://www.nabble.com/Newbiew-need-help-optimizing-code-and-understanding-Octave-tp16397209p16418006.html
Sent from the Octave - General mailing list archive at Nabble.com.
- Re: Newbiew need help optimizing code and understanding Octave,
malhotrag <=