fem_data.py

import numpy as np


class Setup(object):
    def __init__(self):
        self.nodes = np.array([
            [4.0000000e+001, 0.0000000e+000],
            [3.8000000e+001, 0.0000000e+000],
            [4.0000000e+001, 2.0000000e+000],
            [3.6000000e+001, 0.0000000e+000],
            [3.7283235e+001, 1.7316509e+000],
            [4.0000000e+001, 4.0000000e+000],
            [3.8384202e+001, 3.2723653e+000],
            [3.4000000e+001, 0.0000000e+000],
            [3.5061637e+001, 1.7597298e+000],
            [3.6361846e+001, 3.7748572e+000],
            [4.0000000e+001, 6.0000000e+000],
            [3.8258316e+001, 5.0078704e+000],
            [3.2000000e+001, 0.0000000e+000],
            [3.3000000e+001, 1.7320508e+000],
            [3.4070581e+001, 3.5205074e+000],
            [3.6537709e+001, 5.9966986e+000],
            [3.4983944e+001, 5.1019165e+000],
            [4.0000000e+001, 8.0000000e+000],
            [3.8266344e+001, 7.0013117e+000],
            [3.0000000e+001, 0.0000000e+000],
            [3.1000000e+001, 1.7320508e+000],
            [3.2000000e+001, 3.4641016e+000],
            [3.3145738e+001, 5.2175195e+000],
            [3.6535888e+001, 7.9997048e+000],
            [3.4439846e+001, 6.8920061e+000],
            [4.0000000e+001, 1.0000000e+001],
            [3.8267682e+001, 9.0002186e+000],
            [2.8000000e+001, 0.0000000e+000],
            [2.9000000e+001, 1.7320508e+000],
            [3.0000000e+001, 3.4641016e+000],
            [3.1000000e+001, 5.1961524e+000],
            [3.2153732e+001, 6.8375822e+000],
            [3.6535845e+001, 9.9999933e+000],
            [3.4743170e+001, 8.9819507e+000],
            [3.3127140e+001, 8.2968771e+000],
            [4.0000000e+001, 1.2000000e+001],
            [3.8267905e+001, 1.1000036e+001],
            [2.6000000e+001, 0.0000000e+000],
            [2.7000000e+001, 1.7320508e+000],
            [2.8000000e+001, 3.4641016e+000],
            [2.9000000e+001, 5.1961524e+000],
            [3.0037859e+001, 6.9469686e+000],
            [3.1231097e+001, 8.7966472e+000],
            [3.6535881e+001, 1.2000006e+001],
            [3.4793723e+001, 1.0996992e+001],
            [3.3032347e+001, 1.0049937e+001],
            [4.0000000e+001, 1.4000000e+001],
            [3.8267942e+001, 1.3000006e+001],
            [2.4000000e+001, 0.0000000e+000],
            [2.5000000e+001, 1.7320508e+000],
            [2.6000000e+001, 3.4641016e+000],
            [2.7000000e+001, 5.1961524e+000],
            [2.8000000e+001, 6.9282032e+000],
            [2.9044826e+001, 8.3773554e+000],
            [2.9713233e+001, 9.8705404e+000],
            [3.1224572e+001, 1.1109035e+001],
            [3.6479488e+001, 1.4070583e+001],
            [3.4397963e+001, 1.3137978e+001],
            [3.3116759e+001, 1.1718142e+001],
            [4.0000000e+001, 1.6000000e+001],
            [3.8240269e+001, 1.5061638e+001],
            [2.2000000e+001, 0.0000000e+000],
            [2.3000000e+001, 1.7320508e+000],
            [2.4000000e+001, 3.4641016e+000],
            [2.5000000e+001, 5.1961524e+000],
            [2.6000000e+001, 6.9282032e+000],
            [2.7000000e+001, 8.3769228e+000],
            [2.8017108e+001, 9.7676951e+000],
            [2.9048562e+001, 1.1301969e+001],
            [3.2132607e+001, 1.3178653e+001],
            [3.0059197e+001, 1.3072426e+001],
            [3.6225143e+001, 1.6361846e+001],
            [3.4841257e+001, 1.4987152e+001],
            [3.3101004e+001, 1.4826466e+001],
            [4.0000000e+001, 1.8000000e+001],
            [3.8268349e+001, 1.7283235e+001],
            [2.0000000e+001, 0.0000000e+000],
            [2.1000000e+001, 1.7320508e+000],
            [2.2000000e+001, 3.4641017e+000],
            [2.3000000e+001, 5.1961525e+000],
            [2.4000000e+001, 6.9282033e+000],
            [2.4999999e+001, 8.5294860e+000],
            [2.6004572e+001, 9.9786045e+000],
            [2.7009259e+001, 1.1334533e+001],
            [2.8007604e+001, 1.3075217e+001],
            [3.0997726e+001, 1.4808522e+001],
            [2.8997695e+001, 1.4808686e+001],
            [3.6727635e+001, 1.8384202e+001],
            [3.4992130e+001, 1.8258316e+001],
            [3.4002687e+001, 1.6538889e+001],
            [3.1997913e+001, 1.6539588e+001],
            [4.0000000e+001, 2.0000000e+001],
            [3.8000000e+001, 2.0000000e+001],
            [1.8000000e+001, 0.0000000e+000],
            [1.9000000e+001, 1.7320509e+000],
            [1.9999999e+001, 3.4641019e+000],
            [2.0999999e+001, 5.1961529e+000],
            [2.1999997e+001, 6.9282037e+000],
            [2.2999996e+001, 8.6602551e+000],
            [2.4000934e+001, 1.0106964e+001],
            [2.5001063e+001, 1.1470685e+001],
            [2.5999390e+001, 1.3074876e+001],
            [2.6988735e+001, 1.4818904e+001],
            [2.9988444e+001, 1.6560255e+001],
            [2.7934714e+001, 1.6621398e+001],
            [3.6000000e+001, 2.0000000e+001],
            [3.4000000e+001, 2.0000000e+001],
            [3.2997087e+001, 1.8269735e+001],
            [3.0990393e+001, 1.8288299e+001],
            [1.6000000e+001, 0.0000000e+000],
            [1.6999998e+001, 1.7320514e+000],
            [1.7999996e+001, 3.4641035e+000],
            [1.8999993e+001, 5.1961548e+000],
            [1.9999986e+001, 6.9282055e+000],
            [2.0999981e+001, 8.6602601e+000],
            [2.1999976e+001, 1.0663005e+001],
            [2.3399941e+001, 1.1677481e+001],
            [2.3997732e+001, 1.3074954e+001],
            [2.4986999e+001, 1.4819711e+001],
            [2.5933872e+001, 1.6612941e+001],
            [2.8942360e+001, 1.8390048e+001],
            [2.6679783e+001, 1.8629651e+001],
            [3.2000000e+001, 2.0000000e+001],
            [3.0000000e+001, 2.0000000e+001],
            [1.4000000e+001, 0.0000000e+000],
            [1.4999986e+001, 1.7320543e+000],
            [1.5999976e+001, 3.4641121e+000],
            [1.6999962e+001, 5.1961643e+000],
            [1.7999922e+001, 6.9282138e+000],
            [1.8999904e+001, 8.6602875e+000],
            [1.9999854e+001, 1.0392318e+001],
            [2.2398825e+001, 1.2859237e+001],
            [2.0742584e+001, 1.2295881e+001],
            [2.2995300e+001, 1.4811780e+001],
            [2.3987069e+001, 1.6564472e+001],
            [2.4839225e+001, 1.7975593e+001],
            [2.8049097e+001, 2.0192147e+001],
            [2.6173166e+001, 2.0761205e+001],
            [2.4593860e+001, 1.9724412e+001],
            [1.2000000e+001, 0.0000000e+000],
            [1.2999917e+001, 1.7320715e+000],
            [1.3999873e+001, 3.4641602e+000],
            [1.4999802e+001, 5.1962114e+000],
            [1.5999577e+001, 6.9282527e+000],
            [1.6999516e+001, 8.6604418e+000],
            [1.7999237e+001, 1.0392346e+001],
            [1.8997956e+001, 1.2396475e+001],
            [2.1058897e+001, 1.4392453e+001],
            [2.0036097e+001, 1.3384581e+001],
            [2.2018152e+001, 1.6540254e+001],
            [2.3126540e+001, 1.8257576e+001],
            [2.4444298e+001, 2.1685304e+001],
            [2.2339087e+001, 2.0477129e+001],
            [1.0000000e+001, 0.0000000e+000],
            [1.0999504e+001, 1.7321751e+000],
            [1.1999335e+001, 3.4644292e+000],
            [1.2998964e+001, 5.1964371e+000],
            [1.3997698e+001, 6.9284290e+000],
            [1.4997597e+001, 8.6613204e+000],
            [1.5996004e+001, 1.0392331e+001],
            [1.6986600e+001, 1.2134254e+001],
            [1.9344949e+001, 1.4453515e+001],
            [1.7728637e+001, 1.4112246e+001],
            [2.0012841e+001, 1.6468058e+001],
            [2.1082770e+001, 1.8335161e+001],
            [2.2928932e+001, 2.2928932e+001],
            [2.0511916e+001, 2.2425470e+001],
            [2.0217714e+001, 2.0159060e+001],
            [8.0000000e+000, 0.0000000e+000],
            [8.9970241e+000, 1.7327965e+000],
            [9.9965889e+000, 3.4659221e+000],
            [1.0994578e+001, 5.1974743e+000],
            [1.1987422e+001, 6.9292141e+000],
            [1.2988308e+001, 8.6663772e+000],
            [1.3978910e+001, 1.0391205e+001],
            [1.4924361e+001, 1.2183676e+001],
            [1.5781491e+001, 1.3939983e+001],
            [1.8260829e+001, 1.6160791e+001],
            [1.6352726e+001, 1.5785342e+001],
            [1.8671285e+001, 1.8598308e+001],
            [2.1685304e+001, 2.4444298e+001],
            [1.9774094e+001, 2.4644026e+001],
            [1.8273046e+001, 2.3150230e+001],
            [1.8483512e+001, 2.0959233e+001],
            [6.0000000e+000, 0.0000000e+000],
            [6.9821448e+000, 1.7365247e+000],
            [7.9830052e+000, 3.4741548e+000],
            [8.9715440e+000, 5.2019355e+000],
            [9.9309891e+000, 6.9326618e+000],
            [1.0944727e+001, 8.6957562e+000],
            [1.1887557e+001, 1.0378516e+001],
            [1.2571254e+001, 1.2481353e+001],
            [1.3921575e+001, 1.3877106e+001],
            [1.4409391e+001, 1.5705652e+001],
            [1.6347552e+001, 1.7793380e+001],
            [1.5313208e+001, 1.6740333e+001],
            [1.6583970e+001, 1.9889804e+001],
            [2.0761205e+001, 2.6173166e+001],
            [1.8624129e+001, 2.6664506e+001],
            [1.7984162e+001, 2.4845465e+001],
            [1.6564255e+001, 2.3986418e+001],
            [1.6532138e+001, 2.2021948e+001],
            [4.0000000e+000, 0.0000000e+000],
            [4.8928689e+000, 1.7588943e+000],
            [5.9188625e+000, 3.5192010e+000],
            [6.8496700e+000, 5.2189773e+000],
            [7.8469662e+000, 6.8352573e+000],
            [8.7499531e+000, 8.8677978e+000],
            [1.0117498e+001, 1.0195822e+001],
            [1.0469459e+001, 1.1764961e+001],
            [1.0538494e+001, 1.3854029e+001],
            [1.2255853e+001, 1.4864095e+001],
            [1.2175201e+001, 1.6992966e+001],
            [1.4143162e+001, 1.7676956e+001],
            [1.4471916e+001, 1.9316303e+001],
            [1.4363988e+001, 2.1048926e+001],
            [2.0192147e+001, 2.8049097e+001],
            [1.8390048e+001, 2.8942360e+001],
            [1.6620487e+001, 2.7932176e+001],
            [1.6611844e+001, 2.5930833e+001],
            [1.4819504e+001, 2.4986131e+001],
            [1.4747832e+001, 2.2995038e+001],
            [2.0000000e+000, 0.0000000e+000],
            [2.7270378e+000, 1.7270378e+000],
            [3.6381610e+000, 3.7633806e+000],
            [5.0135757e+000, 5.0969083e+000],
            [5.5595337e+000, 6.8895816e+000],
            [6.8644935e+000, 8.2919653e+000],
            [6.9469686e+000, 9.9621411e+000],
            [8.6129481e+000, 1.1037041e+001],
            [8.6602589e+000, 1.3000000e+001],
            [8.6602831e+000, 1.5000000e+001],
            [1.0393761e+001, 1.5999816e+001],
            [1.0400837e+001, 1.7998896e+001],
            [1.2419438e+001, 1.8959087e+001],
            [1.3372973e+001, 2.0011413e+001],
            [1.2732049e+001, 2.2325493e+001],
            [1.2236548e+001, 2.0721339e+001],
            [2.0000000e+001, 3.0000000e+001],
            [1.8288299e+001, 3.0990393e+001],
            [1.6560299e+001, 2.9988481e+001],
            [1.4808702e+001, 2.8997708e+001],
            [1.4818755e+001, 2.6988314e+001],
            [1.3074880e+001, 2.5997641e+001],
            [1.2691689e+001, 2.3997682e+001],
            [0.0000000e+000, 0.0000000e+000],
            [0.0000000e+000, 2.0000000e+000],
            [1.6233887e+000, 3.2591152e+000],
            [1.7597298e+000, 4.9383626e+000],
            [3.4650326e+000, 5.9829826e+000],
            [3.4641016e+000, 8.0000000e+000],
            [5.1961524e+000, 9.0000000e+000],
            [5.1961524e+000, 1.1000000e+001],
            [6.9282032e+000, 1.2000000e+001],
            [6.9282032e+000, 1.4000000e+001],
            [6.9282032e+000, 1.6000000e+001],
            [8.6604283e+000, 1.7000000e+001],
            [8.6612997e+000, 1.9000000e+001],
            [1.0442277e+001, 1.9993373e+001],
            [1.1652419e+001, 2.3001274e+001],
            [1.0684821e+001, 2.1960240e+001],
            [2.0000000e+001, 3.2000000e+001],
            [1.8271341e+001, 3.2998399e+001],
            [1.6539855e+001, 3.1998132e+001],
            [1.4808574e+001, 3.0997769e+001],
            [1.3073571e+001, 2.9999162e+001],
            [1.3075386e+001, 2.7997531e+001],
            [1.1341020e+001, 2.6999097e+001],
            [1.0998951e+001, 2.5332403e+001],
            [1.0501118e+001, 2.3721681e+001],
            [0.0000000e+000, 4.0000000e+000],
            [0.0000000e+000, 6.0000000e+000],
            [1.7320508e+000, 7.0000000e+000],
            [1.7320508e+000, 9.0000000e+000],
            [3.4641016e+000, 1.0000000e+001],
            [3.4641016e+000, 1.2000000e+001],
            [5.1961524e+000, 1.3000000e+001],
            [5.1961524e+000, 1.5000000e+001],
            [5.1961524e+000, 1.7000000e+001],
            [6.9282032e+000, 1.8000000e+001],
            [6.9282032e+000, 2.0000000e+001],
            [8.6665280e+000, 2.1000000e+001],
            [8.6978979e+000, 2.3000000e+001],
            [2.0000000e+001, 3.4000000e+001],
            [1.8268514e+001, 3.4999733e+001],
            [1.6537217e+001, 3.3999377e+001],
            [1.4805515e+001, 3.2999213e+001],
            [1.3072862e+001, 3.1999497e+001],
            [1.1340219e+001, 3.0999777e+001],
            [1.1340719e+001, 2.8999412e+001],
            [9.7955554e+000, 2.7999728e+001],
            [9.8736113e+000, 2.6466246e+001],
            [8.8861171e+000, 2.5000000e+001],
            [0.0000000e+000, 8.0000000e+000],
            [0.0000000e+000, 1.0000000e+001],
            [1.7320508e+000, 1.1000000e+001],
            [1.7320508e+000, 1.3000000e+001],
            [3.4641016e+000, 1.4000000e+001],
            [3.4641016e+000, 1.6000000e+001],
            [3.4641016e+000, 1.8000000e+001],
            [5.1961524e+000, 1.9000000e+001],
            [5.1961524e+000, 2.1000000e+001],
            [6.9282032e+000, 2.2000000e+001],
            [6.9282033e+000, 2.4000000e+001],
            [2.0000000e+001, 3.6000000e+001],
            [1.8268043e+001, 3.6999956e+001],
            [1.6536228e+001, 3.5999844e+001],
            [1.4804400e+001, 3.4999739e+001],
            [1.3072344e+001, 3.3999742e+001],
            [1.1340094e+001, 3.2999836e+001],
            [9.7648297e+000, 3.1999950e+001],
            [9.7910616e+000, 2.9999859e+001],
            [8.3464240e+000, 2.9000002e+001],
            [8.3464128e+000, 2.7000000e+001],
            [6.9282036e+000, 2.6000000e+001],
            [0.0000000e+000, 1.2000000e+001],
            [0.0000000e+000, 1.4000000e+001],
            [1.7320508e+000, 1.5000000e+001],
            [1.7320508e+000, 1.7000000e+001],
            [1.7320508e+000, 1.9000000e+001],
            [3.4641016e+000, 2.0000000e+001],
            [3.4641016e+000, 2.2000000e+001],
            [5.1961524e+000, 2.3000000e+001],
            [5.1961525e+000, 2.5000000e+001],
            [2.0000000e+001, 3.8000000e+001],
            [1.8267965e+001, 3.8999993e+001],
            [1.6535972e+001, 3.7999965e+001],
            [1.4804007e+001, 3.6999925e+001],
            [1.3072007e+001, 3.5999901e+001],
            [1.1339930e+001, 3.4999913e+001],
            [9.7912980e+000, 3.4000044e+001],
            [8.3467920e+000, 3.3000074e+001],
            [8.3464827e+000, 3.1000012e+001],
            [6.9282121e+000, 3.0000000e+001],
            [6.9282050e+000, 2.8000000e+001],
            [5.1961526e+000, 2.7000000e+001],
            [0.0000000e+000, 1.6000000e+001],
            [0.0000000e+000, 1.8000000e+001],
            [0.0000000e+000, 2.0000000e+001],
            [1.7320508e+000, 2.1000000e+001],
            [1.7320508e+000, 2.3000000e+001],
            [3.4641016e+000, 2.4000000e+001],
            [3.4641016e+000, 2.6000000e+001],
            [2.0000000e+001, 4.0000000e+001],
            [1.8267952e+001, 4.0999999e+001],
            [1.6535914e+001, 3.9999993e+001],
            [1.4803889e+001, 3.8999981e+001],
            [1.3071865e+001, 3.7999968e+001],
            [1.1339823e+001, 3.6999964e+001],
            [9.7974474e+000, 3.6000625e+001],
            [8.3484470e+000, 3.5000492e+001],
            [6.9284133e+000, 3.3999958e+001],
            [6.9282475e+000, 3.1999996e+001],
            [5.1961604e+000, 3.1000002e+001],
            [5.1961538e+000, 2.9000000e+001],
            [3.4641016e+000, 2.8000000e+001],
            [0.0000000e+000, 2.2000000e+001],
            [0.0000000e+000, 2.4000000e+001],
            [1.7320508e+000, 2.5000000e+001],
            [1.7320508e+000, 2.7000000e+001],
            [2.0000000e+001, 4.2000000e+001],
            [1.8267950e+001, 4.3000000e+001],
            [1.6535901e+001, 4.1999999e+001],
            [1.4803858e+001, 4.0999995e+001],
            [1.3053051e+001, 4.0037850e+001],
            [1.1295233e+001, 3.9044813e+001],
            [9.8010076e+000, 3.8011963e+001],
            [8.3574909e+000, 3.7003339e+001],
            [6.9291337e+000, 3.5999654e+001],
            [5.1964352e+000, 3.5000083e+001],
            [5.1961998e+000, 3.3000014e+001],
            [3.4641018e+000, 3.2000000e+001],
            [3.4641017e+000, 3.0000000e+001],
            [1.7320508e+000, 2.9000000e+001],
            [0.0000000e+000, 2.6000000e+001],
            [0.0000000e+000, 2.8000000e+001],
            [2.0000000e+001, 4.4000000e+001],
            [1.8240270e+001, 4.5061637e+001],
            [1.6479493e+001, 4.4070580e+001],
            [1.4782483e+001, 4.3145737e+001],
            [1.3162423e+001, 4.2153730e+001],
            [1.1091269e+001, 4.1231096e+001],
            [9.8601024e+000, 3.9694249e+001],
            [8.4083559e+000, 3.9023038e+001],
            [6.9318114e+000, 3.7997341e+001],
            [5.1978440e+000, 3.7000499e+001],
            [3.4641060e+000, 3.5999998e+001],
            [3.4641026e+000, 3.3999999e+001],
            [1.7320510e+000, 3.3000000e+001],
            [1.7320508e+000, 3.1000000e+001],
            [0.0000000e+000, 3.0000000e+001],
            [2.0000000e+001, 4.6000000e+001],
            [1.8268349e+001, 4.7283235e+001],
            [1.6225143e+001, 4.6361846e+001],
            [1.4898084e+001, 4.4983944e+001],
            [1.3107994e+001, 4.4439846e+001],
            [1.1703124e+001, 4.3127139e+001],
            [9.9313829e+000, 4.3032346e+001],
            [8.7046457e+000, 4.1160337e+001],
            [6.9380336e+000, 3.9980550e+001],
            [5.2062796e+000, 3.9003000e+001],
            [3.4641193e+000, 3.7999994e+001],
            [1.7320582e+000, 3.6999994e+001],
            [1.7320520e+000, 3.4999999e+001],
            [0.0000000e+000, 3.4000000e+001],
            [0.0000000e+000, 3.2000000e+001],
            [2.0000000e+001, 4.8000000e+001],
            [1.8000000e+001, 5.0000000e+001],
            [1.6727635e+001, 4.8384202e+001],
            [1.4992130e+001, 4.8258316e+001],
            [1.4003301e+001, 4.6537709e+001],
            [1.2000295e+001, 4.6535888e+001],
            [1.1018049e+001, 4.4743170e+001],
            [9.0030083e+000, 4.4793723e+001],
            [7.8417404e+000, 4.3156953e+001],
            [6.8719748e+000, 4.1701962e+001],
            [5.2568443e+000, 4.1018001e+001],
            [3.4641554e+000, 4.0000007e+001],
            [1.7320954e+000, 3.8999964e+001],
            [0.0000000e+000, 3.8000000e+001],
            [0.0000000e+000, 3.6000000e+001],
            [2.0000000e+001, 5.0000000e+001],
            [1.6000000e+001, 5.0000000e+001],
            [1.4000000e+001, 5.0000000e+001],
            [1.2998688e+001, 4.8266344e+001],
            [1.0999781e+001, 4.8267682e+001],
            [1.0000007e+001, 4.6535845e+001],
            [7.9999940e+000, 4.6535881e+001],
            [6.8547620e+000, 4.4780789e+001],
            [5.5602395e+000, 4.3107704e+001],
            [3.4641121e+000, 4.2000295e+001],
            [1.7323184e+000, 4.0999781e+001],
            [0.0000000e+000, 4.0000000e+001],
            [1.2000000e+001, 5.0000000e+001],
            [1.0000000e+001, 5.0000000e+001],
            [8.9999636e+000, 4.8267905e+001],
            [6.9999939e+000, 4.8267942e+001],
            [5.9294172e+000, 4.6479488e+001],
            [5.0161554e+000, 4.4897744e+001],
            [3.4622914e+000, 4.4003301e+001],
            [1.7336564e+000, 4.2998688e+001],
            [0.0000000e+000, 4.2000000e+001],
            [8.0000000e+000, 5.0000000e+001],
            [6.0000000e+000, 5.0000000e+001],
            [4.9383616e+000, 4.8240269e+001],
            [3.6381539e+000, 4.6225143e+001],
            [1.7416841e+000, 4.4992130e+001],
            [0.0000000e+000, 4.4000000e+001],
            [4.0000000e+000, 5.0000000e+001],
            [2.7167647e+000, 4.8268349e+001],
            [1.6157983e+000, 4.6727635e+001],
            [0.0000000e+000, 4.6000000e+001],
            [2.0000000e+000, 5.0000000e+001],
            [0.0000000e+000, 4.8000000e+001],
            [0.0000000e+000, 5.0000000e+001]
        ])
    
        elements = np.array([
            [124, 137, 121],
            [217, 239, 218],
            [198, 217, 199],
            [181, 198, 182],
            [166, 181, 167],
            [152, 166, 153],
            [138, 152, 139],
            [137, 138, 122],
            [246, 223, 247],
            [2, 1, 3],
            [75, 92, 93],
            [123, 124, 109],
            [124, 121, 109],
            [239, 262, 240],
            [218, 239, 240],
            [407, 422, 408],
            [453, 455, 454],
            [223, 203, 224],
            [223, 224, 247],
            [203, 185, 204],
            [224, 203, 204],
            [185, 169, 186],
            [204, 185, 186],
            [169, 154, 170],
            [186, 169, 170],
            [154, 140, 155],
            [170, 154, 155],
            [140, 125, 141],
            [155, 140, 141],
            [125, 110, 126],
            [141, 125, 126],
            [110, 94, 111],
            [126, 110, 111],
            [94, 77, 95],
            [111, 94, 95],
            [77, 62, 78],
            [95, 77, 78],
            [62, 49, 63],
            [78, 62, 63],
            [49, 38, 50],
            [63, 49, 50],
            [38, 28, 39],
            [50, 38, 39],
            [28, 20, 29],
            [39, 28, 29],
            [20, 13, 21],
            [29, 20, 21],
            [13, 8, 14],
            [21, 13, 14],
            [8, 4, 9],
            [14, 8, 9],
            [4, 2, 5],
            [9, 4, 5],
            [5, 2, 3],
            [3, 6, 7],
            [3, 7, 5],
            [5, 7, 10],
            [6, 11, 12],
            [7, 6, 12],
            [11, 18, 19],
            [12, 11, 19],
            [18, 26, 27],
            [19, 18, 27],
            [26, 36, 37],
            [27, 26, 37],
            [36, 47, 48],
            [37, 36, 48],
            [47, 60, 61],
            [48, 47, 61],
            [60, 75, 76],
            [61, 60, 76],
            [76, 75, 93],
            [93, 106, 88],
            [93, 88, 76],
            [76, 88, 72],
            [106, 107, 89],
            [88, 106, 89],
            [107, 123, 108],
            [89, 107, 108],
            [123, 109, 108],
            [217, 218, 199],
            [198, 199, 182],
            [181, 182, 167],
            [166, 167, 153],
            [152, 153, 139],
            [138, 139, 122],
            [121, 137, 122],
            [262, 284, 263],
            [240, 262, 263],
            [284, 305, 285],
            [263, 284, 285],
            [305, 325, 306],
            [285, 305, 306],
            [325, 344, 326],
            [306, 325, 326],
            [344, 361, 345],
            [326, 344, 345],
            [361, 377, 362],
            [345, 361, 362],
            [377, 392, 378],
            [362, 377, 378],
            [392, 407, 393],
            [378, 392, 393],
            [393, 407, 408],
            [408, 423, 409],
            [408, 409, 393],
            [393, 409, 394],
            [423, 424, 410],
            [409, 423, 410],
            [424, 434, 425],
            [410, 424, 425],
            [434, 435, 426],
            [425, 434, 426],
            [435, 443, 436],
            [426, 435, 436],
            [443, 444, 437],
            [436, 443, 437],
            [444, 449, 445],
            [437, 444, 445],
            [449, 453, 450],
            [445, 449, 450],
            [450, 453, 454],
            [454, 452, 451],
            [454, 451, 450],
            [450, 451, 446],
            [452, 448, 447],
            [451, 452, 447],
            [448, 442, 441],
            [447, 448, 441],
            [442, 433, 432],
            [441, 442, 432],
            [433, 420, 419],
            [432, 433, 419],
            [420, 421, 403],
            [419, 420, 403],
            [421, 405, 404],
            [403, 421, 404],
            [405, 406, 389],
            [404, 405, 389],
            [406, 391, 390],
            [389, 406, 390],
            [391, 376, 374],
            [390, 391, 374],
            [376, 375, 360],
            [374, 376, 360],
            [375, 358, 359],
            [360, 375, 359],
            [358, 357, 341],
            [359, 358, 341],
            [357, 339, 340],
            [341, 357, 340],
            [339, 338, 320],
            [340, 339, 320],
            [338, 337, 319],
            [320, 338, 319],
            [337, 317, 318],
            [319, 337, 318],
            [317, 316, 297],
            [318, 317, 297],
            [316, 295, 296],
            [297, 316, 296],
            [295, 294, 274],
            [296, 295, 274],
            [294, 272, 273],
            [274, 294, 273],
            [272, 271, 249],
            [273, 272, 249],
            [271, 247, 248],
            [249, 271, 248],
            [248, 247, 224],
            [248, 224, 225],
            [108, 109, 91],
            [91, 109, 104],
            [104, 109, 121],
            [240, 263, 264],
            [240, 264, 241],
            [240, 241, 218],
            [224, 204, 225],
            [204, 186, 205],
            [204, 205, 225],
            [186, 170, 187],
            [205, 186, 187],
            [170, 155, 171],
            [187, 170, 171],
            [155, 141, 156],
            [171, 155, 156],
            [141, 126, 142],
            [156, 141, 142],
            [126, 111, 127],
            [142, 126, 127],
            [111, 95, 112],
            [127, 111, 112],
            [95, 78, 96],
            [112, 95, 96],
            [78, 63, 79],
            [96, 78, 79],
            [63, 50, 64],
            [79, 63, 64],
            [50, 39, 51],
            [64, 50, 51],
            [39, 29, 40],
            [51, 39, 40],
            [29, 21, 30],
            [40, 29, 30],
            [21, 14, 22],
            [30, 21, 22],
            [14, 9, 15],
            [22, 14, 15],
            [9, 5, 10],
            [15, 9, 10],
            [10, 7, 12],
            [12, 19, 16],
            [12, 16, 10],
            [19, 27, 24],
            [16, 19, 24],
            [27, 37, 33],
            [24, 27, 33],
            [37, 48, 44],
            [33, 37, 44],
            [48, 61, 57],
            [44, 48, 57],
            [61, 76, 72],
            [57, 61, 72],
            [72, 88, 89],
            [89, 108, 90],
            [89, 90, 72],
            [108, 91, 90],
            [263, 285, 286],
            [264, 263, 286],
            [285, 306, 307],
            [286, 285, 307],
            [306, 326, 327],
            [307, 306, 327],
            [326, 345, 346],
            [327, 326, 346],
            [345, 362, 363],
            [346, 345, 363],
            [362, 378, 379],
            [363, 362, 379],
            [378, 393, 394],
            [379, 378, 394],
            [394, 409, 410],
            [410, 425, 411],
            [410, 411, 394],
            [425, 426, 412],
            [411, 425, 412],
            [426, 436, 427],
            [412, 426, 427],
            [436, 437, 428],
            [427, 436, 428],
            [437, 445, 438],
            [428, 437, 438],
            [445, 450, 446],
            [438, 445, 446],
            [446, 451, 447],
            [447, 441, 440],
            [447, 440, 446],
            [441, 432, 431],
            [440, 441, 431],
            [432, 419, 418],
            [431, 432, 418],
            [419, 403, 402],
            [418, 419, 402],
            [403, 404, 387],
            [402, 403, 387],
            [404, 389, 388],
            [387, 404, 388],
            [389, 390, 372],
            [388, 389, 372],
            [390, 374, 373],
            [372, 390, 373],
            [374, 360, 356],
            [373, 374, 356],
            [360, 359, 343],
            [356, 360, 343],
            [359, 341, 342],
            [343, 359, 342],
            [341, 340, 322],
            [342, 341, 322],
            [340, 320, 321],
            [322, 340, 321],
            [320, 319, 300],
            [321, 320, 300],
            [319, 318, 299],
            [300, 319, 299],
            [318, 297, 298],
            [299, 318, 298],
            [297, 296, 276],
            [298, 297, 276],
            [296, 274, 275],
            [276, 296, 275],
            [274, 273, 251],
            [275, 274, 251],
            [273, 249, 250],
            [251, 273, 250],
            [249, 248, 225],
            [250, 249, 225],
            [91, 104, 86],
            [90, 91, 74],
            [91, 86, 74],
            [86, 104, 87],
            [87, 104, 105],
            [105, 104, 121],
            [241, 264, 265],
            [264, 286, 287],
            [265, 264, 287],
            [241, 265, 242],
            [241, 242, 219],
            [241, 219, 218],
            [205, 187, 206],
            [205, 206, 226],
            [205, 226, 225],
            [187, 171, 188],
            [206, 187, 188],
            [171, 156, 172],
            [188, 171, 172],
            [156, 142, 157],
            [172, 156, 157],
            [142, 127, 143],
            [157, 142, 143],
            [127, 112, 128],
            [143, 127, 128],
            [112, 96, 113],
            [128, 112, 113],
            [96, 79, 97],
            [113, 96, 97],
            [79, 64, 80],
            [97, 79, 80],
            [64, 51, 65],
            [80, 64, 65],
            [51, 40, 52],
            [65, 51, 52],
            [40, 30, 41],
            [52, 40, 41],
            [30, 22, 31],
            [41, 30, 31],
            [22, 15, 23],
            [31, 22, 23],
            [23, 15, 17],
            [17, 15, 10],
            [16, 24, 25],
            [16, 25, 17],
            [16, 17, 10],
            [17, 25, 23],
            [24, 33, 34],
            [25, 24, 34],
            [33, 44, 45],
            [34, 33, 45],
            [44, 57, 58],
            [45, 44, 58],
            [58, 57, 73],
            [73, 57, 72],
            [90, 74, 73],
            [90, 73, 72],
            [73, 74, 58],
            [286, 307, 308],
            [287, 286, 308],
            [307, 327, 328],
            [308, 307, 328],
            [327, 346, 347],
            [328, 327, 347],
            [346, 363, 364],
            [347, 346, 364],
            [363, 379, 380],
            [364, 363, 380],
            [380, 379, 395],
            [395, 379, 394],
            [411, 412, 396],
            [411, 396, 395],
            [411, 395, 394],
            [395, 396, 380],
            [412, 427, 413],
            [396, 412, 413],
            [427, 428, 414],
            [413, 427, 414],
            [428, 438, 429],
            [414, 428, 429],
            [429, 438, 439],
            [439, 438, 446],
            [440, 431, 430],
            [440, 430, 439],
            [440, 439, 446],
            [439, 430, 429],
            [431, 418, 417],
            [430, 431, 417],
            [418, 402, 401],
            [417, 418, 401],
            [402, 387, 386],
            [401, 402, 386],
            [387, 388, 370],
            [386, 387, 370],
            [388, 372, 371],
            [370, 388, 371],
            [372, 373, 354],
            [371, 372, 354],
            [373, 356, 355],
            [354, 373, 355],
            [356, 343, 336],
            [355, 356, 336],
            [343, 342, 324],
            [336, 343, 324],
            [342, 322, 323],
            [324, 342, 323],
            [322, 321, 302],
            [323, 322, 302],
            [321, 300, 301],
            [302, 321, 301],
            [300, 299, 279],
            [301, 300, 279],
            [299, 298, 278],
            [279, 299, 278],
            [298, 276, 277],
            [278, 298, 277],
            [276, 275, 253],
            [277, 276, 253],
            [275, 251, 252],
            [253, 275, 252],
            [251, 250, 227],
            [252, 251, 227],
            [227, 250, 226],
            [226, 250, 225],
            [227, 226, 206],
            [74, 86, 70],
            [86, 87, 71],
            [70, 86, 71],
            [74, 70, 58],
            [87, 105, 103],
            [71, 87, 85],
            [87, 103, 85],
            [103, 105, 120],
            [120, 105, 122],
            [265, 287, 288],
            [242, 265, 266],
            [265, 288, 266],
            [287, 308, 309],
            [288, 287, 309],
            [219, 242, 243],
            [242, 266, 267],
            [243, 242, 267],
            [219, 243, 220],
            [219, 220, 199],
            [206, 188, 207],
            [206, 207, 227],
            [188, 172, 189],
            [207, 188, 189],
            [172, 157, 173],
            [189, 172, 173],
            [157, 143, 158],
            [173, 157, 158],
            [143, 128, 144],
            [158, 143, 144],
            [128, 113, 129],
            [144, 128, 129],
            [113, 97, 114],
            [129, 113, 114],
            [97, 80, 98],
            [114, 97, 98],
            [80, 65, 81],
            [98, 80, 81],
            [65, 52, 66],
            [81, 65, 66],
            [52, 41, 53],
            [66, 52, 53],
            [41, 31, 42],
            [53, 41, 42],
            [31, 23, 32],
            [42, 31, 32],
            [32, 23, 25],
            [25, 34, 35],
            [25, 35, 32],
            [32, 35, 43],
            [34, 45, 46],
            [35, 34, 46],
            [45, 58, 59],
            [46, 45, 59],
            [59, 58, 70],
            [59, 70, 56],
            [308, 328, 329],
            [309, 308, 329],
            [328, 347, 348],
            [329, 328, 348],
            [347, 364, 365],
            [348, 347, 365],
            [364, 380, 381],
            [365, 364, 381],
            [381, 380, 396],
            [396, 413, 397],
            [396, 397, 381],
            [381, 397, 382],
            [413, 414, 398],
            [397, 413, 398],
            [414, 429, 415],
            [398, 414, 415],
            [415, 429, 430],
            [430, 417, 416],
            [430, 416, 415],
            [415, 416, 399],
            [417, 401, 400],
            [416, 417, 400],
            [401, 386, 385],
            [400, 401, 385],
            [386, 370, 369],
            [385, 386, 369],
            [370, 371, 352],
            [369, 370, 352],
            [371, 354, 353],
            [352, 371, 353],
            [354, 355, 334],
            [353, 354, 334],
            [355, 336, 335],
            [334, 355, 335],
            [336, 324, 315],
            [335, 336, 315],
            [324, 323, 304],
            [315, 324, 304],
            [323, 302, 303],
            [304, 323, 303],
            [302, 301, 281],
            [303, 302, 281],
            [301, 279, 280],
            [281, 301, 280],
            [279, 278, 256],
            [280, 279, 256],
            [278, 277, 255],
            [256, 278, 255],
            [277, 253, 254],
            [255, 277, 254],
            [253, 252, 229],
            [254, 253, 229],
            [252, 227, 228],
            [229, 252, 228],
            [228, 227, 207],
            [228, 207, 208],
            [70, 71, 56],
            [71, 85, 69],
            [71, 69, 56],
            [85, 103, 102],
            [103, 120, 119],
            [102, 103, 119],
            [69, 85, 84],
            [85, 102, 84],
            [119, 120, 135],
            [135, 120, 136],
            [136, 120, 122],
            [136, 122, 139],
            [136, 139, 151],
            [151, 139, 153],
            [266, 288, 289],
            [288, 309, 310],
            [289, 288, 310],
            [267, 266, 290],
            [266, 289, 290],
            [309, 329, 330],
            [310, 309, 330],
            [243, 267, 244],
            [220, 243, 221],
            [243, 244, 221],
            [267, 290, 268],
            [244, 267, 268],
            [220, 221, 201],
            [220, 201, 200],
            [220, 200, 199],
            [199, 200, 182],
            [182, 200, 183],
            [182, 183, 167],
            [207, 189, 208],
            [189, 173, 190],
            [189, 190, 208],
            [173, 158, 174],
            [190, 173, 174],
            [158, 144, 159],
            [174, 158, 159],
            [144, 129, 145],
            [159, 144, 145],
            [129, 114, 130],
            [145, 129, 130],
            [114, 98, 115],
            [130, 114, 115],
            [98, 81, 99],
            [115, 98, 99],
            [81, 66, 82],
            [99, 81, 82],
            [66, 53, 67],
            [82, 66, 67],
            [53, 42, 54],
            [67, 53, 54],
            [42, 32, 43],
            [54, 42, 43],
            [43, 35, 46],
            [46, 59, 56],
            [329, 348, 349],
            [330, 329, 349],
            [348, 365, 366],
            [349, 348, 366],
            [365, 381, 382],
            [366, 365, 382],
            [382, 397, 398],
            [398, 415, 399],
            [399, 416, 400],
            [400, 385, 384],
            [400, 384, 399],
            [385, 369, 368],
            [384, 385, 368],
            [369, 352, 351],
            [368, 369, 351],
            [352, 353, 332],
            [351, 352, 332],
            [353, 334, 333],
            [332, 353, 333],
            [334, 335, 313],
            [333, 334, 313],
            [335, 315, 314],
            [313, 335, 314],
            [315, 304, 293],
            [314, 315, 293],
            [304, 303, 283],
            [293, 304, 283],
            [303, 281, 282],
            [283, 303, 282],
            [281, 280, 258],
            [282, 281, 258],
            [280, 256, 257],
            [258, 280, 257],
            [256, 255, 232],
            [257, 256, 232],
            [255, 254, 231],
            [232, 255, 231],
            [254, 229, 230],
            [231, 254, 230],
            [229, 228, 208],
            [230, 229, 208],
            [69, 84, 68],
            [69, 68, 55],
            [69, 55, 56],
            [102, 119, 118],
            [84, 102, 101],
            [102, 118, 101],
            [119, 135, 134],
            [118, 119, 134],
            [68, 84, 83],
            [84, 101, 83],
            [135, 136, 151],
            [134, 135, 150],
            [150, 135, 151],
            [289, 310, 311],
            [290, 289, 312],
            [289, 311, 312],
            [310, 330, 331],
            [311, 310, 331],
            [268, 290, 291],
            [290, 312, 291],
            [330, 349, 350],
            [331, 330, 350],
            [221, 244, 245],
            [244, 268, 269],
            [245, 244, 269],
            [201, 221, 222],
            [221, 245, 222],
            [268, 291, 292],
            [269, 268, 292],
            [200, 201, 183],
            [201, 222, 202],
            [201, 202, 183],
            [190, 174, 191],
            [190, 191, 209],
            [190, 209, 208],
            [174, 159, 175],
            [191, 174, 175],
            [159, 145, 160],
            [175, 159, 160],
            [145, 130, 146],
            [160, 145, 146],
            [130, 115, 131],
            [146, 130, 131],
            [115, 99, 116],
            [131, 115, 116],
            [99, 82, 100],
            [116, 99, 100],
            [82, 67, 83],
            [82, 83, 100],
            [83, 67, 68],
            [68, 67, 54],
            [68, 54, 55],
            [55, 54, 43],
            [55, 43, 56],
            [349, 366, 367],
            [350, 349, 367],
            [367, 366, 383],
            [383, 366, 382],
            [384, 368, 367],
            [384, 367, 383],
            [384, 383, 399],
            [367, 368, 350],
            [350, 368, 351],
            [350, 351, 331],
            [331, 351, 332],
            [331, 332, 311],
            [311, 332, 333],
            [311, 333, 312],
            [312, 333, 313],
            [312, 313, 291],
            [291, 313, 314],
            [291, 314, 292],
            [292, 314, 293],
            [292, 293, 269],
            [399, 383, 382],
            [293, 283, 270],
            [293, 270, 269],
            [269, 270, 245],
            [283, 282, 261],
            [270, 283, 261],
            [282, 258, 259],
            [261, 282, 259],
            [258, 257, 234],
            [259, 258, 234],
            [257, 232, 233],
            [234, 257, 233],
            [232, 231, 211],
            [233, 232, 211],
            [231, 230, 210],
            [211, 231, 210],
            [210, 230, 209],
            [209, 230, 208],
            [210, 209, 191],
            [101, 118, 117],
            [118, 134, 132],
            [117, 118, 132],
            [83, 101, 100],
            [100, 101, 117],
            [100, 117, 116],
            [116, 117, 132],
            [134, 150, 148],
            [132, 134, 148],
            [148, 150, 164],
            [164, 150, 165],
            [165, 150, 151],
            [222, 245, 237],
            [237, 245, 260],
            [260, 245, 270],
            [260, 270, 261],
            [260, 261, 237],
            [202, 222, 216],
            [222, 237, 216],
            [202, 216, 197],
            [202, 197, 184],
            [202, 184, 183],
            [191, 175, 192],
            [191, 192, 210],
            [175, 160, 176],
            [192, 175, 176],
            [160, 146, 161],
            [176, 160, 161],
            [146, 131, 147],
            [161, 146, 147],
            [131, 116, 133],
            [147, 131, 133],
            [133, 116, 132],
            [133, 132, 148],
            [261, 259, 238],
            [261, 238, 237],
            [237, 238, 216],
            [259, 234, 235],
            [238, 259, 235],
            [234, 233, 213],
            [235, 234, 213],
            [233, 211, 212],
            [213, 233, 212],
            [211, 210, 192],
            [212, 211, 192],
            [148, 164, 162],
            [148, 162, 149],
            [148, 149, 133],
            [133, 149, 147],
            [147, 149, 162],
            [164, 165, 180],
            [162, 164, 178],
            [164, 180, 178],
            [180, 165, 168],
            [168, 165, 153],
            [197, 216, 215],
            [215, 216, 236],
            [236, 216, 238],
            [236, 238, 235],
            [236, 235, 215],
            [184, 197, 180],
            [197, 215, 195],
            [197, 195, 180],
            [192, 176, 193],
            [192, 193, 212],
            [176, 161, 177],
            [193, 176, 177],
            [161, 147, 163],
            [177, 161, 163],
            [163, 147, 162],
            [163, 162, 178],
            [235, 213, 214],
            [235, 214, 215],
            [215, 214, 195],
            [213, 212, 194],
            [214, 213, 194],
            [194, 212, 193],
            [178, 180, 195],
            [180, 168, 184],
            [184, 168, 167],
            [193, 177, 194],
            [177, 163, 179],
            [177, 179, 194],
            [179, 163, 178],
            [179, 178, 195],
            [214, 194, 196],
            [214, 196, 195],
            [195, 196, 179],
            [179, 196, 194],
            [165, 151, 153],
            [183, 184, 167],
            [105, 121, 122],
            [218, 219, 199],
            [43, 46, 56],
            [382, 398, 399],
            [153, 167, 168]
        ])
    
        # There are 2 hard problems in computer science:
        # cache invalidation, naming things, and off-by-1 errors.
        self.elements = elements - 1
        
        # Dirichlet BC
        # [node, direction, value]
        neumann = np.array([
            [408, 1, 100.000000],
            [422, 1, 50.000000],
            [423, 1, 100.000000],
            [424, 1, 100.000000],
            [434, 1, 100.000000],
            [435, 1, 100.000000],
            [443, 1, 100.000000],
            [444, 1, 100.000000],
            [449, 1, 100.000000],
            [453, 1, 100.000000],
            [455, 1, 50.000000]
        ])
        # There are 2 hard problems in computer science:
        # cache invalidation, naming things, and off-by-1 errors.
        neumann[:, 0] = neumann[:, 0] - 1
        self.neumann = neumann
    
        # Dirichlet BC
        # [node, direction, value]
        dirichlet = np.array([
            [246, 1, 0],
            [223, 1, 0],
            [203, 1, 0],
            [185, 1, 0],
            [169, 1, 0],
            [154, 1, 0],
            [140, 1, 0],
            [125, 1, 0],
            [110, 1, 0],
            [94, 1, 0],
            [77, 1, 0],
            [62, 1, 0],
            [49, 1, 0],
            [38, 1, 0],
            [28, 1, 0],
            [20, 1, 0],
            [13, 1, 0],
            [8, 1, 0],
            [4, 1, 0],
            [2, 1, 0],
            [1, 1, 0],
            [246, 0, 0],
            [223, 0, 0],
            [203, 0, 0],
            [185, 0, 0],
            [169, 0, 0],
            [154, 0, 0],
            [140, 0, 0],
            [125, 0, 0],
            [110, 0, 0],
            [94, 0, 0],
            [77, 0, 0],
            [62, 0, 0],
            [49, 0, 0],
            [38, 0, 0],
            [28, 0, 0],
            [20, 0, 0],
            [13, 0, 0],
            [8, 0, 0],
            [4, 0, 0],
            [2, 0, 0],
            [1, 0, 0],
        ])
    
        # There are 2 hard problems in computer science:
        # cache invalidation, naming things, and off-by-1 errors.
        dirichlet[:, 0] = dirichlet[:, 0] - 1
        self.dirichlet = dirichlet
    
        # internal forces
        self.force = np.array([0.0, 0.0])

        # material properties
        self.young_modulus = 2.1e5
        self.poisson_ratio = 0.3



class TestSetup:
    def __init__(self):
        nodes = [[0.0, 0.0],
                 [0.0, 1.0],
                 [1.0, 0.0],
                 [1.0, 1.0],
                 [2.0, 0.0],
                 [2.0, 1.0]]

        self.nodes = np.array(nodes)

        elements = [[0, 3, 1],
                    [0, 2, 3],
                    [2, 5, 3],
                    [2, 4, 5]]
        self.elements = np.array(elements)

        self.material_properties = np.array([2000, 0.25])

        # Dirichlet BC
        # [node, direction, value]
        dirichlet = [[0, 0, 0],
                     [0, 1, 0],
                     [1, 0, 0]]
        self.dirichlet = np.array(dirichlet)

        # Neumann BC
        # [node, direction, value]
        neumann = [[4, 0, 1],
                   [5, 0, 1]]
        self.neumann = np.array(neumann)

        # internal forces
        self.force = np.array([0.0, -0.2])

        # material properties
        self.young_modulus = 2000
        self.poisson_ratio = 0.25