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

    +147

    1. 01
    2. 02
    3. 03
    4. 04
    5. 05
    6. 06
    7. 07
    8. 08
    9. 09
    10. 10
    ...
char det(char a[4][4])
{ float det;
det=a[0][0]*(a[1][1]*(a[2][2]*a[3][3]-a[3][2]*a[2][3])-a[1][2]*(a[2][1]*a[3][3]-a[2][3]*a[3][1])+a[1][3]*(a[2][1]*a[3][2]-a[2][2]*a[3][1]))
-a[0][1]*(-a[1][0]*(a[2][2]*a[3][3]-a[3][2]*a[2][3])-a[1][2]*(a[2][0]*a[3][3]-a[3][0]*a[2][3])+a[1][3]*(a[2][0]*a[3][2]-a[3][0]*a[2][2]))
+a[0][2]*(-a[1][0]*(a[2][1]*a[3][3]-a[3][1]*a[2][3])+a[1][1]*(a[2][0]*a[3][3]-a[3][0]*a[2][3])+a[1][3]*(a[2][0]*a[3][1]-a[2][1]*a[3][0]))
-a[0][3]*(-a[1][0]*(a[2][1]*a[3][2]-a[3][1]*a[2][2])+a[1][1]*(a[2][0]*a[3][2]-a[3][0]*a[2][2])-a[1][2]*(a[2][0]*a[3][1]-a[3][0]*a[2][1]));
return(det);
}; 
...

    http://otvet.mail.ru/question/59918103/

    Запостил: 1_and_0, 23 Июня 2011

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

    • ava Dummy00001 23.06.2011 18:50 # −5

      показать все, что скрытонормальная такая себе оптимизация вычисления определителя.

      если это говно, то повально все реализации алгоритмов кодирования и хэширования ( например http://people.csail.mit.edu/rivest/Md5.c ) тоже говно.
      Ответить

      • ava eth0 23.06.2011 18:59 # 0

        А они, фактически, и есть.
        Ответить

        • ava Dummy00001 23.06.2011 19:21 # +2

          для алгоритма которым пользуются может быть миллионы раз в секунду было бы говном *НЕ* воспользоватся даже самой мелкой доступной оптимизацией.
          Ответить

          • ava eth0 23.06.2011 20:03 # 0

            Алгоритмам хэширования в частных случаях оптимизация противопоказана.
            Ответить

      • ava macGovno 23.06.2011 19:14 # 0

        Если минор переписать через макрос, то будет уже лучше.
        Ответить

        • ava guest 23.06.2011 22:59 # 0

          C++ macros are evil.
          Ответить

          • ava macGovno 23.06.2011 23:29 # +1

            Где тут С++ ?
            Тем более что речь об оптимизации любой ценой.
            Ответить

            • ava guest 23.06.2011 23:36 # 0

              >Где тут С++ ?
              Ох, мне мерещится... :(
              Ответить

            • ava guest 25.06.2011 12:28 # 0

              А тут ни слова не сказано об оптимизации любой ценой.

              Просто школяр не знал циклов.
              Ответить

      • ava wvxvw 23.06.2011 23:41 # 0

        На лицо отсутсвие оптимизации, т.как куча повторных вычислений, которые, компилятор, наверняка сам соптимизирует, но, как ни странно, получится более читаемый код. Т.е. ну не будет он по нескольку раз пересчитывать a[1][1], а один раз посчитает и запишет. Человеку тоже было бы так удобнее прочитать, но, если лень подумать, и действительно соптимизировать, то получается такое.
        Ответить

        • ava Dummy00001 24.06.2011 00:18 # 0

          a[i][j] в случае матрицы 4х4 - 16. и куда по твоему компилятор эти 16 значений сможет записать? если например регистров хорошо если 8 свободных есть?

          современный компилятор при виде кода выше навярняка прикрутит какие CPU-специфичные SIMD инструкции. соптимизировать руками то можно - но оно будет непортабельно и я сомневаюсь что быстрее того что компилятор сам может сделать.
          Ответить

          • ava wvxvw 24.06.2011 00:57 # 0

            А куда он будет записывать индексы массива? Какой-то странный вопрос... Нам же не нужны все значения из матрицы одновременно...
            Смысл ответа заключался не в том, как компилятор будет оптимизировать (откуда вы знаете, может он RTL сначала создает, и его, или похожий, и его оптимизирует - техник же разных море, но не об том речь), а в том, что код не оптимизирован / оптимизация явно не являлась тут приоритетом.
            Ответить

            • ava Dummy00001 24.06.2011 01:06 # 0

              > код не оптимизирован

              ёпрст. а цыклы сами собой развернулись??
              Ответить

              • ava wvxvw 24.06.2011 01:24 # 0

                Оптимизация, это когда после переделывания кода он начинает делать лучше то, для чего его оптимизировали. Я не знаю, как нужно было хуже написать тот же код, чтобы это было его оптимизированой формой :S
                В смысле, само по себе разворачивание цикла может ничего и не дать. Ну, опять же, код не дает никаких оснований предполагать, что автор что-то замерял / проверял, и пришел к такому решению потому, что оно было оптимальнее какого-то другого. В коде нет даже намека на такие действия. Просто человек написал первое, что ему пришло в голову, и нет тут никакого скрытого смысла...
                Ответить

                • ava SmackMyBitchUp 24.06.2011 15:10 # 0

                  >Оптимизация, это когда после переделывания кода он начинает делать лучше то, для чего его оптимизировали.
                  %\
                  Ответить

                  • ava wvxvw 24.06.2011 15:20 # 0

                    Ну можно оптимизировать длину кода, а можно время выполнения, а можно объем занимаемой памяти. Т.е. если вы оптимизируете какой-то код, то вы его не просто оптимизиуете, а для чего-то. Ну это если в этом был вопрос...
                    Ответить

        • ava Esper 24.06.2011 12:04 # +1

          >Т.е. ну не будет он по нескольку раз пересчитывать a[1][1], а один раз посчитает и запишет.
          Между делом. В этом коде каждый минор второго порядка считается 2 раза, многие из них - неявно, нужно пораскрывать скобочки, перегруппировать. Просто поищите глазами выражение "(a[2][2]*a[3][3]-a[3][2]*a[2][3])".
          А если оптимизировать, то я бы взял Лапласа и SSE-инструкции.
          Ответить

    • ava Esper 23.06.2011 19:20 # +5

      Т.е. он берет матрицу из char'ов (байт со знаком в наиболее распространенном случае). Потом считает определитель (кстати произведение char'ов - это по умолчанию char или уже short или int?), преобразует его во float (single float, 32-битный с 24-битной мантиссой в наиболее распространенном случае), а потом - обратно к char'у. Но произведение 4 знаковых char'ов требует 4*(8-1)=28 бит. Т.е. из-за преобразования к типу с плавающей точкой потеряна точность.
      Это победа.
      Ответить

      • ava guest 23.06.2011 23:01 # 0

        Хотите сказать, double решил бы проблему?
        Ответить

      • ava gegMOPO4 25.06.2011 14:17 # +1

        char промоутится до int, так что произведение-то наверняка посчитает (если где char 9-битный, то там и int 36-битный). Произведение 4 знаковых 8-битных чисел требует 4*(8-1)+1 бит на знак. Да ещё при сложении/вычитании 24 чисел нужно до 5 бит, чтобы не было переполнения. Итого — на результат 34 бита.
        Ответить

        • ava Esper 28.06.2011 20:47 # 0

          Справедливо от и до.
          My bad: я увидел только часть ИСТИНЫ.
          Ответить

    • ava Esper 24.06.2011 12:10 # +6

      Тут ко всему еще косяки со знаками, если я не забыл линейную алгебру, то должно быть:
      det=a[0][0]*(a[1][1]*(a[2][2]*a[3][3]-a[3][2]*a[2][3])-a[1][2]*(a[2][1]*a[3][3]-a[2][3]*a[3][1])+a[1][3]*(a[2][1]*a[3][2]-a[2][2]*a[3][1]))
-a[0][1]*(a[1][0]*(a[2][2]*a[3][3]-a[3][2]*a[2][3])-a[1][2]*(a[2][0]*a[3][3]-a[3][0]*a[2][3])+a[1][3]*(a[2][0]*a[3][2]-a[3][0]*a[2][2]))
+a[0][2]*(a[1][0]*(a[2][1]*a[3][3]-a[3][1]*a[2][3])-a[1][1]*(a[2][0]*a[3][3]-a[3][0]*a[2][3])+a[1][3]*(a[2][0]*a[3][1]-a[2][1]*a[3][0]))
-a[0][3]*(a[1][0]*(a[2][1]*a[3][2]-a[3][1]*a[2][2])-a[1][1]*(a[2][0]*a[3][2]-a[3][0]*a[2][2])+a[1][2]*(a[2][0]*a[3][1]-a[3][0]*a[2][1]));

      В общем: полное беспросветное говно, не имеющее никаких оправданий.
      Ответить

    • ava Lure Of Chaos 24.06.2011 17:09 # +2

      а вот мне нужен детерминант матрицы 3х3 и 10х10. неуниверсально - печалько...
      Ответить

    • ava TheHamstertamer 25.06.2011 10:57 # 0

      сделайте меня развидеть это
      Ответить

    • ava guest 26.06.2011 22:45 # 0

      google> site:otvet.mail.ru 1_and_0
      Ответить

    • ava guest8 08.04.2019 20:59 # −999

      показать все, что скрытоvanished
      Ответить

    • ava guest8 09.04.2019 11:02 # −999

      показать все, что скрытоvanished
      Ответить

    • ava guest8 12.04.2019 14:25 # −999

      показать все, что скрытоvanished
      Ответить

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