<commit_before>#ifndef DOUBLECONTRACT_H
#define DOUBLECONTRACT_H

#include "Fastor/commons/commons.h"
#include "Fastor/simd_vector/SIMDVector.h"

namespace Fastor {

template<typename T, size_t M, size_t N>
FASTOR_INLINE T _doublecontract(const T* FASTOR_RESTRICT a, const T* FASTOR_RESTRICT b) {

    using V = SIMDVector<T,DEFAULT_ABI>;
    constexpr int size = M*N;
    constexpr int stride = V::Size;
    int i = 0;

    V vec_a, vec_b, vec_out;
    for (; i< ROUND_DOWN(size,stride); i+=stride) {
        vec_a.load(a+i);
        vec_b.load(b+i);
#ifndef __FMA__
        vec_out += vec_a*vec_b;
#else
        vec_out = fmadd(vec_a,vec_b,vec_out);
#endif
    }
    T scalar = static_cast<T>(0);<commit_msg>Minor change<commit_after>#ifndef DOUBLECONTRACT_H
#define DOUBLECONTRACT_H

#include "Fastor/commons/commons.h"
#include "Fastor/simd_vector/SIMDVector.h"

namespace Fastor {

template<typename T, size_t M, size_t N>
FASTOR_INLINE T _doublecontract(const T* FASTOR_RESTRICT a, const T* FASTOR_RESTRICT b) {

    using V = SIMDVector<T,DEFAULT_ABI>;
    constexpr int size = M*N;
    constexpr int stride = V::Size;
    int i = 0;

    V vec_a, vec_b, vec_out;
    for (; i< ROUND_DOWN(size,stride); i+=stride) {
        vec_a.load(a+i,false);
        vec_b.load(b+i,false);
#ifndef __FMA__
        vec_out += vec_a*vec_b;
#else
        vec_out = fmadd(vec_a,vec_b,vec_out);
#endif
    }
    T scalar = static_cast<T>(0);
{'fraction_non_alphanumeric': 0.12159329140461216, 'fraction_numerical': 0.002795248078266946, 'mean_word_length': 3.7733333333333334, 'pattern_counts': {'":': 0, '<': 11, '<?xml version=': 0, '>': 9, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'dir': 'git-commits-cleaned', 'id': '3198126', 'n_tokens_mistral': 593, 'n_tokens_neox': 570, 'n_words': 122}