Statistics
| Branch: | Revision:

root / target-s390x / fpu_helper.c @ 5b50e790

History | View | Annotate | Download (20.8 kB)

1
/*
2
 *  S/390 FPU helper routines
3
 *
4
 *  Copyright (c) 2009 Ulrich Hecht
5
 *  Copyright (c) 2009 Alexander Graf
6
 *
7
 * This library is free software; you can redistribute it and/or
8
 * modify it under the terms of the GNU Lesser General Public
9
 * License as published by the Free Software Foundation; either
10
 * version 2 of the License, or (at your option) any later version.
11
 *
12
 * This library is distributed in the hope that it will be useful,
13
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
14
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
15
 * Lesser General Public License for more details.
16
 *
17
 * You should have received a copy of the GNU Lesser General Public
18
 * License along with this library; if not, see <http://www.gnu.org/licenses/>.
19
 */
20

    
21
#include "cpu.h"
22
#include "helper.h"
23

    
24
#if !defined(CONFIG_USER_ONLY)
25
#include "exec/softmmu_exec.h"
26
#endif
27

    
28
/* #define DEBUG_HELPER */
29
#ifdef DEBUG_HELPER
30
#define HELPER_LOG(x...) qemu_log(x)
31
#else
32
#define HELPER_LOG(x...)
33
#endif
34

    
35
#define RET128(F) (env->retxl = F.low, F.high)
36

    
37
#define convert_bit(mask, from, to) \
38
    (to < from                      \
39
     ? (mask / (from / to)) & to    \
40
     : (mask & from) * (to / from))
41

    
42
static void ieee_exception(CPUS390XState *env, uint32_t dxc, uintptr_t retaddr)
43
{
44
    /* Install the DXC code.  */
45
    env->fpc = (env->fpc & ~0xff00) | (dxc << 8);
46
    /* Trap.  */
47
    runtime_exception(env, PGM_DATA, retaddr);
48
}
49

    
50
/* Should be called after any operation that may raise IEEE exceptions.  */
51
static void handle_exceptions(CPUS390XState *env, uintptr_t retaddr)
52
{
53
    unsigned s390_exc, qemu_exc;
54

    
55
    /* Get the exceptions raised by the current operation.  Reset the
56
       fpu_status contents so that the next operation has a clean slate.  */
57
    qemu_exc = env->fpu_status.float_exception_flags;
58
    if (qemu_exc == 0) {
59
        return;
60
    }
61
    env->fpu_status.float_exception_flags = 0;
62

    
63
    /* Convert softfloat exception bits to s390 exception bits.  */
64
    s390_exc = 0;
65
    s390_exc |= convert_bit(qemu_exc, float_flag_invalid, 0x80);
66
    s390_exc |= convert_bit(qemu_exc, float_flag_divbyzero, 0x40);
67
    s390_exc |= convert_bit(qemu_exc, float_flag_overflow, 0x20);
68
    s390_exc |= convert_bit(qemu_exc, float_flag_underflow, 0x10);
69
    s390_exc |= convert_bit(qemu_exc, float_flag_inexact, 0x08);
70

    
71
    /* Install the exceptions that we raised.  */
72
    env->fpc |= s390_exc << 16;
73

    
74
    /* Send signals for enabled exceptions.  */
75
    s390_exc &= env->fpc >> 24;
76
    if (s390_exc) {
77
        ieee_exception(env, s390_exc, retaddr);
78
    }
79
}
80

    
81
static inline int float_comp_to_cc(CPUS390XState *env, int float_compare)
82
{
83
    switch (float_compare) {
84
    case float_relation_equal:
85
        return 0;
86
    case float_relation_less:
87
        return 1;
88
    case float_relation_greater:
89
        return 2;
90
    case float_relation_unordered:
91
        return 3;
92
    default:
93
        cpu_abort(env, "unknown return value for float compare\n");
94
    }
95
}
96

    
97
/* condition codes for unary FP ops */
98
uint32_t set_cc_nz_f32(float32 v)
99
{
100
    if (float32_is_any_nan(v)) {
101
        return 3;
102
    } else if (float32_is_zero(v)) {
103
        return 0;
104
    } else if (float32_is_neg(v)) {
105
        return 1;
106
    } else {
107
        return 2;
108
    }
109
}
110

    
111
uint32_t set_cc_nz_f64(float64 v)
112
{
113
    if (float64_is_any_nan(v)) {
114
        return 3;
115
    } else if (float64_is_zero(v)) {
116
        return 0;
117
    } else if (float64_is_neg(v)) {
118
        return 1;
119
    } else {
120
        return 2;
121
    }
122
}
123

    
124
uint32_t set_cc_nz_f128(float128 v)
125
{
126
    if (float128_is_any_nan(v)) {
127
        return 3;
128
    } else if (float128_is_zero(v)) {
129
        return 0;
130
    } else if (float128_is_neg(v)) {
131
        return 1;
132
    } else {
133
        return 2;
134
    }
135
}
136

    
137
/* 32-bit FP addition */
138
uint64_t HELPER(aeb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
139
{
140
    float32 ret = float32_add(f1, f2, &env->fpu_status);
141
    handle_exceptions(env, GETPC());
142
    return ret;
143
}
144

    
145
/* 64-bit FP addition */
146
uint64_t HELPER(adb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
147
{
148
    float64 ret = float64_add(f1, f2, &env->fpu_status);
149
    handle_exceptions(env, GETPC());
150
    return ret;
151
}
152

    
153
/* 128-bit FP addition */
154
uint64_t HELPER(axb)(CPUS390XState *env, uint64_t ah, uint64_t al,
155
                     uint64_t bh, uint64_t bl)
156
{
157
    float128 ret = float128_add(make_float128(ah, al),
158
                                make_float128(bh, bl),
159
                                &env->fpu_status);
160
    handle_exceptions(env, GETPC());
161
    return RET128(ret);
162
}
163

    
164
/* 32-bit FP subtraction */
165
uint64_t HELPER(seb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
166
{
167
    float32 ret = float32_sub(f1, f2, &env->fpu_status);
168
    handle_exceptions(env, GETPC());
169
    return ret;
170
}
171

    
172
/* 64-bit FP subtraction */
173
uint64_t HELPER(sdb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
174
{
175
    float64 ret = float64_sub(f1, f2, &env->fpu_status);
176
    handle_exceptions(env, GETPC());
177
    return ret;
178
}
179

    
180
/* 128-bit FP subtraction */
181
uint64_t HELPER(sxb)(CPUS390XState *env, uint64_t ah, uint64_t al,
182
                     uint64_t bh, uint64_t bl)
183
{
184
    float128 ret = float128_sub(make_float128(ah, al),
185
                                make_float128(bh, bl),
186
                                &env->fpu_status);
187
    handle_exceptions(env, GETPC());
188
    return RET128(ret);
189
}
190

    
191
/* 32-bit FP division */
192
uint64_t HELPER(deb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
193
{
194
    float32 ret = float32_div(f1, f2, &env->fpu_status);
195
    handle_exceptions(env, GETPC());
196
    return ret;
197
}
198

    
199
/* 64-bit FP division */
200
uint64_t HELPER(ddb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
201
{
202
    float64 ret = float64_div(f1, f2, &env->fpu_status);
203
    handle_exceptions(env, GETPC());
204
    return ret;
205
}
206

    
207
/* 128-bit FP division */
208
uint64_t HELPER(dxb)(CPUS390XState *env, uint64_t ah, uint64_t al,
209
                     uint64_t bh, uint64_t bl)
210
{
211
    float128 ret = float128_div(make_float128(ah, al),
212
                                make_float128(bh, bl),
213
                                &env->fpu_status);
214
    handle_exceptions(env, GETPC());
215
    return RET128(ret);
216
}
217

    
218
/* 32-bit FP multiplication */
219
uint64_t HELPER(meeb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
220
{
221
    float32 ret = float32_mul(f1, f2, &env->fpu_status);
222
    handle_exceptions(env, GETPC());
223
    return ret;
224
}
225

    
226
/* 64-bit FP multiplication */
227
uint64_t HELPER(mdb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
228
{
229
    float64 ret = float64_mul(f1, f2, &env->fpu_status);
230
    handle_exceptions(env, GETPC());
231
    return ret;
232
}
233

    
234
/* 64/32-bit FP multiplication */
235
uint64_t HELPER(mdeb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
236
{
237
    float64 ret = float32_to_float64(f2, &env->fpu_status);
238
    ret = float64_mul(f1, ret, &env->fpu_status);
239
    handle_exceptions(env, GETPC());
240
    return ret;
241
}
242

    
243
/* 128-bit FP multiplication */
244
uint64_t HELPER(mxb)(CPUS390XState *env, uint64_t ah, uint64_t al,
245
                     uint64_t bh, uint64_t bl)
246
{
247
    float128 ret = float128_mul(make_float128(ah, al),
248
                                make_float128(bh, bl),
249
                                &env->fpu_status);
250
    handle_exceptions(env, GETPC());
251
    return RET128(ret);
252
}
253

    
254
/* 128/64-bit FP multiplication */
255
uint64_t HELPER(mxdb)(CPUS390XState *env, uint64_t ah, uint64_t al,
256
                      uint64_t f2)
257
{
258
    float128 ret = float64_to_float128(f2, &env->fpu_status);
259
    ret = float128_mul(make_float128(ah, al), ret, &env->fpu_status);
260
    handle_exceptions(env, GETPC());
261
    return RET128(ret);
262
}
263

    
264
/* convert 32-bit float to 64-bit float */
265
uint64_t HELPER(ldeb)(CPUS390XState *env, uint64_t f2)
266
{
267
    float64 ret = float32_to_float64(f2, &env->fpu_status);
268
    handle_exceptions(env, GETPC());
269
    return ret;
270
}
271

    
272
/* convert 128-bit float to 64-bit float */
273
uint64_t HELPER(ldxb)(CPUS390XState *env, uint64_t ah, uint64_t al)
274
{
275
    float64 ret = float128_to_float64(make_float128(ah, al), &env->fpu_status);
276
    handle_exceptions(env, GETPC());
277
    return ret;
278
}
279

    
280
/* convert 64-bit float to 128-bit float */
281
uint64_t HELPER(lxdb)(CPUS390XState *env, uint64_t f2)
282
{
283
    float128 ret = float64_to_float128(f2, &env->fpu_status);
284
    handle_exceptions(env, GETPC());
285
    return RET128(ret);
286
}
287

    
288
/* convert 32-bit float to 128-bit float */
289
uint64_t HELPER(lxeb)(CPUS390XState *env, uint64_t f2)
290
{
291
    float128 ret = float32_to_float128(f2, &env->fpu_status);
292
    handle_exceptions(env, GETPC());
293
    return RET128(ret);
294
}
295

    
296
/* convert 64-bit float to 32-bit float */
297
uint64_t HELPER(ledb)(CPUS390XState *env, uint64_t f2)
298
{
299
    float32 ret = float64_to_float32(f2, &env->fpu_status);
300
    handle_exceptions(env, GETPC());
301
    return ret;
302
}
303

    
304
/* convert 128-bit float to 32-bit float */
305
uint64_t HELPER(lexb)(CPUS390XState *env, uint64_t ah, uint64_t al)
306
{
307
    float32 ret = float128_to_float32(make_float128(ah, al), &env->fpu_status);
308
    handle_exceptions(env, GETPC());
309
    return ret;
310
}
311

    
312
/* 32-bit FP compare */
313
uint32_t HELPER(ceb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
314
{
315
    int cmp = float32_compare_quiet(f1, f2, &env->fpu_status);
316
    handle_exceptions(env, GETPC());
317
    return float_comp_to_cc(env, cmp);
318
}
319

    
320
/* 64-bit FP compare */
321
uint32_t HELPER(cdb)(CPUS390XState *env, uint64_t f1, uint64_t f2)
322
{
323
    int cmp = float64_compare_quiet(f1, f2, &env->fpu_status);
324
    handle_exceptions(env, GETPC());
325
    return float_comp_to_cc(env, cmp);
326
}
327

    
328
/* 128-bit FP compare */
329
uint32_t HELPER(cxb)(CPUS390XState *env, uint64_t ah, uint64_t al,
330
                     uint64_t bh, uint64_t bl)
331
{
332
    int cmp = float128_compare_quiet(make_float128(ah, al),
333
                                     make_float128(bh, bl),
334
                                     &env->fpu_status);
335
    handle_exceptions(env, GETPC());
336
    return float_comp_to_cc(env, cmp);
337
}
338

    
339
static int swap_round_mode(CPUS390XState *env, int m3)
340
{
341
    int ret = env->fpu_status.float_rounding_mode;
342
    switch (m3) {
343
    case 0:
344
        /* current mode */
345
        break;
346
    case 1:
347
        /* biased round no nearest */
348
    case 4:
349
        /* round to nearest */
350
        set_float_rounding_mode(float_round_nearest_even, &env->fpu_status);
351
        break;
352
    case 5:
353
        /* round to zero */
354
        set_float_rounding_mode(float_round_to_zero, &env->fpu_status);
355
        break;
356
    case 6:
357
        /* round to +inf */
358
        set_float_rounding_mode(float_round_up, &env->fpu_status);
359
        break;
360
    case 7:
361
        /* round to -inf */
362
        set_float_rounding_mode(float_round_down, &env->fpu_status);
363
        break;
364
    }
365
    return ret;
366
}
367

    
368
/* convert 64-bit int to 32-bit float */
369
uint64_t HELPER(cegb)(CPUS390XState *env, int64_t v2, uint32_t m3)
370
{
371
    int hold = swap_round_mode(env, m3);
372
    float32 ret = int64_to_float32(v2, &env->fpu_status);
373
    set_float_rounding_mode(hold, &env->fpu_status);
374
    handle_exceptions(env, GETPC());
375
    return ret;
376
}
377

    
378
/* convert 64-bit int to 64-bit float */
379
uint64_t HELPER(cdgb)(CPUS390XState *env, int64_t v2, uint32_t m3)
380
{
381
    int hold = swap_round_mode(env, m3);
382
    float64 ret = int64_to_float64(v2, &env->fpu_status);
383
    set_float_rounding_mode(hold, &env->fpu_status);
384
    handle_exceptions(env, GETPC());
385
    return ret;
386
}
387

    
388
/* convert 64-bit int to 128-bit float */
389
uint64_t HELPER(cxgb)(CPUS390XState *env, int64_t v2, uint32_t m3)
390
{
391
    int hold = swap_round_mode(env, m3);
392
    float128 ret = int64_to_float128(v2, &env->fpu_status);
393
    set_float_rounding_mode(hold, &env->fpu_status);
394
    handle_exceptions(env, GETPC());
395
    return RET128(ret);
396
}
397

    
398
/* convert 64-bit uint to 32-bit float */
399
uint64_t HELPER(celgb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
400
{
401
    int hold = swap_round_mode(env, m3);
402
    float32 ret = uint64_to_float32(v2, &env->fpu_status);
403
    set_float_rounding_mode(hold, &env->fpu_status);
404
    handle_exceptions(env, GETPC());
405
    return ret;
406
}
407

    
408
/* convert 64-bit uint to 64-bit float */
409
uint64_t HELPER(cdlgb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
410
{
411
    int hold = swap_round_mode(env, m3);
412
    float64 ret = uint64_to_float64(v2, &env->fpu_status);
413
    set_float_rounding_mode(hold, &env->fpu_status);
414
    handle_exceptions(env, GETPC());
415
    return ret;
416
}
417

    
418
/* convert 64-bit uint to 128-bit float */
419
uint64_t HELPER(cxlgb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
420
{
421
    int hold = swap_round_mode(env, m3);
422
    float128 ret = uint64_to_float128(v2, &env->fpu_status);
423
    set_float_rounding_mode(hold, &env->fpu_status);
424
    handle_exceptions(env, GETPC());
425
    return RET128(ret);
426
}
427

    
428
/* convert 32-bit float to 64-bit int */
429
uint64_t HELPER(cgeb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
430
{
431
    int hold = swap_round_mode(env, m3);
432
    int64_t ret = float32_to_int64(v2, &env->fpu_status);
433
    set_float_rounding_mode(hold, &env->fpu_status);
434
    handle_exceptions(env, GETPC());
435
    return ret;
436
}
437

    
438
/* convert 64-bit float to 64-bit int */
439
uint64_t HELPER(cgdb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
440
{
441
    int hold = swap_round_mode(env, m3);
442
    int64_t ret = float64_to_int64(v2, &env->fpu_status);
443
    set_float_rounding_mode(hold, &env->fpu_status);
444
    handle_exceptions(env, GETPC());
445
    return ret;
446
}
447

    
448
/* convert 128-bit float to 64-bit int */
449
uint64_t HELPER(cgxb)(CPUS390XState *env, uint64_t h, uint64_t l, uint32_t m3)
450
{
451
    int hold = swap_round_mode(env, m3);
452
    float128 v2 = make_float128(h, l);
453
    int64_t ret = float128_to_int64(v2, &env->fpu_status);
454
    set_float_rounding_mode(hold, &env->fpu_status);
455
    handle_exceptions(env, GETPC());
456
    return ret;
457
}
458

    
459
/* convert 32-bit float to 32-bit int */
460
uint64_t HELPER(cfeb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
461
{
462
    int hold = swap_round_mode(env, m3);
463
    int32_t ret = float32_to_int32(v2, &env->fpu_status);
464
    set_float_rounding_mode(hold, &env->fpu_status);
465
    handle_exceptions(env, GETPC());
466
    return ret;
467
}
468

    
469
/* convert 64-bit float to 32-bit int */
470
uint64_t HELPER(cfdb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
471
{
472
    int hold = swap_round_mode(env, m3);
473
    int32_t ret = float64_to_int32(v2, &env->fpu_status);
474
    set_float_rounding_mode(hold, &env->fpu_status);
475
    handle_exceptions(env, GETPC());
476
    return ret;
477
}
478

    
479
/* convert 128-bit float to 32-bit int */
480
uint64_t HELPER(cfxb)(CPUS390XState *env, uint64_t h, uint64_t l, uint32_t m3)
481
{
482
    int hold = swap_round_mode(env, m3);
483
    float128 v2 = make_float128(h, l);
484
    int32_t ret = float128_to_int32(v2, &env->fpu_status);
485
    set_float_rounding_mode(hold, &env->fpu_status);
486
    handle_exceptions(env, GETPC());
487
    return ret;
488
}
489

    
490
/* convert 32-bit float to 64-bit uint */
491
uint64_t HELPER(clgeb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
492
{
493
    int hold = swap_round_mode(env, m3);
494
    uint64_t ret;
495
    v2 = float32_to_float64(v2, &env->fpu_status);
496
    ret = float64_to_uint64(v2, &env->fpu_status);
497
    set_float_rounding_mode(hold, &env->fpu_status);
498
    handle_exceptions(env, GETPC());
499
    return ret;
500
}
501

    
502
/* convert 64-bit float to 64-bit uint */
503
uint64_t HELPER(clgdb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
504
{
505
    int hold = swap_round_mode(env, m3);
506
    uint64_t ret = float64_to_uint64(v2, &env->fpu_status);
507
    set_float_rounding_mode(hold, &env->fpu_status);
508
    handle_exceptions(env, GETPC());
509
    return ret;
510
}
511

    
512
/* convert 128-bit float to 64-bit uint */
513
uint64_t HELPER(clgxb)(CPUS390XState *env, uint64_t h, uint64_t l, uint32_t m3)
514
{
515
    int hold = swap_round_mode(env, m3);
516
    float128 v2 = make_float128(h, l);
517
    /* ??? Not 100% correct.  */
518
    uint64_t ret = float128_to_int64(v2, &env->fpu_status);
519
    set_float_rounding_mode(hold, &env->fpu_status);
520
    handle_exceptions(env, GETPC());
521
    return ret;
522
}
523

    
524
/* convert 32-bit float to 32-bit uint */
525
uint64_t HELPER(clfeb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
526
{
527
    int hold = swap_round_mode(env, m3);
528
    uint32_t ret = float32_to_uint32(v2, &env->fpu_status);
529
    set_float_rounding_mode(hold, &env->fpu_status);
530
    handle_exceptions(env, GETPC());
531
    return ret;
532
}
533

    
534
/* convert 64-bit float to 32-bit uint */
535
uint64_t HELPER(clfdb)(CPUS390XState *env, uint64_t v2, uint32_t m3)
536
{
537
    int hold = swap_round_mode(env, m3);
538
    uint32_t ret = float64_to_uint32(v2, &env->fpu_status);
539
    set_float_rounding_mode(hold, &env->fpu_status);
540
    handle_exceptions(env, GETPC());
541
    return ret;
542
}
543

    
544
/* convert 128-bit float to 32-bit uint */
545
uint64_t HELPER(clfxb)(CPUS390XState *env, uint64_t h, uint64_t l, uint32_t m3)
546
{
547
    int hold = swap_round_mode(env, m3);
548
    float128 v2 = make_float128(h, l);
549
    /* Not 100% correct.  */
550
    uint32_t ret = float128_to_int64(v2, &env->fpu_status);
551
    set_float_rounding_mode(hold, &env->fpu_status);
552
    handle_exceptions(env, GETPC());
553
    return ret;
554
}
555

    
556
/* 32-bit FP multiply and add */
557
uint64_t HELPER(maeb)(CPUS390XState *env, uint64_t f1,
558
                      uint64_t f2, uint64_t f3)
559
{
560
    float32 ret = float32_muladd(f2, f3, f1, 0, &env->fpu_status);
561
    handle_exceptions(env, GETPC());
562
    return ret;
563
}
564

    
565
/* 64-bit FP multiply and add */
566
uint64_t HELPER(madb)(CPUS390XState *env, uint64_t f1,
567
                      uint64_t f2, uint64_t f3)
568
{
569
    float64 ret = float64_muladd(f2, f3, f1, 0, &env->fpu_status);
570
    handle_exceptions(env, GETPC());
571
    return ret;
572
}
573

    
574
/* 32-bit FP multiply and subtract */
575
uint64_t HELPER(mseb)(CPUS390XState *env, uint64_t f1,
576
                      uint64_t f2, uint64_t f3)
577
{
578
    float32 ret = float32_muladd(f2, f3, f1, float_muladd_negate_c,
579
                                 &env->fpu_status);
580
    handle_exceptions(env, GETPC());
581
    return ret;
582
}
583

    
584
/* 64-bit FP multiply and subtract */
585
uint64_t HELPER(msdb)(CPUS390XState *env, uint64_t f1,
586
                      uint64_t f2, uint64_t f3)
587
{
588
    float64 ret = float64_muladd(f2, f3, f1, float_muladd_negate_c,
589
                                 &env->fpu_status);
590
    handle_exceptions(env, GETPC());
591
    return ret;
592
}
593

    
594
/* test data class 32-bit */
595
uint32_t HELPER(tceb)(uint64_t f1, uint64_t m2)
596
{
597
    float32 v1 = f1;
598
    int neg = float32_is_neg(v1);
599
    uint32_t cc = 0;
600

    
601
    if ((float32_is_zero(v1) && (m2 & (1 << (11-neg)))) ||
602
        (float32_is_infinity(v1) && (m2 & (1 << (5-neg)))) ||
603
        (float32_is_any_nan(v1) && (m2 & (1 << (3-neg)))) ||
604
        (float32_is_signaling_nan(v1) && (m2 & (1 << (1-neg))))) {
605
        cc = 1;
606
    } else if (m2 & (1 << (9-neg))) {
607
        /* assume normalized number */
608
        cc = 1;
609
    }
610
    /* FIXME: denormalized? */
611
    return cc;
612
}
613

    
614
/* test data class 64-bit */
615
uint32_t HELPER(tcdb)(uint64_t v1, uint64_t m2)
616
{
617
    int neg = float64_is_neg(v1);
618
    uint32_t cc = 0;
619

    
620
    if ((float64_is_zero(v1) && (m2 & (1 << (11-neg)))) ||
621
        (float64_is_infinity(v1) && (m2 & (1 << (5-neg)))) ||
622
        (float64_is_any_nan(v1) && (m2 & (1 << (3-neg)))) ||
623
        (float64_is_signaling_nan(v1) && (m2 & (1 << (1-neg))))) {
624
        cc = 1;
625
    } else if (m2 & (1 << (9-neg))) {
626
        /* assume normalized number */
627
        cc = 1;
628
    }
629
    /* FIXME: denormalized? */
630
    return cc;
631
}
632

    
633
/* test data class 128-bit */
634
uint32_t HELPER(tcxb)(uint64_t ah, uint64_t al, uint64_t m2)
635
{
636
    float128 v1 = make_float128(ah, al);
637
    int neg = float128_is_neg(v1);
638
    uint32_t cc = 0;
639

    
640
    if ((float128_is_zero(v1) && (m2 & (1 << (11-neg)))) ||
641
        (float128_is_infinity(v1) && (m2 & (1 << (5-neg)))) ||
642
        (float128_is_any_nan(v1) && (m2 & (1 << (3-neg)))) ||
643
        (float128_is_signaling_nan(v1) && (m2 & (1 << (1-neg))))) {
644
        cc = 1;
645
    } else if (m2 & (1 << (9-neg))) {
646
        /* assume normalized number */
647
        cc = 1;
648
    }
649
    /* FIXME: denormalized? */
650
    return cc;
651
}
652

    
653
/* square root 32-bit */
654
uint64_t HELPER(sqeb)(CPUS390XState *env, uint64_t f2)
655
{
656
    float32 ret = float32_sqrt(f2, &env->fpu_status);
657
    handle_exceptions(env, GETPC());
658
    return ret;
659
}
660

    
661
/* square root 64-bit */
662
uint64_t HELPER(sqdb)(CPUS390XState *env, uint64_t f2)
663
{
664
    float64 ret = float64_sqrt(f2, &env->fpu_status);
665
    handle_exceptions(env, GETPC());
666
    return ret;
667
}
668

    
669
/* square root 128-bit */
670
uint64_t HELPER(sqxb)(CPUS390XState *env, uint64_t ah, uint64_t al)
671
{
672
    float128 ret = float128_sqrt(make_float128(ah, al), &env->fpu_status);
673
    handle_exceptions(env, GETPC());
674
    return RET128(ret);
675
}
676

    
677
static const int fpc_to_rnd[4] = {
678
    float_round_nearest_even,
679
    float_round_to_zero,
680
    float_round_up,
681
    float_round_down
682
};
683

    
684
/* set fpc */
685
void HELPER(sfpc)(CPUS390XState *env, uint64_t fpc)
686
{
687
    /* Install everything in the main FPC.  */
688
    env->fpc = fpc;
689

    
690
    /* Install the rounding mode in the shadow fpu_status.  */
691
    set_float_rounding_mode(fpc_to_rnd[fpc & 3], &env->fpu_status);
692
}
693

    
694
/* set fpc and signal */
695
void HELPER(sfas)(CPUS390XState *env, uint64_t val)
696
{
697
    uint32_t signalling = env->fpc;
698
    uint32_t source = val;
699
    uint32_t s390_exc;
700

    
701
    /* The contents of the source operand are placed in the FPC register;
702
       then the flags in the FPC register are set to the logical OR of the
703
       signalling flags and the source flags.  */
704
    env->fpc = source | (signalling & 0x00ff0000);
705
    set_float_rounding_mode(fpc_to_rnd[source & 3], &env->fpu_status);
706

    
707
    /* If any signalling flag is 1 and the corresponding source mask
708
       is also 1, a simulated-iee-exception trap occurs.  */
709
    s390_exc = (signalling >> 16) & (source >> 24);
710
    if (s390_exc) {
711
        ieee_exception(env, s390_exc | 3, GETPC());
712
    }
713
}