X
PrevPrev Go to previous topic
NextNext Go to next topic
Last Post 10 Sep 2013 08:48 AM by  anon
problem with IMSL
 2 Replies
Sort:
You are not authorized to post a reply.
Author Messages

anon



New Member


Posts:
New Member


--
10 Sep 2013 08:48 AM
    When I run the model of IMSL, the problem appears % IMSL Error: IMSL_LINLSQ: Fatal error: MATH_TOO_MANY_ITN: Maximum number of iterations exceeded. % Execution halted at: MESMA 115 E:\PHD-UWM\study\IDL\mesma.pro % READSPECTRA 70 E:\PHD-UWM\study\IDL\readspectra.pro % $MAIN$ IDL> and my codes are b =[0.192377,0.201996,0.221146 ,0.235080,0.322395,0.241512,0.204528,0.0305424] a =[[0.0853720,0.0615450,0.0364250,0.295134,0.123997,0.0387040,-0.146001,0.779606], [0.0862940,0.0711650,0.0418080,0.283938,0.156914,0.0579790,-0.143244,0.742956], [0.0925340,0.0865570,0.0583260,0.316509,0.206240,0.0888180,-0.149263,0.681642], [0.218756,0.241758,0.261851,0.302259,0.340694,0.257684,0.241602,0.0752340]] c4 = [[1.00000], [1.00000], [1.00000], [1.00000]] bc =[1.0000000] contype = [0] xlb4 =[0.000000,0.000000,0.000000,0.000000] ite = 10 xub4 =[1,1,1,1] sol = IMSL_LINLSQ(b, a, c4, bc, bc, contype, Xlb = xlb4,ITMAX=ite, $ ; a is the spectra library, a set of spectra Xub = xub4, Residual = residual1) ; sol is the percentage of each endmember So, could anybody can help me with this problem? Thank you very much.

    Deleted User



    New Member


    Posts:
    New Member


    --
    10 Sep 2013 05:17 PM
    It is possible that for your particular data, you will need to use a =See if this makes sense for your particular data. Maybe you will need to use other constraints as well, I'm not sure about it. I hope this helps. Fernando

    Deleted User



    New Member


    Posts:
    New Member


    --
    11 Sep 2013 02:56 PM
    As I was saying, knowing your data will make the difference to understand why certain constrains fail. Also, notice that there are other parameters that you can play with when using ISML_LINSQL(), like ABS_TOLERANCE or DOUBLE or even: ITMAX Set the maximum number of iterations. Default: ITMAX = 5*max(nra, nca) Cheers, fernando
    You are not authorized to post a reply.

    Forums Products IDL