Subversion Repositories HelenOS

Rev

Rev 91 | Rev 94 | Go to most recent revision | Details | Compare with Previous | Last modification | View Log | RSS feed

Rev Author Line No. Line
86 jermar 1
/*
2
 * Copyright (C) 2005 Jakub Vana
91 jermar 3
 * Copyright (C) 2005 Jakub Jermar
86 jermar 4
 * All rights reserved.
5
 *
6
 * Redistribution and use in source and binary forms, with or without
7
 * modification, are permitted provided that the following conditions
8
 * are met:
9
 *
10
 * - Redistributions of source code must retain the above copyright
11
 *   notice, this list of conditions and the following disclaimer.
12
 * - Redistributions in binary form must reproduce the above copyright
13
 *   notice, this list of conditions and the following disclaimer in the
14
 *   documentation and/or other materials provided with the distribution.
15
 * - The name of the author may not be used to endorse or promote products
16
 *   derived from this software without specific prior written permission.
17
 *
18
 * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
19
 * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
20
 * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
21
 * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
22
 * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
23
 * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
24
 * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
25
 * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
26
 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
27
 * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
28
 */
84 vana 29
 
30
#include <arch/interrupt.h>
31
#include <print.h>
32
#include <debug.h>
33
#include <panic.h>
34
#include <arch/i8259.h>
35
#include <func.h>
36
#include <cpu.h>
37
#include <arch/asm.h>
38
#include <mm/tlb.h>
39
 
40
#include <test.h>
41
#include <arch.h>
42
#include <arch/smp/atomic.h>
43
#include <proc/thread.h>
44
 
91 jermar 45
#define THREADS     150*2
92 jermar 46
#define ATTEMPTS    10
91 jermar 47
 
48
#define E_10e8  271828182
49
#define PI_10e8 314159265
50
 
51
static inline double sqrt(double x) { double v; __asm__ ("fsqrt\n" : "=t" (v) : "0" (x)); return v; }
52
 
53
static volatile int threads_ok;
54
static waitq_t can_start;
55
 
86 jermar 56
static void e(void *data)
84 vana 57
{
92 jermar 58
    int i;
91 jermar 59
    double e,d,le,f;
60
 
61
    waitq_sleep(&can_start);
62
 
92 jermar 63
    for (i = 0; i<ATTEMPTS; i++) {
64
        le=-1;
65
        e=0;
66
        f=1;
67
 
68
        for(d=1;e!=le;d*=f,f+=1) {
69
            le=e;
70
            e=e+1/d;
71
        }
72
 
73
        if((int)(100000000*e)!=E_10e8)
74
            panic("tid%d: e*10e8=%d\n", THREAD->tid, (int) 100000000*e);
86 jermar 75
    }
91 jermar 76
 
92 jermar 77
    atomic_inc((int *) &threads_ok);
84 vana 78
}
79
 
91 jermar 80
static void pi(void *data)
81
{
92 jermar 82
    int i;
83
        double lpi, pi;
84
        double n, ab, ad;
84 vana 85
 
91 jermar 86
    waitq_sleep(&can_start);
84 vana 87
 
91 jermar 88
 
92 jermar 89
    for (i = 0; i<ATTEMPTS; i++) {
90
        lpi = -1;
91
        pi = 0;
91 jermar 92
 
92 jermar 93
            for (n=2, ab = sqrt(2); lpi != pi; n *= 2, ab = ad) {
94
                    double sc, cd;
95
 
96
                    sc = sqrt(1 - (ab*ab/4));
97
                    cd = 1 - sc;
98
                    ad = sqrt(ab*ab/4 + cd*cd);
99
                    lpi = pi;
100
                    pi = 2 * n * ad;
101
            }
102
 
103
        if((int)(100000000*pi)!=PI_10e8)
104
            panic("tid%d: pi*10e8=%d\n", THREAD->tid, (int) 100000000*pi);
91 jermar 105
    }
92 jermar 106
 
107
    atomic_inc((int *) &threads_ok);
91 jermar 108
}
109
 
110
 
84 vana 111
void test(void)
112
{
86 jermar 113
    thread_t *t;
114
    int i;
84 vana 115
 
91 jermar 116
    waitq_initialize(&can_start);
117
 
118
    printf("FPU test #1\n");
119
    printf("Creating %d threads... ", THREADS);
120
 
121
    for (i=0; i<THREADS/2; i++) {  
122
        if (!(t = thread_create(e, NULL, TASK, 0)))
123
            panic("could not create thread\n");
86 jermar 124
        thread_ready(t);
91 jermar 125
        if (!(t = thread_create(pi, NULL, TASK, 0)))
126
            panic("could not create thread\n");
127
        thread_ready(t);
86 jermar 128
    }
91 jermar 129
    printf("ok\n");
90 vana 130
 
91 jermar 131
    thread_sleep(1);
132
    waitq_wakeup(&can_start, WAKEUP_ALL);
84 vana 133
 
91 jermar 134
    while (threads_ok != THREADS)
135
        ;
136
 
137
    printf("Test passed.\n");
84 vana 138
}