Subversion Repositories HelenOS

Rev

Rev 283 | Rev 501 | 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 <print.h>
31
#include <debug.h>
32
#include <panic.h>
33
 
34
#include <test.h>
111 palkovsky 35
#include <arch/atomic.h>
84 vana 36
#include <proc/thread.h>
37
 
94 jermar 38
#include <arch.h>
39
 
91 jermar 40
#define THREADS     150*2
283 palkovsky 41
#define ATTEMPTS    100
91 jermar 42
 
43
#define E_10e8  271828182
44
#define PI_10e8 314159265
45
 
46
static inline double sqrt(double x) { double v; __asm__ ("fsqrt\n" : "=t" (v) : "0" (x)); return v; }
47
 
483 jermar 48
static atomic_t threads_ok;
91 jermar 49
static waitq_t can_start;
50
 
86 jermar 51
static void e(void *data)
84 vana 52
{
92 jermar 53
    int i;
91 jermar 54
    double e,d,le,f;
55
 
56
    waitq_sleep(&can_start);
57
 
92 jermar 58
    for (i = 0; i<ATTEMPTS; i++) {
59
        le=-1;
60
        e=0;
61
        f=1;
62
 
63
        for(d=1;e!=le;d*=f,f+=1) {
64
            le=e;
65
            e=e+1/d;
66
        }
67
 
68
        if((int)(100000000*e)!=E_10e8)
69
            panic("tid%d: e*10e8=%d\n", THREAD->tid, (int) 100000000*e);
86 jermar 70
    }
91 jermar 71
 
119 jermar 72
    atomic_inc(&threads_ok);
84 vana 73
}
74
 
91 jermar 75
static void pi(void *data)
76
{
92 jermar 77
    int i;
125 jermar 78
    double lpi, pi;
79
    double n, ab, ad;
84 vana 80
 
91 jermar 81
    waitq_sleep(&can_start);
84 vana 82
 
91 jermar 83
 
92 jermar 84
    for (i = 0; i<ATTEMPTS; i++) {
85
        lpi = -1;
86
        pi = 0;
91 jermar 87
 
125 jermar 88
        for (n=2, ab = sqrt(2); lpi != pi; n *= 2, ab = ad) {
89
            double sc, cd;
92 jermar 90
 
125 jermar 91
            sc = sqrt(1 - (ab*ab/4));
92
            cd = 1 - sc;
93
            ad = sqrt(ab*ab/4 + cd*cd);
94
            lpi = pi;
95
            pi = 2 * n * ad;
96
        }
92 jermar 97
 
98
        if((int)(100000000*pi)!=PI_10e8)
99
            panic("tid%d: pi*10e8=%d\n", THREAD->tid, (int) 100000000*pi);
91 jermar 100
    }
92 jermar 101
 
119 jermar 102
    atomic_inc(&threads_ok);
91 jermar 103
}
104
 
105
 
84 vana 106
void test(void)
107
{
86 jermar 108
    thread_t *t;
109
    int i;
84 vana 110
 
91 jermar 111
    waitq_initialize(&can_start);
112
 
113
    printf("FPU test #1\n");
114
    printf("Creating %d threads... ", THREADS);
115
 
116
    for (i=0; i<THREADS/2; i++) {  
117
        if (!(t = thread_create(e, NULL, TASK, 0)))
118
            panic("could not create thread\n");
86 jermar 119
        thread_ready(t);
91 jermar 120
        if (!(t = thread_create(pi, NULL, TASK, 0)))
121
            panic("could not create thread\n");
122
        thread_ready(t);
86 jermar 123
    }
91 jermar 124
    printf("ok\n");
90 vana 125
 
91 jermar 126
    thread_sleep(1);
127
    waitq_wakeup(&can_start, WAKEUP_ALL);
84 vana 128
 
91 jermar 129
    while (threads_ok != THREADS)
130
        ;
131
 
132
    printf("Test passed.\n");
84 vana 133
}