[Help-glpk] Re: "No value" problem, not part of the data set - (MathProg)

[Noli Sicad]

I searched the help-glpk archive. I found out that I have this problem before,

http://lists.gnu.org/archive/html/help-glpk/2004-09/msg00006.html

And the answer is just putting a default 0 to the parameter according to Andrew.

http://lists.gnu.org/archive/html/help-glpk/2004-09/msg00008.html

Answer to my present question is, putting default 0;

param Area_Stand{i in I,j in J,k in K,l in L} default 0;

instead of just this,

param Area_Stand{i in I,j in J,k in K,l in L};

Thanks, Noli


On 12/17/09, Noli Sicad <address@hidden> wrote:
> Hi,
>
> I am trying to read a dbf file (*.dbf) from shapefile (i.e. ESRI
> ArcGIS)  or GIS file (*.shp) and writing a dbf table.
>
> I am wondering how can I solve this problem  - "no value" problem as
> you can see in the console output below, it is not part of the data
> set from dbf file.
>
> ~~~~~~~~~~~~~~~
> Reading tab8...
> Display statement at line 30
> Area_Stand[A,R,1,3] = 131.507
> Area_Stand[A,R,2,3] = 122.955
> Area_Stand[A,'',1,1] = 57.899
> Area_Stand[A,'',2,1] = 15.573
> Area_Stand[A,'',3,1] = 78.544
> Area_Stand[A,'',4,1] = 67.906
> Area_Stand[A,'',5,1] = 145.412
> Area_Stand[A,R,1,2] = 155.223
> Area_Stand[A,'',2,2] = 64.358
> Area_Stand[A,'',3,2] = 35.798
> Area_Stand[A,R,4,2] = 124.597
> Area_Stand[B,'',1,3] = 46.048
> Area_Stand[B,'',2,3] = 73.329
> Area_Stand[B,'',3,3] = 40.189
> Area_Stand[B,'',1,1] = 4.184
> Area_Stand[B,'',2,1] = 15.175
> Area_Stand[B,'',3,1] = 52.294
> Area_Stand[B,'',1,2] = 78.506
> Area_Stand[B,'',2,2] = 32.534
> Area_Stand[B,'',3,2] = 22.573
> Area_Stand[B,'',4,2] = 6.566
> Area_Stand[B,'',5,2] = 58.724
> Area_Stand[B,'',4,3] = 69.777
> Area_Stand[C,'',1,3] = 96.944
> Area_Stand[C,'',2,3] = 23.426
> Area_Stand[C,'',3,3] = 43.471
> Area_Stand[C,'',1,1] = 208.631
> Area_Stand[C,'',2,1] = 245.654
> Area_Stand[C,'',3,1] = 43.55
> Area_Stand[C,'',4,2] = 6.934
> Area_Stand[C,'',5,2] = 163.692
> Area_Stand[C,'',6,2] = 72.195
> Area_Stand[C,'',7,2] = 23.371
> Area_Stand[C,'',8,2] = 72.455
> Generating Area...
> ReadingDbfGIS_7.mod:41: no value for Area_Stand[A,R,1,1]
> MathProg model processing error
> ~~~~~~~~~~~~~~~~~~~~~
>
> My problem are in here:
>
> # Problem here
> param A{i in I}:= sum {j in J, k in K, l in L} Area_Stand[i,j,k,l];
>
>
> # Problem here as well.
> subject to Area {i in I}:
>    sum {j in J,k in K,l in L} XForestLand[i,j,k,l] = sum {j in J,k in
> K,l in L} Area_Stand[i,j,k,l];
>
> How can I put conditions on these two equation e.g if statements in
> these to equation.
>
> Please see the mathprog script and console output below.
>
> Thanks in advance, Noli
> ~~~~~~~~~~~~~~~~~~~~~
>
> Mathprog script
> ~~~~~~
> #Reading dbf (*,dbf) from shapefile (*..shp), GIS file i.e. ArcGIS
> file and writing dbf table.
>
> /*
> Noli Sicad <address@hidden>
>
> stands8.shp
> stands8.dbf
> stands8.shx
> stabds8.prj
>
> */
>
> set I;
> table tab IN "xBASE" "vegtype.dbf": I <- [VEG_TYPE];
> display I;
>
> set J;
> table tab2 IN "xBASE" "vegtype2.dbf": J <- [VEG_TYPE2];
> display J;
>
> set K;
> table tab3 IN "xBASE" "climatec.dbf": K <- [CLIMATE_C];
> display K;
>
> set L;
> table tab4 IN "xBASE" "MGT.dbf": L <- [MANAGEMENT];
> display L;
>
> param Area_Stand{i in I,j in J,k in K,l in L};
>
> set S, dimen 4;
> table tab8 IN "xBASE" "stands8.dbf": S <- [VEG_TYPE,
> VEG_TYPE2,CLIMATE_C,MANAGEMENT], Area_Stand ~ACRES;
> display Area_Stand;
>
> var XForestLand{i in I,j in J,k in K,l in L};
>
> # Problem here
> param A{i in I}:= sum {j in J, k in K, l in L} Area_Stand[i,j,k,l];
>
> # Problem here as well.
> subject to Area {i in I}:
>    sum {j in J,k in K,l in L} XForestLand[i,j,k,l] = sum {j in J,k in
> K,l in L} Area_Stand[i,j,k,l];
>
> ~~~~~~~~~~~~~~~
> Console output
>
> ~~~~~~~~~~~~~~~
> Parameter(s) specified in the command line:
>  --cover --clique --gomory --mir -m ReadingDbfGIS_7.mod
> Reading model section from ReadingDbfGIS_7.mod...
> ReadingDbfGIS_7.mod:78: warning: unexpected end of file; missing end
> statement inserted
> 78 lines were read
> Reading tab...
> Display statement at line 12
> I:
>    A
>    B
>    C
>    D
> Reading tab2...
> Display statement at line 16
> J:
>    R
>    ''
> Reading tab3...
> Display statement at line 20
> K:
>    1
>    2
>    3
>    4
>    5
>    6
>    7
>    8
> Reading tab4...
> Display statement at line 24
> L:
>    1
>    2
>    3
> Reading tab8...
> Display statement at line 30
> Area_Stand[A,R,1,3] = 131.507
> Area_Stand[A,R,2,3] = 122.955
> Area_Stand[A,'',1,1] = 57.899
> Area_Stand[A,'',2,1] = 15.573
> Area_Stand[A,'',3,1] = 78.544
> Area_Stand[A,'',4,1] = 67.906
> Area_Stand[A,'',5,1] = 145.412
> Area_Stand[A,R,1,2] = 155.223
> Area_Stand[A,'',2,2] = 64.358
> Area_Stand[A,'',3,2] = 35.798
> Area_Stand[A,R,4,2] = 124.597
> Area_Stand[B,'',1,3] = 46.048
> Area_Stand[B,'',2,3] = 73.329
> Area_Stand[B,'',3,3] = 40.189
> Area_Stand[B,'',1,1] = 4.184
> Area_Stand[B,'',2,1] = 15.175
> Area_Stand[B,'',3,1] = 52.294
> Area_Stand[B,'',1,2] = 78.506
> Area_Stand[B,'',2,2] = 32.534
> Area_Stand[B,'',3,2] = 22.573
> Area_Stand[B,'',4,2] = 6.566
> Area_Stand[B,'',5,2] = 58.724
> Area_Stand[B,'',4,3] = 69.777
> Area_Stand[C,'',1,3] = 96.944
> Area_Stand[C,'',2,3] = 23.426
> Area_Stand[C,'',3,3] = 43.471
> Area_Stand[C,'',1,1] = 208.631
> Area_Stand[C,'',2,1] = 245.654
> Area_Stand[C,'',3,1] = 43.55
> Area_Stand[C,'',4,2] = 6.934
> Area_Stand[C,'',5,2] = 163.692
> Area_Stand[C,'',6,2] = 72.195
> Area_Stand[C,'',7,2] = 23.371
> Area_Stand[C,'',8,2] = 72.455
> Generating Area...
> ReadingDbfGIS_7.mod:41: no value for Area_Stand[A,R,1,1]
> MathProg model processing error
>>Exit code: 1    Time: 0.209
>