Well, you have the "where" instructions.
Then you need to keep track (e.g. with logical vectors) of what's in A,
B, C. For instance, I have the following in my code:
where (.not. genotyped)
sol=0d0
endwhere
There is also the possibility to indicate elements of an array with
another, if I remember correctly:
A(B)=0d0
where, for instance B=(/1,5,9,10/)
Andres
Le 21/03/2012 15:57, Alexander Dimitrov a écrit :
> I have to solve quite complex system of equations. To each equation i have solutions that are arrays with some size. For the different equations that size is different. There are some complex interactions between the equations. Like if in array A elements from some complex pattern did not existed then with the reminding part I would like to make some "whole array" operations on another array B. Those absent elements carry interactions with equation C ets. This must be a common situation solved by Fortran but I do not see the natural way I have to deal with arrays that have different sizes. Do you know how such problem is solved for large arrays?
>
--
Andres Legarra
+33 561285182
INRA, UR 631 SAGA
BP 52627
31326 Castanet Tolosan, France
|