← Back to branch summary

~ubuntu-branches/ubuntu/saucy/lazarus/saucy

« back to all changes in this revision

Viewing changes to components/aggpas/src/agg_simul_eq.pas

  • Committer: Package Import Robot
  • Author(s): Paul Gevers, Abou Al Montacir, Bart Martens, Paul Gevers
  • Date: 2013-06-08 14:12:17 UTC
  • mfrom: (1.1.9)
  • Revision ID: package-import@ubuntu.com-20130608141217-7k0cy9id8ifcnutc
Tags: 1.0.8+dfsg-1
http://bugs.debian.org/708695
[ Abou Al Montacir ]
* New upstream major release and multiple maintenace release offering many
  fixes and new features marking a new milestone for the Lazarus development
  and its stability level.
  - The detailed list of changes can be found here:
    http://wiki.lazarus.freepascal.org/Lazarus_1.0_release_notes
    http://wiki.lazarus.freepascal.org/Lazarus_1.0_fixes_branch
* LCL changes:
  - LCL is now a normal package.
      + Platform independent parts of the LCL are now in the package LCLBase
      + LCL is automatically recompiled when switching the target platform,
        unless pre-compiled binaries for this target are already installed.
      + No impact on existing projects.
      + Linker options needed by LCL are no more added to projects that do
        not use the LCL package.
  - Minor changes in LCL basic classes behaviour
      + TCustomForm.Create raises an exception if a form resource is not
        found.
      + TNotebook and TPage: a new implementation of these classes was added.
      + TDBNavigator: It is now possible to have focusable buttons by setting
        Options = [navFocusableButtons] and TabStop = True, useful for
        accessibility and for devices with neither mouse nor touch screen.
      + Names of TControlBorderSpacing.GetSideSpace and GetSpace were swapped
        and are now consistent. GetSideSpace = Around + GetSpace.
      + TForm.WindowState=wsFullscreen was added
      + TCanvas.TextFitInfo was added to calculate how many characters will
        fit into a specified Width. Useful for word-wrapping calculations.
      + TControl.GetColorResolvingParent and
        TControl.GetRGBColorResolvingParent were added, simplifying the work
        to obtain the final color of the control while resolving clDefault
        and the ParentColor.
      + LCLIntf.GetTextExtentExPoint now has a good default implementation
        which works in any platform not providing a specific implementation.
        However, Widgetset specific implementation is better, when available.
      + TTabControl was reorganized. Now it has the correct class hierarchy
        and inherits from TCustomTabControl as it should.
  - New unit in the LCL:
      + lazdialogs.pas: adds non-native versions of various native dialogs,
        for example TLazOpenDialog, TLazSaveDialog, TLazSelectDirectoryDialog.
        It is used by widgetsets which either do not have a native dialog, or
        do not wish to use it because it is limited. These dialogs can also be
        used by user applications directly.
      + lazdeviceapis.pas: offers an interface to more hardware devices such
        as the accelerometer, GPS, etc. See LazDeviceAPIs
      + lazcanvas.pas: provides a TFPImageCanvas descendent implementing
        drawing in a LCL-compatible way, but 100% in Pascal.
      + lazregions.pas. LazRegions is a wholly Pascal implementation of
        regions for canvas clipping, event clipping, finding in which control
        of a region tree one an event should reach, for drawing polygons, etc.
      + customdrawncontrols.pas, customdrawndrawers.pas,
        customdrawn_common.pas, customdrawn_android.pas and
        customdrawn_winxp.pas: are the Lazarus Custom Drawn Controls -controls
        which imitate the standard LCL ones, but with the difference that they
        are non-native and support skinning.
  - New APIs added to the LCL to improve support of accessibility software
    such as screen readers.
* IDE changes:
  - Many improvments.
  - The detailed list of changes can be found here:
    http://wiki.lazarus.freepascal.org/New_IDE_features_since#v1.0_.282012-08-29.29
    http://wiki.lazarus.freepascal.org/Lazarus_1.0_release_notes#IDE_Changes
* Debugger / Editor changes:
  - Added pascal sources and breakpoints to the disassembler
  - Added threads dialog.
* Components changes:
  - TAChart: many fixes and new features
  - CodeTool: support Delphi style generics and new syntax extensions.
  - AggPas: removed to honor free licencing. (Closes: Bug#708695)
[Bart Martens]
* New debian/watch file fixing issues with upstream RC release.
[Abou Al Montacir]
* Avoid changing files in .pc hidden directory, these are used by quilt for
  internal purpose and could lead to surprises during build.
[Paul Gevers]
* Updated get-orig-source target and it compinion script orig-tar.sh so that they
  repack the source file, allowing bug 708695 to be fixed.

Show diffs side-by-side

added added

removed removed

Lines of Context:
components/aggpas/src/agg_simul_eq.pas
1
 
//----------------------------------------------------------------------------
2
 
// Anti-Grain Geometry - Version 2.4 (Public License)
3
 
// Copyright (C) 2002-2005 Maxim Shemanarev (http://www.antigrain.com)
4
 
//
5
 
// Anti-Grain Geometry - Version 2.4 Release Milano 3 (AggPas 2.4 RM3)
6
 
// Pascal Port By: Milan Marusinec alias Milano
7
 
//                 milan@marusinec.sk
8
 
//                 http://www.aggpas.org
9
 
// Copyright (c) 2005-2006
10
 
//
11
 
// Permission to copy, use, modify, sell and distribute this software
12
 
// is granted provided this copyright notice appears in all copies.
13
 
// This software is provided "as is" without express or implied
14
 
// warranty, and with no claim as to its suitability for any purpose.
15
 
//
16
 
//----------------------------------------------------------------------------
17
 
// Contact: mcseem@antigrain.com
18
 
//          mcseemagg@yahoo.com
19
 
//          http://www.antigrain.com
20
 
//
21
 
// [Pascal Port History] -----------------------------------------------------
22
 
//
23
 
// 23.06.2006-Milano: ptrcomp adjustments
24
 
// 06.02.2006-Milano: Unit port establishment
25
 
//
26
 
{ agg_simul_eq.pas }
27
 
unit
28
 
 agg_simul_eq ;
29
 
 
30
 
INTERFACE
31
 
 
32
 
{$I agg_mode.inc }
33
 
 
34
 
uses
35
 
 agg_basics ;
36
 
 
37
 
{ GLOBAL PROCEDURES }
38
 
 procedure swap_arrays   (a1 ,a2 : double_ptr; n : unsigned );
39
 
 function  matrix_pivot  (m : double_ptr; row ,Rows ,Cols : unsigned ) : int;
40
 
 function  simul_eq_solve(left ,right ,result_ : double_ptr; Size ,RightCols : unsigned ) : boolean;
41
 
 
42
 
 
43
 
IMPLEMENTATION
44
 
{ LOCAL VARIABLES & CONSTANTS }
45
 
{ UNIT IMPLEMENTATION }
46
 
{ SWAP_ARRAYS }
47
 
procedure swap_arrays;
48
 
var
49
 
 i   : unsigned;
50
 
 tmp : double;
51
 
 
52
 
begin
53
 
 i:=0;
54
 
 
55
 
 while i < n do
56
 
  begin
57
 
   tmp:=a1^;
58
 
   a1^:=a2^;
59
 
   a2^:=tmp;
60
 
 
61
 
   inc(ptrcomp(a1 ) ,sizeof(double ) );
62
 
   inc(ptrcomp(a2 ) ,sizeof(double ) );
63
 
   inc(i );
64
 
 
65
 
  end;
66
 
 
67
 
end;
68
 
 
69
 
{ MATRIX_PIVOT }
70
 
function matrix_pivot;
71
 
var
72
 
 i : unsigned;
73
 
 k : int;
74
 
 
75
 
 max_val ,tmp : double;
76
 
 
77
 
begin
78
 
 k:=row;
79
 
 
80
 
 max_val:=-1.0;
81
 
 
82
 
 i:=row;
83
 
 
84
 
 while i < Rows do
85
 
  begin
86
 
   tmp:=Abs(double_ptr(ptrcomp(m ) + (i * Cols + row ) * sizeof(double ) )^ );
87
 
 
88
 
   if (tmp > max_val ) and
89
 
      (tmp <> 0.0 ) then
90
 
    begin
91
 
     max_val:=tmp;
92
 
 
93
 
     k:=i;
94
 
 
95
 
    end;
96
 
 
97
 
   inc(i );
98
 
 
99
 
  end;
100
 
 
101
 
 if double_ptr(ptrcomp(m ) + (k * Cols + row ) * sizeof(double ) )^ = 0.0 then
102
 
  begin
103
 
   result:=-1;
104
 
 
105
 
   exit;
106
 
 
107
 
  end;
108
 
 
109
 
 if k <> row then
110
 
  begin
111
 
   swap_arrays(
112
 
    double_ptr(ptrcomp(m ) + k * Cols * sizeof(double ) ) ,
113
 
    double_ptr(ptrcomp(m ) + row * Cols * sizeof(double ) ) ,
114
 
    Cols );
115
 
 
116
 
   result:=k;
117
 
 
118
 
   exit;
119
 
 
120
 
  end;
121
 
 
122
 
 result:=0;
123
 
 
124
 
end;
125
 
 
126
 
{ SIMUL_EQ_SOLVE }
127
 
function simul_eq_solve;
128
 
var
129
 
 m : int;
130
 
 
131
 
 i , j , k ,adx : unsigned;
132
 
 
133
 
 a1  : double;
134
 
 tmp : double_ptr;
135
 
 
136
 
label
137
 
 return_false ,free ;
138
 
 
139
 
begin
140
 
// Alloc
141
 
 adx:=Size + RightCols;
142
 
 
143
 
 agg_getmem(pointer(tmp ) ,Size * adx * sizeof(double ) );
144
 
 
145
 
 for i:=0 to Size - 1 do
146
 
  begin
147
 
   for j:=0 to Size - 1 do
148
 
    double_ptr(ptrcomp(tmp ) + (i * adx + j ) * sizeof(double ) )^:=
149
 
     double_ptr(ptrcomp(left ) + (i * Size + j ) * sizeof(double ) )^;
150
 
 
151
 
   for j:=0 to RightCols - 1 do
152
 
    double_ptr(ptrcomp(tmp ) + (i * adx + Size + j ) * sizeof(double ) )^:=
153
 
     double_ptr(ptrcomp(right ) + (i * RightCols + j ) * sizeof(double ) )^;
154
 
 
155
 
  end;
156
 
 
157
 
 for k:=0 to Size - 1 do
158
 
  begin
159
 
   if matrix_pivot(tmp ,k ,Size ,Size + RightCols ) < 0 then
160
 
    goto return_false; // Singularity....
161
 
 
162
 
   a1:=double_ptr(ptrcomp(tmp ) + (k * adx + k ) * sizeof(double ) )^;
163
 
   j :=k;
164
 
 
165
 
   while j < Size + RightCols do
166
 
    begin
167
 
     double_ptr(ptrcomp(tmp ) + (k * adx + j ) * sizeof(double ) )^:=
168
 
      double_ptr(ptrcomp(tmp ) + (k * adx + j ) * sizeof(double ) )^ / a1;
169
 
 
170
 
     inc(j );
171
 
 
172
 
    end;
173
 
 
174
 
   i:=k + 1;
175
 
 
176
 
   while i < Size do
177
 
    begin
178
 
     a1:=double_ptr(ptrcomp(tmp ) + (i * adx + k ) * sizeof(double ) )^;
179
 
     j :=k;
180
 
 
181
 
     while j < Size + RightCols do
182
 
      begin
183
 
       double_ptr(ptrcomp(tmp ) + (i * adx + j ) * sizeof(double ) )^:=
184
 
        double_ptr(ptrcomp(tmp ) + (i * adx + j ) * sizeof(double ) )^ -
185
 
        a1 * double_ptr(ptrcomp(tmp ) + (k * adx + j ) * sizeof(double ) )^;
186
 
 
187
 
       inc(j );
188
 
 
189
 
      end;
190
 
 
191
 
     inc(i );
192
 
 
193
 
    end;
194
 
 
195
 
  end;
196
 
 
197
 
 for k:=0 to RightCols - 1 do
198
 
  begin
199
 
   m:=int(Size - 1 );
200
 
 
201
 
   while m >= 0 do
202
 
    begin
203
 
     double_ptr(ptrcomp(result_ ) + (m * RightCols + k ) * sizeof(double ) )^:=
204
 
      double_ptr(ptrcomp(tmp ) + (m * adx + Size + k ) * sizeof(double ) )^;
205
 
 
206
 
     j:=m + 1;
207
 
 
208
 
     while j < Size do
209
 
      begin
210
 
       double_ptr(ptrcomp(result_ ) + (m * RightCols + k ) * sizeof(double ) )^:=
211
 
        double_ptr(ptrcomp(result_ ) + (m * RightCols + k ) * sizeof(double ) )^ -
212
 
        (double_ptr(ptrcomp(tmp ) + (m * adx + j ) * sizeof(double ) )^ *
213
 
         double_ptr(ptrcomp(result_ ) + (j * RightCols + k ) * sizeof(double ) )^ );
214
 
 
215
 
       inc(j );
216
 
 
217
 
      end;
218
 
 
219
 
     dec(m );
220
 
 
221
 
    end;
222
 
 
223
 
  end;
224
 
 
225
 
// Done
226
 
 result:=true;
227
 
 
228
 
 goto free;
229
 
 
230
 
return_false:
231
 
 result:=false;
232
 
 
233
 
// Free
234
 
free:
235
 
 agg_freemem(pointer(tmp ) ,Size * adx * sizeof(double ) );
236
 
 
237
 
end;
238
 
 
239
 
END.
240