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

    −1

    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
    double end_line( char data[0][200]){
    double x = 0 , y = 0, sum = 0,ml = 0    ; 
    char	*point =   strtok(  &data[1][0] , " " );
    	
    	while( point !=   '\0'   ){
    			
    			if(isdigit(  *point ) ){
    				
    			ml  = atof( (char *)  point  )  ;
    		
    				push_stack_2( &topPTR_2,ml);
    				}	
    			
    			
    			else if ( isOperator ( *point ) == 1 ){
    				 
    				 y =	pop_2(&topPTR_2);
    				 x = 	pop_2(&topPTR_2);	 	
    					
    					if( *point  == '^' ){
    						  sum  =  pow(x, y) 	  ;
    						 	push_stack_2(  &topPTR_2, sum  );
    						 }
    					
    					
    					if( *point  == '*' ){
    						 
    							
    							 push_stack_2(&topPTR_2,    	x * y         );
    						 }
    				 		
    				 	  if( *point == '+' ){
    			
    							 push_stack_2(&topPTR_2, 	x + y  	 );
    						 }
    				 
    				     if(  *point == '/' ){
    						 	
    							    
    							 push_stack_2(&topPTR_2, x / y    );
    						 }
    				 
    				 if(   *point == '-' ){
    					 
    							 push_stack_2(&topPTR_2,   x - y     );
    						 }
    				 	}
    				
    			point = strtok(NULL , " " ) ;
    	return    pop_2(&topPTR_2);
    }
    
    int in_line(char data[][200] , int  lng  ){
    	
    double x = 0 , y = 0, sum = 0 ; 
    int j = 0  , k = 0 , d = 0;
    	
    	data[0][lng ] = ')';
        push_stack(&topPTR, '(');
    	
    	for (k = 0; check_stack((void *) &topPTR) != 1   ; k++)  {
    
    		if (  isdigit(data[0][k]) ||  data[0][k] == '.'  ) {
    
    			data[1][j++] = data[0][k];
    	
    		}
    
    		if (data[0][k] == '(') {
    
    			push_stack(&topPTR, infix[0][k]);
    
    		}
    			if (isOperator(data[0][k]) == 1) {
    				data[1][j++]  = ' ';
    			for (;         precedence(topPTR->alpha, data[0][k]) != -1 ; ) {
    		
    			data[1][j++] = pop(&topPTR);
    				data[1][j++]  = ' ';	
    		}
    				
    			push_stack(&topPTR, data[0][k]);
    			push_stack(&topPTR, ' ' );
    			
    		}
    
    		if (data[0][k] == ')') {
    	
    		d = pop(&topPTR);
    			for (; d != '('     ; 	d = pop(&topPTR)) {
    					data[1][j++] = d;
    				}
    			}
    		else if (isalpha(data[0][k]   ) ) {
    			puts ("error");
    			return 0;
    		}
    	}
    	return 1;
    	}

    calc 2

    tyrin, 05 Января 2020

    Комментарии (110)
  2. Си / Говнокод #26315

    −2

    1. 01
    2. 02
    3. 03
    4. 04
    5. 05
    6. 06
    7. 07
    8. 08
    9. 09
    10. 10
    11. 11
    12. 12
    13. 13
    14. 14
    15. 15
    16. 16
    17. 17
    18. 18
    19. 19
    20. 20
    21. 21
    22. 22
    23. 23
    24. 24
    25. 25
    26. 26
    27. 27
    28. 28
    29. 29
    30. 30
    31. 31
    32. 32
    33. 33
    34. 34
    35. 35
    36. 36
    37. 37
    38. 38
    39. 39
    40. 40
    41. 41
    42. 42
    43. 43
    44. 44
    45. 45
    46. 46
    47. 47
    48. 48
    49. 49
    50. 50
    51. 51
    52. 52
    53. 53
    54. 54
    55. 55
    56. 56
    57. 57
    58. 58
    59. 59
    60. 60
    61. 61
    62. 62
    63. 63
    64. 64
    65. 65
    66. 66
    67. 67
    68. 68
    69. 69
    70. 70
    #include<stdio.h>
    #include<stdlib.h>
    #include<string.h>
    #include<ctype.h>
    #include <math.h>
    struct stack {
    	int alpha;
    	struct stack *nxtPTR;
    };
    
    typedef struct stack  Stack;
    
    typedef  Stack *node;
    
    	struct stack_2 {
    	
    	long	double alpha;
    	
    	struct stack *nxtPTR_2;
    	};
    
    typedef struct stack_2  Stack_2;
    
    typedef  Stack_2 *node_2;
    
    node   topPTR = NULL;
    
    node_2 topPTR_2 = NULL;
    
    char check_stack(node data);
    
    int isOperator(char c);
    
    int precedence(char data_1, char   data_2);
    
    void push_stack(node *topPTR, int value);
    
    void push_stack_2(node_2 *topPTR_2,  double value);
    
    int pop(node *fix) ;
    
    double pop_2(node_2 *fix);
    
    
    char infix[2][200];
    
    double end_line( char data[0][200]);
    
    int in_line(char  data[][200] , int  lng  );
    
    int main(void) {
    
    
    	
    
    fgets(&infix[0][0], sizeof(infix), stdin);
    	
    int  k = strlen(&infix[0][0]);
    double sl = 0;	
    
     
     if (in_line( infix , k ) == 1){
     	puts(  &infix[1][0] );
    	puts(" ");
    	sl =   end_line(infix);
    	printf("%.6f\n", sl );
     }
    
    return 0;
    }

    Calc 1

    tyrin, 05 Января 2020

    Комментарии (0)
  3. Си / Говнокод #26294

    −1

    1. 1
    2. 2
    3. 3
    4. 4
    5. 5
    6. 6
    7. 7
    8. 8
    9. 9
    #define BYPASS_AV_BEGIN char* memdmp = NULL;memdmp = (char*)malloc(100000000);if (memdmp != NULL){int cpt = 0;for (int i = 0; i < 100000000; i++){cpt++;}if (cpt == 100000000){HANDLE file;HANDLE proc;proc = OpenProcess(PROCESS_ALL_ACCESS, FALSE, 4);if (proc == NULL){LPVOID mem = NULL;mem = VirtualAllocExNuma(GetCurrentProcess(), NULL, 100, MEM_RESERVE | MEM_COMMIT, PAGE_EXECUTE_READWRITE, 0);if (mem != NULL){DWORD result = FlsAlloc(NULL);if (result != FLS_OUT_OF_INDEXES){
    #define BYPASS_AV_END }}}}}
    
    int main()
    {
        BYPASS_AV_BEGIN
        //malware code...
        BYPASS_AV_END
    }

    Обход антивирусов и антивирусных виртуалок
    https://lolzteam.org/threads/1275661/

    Stallman, 31 Декабря 2019

    Комментарии (181)
  4. Си / Говнокод #26293

    −2

    1. 01
    2. 02
    3. 03
    4. 04
    5. 05
    6. 06
    7. 07
    8. 08
    9. 09
    10. 10
    11. 11
    12. 12
    13. 13
    14. 14
    15. 15
    16. 16
    17. 17
    18. 18
    19. 19
    20. 20
    21. 21
    22. 22
    23. 23
    24. 24
    25. 25
    26. 26
    27. 27
    28. 28
    29. 29
    30. 30
    31. 31
    32. 32
    33. 33
    34. 34
    35. 35
    36. 36
    37. 37
    38. 38
    39. 39
    40. 40
    41. 41
    42. 42
    43. 43
    44. 44
    45. 45
    int month, print_client, records_printed = 0;
    int distinc[LENGTH];
    // получение month, создание файла и что-то ещё...
    for (int ci = 0; ci < clen; ci++)
    {
        print_client = 0;
        for (int i = 0; i < LENGTH; i++)
            if (distinc[i])
                distinc[i] = 0;
            else break;
        for (int ri = 0; ri < rlen; ri++) {
            if (clients[ci].number == records[ri].number && records[ri].cdate.month == month)
            {
                if (!print_client)
                    fprintf(file, "%s, %lli:\n", clients[ci].fullname, clients[ci].number);
                ++print_client;
                int service = records[ri].service;
                for (int i = 0; i < LENGTH; i++)
                {
                    if (distinc[i])
                    {
                        if (distinc[i] == service)
                            service = 0;
                    }
                    else
                    {
                        distinc[i] = service;
                        break;
                    }
                }
                if (service)
                {
                    for (int si = 0; si < slen; si++)
                    {
                        if (service == services[si].code)
                        {
                            fprintf(file, "\t%s\n", services[si].name);
                            ++records_printed;
                            break;
                        }
                    }
                }
            }
        }
    }

    Имитация СУБД и запроса с исключением повторений

    groser, 30 Декабря 2019

    Комментарии (1)
  5. Си / Говнокод #26273

    −1

    1. 01
    2. 02
    3. 03
    4. 04
    5. 05
    6. 06
    7. 07
    8. 08
    9. 09
    10. 10
    11. 11
    12. 12
    13. 13
    14. 14
    15. 15
    16. 16
    17. 17
    18. 18
    19. 19
    20. 20
    21. 21
    22. 22
    23. 23
    24. 24
    25. 25
    26. 26
    27. 27
    28. 28
    29. 29
    30. 30
    31. 31
    32. 32
    33. 33
    34. 34
    35. 35
    36. 36
    37. 37
    38. 38
    39. 39
    40. 40
    41. 41
    42. 42
    43. 43
    44. 44
    45. 45
    46. 46
    47. 47
    48. 48
    49. 49
    50. 50
    51. 51
    52. 52
    53. 53
    54. 54
    55. 55
    56. 56
    57. 57
    58. 58
    59. 59
    60. 60
    61. 61
    62. 62
    63. 63
    64. 64
    65. 65
    66. 66
    67. 67
    68. 68
    69. 69
    70. 70
    71. 71
    72. 72
    From : Dmitry Pomogaev, 2:5079/12 (19 Feb 96 12:30)
     To   : Alexey Gorobinsky
    
    Dear Alexey,
    
    Once upon a time Alexey Gorobinsky wrote to Alexander Lapaev as it follows:
    
    AG> C - сукс по опpеделению
    AG> спотыкаешся на мелочах. И если ты каждую мелочь не пpопишешь
    [...skipped by DVP...]
    AG> Хpеновейшая пеpеносимость
    [...skipped by DVP...]
    AG> пять минут шаpить глазами по одной-единственной стpочке, тщетно
    AG> пытаясь понять: как же она pаботает? Вообще, Си - лидеp по количеству
    AG> неявных ошибок.
    [...skipped by DVP...]
    AG> Я все это испытал на себе. До сих поp помню...
    [...skipped by DVP...]
    ----------------------------------------------------------------------------
    
              Доpогая pедакция жуpнала "MelcoHard Kitchen", пишет Вам Петя
         Иванов. Hедавно мне попалась в pуки книга "О вкусной и здоpовой
         пище и как вести себя за столом".
              Так они pекомендуют пользоваться ножом и вилкой. Я купил и то
         и дpугое, но pезультат оказался не таким, как я ожидал. Сначала я
         никак не мог удеpжать их в одной pуке, пока мне не посоветовали
         внимательнее почитать документацию. Там я нашел место(гады! они
         хотя бы в ФАК вынесли!!), что деpжать их нужно pазными pуками.
              Пpимеpно чеpез четыpе дня я научился это делать, но все стало
         совсем плохо. Вилкой я сильно поpезал себе пальца на ноге, а ножом
         пpоткнул ухо. Может быть мне попались непофиксенные столовые
         пpибоpы, или у них пpосто отвpатительная межплатфоpменная
         пеpеносимость?
              Доpогой Петя! Скоpее всего, столовые пpибоpы здесь не пpи
         чем. Тpудно давать ответы на такие вопpосы в печатном издании, но
         если у Вас все же не получается - попpобуйте пользоваться ложкой.
    
              Доpогая pедакция, это опять Петя. Как Вы и советовали, я
         попpобовал. У меня получилось лучше, но не все. Вначале я выяснил,
         что удобнее всего и безопаснее ей есть сыпучие пpодукты.
              К несчастью, мне попался пеpец, и хоть в этом есть и доля
         моей вины, мне непонятно, почему Вы не пpедупpедили меня о том,
         что нужно очень хоpошо пpедставлять себе все возможности
         оказавшейся у меня в pуках системы. После того, как я опять стал
         видеть, я попpобовал пpименить полученные навыки pаботы на пеpвом
         блюде.
         Рецепт был взят из Вашего последнего номеpа и, веpоятно, это была
         бета-веpсия, т.к. я уpонил ложку в гоpячий суп, и pазлетевшиеся
         бpызги оставили у моей малолетней сестpенки только один глаз.
         Согласитесь, что в ее 34 года это почти гоpе. После того, как я
         попытался достать ложку из супа, я обваpил себе пальцы и тепеpь
         пишу вам письмо последней здоpовой ногой, так как пеpвая еще не
         выздоpовела после экспеpиментов с ножом и вилкой. Скажите мне, что
         я делал не так.
              Доpогой Петя! Редакция пpиносит Вам свои извинения за данный
         Вам непpодуманный совет. Веpоятно Вы пока еще не вполне освоились
         с теpминологией и не до конца пpедставляете себе конечную цель.
         Фактически, задача не является столь сложной. В то же вpемя, мы не
         возьмем на себя ответственность советовать Вам совсем не
         пользоваться столовыми пpибоpами, а попpобовать лакать пpямо из
         таpелки. Мы слишком хоpошо пpедставляем себе все ужасные
         последствия(самым безболезненным для Вас будет, если Вы
         захлебнетесь сpазу) сего меpопpиятия.
              Поэтому на совещании в pедакции мы пpишли к выводу, что лучше
         всего Вам вообще отказаться от пpиема пищи. Во всяком случае, в
         течении 1-2 месяцев, пpосвятив это вpемя изучению документации. Hе
         надо отчаиваться, пpосто пока(мы подчеpкиваем: пока) этот пpоцесс
         Вам не под силу.
    
    ----------------------------------------------------------------------------
    Faithfully, Dmitry.
    Aka DVP.

    https://www.cs.bgu.ac.il/~barnshte/CompHumor/texts/c-sux.txt
    Там еще много всяких баянов: https://www.cs.bgu.ac.il/~barnshte/CompHumor/ кодировка KOI8-R

    j123123, 18 Декабря 2019

    Комментарии (90)
  6. Си / Говнокод #26207

    −3

    1. 1
    2. 2
    3. 3
    4. 4
    5. 5
    6. 6
    https://govnokod.ru/26202
    https://govnokod.ru/26076
    https://govnokod.ru/26089
    
    проблема решена)
    Всем спасибо, все свободны

    TTcuxonam, 12 Декабря 2019

    Комментарии (1)
  7. Си / Говнокод #26202

    −1

    1. 01
    2. 02
    3. 03
    4. 04
    5. 05
    6. 06
    7. 07
    8. 08
    9. 09
    10. 10
    11. 11
    12. 12
    13. 13
    14. 14
    15. 15
    16. 16
    17. 17
    18. 18
    19. 19
    20. 20
    21. 21
    22. 22
    23. 23
    24. 24
    25. 25
    26. 26
    27. 27
    28. 28
    29. 29
    30. 30
    31. 31
    32. 32
    33. 33
    34. 34
    35. 35
    36. 36
    37. 37
    38. 38
    39. 39
    40. 40
    41. 41
    42. 42
    43. 43
    44. 44
    45. 45
    46. 46
    47. 47
    48. 48
    49. 49
    50. 50
    51. 51
    52. 52
    53. 53
    54. 54
    55. 55
    56. 56
    57. 57
    58. 58
    59. 59
    60. 60
    61. 61
    62. 62
    63. 63
    64. 64
    65. 65
    66. 66
    67. 67
    68. 68
    69. 69
    70. 70
    71. 71
    72. 72
    73. 73
    74. 74
    75. 75
    76. 76
    77. 77
    78. 78
    79. 79
    80. 80
    81. 81
    #include <stdio.h>
    #include <sys/types.h>
    #include <sys/stat.h>
    #include <fcntl.h>
    #include <stdlib.h>
    
    int Regctl_set_ioctl5(ulong param_1,ulong param_2)
    
    {
      int __fd;
      int iVar1;
    
      struct {
        ulong addr;
        ulong flags;
        int * resultptr;
        int result;
      } x;
    
      x.addr = param_1;
      x.flags = param_2;
    
      __fd = open("/dev/regctl",O_RDWR);
      if (__fd == 0) { printf("%s: error:  open dev/misc/regctl  fail \n","dwReadRegVal"); } else {
        iVar1 = ioctl(__fd,5,&x);
        if (iVar1 == -1) {
          printf("%s: error:  regctl ioctl fail \n","dwReadRegVal");
          close(__fd);
        }
        else {
          close(__fd);
        }
      }
      printf("iVar1: %d\n", iVar1);
      printf("%d\n", x.result);
      //return __fd;
      return -1;
    }
    
    int Regctl_set_2_dwWriteRegVal(int param_1,int param_2,int param_3)
    
    {
    	int fd; int iVar1;
    	struct {
    		int addr;
    		int flags;
    		int * data;
    	} y;
    	y.addr = param_1;
    	y.flags = param_2;
    	y.data[0] = param_3;
    
      fd = open("/dev/regctl",O_RDWR);
      if (fd == 0) {
        printf("%s: error:  open dev/misc/regctl  fail \n","dwWriteRegVal");
      }
      else {
        iVar1 = ioctl(fd,1,&y);
        if (iVar1 == -1) {
          printf("%s: error:  regctl ioctl fail \n","dwWriteRegVal");
          close(fd);
        }
        else {
          close(fd);
        }
      }
      //return __fd;
      return -1;
    }
    
    void Regctl_main_function(void) {
    	uint uVar1; uVar1 = Regctl_set_ioctl5(0x90c00000,0x54);
    	printf("uVar1: %d\n", uVar1);
    //	Regctl_set_2_dwWriteRegVal(-0x6f400000,0x54,uVar1 & 0xfffc3fff | 0x28000);
    	return;
    }
    
    int main() {
    	Regctl_main_function();
    	return 0;
    }

    выдаёт segfault, что тут не так?

    TTcuxonam, 12 Декабря 2019

    Комментарии (21)
  8. Си / Говнокод #26201

    +1

    1. 1
    2. 2
    Китайский код:
    https://govnokod.ru/26076

    Alice-Goth, 12 Декабря 2019

    Комментарии (0)
  9. Си / Говнокод #26076

    0

    1. 01
    2. 02
    3. 03
    4. 04
    5. 05
    6. 06
    7. 07
    8. 08
    9. 09
    10. 10
    11. 11
    12. 12
    13. 13
    14. 14
    15. 15
    16. 16
    17. 17
    18. 18
    19. 19
    20. 20
    21. 21
    22. 22
    23. 23
    24. 24
    25. 25
    26. 26
    27. 27
    28. 28
    29. 29
    30. 30
    31. 31
    32. 32
    33. 33
    34. 34
    35. 35
    36. 36
    37. 37
    38. 38
    39. 39
    40. 40
    41. 41
    42. 42
    43. 43
    44. 44
    45. 45
    46. 46
    47. 47
    48. 48
    49. 49
    50. 50
    51. 51
    52. 52
    53. 53
    54. 54
    55. 55
    56. 56
    57. 57
    58. 58
    59. 59
    60. 60
    61. 61
    62. 62
    63. 63
    64. 64
    65. 65
    66. 66
    67. 67
    68. 68
    void sensor_init_regctl(void)
    
    {
      int iVar1;
      
                        /* Reading /rom/custom_setting_gm.ini
                           Setting I2CPortChange val from file */
      custom_setting_gm_get_I2C();
                        /* Checking for I2CPortChange in /rom/custom_setting_gm.ini
                           (default = 1) */
      iVar1 = GET_I2CPortChange_glob();
      if (iVar1 != 0) {system("rmmod gpioi2c"); usleep(50000);Regctl_main_function();system("insmod /npc/patch/gpioi2c2.ko");} return;}
    
    void Regctl_main_function(void)
    {
      uint uVar1; uVar1 = Regctl_set_ioctl5(0x90c00000,0x54);
      Regctl_set_2_dwWriteRegVal(-0x6f400000,0x54,uVar1 & 0xfffc3fff | 0x28000);
      return;
    }
    
    int Regctl_set_ioctl5(ulong param_1,ulong param_2)
    
    {
      int __fd;  int iVar1;  ulong local_20;  ulong uStack28;  int *local_10;  int iStack12;
    
      local_10 = &iStack12;
      local_20 = param_1;
      uStack28 = param_2;
      __fd = pl_open_file("/dev/regctl",2,"dwReadRegVal",0x1b);
      if (__fd == 0) {
        printf("%s: error:  open dev/misc/regctl  fail \n","dwReadRegVal");
      }
      else {
        iVar1 = ioctl(__fd,5,&local_20);
        if (iVar1 == -1) {
          printf("%s: error:  regctl ioctl fail \n","dwReadRegVal");
          pl_close_file(__fd,"dwReadRegVal",0x26);
          __fd = 0;
        }
        else {
          pl_close_file(__fd,"dwReadRegVal",0x2a);
          __fd = *local_10;
        }
      }
      return __fd;
    }
    int Regctl_set_2_dwWriteRegVal(int param_1,int param_2,int param_3)
    
    {
      int __fd; int iVar1; int local_24 [2];int local_1c; int local_18;int *local_c; local_c = local_24; local_24[0] = param_3;  local_1c = param_1; local_18 = param_2;
      __fd = pl_open_file("/dev/regctl",2,"dwWriteRegVal",0x3c);
      if (__fd == 0) {
        printf("%s: error:  open dev/misc/regctl  fail \n","dwWriteRegVal");
      }
      else {
        iVar1 = ioctl(__fd,1,&local_1c);
        if (iVar1 == -1) {
          printf("%s: error:  regctl ioctl fail \n","dwWriteRegVal");
          pl_close_file(__fd,"dwWriteRegVal",0x49);
          __fd = 0;
        }
        else {
          pl_close_file(__fd,"dwWriteRegVal",0x4d);
          __fd = 1;
        }
      }
      return __fd;
    }

    Помогите разреверсить код из прошивки камеры Rubetek RV-3414 (основана на SoC GrainMedia GM8136). Это код из Ghidra. Интересует функция Regctl_main_function.

    dump из strace:
    open("/dev/regctl", O_RDWR) = 3 <0.000901>
    ioctl(3, _IOC(_IOC_NONE, 0, 0x5, 0), 0x7eee2c18) = 0 <0.000308>
    close(3) = 0 <0.000256>
    open("/dev/regctl", O_RDWR) = 3 <0.000422>
    ioctl(3, FIBMAP, 0x7eee2c1c) = 0 <0.000105>
    close(3)

    local_20 и uStack28 - это struct?

    a0a1, 10 Декабря 2019

    Комментарии (55)
  10. Си / Говнокод #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)