1. Си / Говнокод #26067

    +2

    1. 001
    2. 002
    3. 003
    4. 004
    5. 005
    6. 006
    7. 007
    8. 008
    9. 009
    10. 010
    11. 011
    12. 012
    13. 013
    14. 014
    15. 015
    16. 016
    17. 017
    18. 018
    19. 019
    20. 020
    21. 021
    22. 022
    23. 023
    24. 024
    25. 025
    26. 026
    27. 027
    28. 028
    29. 029
    30. 030
    31. 031
    32. 032
    33. 033
    34. 034
    35. 035
    36. 036
    37. 037
    38. 038
    39. 039
    40. 040
    41. 041
    42. 042
    43. 043
    44. 044
    45. 045
    46. 046
    47. 047
    48. 048
    49. 049
    50. 050
    51. 051
    52. 052
    53. 053
    54. 054
    55. 055
    56. 056
    57. 057
    58. 058
    59. 059
    60. 060
    61. 061
    62. 062
    63. 063
    64. 064
    65. 065
    66. 066
    67. 067
    68. 068
    69. 069
    70. 070
    71. 071
    72. 072
    73. 073
    74. 074
    75. 075
    76. 076
    77. 077
    78. 078
    79. 079
    80. 080
    81. 081
    82. 082
    83. 083
    84. 084
    85. 085
    86. 086
    87. 087
    88. 088
    89. 089
    90. 090
    91. 091
    92. 092
    93. 093
    94. 094
    95. 095
    96. 096
    97. 097
    98. 098
    99. 099
    100. 100
    /* ../COBYLA2/cobyla2.f -- translated by f2c (version 20100827) */
    ...
    /* Subroutine */ int cobyla_(integer *n, integer *m, doublereal *x, 
    	doublereal *rhobeg, doublereal *rhoend, integer *iprint, integer *
    	maxfun, doublereal *w, integer *iact)
    {
        static integer ia, idx, mpp, icon, isim, isigb, idatm, iveta, isimi, 
    	    ivsig, iwork;
        extern /* Subroutine */ int cobylb_(integer *, integer *, integer *, 
    	    doublereal *, doublereal *, doublereal *, integer *, integer *, 
    	    doublereal *, doublereal *, doublereal *, doublereal *, 
    	    doublereal *, doublereal *, doublereal *, doublereal *, 
    	    doublereal *, doublereal *, integer *);
        /* Parameter adjustments */
        --iact;
        --w;
        --x;
    
        /* Function Body */
        mpp = *m + 2;
        icon = 1;
        isim = icon + mpp;
        isimi = isim + *n * *n + *n;
        idatm = isimi + *n * *n;
        ia = idatm + *n * mpp + mpp;
        ivsig = ia + *m * *n + *n;
        iveta = ivsig + *n;
        isigb = iveta + *n;
        idx = isigb + *n;
        iwork = idx + *n;
        cobylb_(n, m, &mpp, &x[1], rhobeg, rhoend, iprint, maxfun, &w[icon], &w[
    	    isim], &w[isimi], &w[idatm], &w[ia], &w[ivsig], &w[iveta], &w[
    	    isigb], &w[idx], &w[iwork], &iact[1]);
        return 0;
    } /* cobyla_ */
    ...
    /* Subroutine */ int cobylb_(integer *n, integer *m, integer *mpp, doublereal 
    	*x, doublereal *rhobeg, doublereal *rhoend, integer *iprint, integer *
    	maxfun, doublereal *con, doublereal *sim, doublereal *simi, 
    	doublereal *datmat, doublereal *a, doublereal *vsig, doublereal *veta,
    	 doublereal *sigbar, doublereal *dx, doublereal *w, integer *iact)
    {
        /* Format strings */
        static char fmt_10[] = "(/3x,\002The initial value of RHO is\002,1pe13.6"
    	    ",2x,\002and PARMU is set to zero.\002)";
        static char fmt_50[] = "(/3x,\002Return from subroutine COBYLA because t"
    	    "he \002,\002MAXFUN limit has been reached.\002)";
        static char fmt_70[] = "(/3x,\002NFVALS =\002,i5,3x,\002F =\002,1pe13.6,"
    	    "4x,\002MAXCV =\002,1pe13.6/3x,\002X =\002,1pe13.6,1p4e15.6)";
        static char fmt_80[] = "(1pe19.6,1p4e15.6)";
        static char fmt_210[] = "(/3x,\002Return from subroutine COBYLA because"
    	    " \002,\002rounding errors are becoming damaging.\002)";
        static char fmt_410[] = "(/3x,\002Increase in PARMU to\002,1pe13.6)";
        static char fmt_580[] = "(/3x,\002Reduction in RHO to\002,1pe13.6,\002  "
    	    "and PARMU =\002,1pe13.6)";
        static char fmt_590[] = "(/3x,\002Normal return from subroutine COBYL"
    	    "A\002)";
    
        /* System generated locals */
        integer sim_dim1, sim_offset, simi_dim1, simi_offset, datmat_dim1, 
    	    datmat_offset, a_dim1, a_offset, i__1, i__2, i__3;
        doublereal d__1, d__2;
    
        /* Builtin functions */
        integer s_wsfe(cilist *), do_fio(integer *, char *, ftnlen), e_wsfe(void);
        double sqrt(doublereal);
    
        /* Local variables */
        static doublereal f;
        static integer i__, j, k, l, mp, np, iz;
        static doublereal phi, rho, sum, beta, cmin, cmax;
        static integer ivmc;
        static doublereal weta;
        static integer ivmd;
        static doublereal temp, wsig, gamma;
        static integer iflag;
        static doublereal alpha, delta, denom, tempa, barmu;
        static integer nbest, ifull, iptem, jdrop;
        static doublereal ratio, vmold, parmu, error, vmnew;
        extern /* Subroutine */ int calcfc_(integer *, integer *, doublereal *, 
    	    doublereal *, doublereal *);
        static integer ibrnch;
        static doublereal edgmax, pareta, prerec, phimin, parsig;
        static integer isdirn, nfvals, izdota;
        static doublereal cvmaxm, dxsign, prerem;
        static integer iptemp;
        static doublereal resmax, cvmaxp;
        static integer idxnew;
        static doublereal resnew, trured;
        extern /* Subroutine */ int trstlp_(integer *, integer *, doublereal *, 
    	    doublereal *, doublereal *, doublereal *, integer *, integer *, 
    	    doublereal *, doublereal *, doublereal *, doublereal *, 
    	    doublereal *, doublereal *);
    
        /* Fortran I/O blocks */
        static cilist io___22 = { 0, 6, 0, fmt_10, 0 };
        static cilist io___29 = { 0, 6, 0, fmt_50, 0 };
        static cilist io___33 = { 0, 6, 0, fmt_70, 0 };
        static cilist io___34 = { 0, 6, 0, fmt_80, 0 };
        static cilist io___39 = { 0, 6, 0, fmt_210, 0 };

    Найдено в зависимостях «SlangTNG». Библиотека «COBYLA», автоматически переведённая с «Фортрана».

    /* This subroutine minimizes an objective function F(X) subject to M */
    /* inequality constraints on X, where X is a vector of variables that has */
    /* N components. The algorithm employs linear approximations to the */
    /* objective and constraint functions, the approximations being formed by */
    /* linear interpolation at N+1 points in the space of the variables. */

    https://tng.tuxfamily.org/index.php?title=Main_Page

    Запостил: XYPO3BO3, 05 Декабря 2019

    Комментарии (10) RSS

    • SEO-пост.

      #COBYLA #fortran #f2c
      Ответить
    • Переведи на "PHP".
      Ответить
      • <?php
        function cobyla_($n,$m,$x,$rhobeg,$rhoend,$iprint,$maxfun,$w,$iact)
        {
            /* Parameter adjustments */
            --$iact;
            --$w;
            --$x;
            /* Function Body */
            $mpp = $m[0] + 2;
            $icon = 1;
            $isim = $icon + $mpp;
            $isimi = $isim + $n[0] * $n[0] + $n[0];
            $idatm = $isimi + $n[0] * $n[0];
            $ia = $idatm + $n[0] * $mpp + $mpp;
            $ivsig = $ia + $m[0] * $n[0] + $n[0];
            $iveta = $ivsig + $n[0];
            $isigb = $iveta + $n[0];
            $idx = $isigb + $n[0];
            $iwork = $idx + $n[0];
            cobylb_($n, $m, $mpp, $x[1], $rhobeg, $rhoend, $iprint, $maxfun, $w[$icon], $w[
        	    $isim], $w[$isimi], $w[$idatm], $w[$ia], $w[$ivsig], $w[$iveta], $w[
        	    $isigb], $w[$idx], $w[$iwork], $iact[1]);
            return 0;
        } /* cobyla_ */
        ...
        function cobylb_($n,$m,$mpp,$x,$rhobeg,$rhoend,$iprint,$maxfun,$con,$sim,$simi,$datmat,$a,$vsig,$veta,$sigbar,$dx,$w,$iact)
        {
            /* Format strings */
            $fmt_10 = "(/3x,\002The initial value of RHO is\002,1pe13.6"
        	    ",2x,\002and PARMU is set to zero.\002)";
            $fmt_50 = "(/3x,\002Return from subroutine COBYLA because t"
        	    "he \002,\002MAXFUN limit has been reached.\002)";
            $fmt_70 = "(/3x,\002NFVALS =\002,i5,3x,\002F =\002,1pe13.6,"
        	    "4x,\002MAXCV =\002,1pe13.6/3x,\002X =\002,1pe13.6,1p4e15.6)";
            $fmt_80 = "(1pe19.6,1p4e15.6)";
            $fmt_210 = "(/3x,\002Return from subroutine COBYLA because"
        	    " \002,\002rounding errors are becoming damaging.\002)";
            $fmt_410 = "(/3x,\002Increase in PARMU to\002,1pe13.6)";
            $fmt_580 = "(/3x,\002Reduction in RHO to\002,1pe13.6,\002  "
        	    "and PARMU =\002,1pe13.6)";
            $fmt_590 = "(/3x,\002Normal return from subroutine COBYL"
        	    "A\002)";
        
            /* Fortran I/O blocks */
            $cilist io___22 = [0, 6, 0, $fmt_10, 0];
            $cilist io___29 = [0, 6, 0, $fmt_50, 0];
            $cilist io___33 = [0, 6, 0, $fmt_70, 0];
            $cilist io___34 = [0, 6, 0, $fmt_80, 0];
            $cilist io___39 = [0, 6, 0, $fmt_210, 0];
        Ответить
    • Это какой-то матан для пидорашек.
      Ответить
      • Хуле тебе надо, пидоркшка?
        Ответить
      • ААААА!

        https://en.wikipedia.org/wiki/COBYLA

        «Constrained optimization by linear approximation (COBYLA) is a numerical optimization method for constrained problems where the derivative of the objective function is not known, invented by Michael J. D. Powell.»

        Математический алгоритм под названием «кобыла». Эта фигню ещё добавили в питоновский пакет «scipy»:
        https://stackoverflow.com/questions/12781622/
        Ответить
        • def fmin_cobyla(func, x0, cons, args=(), consargs=None, rhobeg=1.0,
                          rhoend=1e-4, maxfun=1000, disp=None, catol=2e-4):
              """
              Minimize a function using the Constrained Optimization BY Linear
              Approximation (COBYLA) method. This method wraps a FORTRAN
              implementation of the algorithm.
          ...
          
          def _minimize_cobyla(fun, x0, args=(), constraints=(),
                               rhobeg=1.0, tol=1e-4, maxiter=1000,
                               disp=False, catol=2e-4, **unknown_options):
              """
              Minimize a scalar function of one or more variables using the
              Constrained Optimization BY Linear Approximation (COBYLA) algorithm.
          Ответить

    Добавить комментарий