1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
#include "scicos_block.h"
#include <math.h>
void minmax(scicos_block *block,int flag)
{ /*ipar[0]=1 -> min, ipar[0]=2 -> max */
int i,phase;
double maxmin;
phase=get_phase_simulation();
if (flag == 1) {
if(block->nin==1){
if((block->ng==0)||(phase==1)){
maxmin=block->inptr[0][0];
for (i=1;i<block->insz[0];++i){
if(block->ipar[0]==1){
if(block->inptr[0][i]<maxmin) maxmin=block->inptr[0][i];
} else {
if(block->inptr[0][i]>maxmin) maxmin=block->inptr[0][i];
}
}
}else{
maxmin=block->inptr[0][block->mode[0]-1];
}
block->outptr[0][0]=maxmin;
}else if (block->nin==2){
for (i=0;i<block->insz[0];++i){
if((block->ng==0)||(phase==1)){
if(block->ipar[0]==1){
block->outptr[0][i]=min(block->inptr[0][i],block->inptr[1][i]);
} else {
block->outptr[0][i]=max(block->inptr[0][i],block->inptr[1][i]);
}
}else{
block->outptr[0][i]=block->inptr[block->mode[0]-1][i];
}
}
}
} else if(flag == 9){
if(block->nin==1){
if(block->nin==1){
if (phase==2){
for (i=0;i<block->insz[0];++i){
if(i!=block->mode[0]-1){
block->g[i]=block->inptr[0][i]-block->inptr[0][block->mode[0]-1];
}else{
block->g[i]=1.0;
}
}
} else if (phase==1){
maxmin=block->inptr[0][0];
for (i=1;i<block->insz[0];++i){
block->mode[0]=1;
if(block->ipar[0]==1){
if(block->inptr[0][i]<maxmin) {
maxmin=block->inptr[0][i];
block->mode[0]=i+1;
}
} else {
if(block->inptr[0][i]>maxmin) {
maxmin=block->inptr[0][i];
block->mode[0]=i+1;
}
}
}
}
}
}else if(block->nin==2){
for (i=0;i<block->insz[0];++i){
block->g[i]=block->inptr[0][i]-block->inptr[1][i];
if(phase==1){
if(block->ipar[0]==1){
if(block->g[i]>0) {
block->mode[i]=2;
}else{
block->mode[i]=1;
}
}else{
if(block->g[i]<0) {
block->mode[i]=2;
}else{
block->mode[i]=1;
}
}
}
}
}
}
}
|
