Re: GLPK and OPENMP - GLP_SIMPLEX running concurrently

2020-05-13 Thread Domingo Alvarez Duarte
Hello Sierra ! My 2 cents, have you tried to compile with -fsanitize=thread to see if it catch the problem ? https://github.com/google/sanitizers/wiki/ThreadSanitizerCppManual Cheers ! On 13/5/20 14:50, Sierra Ansuas, Juan Pablo wrote: Hello, and sorry for insisting about this issue.

Re: Adding if/then/else statement to GMPL

2020-08-30 Thread Domingo Alvarez Duarte
https://ampl.com/BOOK/EXAMPLES/EXAMPLES2/cut2.run and the data file is at: https://ampl.com/BOOK/EXAMPLES/EXAMPLES2/cut.dat -Original Message- From: Domingo Alvarez Duarte Sent: Monday, August 24, 2020 10:34 AM To: Meketon, Marc ; Andrew Makhorin ; help-glpk@gnu.org Subject: Re: Adding if/t

Try experimental GMPL extensions online

2020-09-05 Thread Domingo Alvarez Duarte
Hello ! Try it online here https://meimporta.eu/myglpk-ext After looking around and with the help of the author of https://www.strandmark.net/wasm/glpk.html I managed to compile https://github.com/mingodad/GLPK/tree/local-set-param to wasm and make it available here

Added sqlite3 driver fro tables in GMPL

2020-10-08 Thread Domingo Alvarez Duarte
Hello ! I just added a sqlite3 driver for tables in GMPL with this commit https://github.com/mingodad/GLPK/commit/66a2ff9ec7fd5e2506cfebceca4bab95203b9739 , it also include the updated examples in examples/sql . And in doing so I discovered a possible use after free in src/mpl/mpsql.c for

GLPK/glpsol gives different objective values for mip/nomip

2020-10-06 Thread Domingo Alvarez Duarte
Hello ! Testing GLPK with miplib2017 I noticed that GLPK/glpsol gives different values for LP objective funcion value when a MIP problem is given with the option --nomip , is this the expected behavior ? Solving MIP -> *  6593: obj =   4.383346402e+07 inf = 1.946e-09 (0) 29 Soling "--nomip"

How hard would be to add "Indicator Constraints"

2020-10-05 Thread Domingo Alvarez Duarte
Hello ! While testing GLPK with miplib2017 I noticed that several problems can not be handled by GLPK due to having "Indicator Constraints" in the mps file. Looking here https://www.gams.com/latest/docs/UG_LanguageFeatures.html#UG_LanguageFeatures_IndicatorConstraints for a explanation of

Adding if/then/else statement to GMPL

2020-08-23 Thread Domingo Alvarez Duarte
Hello ! I just added the if/then/else statement to GMPL here https://github.com/mingodad/GLPK/commit/8984d5db70ce0cc91189911164f5448390026a4a and hope it could be added to a future release of GLPK/GMPL. Any comment/suggestion is welcome ! param p symbolic := "dad"; check length(p) ==

Re: Adding if/then/else statement to GMPL

2020-08-23 Thread Domingo Alvarez Duarte
Hello ! Also I've added the break/continue statements here https://github.com/mingodad/GLPK/commit/9d70a37b16bd377722eeb3880fcf86bb3b812118 Again any comment/suggestion is welcome ! Cheers !

Re: Adding if/then/else statement to GMPL

2020-08-25 Thread Domingo Alvarez Duarte
m,n) in B:     row >= k and row < m and (col = l or col = n) and     x[i,j,k,l,m,n] = 1}) > 0} "|";     printf{0..0: card({(i,j,k,l,m,n) in B:     row >= k and row < m and (col = l or col = n) and         x[i,j,k,l,m,n] = 1}) = 0} " ";  

Adding scoped local set/param declaration to GMPL

2020-08-25 Thread Domingo Alvarez Duarte
Hello ! I'm experimenting with allow local set/param declarations inside scoped blocks (for/if/then/else), it's the base to possibly experiment further with problem/callback/function declarations you can see/experiment with it in this branch

Multi problems in GMPL

2020-08-26 Thread Domingo Alvarez Duarte
Hello ! Thinking in how to have multiple problem in GMPL I'm asking for other user opinions about how that could be expressed, one example can be like the one bellow where we have common elements and namespaced problems. If you have any other idea that you want to share let me know !

GMPL/GLPK display objective function value

2020-08-28 Thread Domingo Alvarez Duarte
Hello ! While trying to implement multi solve statements I found that in GMPL the display of an objective function after solving do not show the optimal value. GLPSOL warnings: knapsack-3.glpsol:38: warning: unexpected end of file; missing end statement inserted But it shows conflicting

Re: GMPL/GLPK display objective function value

2020-08-29 Thread Domingo Alvarez Duarte
Hello Andrew ! Thanks for reply ! The problem is that the GMPL is not informed/knows anything about the warning/divergence in the results and goes blindly. Te exit code of glpsol also doesn't reflect anything that could be detected, only by visual inspection of the output or redirect the

Re: GMPL/GLPK display objective function value

2020-08-29 Thread Domingo Alvarez Duarte
Hello again ! My bad, I've got really confused, the second result is correct and expected with the removed "shift" term. Sorry by the noise ! Cheers ! On 29/8/20 12:25, Domingo Alvarez Duarte wrote: Hello again ! Doing more tests with this now I've got confused again, I comm

Re: Adding if/then/else statement to GMPL

2020-08-29 Thread Domingo Alvarez Duarte
://ampl.com/BOOK/EXAMPLES/EXAMPLES2/cut.dat -Original Message- From: Domingo Alvarez Duarte Sent: Monday, August 24, 2020 10:34 AM To: Meketon, Marc ; Andrew Makhorin ; help-glpk@gnu.org Subject: Re: Adding if/then/else statement to GMPL Hello Meketon ! Could you share your view of how

Re: GMPL/GLPK display objective function value

2020-08-29 Thread Domingo Alvarez Duarte
on row 0     High quality KKT.DE: max.abs.err = 0.00e+00 on column 0     max.rel.err = 0.00e+00 on column 0     High quality KKT.DB: max.abs.err = 0.00e+00 on row 0     max.rel.err = 0.00e+00 on row 0     High quality End of output ==== Cheers ! On 29/8/20 11:32, Domingo Alvare

Re: GMPL/GLPK display objective function value

2020-08-29 Thread Domingo Alvarez Duarte
term %.12g ig"    "nored\n",    mpl_get_row_name(tran, i), mpl_get_row_c0(tran, i));             ret = WARNING_CONST_TERM_IGNORED; } ... return ret; } Cheers ! On 29/8/20 11:01, Domingo Alvarez Duarte wrote: Hello Andrew ! Thanks for

Re: GMPL/GLPK display objective function value

2020-08-29 Thread Domingo Alvarez Duarte
. Cheers ! On 28/8/20 20:40, Domingo Alvarez Duarte wrote: Hello ! While trying to implement multi solve statements I found that in GMPL the display of an objective function after solving do not show the optimal value. GLPSOL warnings: knapsack-3.glpsol:38: warning: unexpected end of file; m

Re: Adding if/then/else statement to GMPL

2020-08-27 Thread Domingo Alvarez Duarte
AMPLES/EXAMPLES2/cut2.mod https://ampl.com/BOOK/EXAMPLES/EXAMPLES2/cut2.run and the data file is at: https://ampl.com/BOOK/EXAMPLES/EXAMPLES2/cut.dat -Original Message- From: Domingo Alvarez Duarte Sent: Monday, August 24, 2020 10:34 AM To: Meketon, Marc ; Andrew Makhorin ; help-glpk@g

Re: Adding if/then/else statement to GMPL

2020-08-27 Thread Domingo Alvarez Duarte
_intopt(csa->prob, >iocp);          ...          ret = glp_mpl_postsolve(csa->tran, csa->prob, GLP_SOL);          //resume the execution      } = Cheers ! On 27/8/20 11:25, Andrew Makhorin wrote: On Thu, 2020-08-27 at 11:01 +0200, Domingo Alvarez Duarte wrote: Hello ! I just finis

Re: Adding if/then/else statement to GMPL

2020-08-27 Thread Domingo Alvarez Duarte
Hello Andrew ! After replying to your last email I came with this possible solution to allow multiple solve statements in GMPL: https://github.com/mingodad/GLPK/commit/b80bcd86f1a78eba753e919a56fb9b976934b237 Cheers ! On 27/8/20 11:44, Domingo Alvarez Duarte wrote: Hello Andrew ! Thanks

Re: Multi problems in GMPL

2020-08-26 Thread Domingo Alvarez Duarte
ent at line 39 problem Cutting_Opt: Cut, Number, Fill; problem Pattern_Gen: Use, Reduced_Cost, Width_Limit; problem Mix: Cut, Reduced_Cost; Model has been successfully generated >Exit code: 0 Cheers ! On 26/8/20 19:57, Domingo Alvarez Duarte wrote: Hello ! Thinking in how to have mu

Re: Adding if/then/else statement to GMPL

2020-08-24 Thread Domingo Alvarez Duarte
i,jj) in S: (i-ii,j-jj) in D[c] and x[c,ii,jj]} printf" %s", substr(c,1,1); } else printf" ."; } printf"\n"; } Cheers ! On 24/8/20 16:59, Heinrich Schuchardt wrote: On 24.08.20 16:33, Domingo Alvarez Duarte wrote: Hello Meketon ! Could you share your view of how

Re: Adding if/then/else statement to GMPL

2020-08-24 Thread Domingo Alvarez Duarte
Makhorin Sent: Sunday, August 23, 2020 2:56 PM To: Domingo Alvarez Duarte ; help-glpk@gnu.org Subject: Re: Adding if/then/else statement to GMPL On Sun, 2020-08-23 at 15:36 +0200, Domingo Alvarez Duarte wrote: Hello ! Also I've added the break/continue statements here https://github.com/mingodad/GLPK

Re: GLPSOL in webassemby faster than native ?

2020-09-23 Thread Domingo Alvarez Duarte
e glp_double double     #define GLP_DBL_EPSILON DBL_EPSILON     #define GLP_DBL_FMT_G "g"     #define GLP_DBL_FMT_F "f"     #define GLP_DBL_FMT_E "E"     #define GLP_DBL_FMT_e "e"     #endif #endif Cheers ! On 22/9/20 21:49, Mich

Re: GLPSOL in webassemby faster than native ?

2020-09-22 Thread Domingo Alvarez Duarte
at 16:09 +0200, Domingo Alvarez Duarte wrote: Hello Andrew ! Are you saying that floating point calculations are more efficient/precise in webassembly ? No. I meant that due to floating-point computations running the same computer program with the same data as a rule produces different results

Re: GLPSOL in webassemby faster than native ?

2020-09-22 Thread Domingo Alvarez Duarte
generations selected ? Isn't strange that on wasm it's been faster than native ? Doesn't this difference gives insight to select the parameter differently ? Cheers ! On 22/9/20 17:56, Andrew Makhorin wrote: On Tue, 2020-09-22 at 15:53 +0200, Domingo Alvarez Duarte wrote: Hello again ! On an And

Re: GLPSOL in webassemby faster than native ?

2020-09-24 Thread Domingo Alvarez Duarte
22/9/20 21:49, Michael Hennebry wrote: On Tue, 22 Sep 2020, Domingo Alvarez Duarte wrote: Due to the big difference in solver time how could we figure out what's is it in order to use this knowledge to improve the native solver time ? I mean what debug/verbose options could help us have a

Re: GLPSOL in webassemby faster than native ?

2020-09-25 Thread Domingo Alvarez Duarte
; are equal). But also leave the possibility of have it working (somehow) with float32, float64, float80, float128, ... everywhere. Again help/suggestions/comments are welcome ! Cheers ! On 24/9/20 21:18, Michael Hennebry wrote: On Thu, 24 Sep 2020, Domingo Alvarez Duarte wrote: I just got glp

Re: GLPSOL in webassemby faster than native ?

2020-09-25 Thread Domingo Alvarez Duarte
erformance decreased, hashi.mod takes too much time, but defining glp_long_double as glp_double we still have the same results/output as the original GLPK === Cheers ! On 24/9/20 21:18, Michael Hennebry wrote: On Thu, 24 Sep 2020, Domingo Alvarez Duarte wrote: I just got glpsol with "long

Re: GLPSOL in webassemby faster than native ?

2020-09-26 Thread Domingo Alvarez Duarte
| | | " | |   | | | |   | " |   |   " 1---3 | | " 1 |   1 | | 4---3 " |   |   "     | | | "   | | | "   | " |   4===5-2   | | 2-6===6-3 | "   | 3 |   | |   | |   |   " | | "   | | 2  

Re: Solver performance solving examples/life_goe.mod

2020-09-26 Thread Domingo Alvarez Duarte
I'm not very happy with it, so it isn't committed to my tree. Best Regards, Chris Matrakids On Sat, 26 Sep 2020 at 17:55, Domingo Alvarez Duarte mailto:mingo...@gmail.com>> wrote: Hello Chris ! Thank you for reply ! Can you describe in a few lines what your improvements achieve

Re: Solver performance solving examples/life_goe.mod

2020-09-26 Thread Domingo Alvarez Duarte
improving patches a few years ago. You can find them in: https://github.com/cmatraki/GLPK Best Regards, Chris Matrakidis On Sat, 26 Sep 2020 at 14:54, Domingo Alvarez Duarte mailto:mingo...@gmail.com>> wrote: Hello ! Testing GLPK I left it solving examples/lie_goe.mod

Re: GLPSOL in webassemby faster than native ?

2020-09-26 Thread Domingo Alvarez Duarte
   /* compute new d[q], which is the reduced cost of xB[p] in the    * adjacent basis */   d[q] = (dq /= tcol_vec[p]);   /* compute new d[j] for all j != q */   for (k = 1; k <= trow_nnz; k++)   {  j = trow_ind[k]; if (j != q)         d[j] -= trow_vec[j] * dq;    

Re: GLPSOL in webassemby faster than native ?

2020-09-25 Thread Domingo Alvarez Duarte
ct graph... Conflict graph has 530 + 20 = 550 vertices +   739: mip = not found yet >=  -inf (1; 0) Cuts on level 0: gmi = 10; mir = 28; cov = 33; clq = 3; !* not the same = Cheers ! On 22/9/20 17:56, Andrew Makhorin wrote: On Tue, 2020-09-22 at 15:53 +0200,

Re: GLPSOL in webassemby faster than native ?

2020-09-22 Thread Domingo Alvarez Duarte
s could give hints to improve GLPK solver performance by inspecting the decision's criteria and eventually find a better ones. Anyone can give any idea with this data ? Cheers ! On 21/9/20 17:11, Andrew Makhorin wrote: On Mon, 2020-09-21 at 16:09 +0200, Domingo Alvarez Duarte wrote: Hello Andrew !

Re: GLPSOL in webassemby faster than native ?

2020-09-22 Thread Domingo Alvarez Duarte
time to solve the MIP. On Tue, 22 Sep 2020, Domingo Alvarez Duarte wrote: Hello Andrew ! Due to the big difference in solver time how could we figure out what's is it in order to use this knowledge to improve the native solver time ? I mean what debug/verbose options could help us have a clue

Re: GLPSOL in webassemby faster than native ?

2020-10-02 Thread Domingo Alvarez Duarte
ase28i.mps.gz    0.0009    0:00.09    9384 === Cheers ! On 2/10/20 17:59, Michael Hennebry wrote: On Thu, 1 Oct 2020, Domingo Alvarez Duarte wrote: But in reality it seems that musl/qsort results in "stable sort" but actually for hashi.mod and tiling.mod the order of the "tie-brea

Re: Try experimental GMPL extensions online

2020-09-21 Thread Domingo Alvarez Duarte
2:55, Karl Ove Hufthammer wrote: Domingo Alvarez Duarte skreiv 05.09.2020 13:24: After looking around and with the help of the author of https://www.strandmark.net/wasm/glpk.html I managed to compile https://github.com/mingodad/GLPK/tree/local-set-param to wasm and make it available

GLPSOL in webassemby faster than native ?

2020-09-21 Thread Domingo Alvarez Duarte
Hello ! I'm a bit lost here, I did some changes to GLPK/GMPL here https://github.com/mingodad/GLPK/tree/local-set-param and compiled it to webassembly and made it availlable here https://meimporta.eu/myglpk-ext/ then a user (Karl) pointed out that it was segfaulting when using the "--cuts"

Re: GLPSOL in webassemby faster than native ?

2020-09-21 Thread Domingo Alvarez Duarte
Hello Andrew ! Are you saying that floating point calculations are more efficient/precise in webassembly ? Cheers ! On 21/9/20 15:08, Andrew Makhorin wrote: Does someone can give a possible explanation ? floating-point computations

Re: Is this code correct in the GMPL implementation ?

2020-08-01 Thread Domingo Alvarez Duarte
, Domingo Alvarez Duarte wrote: Hello ! Trying to change the actual TUPLE structure by: #ifndef TUPLE_SIZE_T #define TUPLE_SIZE_T int #endif struct TUPLE {    TUPLE_SIZE_T size;    TUPLE_SIZE_T refcount;    SYMBOL sym[1];    /* symbol, which the component refers to; cannot

Re: How to retrieve intermediate solutions of MIP problems without interrupting the solver?

2020-08-10 Thread Domingo Alvarez Duarte
Hello Yuri ! Following the code again I found this lines in examples/glpsol.c: #ifdef GLP_CB_FUNC /* 05/IV-2016 */ {  extern void GLP_CB_FUNC(glp_tree *, void *);     csa->iocp.cb_func = GLP_CB_FUNC;     csa->iocp.cb_info = NULL; } #endif It seems

Re: How to retrieve intermediate solutions of MIP problems without interrupting the solver?

2020-08-10 Thread Domingo Alvarez Duarte
Hello Heinrich ! Thank you for reply ! My bad, I did followed the code instead of reading the manual ! There is a old saying: When everything else fails, read the  manual ! Cheers ! On 10/8/20 11:56, Heinrich Schuchardt wrote: On 10.08.20 11:47, Domingo Alvarez Duarte wrote: Hello Yuri

Re: How to retrieve intermediate solutions of MIP problems without interrupting the solver?

2020-08-10 Thread Domingo Alvarez Duarte
Hello Yuri ! After my last reply I followed the code again and found that it seems to already have a way to do it using src/proxy/proxy.h but I couldn't find a way to set our own proxy function/wrapper for it. Cheers ! On 9/8/20 20:04, Yuri wrote: For my problem Glpk quickly finds some

Is this code correct in the GMPL implementation ?

2020-07-25 Thread Domingo Alvarez Duarte
Hello ! Trying to change the actual TUPLE structure by: #ifndef TUPLE_SIZE_T #define TUPLE_SIZE_T int #endif struct TUPLE {   TUPLE_SIZE_T size;   TUPLE_SIZE_T refcount;   SYMBOL sym[1];   /* symbol, which the component refers to; cannot be NULL */ }; I found the

Re: Excessive copies of set elements in GMPL

2020-07-17 Thread Domingo Alvarez Duarte
Hello Andrew ! Thank you again for reply ! Running all GMPL models in the examples folder and comparing it's output I found that the ones that have "SET UNION" operations was failing, reviewing the code I found that only "src/mpl/mpl3.c::set_union" was not making a working copy to return,

Re: Excessive copies of set elements in GMPL

2020-07-16 Thread Domingo Alvarez Duarte
Hello Andrew ! Thank you for reply ! After doing some experimentation I think I found an easy way to eliminate unnecessary set copies here is a commit on https://github.com/mingodad/GLPK/commit/6f3da6ab31ca8d710f706f60635e5b17cf2ded40, also see bellow, glpsol now uses 1/3 of the memory and

Re: Excessive copies of set elements in GMPL

2020-07-16 Thread Domingo Alvarez Duarte
#with refcount & splaytree GLPSOL: GLPK LP/MIP Solver, v4.65 ... Memory used: 169.0 Mb (177220286 bytes) 1.76user 0.09system 0:01.85elapsed 100%CPU (0avgtext+0avgdata 176188maxresident)k On 16/7/20 18:40, Domingo Alvarez Duarte wrote: Hello Andrew ! Thank you for reply ! After doing

Re: Excessive copies of set elements in GMPL

2020-07-21 Thread Domingo Alvarez Duarte
Hello Michael ! Do you have any example doing what you are talking about ? In theory it should work. Right now it process all models in the "examples" folder and the output is identical to the original GLPK/GMPL except that uses less memory and is slightly faster. Cheers ! On 21/7/20

Re: Excessive copies of set elements in GMPL

2020-07-24 Thread Domingo Alvarez Duarte
0)    (8,9,11)    (9,10,12)    (10,11,13) Cheers ! On 23/7/20 18:40, Michael Hennebry wrote: On Tue, 21 Jul 2020, Domingo Alvarez Duarte wrote: Do you have any example doing what you are talking about ? In theory it should work. Right now it process all models in the "

Re: Excessive copies of set elements in GMPL

2020-07-24 Thread Domingo Alvarez Duarte
uot;, card(C); printf "D count: %d\n", card(D); display D; = On 24/7/20 9:05, Domingo Alvarez Duarte wrote: Here is the output of the script from the previous message: Memory usage: ubuntu glpsol package : 4,544 KB glpk with my changes: 3,268 KB /usr/bin/time glpsol -m test-nes

Re: Excessive copies of set elements in GMPL

2020-07-24 Thread Domingo Alvarez Duarte
100%CPU (0avgtext+0avgdata 3268maxresident)k 0inputs+0outputs (0major+214minor)pagefaults 0swaps On 24/7/20 8:58, Domingo Alvarez Duarte wrote: Hello Michael ! What kind of problem do you expect ? Other than memory size requirements ? Do you mean something like this: set

Re: Excessive copies of set elements in GMPL

2020-07-16 Thread Domingo Alvarez Duarte
Hello Andrew ! Thank you for reply ! I'm aware of it and that's why I've added an assert in "add_tuple" and so far it was not triggered running several models including "osemosys" where the output file "SelectedResults.csv" are identical. Do you have any tests or know someone that has then

Huge memory usage difference

2020-08-16 Thread Domingo Alvarez Duarte
Hello ! I just found a class of models that can exhibit a huge memory usage difference depending on if we add or not add parameter missing values when there is a default value declaration for then. See this model for example

Re: Huge memory usage difference

2020-08-16 Thread Domingo Alvarez Duarte
Hello Andrew ! Thank you for reply ! So testing to detect simple expressions can make it usable in those cases ? Can you give one such example where it'll fail ? Cheers ! On 16/8/20 16:44, Andrew Makhorin wrote: On Sun, 2020-08-16 at 11:53 +0200, Domingo Alvarez Duarte wrote: Hello ! I

Re: Change Path Input/Output file

2020-11-25 Thread Domingo Alvarez Duarte
Hello Manuel ! Something like this should work. table tab_plant IN "CSV" "different_diretory/plants.csv" : I <- [plant], a ~ capacity; Cheers ! On 24/11/20 20:11, Manuel Castro wrote: Hi there, I am currently using gusek as an IDE to model and run GLPK. I use the following statement to

Re: glpk 5.0 release information

2020-12-16 Thread Domingo Alvarez Duarte
I've got a the mailing list link wrongly. Cheers ! On 16/12/20 18:48, Andrew Makhorin wrote: On Wed, 2020-12-16 at 16:59 +0100, Domingo Alvarez Duarte wrote: It seems that the problems reported on the mailing list and it's fixes are not in this release: https://github.com/mingodad/GLPK/commit

Re: glpk 5.0 release information

2020-12-16 Thread Domingo Alvarez Duarte
It seems that the problems reported on the mailing list and it's fixes are not in this release: https://github.com/mingodad/GLPK/commit/6c9c1de01e10a99de89a94746bad656dcc4f9838 https://github.com/mingodad/GLPK/commit/356faa9710ec7375e49615bd794258c914b9ae63

Re: [Fwd: Portable version]

2020-11-02 Thread Domingo Alvarez Duarte
Hello Benhard ! Probably the easy way is to static build/link glpsol. I normally look at the end of the make output and manually replace it: examples$ make gcc -DHAVE_CONFIG_H -I. -I..  -I./../src   -g -O3 -march=native -ffp-contract=off -DWITH_SPLAYTREE -flto -MT glpsol.o -MD -MP -MF

Re: [Fwd: What is this format?]

2021-01-26 Thread Domingo Alvarez Duarte
It seems to be a cplex OPL format. After manually convert it to GMPL (see bellow) I'm getting a segfault with ubuntu 18.04 glpk distribution: glpsol -m cplex2mpl.mod GLPSOL: GLPK LP/MIP Solver, v4.65 Parameter(s) specified in the command line:  -m cplex2mpl.mod Reading model section

Re: [Fwd: What is this format?]

2021-01-26 Thread Domingo Alvarez Duarte
Where is GLPK 6.0 available ? It's the first time I hear of it, and it's not available at https://ftp.gnu.org/gnu/glpk/ . On 26/1/21 15:57, Andrew Makhorin wrote: On Tue, 2021-01-26 at 13:33 +0100, Domingo Alvarez Duarte wrote: It seems to be a cplex OPL format. After manually convert

Re: [Fwd: Constraint not respected]

2021-06-22 Thread Domingo Alvarez Duarte
Without a full model and data don't expect much help ! On 22/6/21 11:20, Andrew Makhorin wrote: Forwarded Message *Date*: Tue, 22 Jun 2021 09:25:01 +0200 *Subject*: Constraint not respected *To*: help-glpk@gnu.org *From*: Gioia Lorusso

Re: Multicore LP solvers

2021-03-27 Thread Domingo Alvarez Duarte
Hello Reginald ! Have looked at https://github.com/ERGO-Code/HiGHS they seem to be doing a parallel LP solver. Cheers ! On 26/3/21 21:26, Reginald Beardsley wrote: I haven't fooled around with GLPK and LP problems in general for several years now. The appearance of off lease machines with

Re: Solver delivers wrong answer when 2 constraints are close

2021-03-04 Thread Domingo Alvarez Duarte
Testing this problem I discover that if we change the order of constraint declarations it seems to give the expected answer as stated by Thiago (what I think could be another bug). /param min_bound default 0;/ /var x >= 0; minimize y: x;/ * */*s.t. PART_MIN_X: x >= 1 + min_bound;* /

Re: Solver delivers wrong answer when 2 constraints are close

2021-03-04 Thread Domingo Alvarez Duarte
Hello Antti ! Thanks for the explanation ! What's your opinion of a possible solution for this "design flaw" ? Cheers ! On 4/3/21 16:44, Antti Lehtila wrote: I guess it may be a considered a design flaw, but I think it should not be called a bug, as it is working as designed and documented.

Re: Solver delivers wrong answer when 2 constraints are close

2021-03-04 Thread Domingo Alvarez Duarte
and documented.  Besides, I think one should use the Bounds section for bounds, instead of using multiple constraints for defining a single lower bound. Best, Antti ___ On 04.03.2021 11:38, Domingo Alvarez Duarte wrote: > > Testing

Re: Solver delivers wrong answer when 2 constraints are close

2021-03-04 Thread Domingo Alvarez Duarte
When I test the provided sample with AMPL I get the stated expected result independent of the order of the constraints. Doesn't seem reasonable to somehow fix the behavior of GLPK/GLPSOL ? Cheers ! On 4/3/21 18:19, Andrew Makhorin wrote: Both 1e-3 and 1e-5 are good enough. If the result is

Re: Solver delivers wrong answer when 2 constraints are close

2021-03-05 Thread Domingo Alvarez Duarte
bug, as it is working as designed and documented. Besides, I think one should use the Bounds section for bounds, instead of using multiple constraints for defining a single lower bound. Best, Antti ___ On 04.03.2021 11:38, Domingo Alvarez Duarte wrote: > > Testing this probl

Re: Solver delivers wrong answer when 2 constraints are close

2021-03-05 Thread Domingo Alvarez Duarte
tion for bounds, instead of using multiple constraints for defining a single lower bound. Best, Antti ___ On 04.03.2021 11:38, Domingo Alvarez Duarte wrote: > > Testing this problem I discover that if we change the order of > constraint declarations it seems to give

Re: Solver delivers wrong answer when 2 constraints are close

2021-03-04 Thread Domingo Alvarez Duarte
. Besides, I think one should use the Bounds section for bounds, instead of using multiple constraints for defining a single lower bound. Best, Antti ___ On 04.03.2021 11:38, Domingo Alvarez Duarte wrote: > > Testing this problem I discover that if we change the order of > constraint decla

Re: Fixing integer variables after MIP solution found

2021-08-27 Thread Domingo Alvarez Duarte
Hello Juan ! Can you provide an example of a problem you are trying to solve ? I did some extensions to GMPL/GLPK on this repository https://github.com/mingodad/GLPK/tree/local-set-param and you can see an online example here https://meimporta.eu/myglpk-ext/ select from "extended-samples" ->

Re: [Fwd: help for gusek]

2021-07-15 Thread Domingo Alvarez Duarte
There is a menu on "View->Paremeters" and there you can add "--tmlim" see "glsplo --help" about it. On 15/7/21 14:26, Andrew Makhorin wrote: Forwarded Message *Date*: Thu, 15 Jul 2021 13:27:24 +0200 *Subject*: help for gusek *To*: help-glpk@gnu.org

Re: [Fwd: gmpl question]

2021-12-17 Thread Domingo Alvarez Duarte
Hello Davor ! Based on your description and hints from Fouad maybe this does what you expect (you can try it online at https://meimporta.eu/myglpk-ext/): /* ASSIGN, Assignment Problem */ /* Written in GNU MathProg by Andrew Makhorin */ /* The assignment problem is one of the

Re: Defining a variable of mixed datatype in Mathprog

2022-03-15 Thread Domingo Alvarez Duarte
It would help if you give a problem description and your expected implementation for it because without then your initial question is too open/vague. On 15/3/22 7:26, Anku Pandey wrote: Hi, I have to formulate a problem in mathprog where a variable can be both integer type and continuous

Re: One-off error in the MIP solver?

2022-01-24 Thread Domingo Alvarez Duarte
Hello Harmut ! Testing it on https://meimporta.eu/myglpk-ext/ that's GLPK 4.65 plus several customizations gives this result: === v: 131072 v: 0001 v: 131071 e: 0 === I'm don't see a reason to upgrade right now. Cheers ! On 24/1/22 10:46, Hartmut Henkel wrote: Hi, here

Re: Public VCS

2022-04-26 Thread Domingo Alvarez Duarte
Hello Maxim ! I'm maintaining this https://github.com/mingodad/GLPK with all the versions that I found through the net and my changes to it too and have a online wasm demo here https://meimporta.eu/myglpk-ext . Cheers ! On 26/4/22 18:14, Maxim Akhmedov wrote: Hi, As long as I see there is

Re: GLPK memory problem

2022-05-03 Thread Domingo Alvarez Duarte
Maybe there is a mistake on the previous reply, GLPK is 32/64 bits depending on the compilation/platform/OS on 64 bits operating systems it can allocated more than 4GB of memory. Can you use gdb to see how much memory was attempted to allocate ? On 3/5/22 13:11, Andrew Makhorin wrote: On

Re: GLPK memory problem

2022-05-04 Thread Domingo Alvarez Duarte
:10, Andrew Makhorin wrote: On Tue, 2022-05-03 at 13:29 +0200, Domingo Alvarez Duarte wrote: Maybe there is a mistake on the previous reply, GLPK is 32/64 bits depending on the compilation/platform/OS on 64 bits operating systems it can allocated more than 4GB of memory. I meant that for indexing

Re: Requesting Repo access regd

2023-04-16 Thread Domingo Alvarez Duarte
I've done a collection of all versions available at https://ftp.gnu.org/gnu/glpk and created this repository https://github.com/mingodad/GLPK , that also has several changes made by me. On 13/4/23 9:42, Heinrich Schuchardt wrote: Am 13. April 2023 09:25:54 MESZ schrieb PRANJAL PRAGYA VERMA